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.展开更多
The harmonic balance method(HBM)is one of the most widely used methods in solving nonlinear vibration problems,and its accuracy and computational efficiency largely depend on the number of the harmonics selected.The a...The harmonic balance method(HBM)is one of the most widely used methods in solving nonlinear vibration problems,and its accuracy and computational efficiency largely depend on the number of the harmonics selected.The adaptive harmonic balance(AHB)method is an improved HBM method.This paper presents a modified AHB method with the asymptotic harmonic selection(AHS)procedure.This new harmonic selection procedure selects harmonics from the frequency spectra of nonlinear terms instead of estimating the contribution of each harmonic to the whole nonlinear response,by which the additional calculation is avoided.A modified continuation method is proposed to deal with the variable size of nonlinear algebraic equations at different values of path parameters,and then all solution branches of the amplitude-frequency response are obtained.Numerical experiments are carried out to verify the performance of the AHB-AHS method.Five typical nonlinear dynamic equations with different types of nonlinearities and excitations are chosen as the illustrative examples.Compared with the classical HBM and Runge-Kutta methods,the proposed AHB-AHS method is of higher accuracy and better convergence.The AHB-AHS method proposed in this paper has the potential to investigate the nonlinear vibrations of complex high-dimensional nonlinear systems.展开更多
The harmonic balance(HB)method is one of the most commonly used methods for solving periodic solutions of both weakly and strongly nonlinear dynamical systems.However,it is confined to low-order approximations due to ...The harmonic balance(HB)method is one of the most commonly used methods for solving periodic solutions of both weakly and strongly nonlinear dynamical systems.However,it is confined to low-order approximations due to complex symbolic operations.Many variants have been developed to improve the HB method,among which the time domain HB-like methods are regarded as crucial improvements because of their fast computation and simple derivation.So far,there are two problems remaining to be addressed.i)A dozen of different versions of HB-like methods,in frequency domain or time domain or in hybrid,have been developed;unfortunately,misclassification pervades among them due to the unclear borderlines of different methods.ii)The time domain HB-like methods suffer from non-physical solutions,which have been shown to be caused by aliasing(mixture of the high-order into the low-order harmonics).Although a series of dealiasing techniques have been developed over the past two decades,the mechanism of aliasing and the final solution to dealiasing are still not well known to the academic community.This paper aims to provide a comprehensive review of the development of HB-like methods and enunciate their principal differences.In particular,the time domain methods are emphasized with the famous aliasing phenomenon clearly addressed.展开更多
A Josephson traveling wave parametric amplifier(JTWPA),which is a quantum-limited amplifier with high gain and large bandwidth,is the core device of large-scale measurement and control systems for quantum computing.A ...A Josephson traveling wave parametric amplifier(JTWPA),which is a quantum-limited amplifier with high gain and large bandwidth,is the core device of large-scale measurement and control systems for quantum computing.A typical JTWPA consists of thousands of Josephson junctions connected in series to form a transmission line and hundreds of shunt LC resonators periodically loaded along the line for phase matching.Because the variation of these capacitors and inductors can be detrimental to their high-frequency characteristics,the fabrication of a JTWPA typically necessitates precise processing equipment.To guide the fabrication process and further improve the design for manufacturability,it is necessary to understand how each electronic component affects the amplifier.In this paper,we use the harmonic balance method to conduct a comprehensive study on the impact of nonuniformity and fabrication yield of the electronic components on the performance of a JTWPA.The results provide insightful and scientific guidance for device design and fabrication processes.展开更多
The dispersed phase in multiphase flows can be modeled by the population balance model(PBM). A typical population balance equation(PBE) contains terms for spatial transport, loss/growth and breakage/coalescence source...The dispersed phase in multiphase flows can be modeled by the population balance model(PBM). A typical population balance equation(PBE) contains terms for spatial transport, loss/growth and breakage/coalescence source terms. The equation is therefore quite complex and difficult to solve analytically or numerically. The quadrature-based moment methods(QBMMs) are a class of methods that solve the PBE by converting the transport equation of the number density function(NDF) into moment transport equations. The unknown source terms are closed by numerical quadrature. Over the years, many QBMMs have been developed for different problems, such as the quadrature method of moments(QMOM), direct quadrature method of moments(DQMOM),extended quadrature method of moments(EQMOM), conditional quadrature method of moments(CQMOM),extended conditional quadrature method of moments(ECQMOM) and hyperbolic quadrature method of moments(Hy QMOM). In this paper, we present a comprehensive algorithm review of these QBMMs. The mathematical equations for spatially homogeneous systems with first-order point processes and second-order point processes are derived in detail. The algorithms are further extended to the inhomogeneous system for multiphase flows, in which the computational fluid dynamics(CFD) can be coupled with the PBE. The physical limitations and the challenging numerical problems of these QBMMs are discussed. Possible solutions are also summarized.展开更多
Comparisons of the common methods for obtaining the periodic responses show that the harmonic balance method with alternating frequency/time(HB-AFT) domain technique has some advantages in dealing with nonlinear probl...Comparisons of the common methods for obtaining the periodic responses show that the harmonic balance method with alternating frequency/time(HB-AFT) domain technique has some advantages in dealing with nonlinear problems of fractional exponential models. By the HB-AFT method, a rigid rotor supported by ball bearings with nonlinearity of Hertz contact and ball passage vibrations is considered. With the aid of the Floquet theory, the movement characteristics of interval stability are deeply studied. Besides, a simple strategy to determine the monodromy matrix is proposed for the stability analysis.展开更多
A nonlinear transformation of the Whitham-Broer-Kaup (WBK) model equations in the shallow water small-amplitude regime is derived by using a simplified homogeneous balance method. The WBK model equations are linearize...A nonlinear transformation of the Whitham-Broer-Kaup (WBK) model equations in the shallow water small-amplitude regime is derived by using a simplified homogeneous balance method. The WBK model equations are linearized under the nonlinear transformation. Various exact solutions of the WBK model equations are obtained via the nonlinear transformation with the aid of solutions for the linear equation.展开更多
In this study, forced nonlinear vibration of a circular micro-plate under two-sided electrostatic, two-sided Casimir and external harmonic forces is investigated analytically. For this purpose, at first, von Karman pl...In this study, forced nonlinear vibration of a circular micro-plate under two-sided electrostatic, two-sided Casimir and external harmonic forces is investigated analytically. For this purpose, at first, von Karman plate theory including geometrical nonlinearity is used to obtain the deflection of the micro-plate. Galerkin decomposition method is then employed, and nonlinear ordinary differential equations (ODEs) of motion are determined. A harmonic balance method (HBM) is applied to equations and analytical relation for nonlineaT frequency response (F-R) curves are derived for two categories (including and neglecting Casimir force) separately. The analytical results for three cases:(1) semi-linear vibration;(2) weakly nonlinear vibration;(3) highly non linear vibration, are validated by comparing with the numerical solutio ns. After validation, the effects of the voltage and Casimir force on the natural frequency of two-sided capacitor system are investigated. It is shown that by assuming Casimir force in small gap distances, reduction of the natural frequency is considerable. The influences of the applied voltage, damping, micro-plate thickness and Casimir force on the frequency response curves have been presented too. The results of this study can be useful for modeling circular parallel-plates in nano /microelectromechanical transducers such as microphones and pressure sensors.展开更多
The dynamical properties of fractional-order Duffing–van der Pol oscillator are studied, and the amplitude–frequency response equation of primary resonance is obtained by the harmonic balance method. The stability c...The dynamical properties of fractional-order Duffing–van der Pol oscillator are studied, and the amplitude–frequency response equation of primary resonance is obtained by the harmonic balance method. The stability condition for steady-state solution is obtained based on Lyapunov theory. The comparison of the approximate analytical results with the numerical results is fulfilled, and the approximations obtained are in good agreement with the numerical solutions. The bifurcations of primary resonance for system parameters are analyzed. The results show that the harmonic balance method is effective and convenient for solving this problem, and it provides a reference for the dynamical analysis of similar nonlinear systems.展开更多
The population balance modeling is regarded as a universally accepted mathematical framework for dynamic simulation of various particulate processes, such as crystallization, granulation and polymerization. This artic...The population balance modeling is regarded as a universally accepted mathematical framework for dynamic simulation of various particulate processes, such as crystallization, granulation and polymerization. This article is concerned with the application of the method of characteristics (MOC) for solving population balance models describing batch crystallization process. The growth and nucleation are considered as dominant phenomena, while the breakage and aggregation are neglected. The numerical solutions of such PBEs require high order accuracy due to the occurrence of steep moving fronts and narrow peaks in the solutions. The MOC has been found to be a very effective technique for resolving sharp discontinuities. Different case studies are carried out to analyze the accuracy of proposed algorithm. For validation, the results of MOC are compared with the available analytical solutions and the results of finite volume schemes. The results of MOC were found to be in good agreement with analytical solutions and superior than those obtained by finite volume schemes.展开更多
Positive results are proved here about the ability of balanced methods to reproduce the mean square stability of the impulsive stochastic differential equations. It is shown that the balanced methods with strong conve...Positive results are proved here about the ability of balanced methods to reproduce the mean square stability of the impulsive stochastic differential equations. It is shown that the balanced methods with strong convergence can preserve the mean square stability with the sufficiently small stepsize. Weak variants and their mean square stability are also considered. Several numerical experiments are given for illustration and show that the fully implicit methods are superior to those of the explicit methods in terms of mean-square stabilities for relatively large stepsizes especially.展开更多
Two types of existing iterative methods for solving the nonlinear balance equation(NBE)are revisited.In the first type,the NBE is rearranged into a linearized equation for a presumably small correction to the initial ...Two types of existing iterative methods for solving the nonlinear balance equation(NBE)are revisited.In the first type,the NBE is rearranged into a linearized equation for a presumably small correction to the initial guess or the subsequent updated solution.In the second type,the NBE is rearranged into a quadratic form of the absolute vorticity with the positive root of this quadratic form used in the form of a Poisson equation to solve NBE iteratively.The two methods are rederived by expanding the solution asymptotically upon a small Rossby number,and a criterion for optimally truncating the asymptotic expansion is proposed to obtain the super-asymptotic approximation of the solution.For each rederived method,two iterative procedures are designed using the integral-form Poisson solver versus the over-relaxation scheme to solve the boundary value problem in each iteration.Upon testing with analytically formulated wavering jet flows on the synoptic,sub-synoptic and meso-αscales,the iterative procedure designed for the first method with the Poisson solver,named M1a,is found to be the most accurate and efficient.For the synoptic wavering jet flow in which the NBE is entirely elliptic,M1a is extremely accurate.For the sub-synoptic wavering jet flow in which the NBE is mostly elliptic,M1a is sufficiently accurate.For the meso-αwavering jet flow in which the NBE is partially hyperbolic so its boundary value problem becomes seriously ill-posed,M1a can effectively reduce the solution error for the cyclonically curved part of the wavering jet flow,but not for the anti-cyclonically curved part.展开更多
We have deduced an iteration scheme in theincremental harmonic balance (IHB) method using the harmonicbalance plus the Newton-Raphson method. Since theconvergence of the iteration is dependent upon the initial valuesi...We have deduced an iteration scheme in theincremental harmonic balance (IHB) method using the harmonicbalance plus the Newton-Raphson method. Since theconvergence of the iteration is dependent upon the initial valuesin the iteration, the convergent region is greatly restrictedfor some cases. In this contribution, in order to enlarge theconvergent region of the IHB method, we constructed thezeroth-order deformation equation using the homotopy analysismethod, in which the IHB method is employed to solvethe deformation equation with an embedding parameter asthe active increment. Taking the Duffing and the van derPol equations as examples, we obtained the highly accuratesolutions. Importantly, the presented approach renders a convenientway to control and adjust the convergence.展开更多
We applied the material balance principle of the denudation volume and sedimentary flux to study the denudation-accumulation system between the Longmen Mountains(Mts.) and the foreland basin. The amount of sediment in...We applied the material balance principle of the denudation volume and sedimentary flux to study the denudation-accumulation system between the Longmen Mountains(Mts.) and the foreland basin. The amount of sediment in each sedimentation stage of the basin was estimated to obtain the denudation volume,erosion thickness and deposit thickness since the Late Triassic Epoch,to enable us to recover the paleoelevation of the provenance and the sedimentary area. The results show the following:(1) Since the Late Triassic Epoch,the elevation of the surface of the Longmen Mts. has uplifted from 0 m to 2751 m,and the crust of the Longmen Mts. has uplifted by 9.8 km. Approximately 72% of the materials introduced have been denuded from the mountains.(2) It is difficult to recover the paleoelevation of each stage of the Longmen Mts. foreland basin quantitatively by the present-day techniques and data.(3) The formation of the Longmen Mts. foreland basin consisted of three stages of thrust belt tectonic load and three stages of thrust belt erosional unload. During tectonic loading stages(Late Triassic Epoch,Late Jurassic–Early Cretaceous,Late Cretaceous–Miocene),the average elevation of Longmen Mts. was lower(approximately 700–1700 m). During erosional unloading stages(Early and Middle Jurassic,Middle Cretaceous and Jiaguan,Late Cenozoic),the average elevation of Longmen Mts. was high at approximately 2000–2800m.展开更多
As the main unconventional natural gas reservoirs,shale gas reservoirs and coalbed methane(CBM)reservoirs belong to adsorptive gas reservoirs,i.e.,gas reservoirs containing adsorbed gas.Shale gas and CBM reservoirs us...As the main unconventional natural gas reservoirs,shale gas reservoirs and coalbed methane(CBM)reservoirs belong to adsorptive gas reservoirs,i.e.,gas reservoirs containing adsorbed gas.Shale gas and CBM reservoirs usually have the characteristics of rich adsorbed gas and obvious dynamic changes of porosity and permeability.A generalized material balance equation and the corresponding reserve evaluation method considering all the mechanisms for both shale gas reservoirs and CBM reservoirs are necessary.In this work,a generalized material balance equation(GMBE)considering the effects of critical desorption pressure,stress sensitivity,matrix shrinkage,water production,water influx,and solubility of natural gas in water is established.Then,by converting the GMBE to a linear relationship between two parameter groups related with known formation/fluid properties and dynamic performance data,the straight-line reserve evaluation method is proposed.By using the slope and the y-intercept of this straight line,the original adsorbed gas in place(OAGIP),original free gas in place(OFGIP),original dissolved gas in place(ODGIP),and the original gas in place(OGIP)can be quickly calculated.Third,two validation cases for shale gas reservoir and CBM reservoir are conducted using commercial reservoir simulator and the coalbed methane dynamic performance analysis software,respectively.Finally,two field studies in the Fuling shale gas field and the Baode CBM field are presented.Results show that the GMBE and the corresponding straight-line reserve evaluation method are rational,accurate,and effective for both shale gas reservoirs and CBM reservoirs.More detailed information about reserves of shale gas and CBM reservoirs can be clarified,and only the straight-line fitting approach is used to determine all kinds of reserves without iteration,proving that the proposed method has great advantages compared with other current methods.展开更多
The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered...The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered in the flutter equations of two-dimensional airfoil. First, the equations were transferred into matrix form, then the vibration process was divided into the persistent incremental processes of vibration moments. And the expression of their solutions could be obtained by using a certain amplitude as control parameter in the harmonic balance process, and then the bifurcation, limit cycle flutter phenomena and the number of harmonic terms were analyzed. Finally, numerical results calculated by the Runge-Kutta method were given to verify the results obtained by the proposed procedure. It has been shown that the incremental harmonic method is effective and precise in the analysis of strongly nonlinear flutter with multiple structural nonlinearities.展开更多
A modified homogeneous balance method is proposed by improving some key steps in the homogeneousbalance method.Bilinear equations of some nonlinear evolution equations are derived by using the modified homogeneousbala...A modified homogeneous balance method is proposed by improving some key steps in the homogeneousbalance method.Bilinear equations of some nonlinear evolution equations are derived by using the modified homogeneousbalance method.Generalized Boussinesq equation,KP equation,and mKdV equation are chosen as examples to illustrateour method.This approach is also applicable to a large variety of nonlinear evolution equations.展开更多
In order to identify the day and night pollution sources of PM10 in ambient air in Longyan City,the authors analyzed the elemental composition of respirable particulate matters in the day and night ambient air samples...In order to identify the day and night pollution sources of PM10 in ambient air in Longyan City,the authors analyzed the elemental composition of respirable particulate matters in the day and night ambient air samples and various pollution sources which were collected in January 2010 in Longyan with inductivity coupled plasma-mass spectrometry(ICP-MS).Then chemical mass balance(CMB) model and factor analysis(FA) method were applied to comparatively study the inorganic components in the sources and receptor samples.The results of factor analysis show that the major sources were road dust,waste incineration and mixed sources which contained automobile exhaust,soil dust/secondary dust and coal dust during the daytime in Longyan City,China.There are two major sources of pollution which are soil dust and mixture sources of automobile exhaust and secondary dust during the night in Longyan.The results of CMB show that the major sources are secondary dust,automobile exhaust and road dust during the daytime in Longyan.The major sources are secondary dust,soil dust and automobile exhaust during the night in Longyan.The results of the two methods are similar to each other and the results will guide us to plan to control the PM10 pollution sources in Longyan.展开更多
A class of general modified split-step balanced methods proposed in the paper can be applied to solve stiff stochastic differential systems with m-dimensional multiplicative noise.Compared to some other already report...A class of general modified split-step balanced methods proposed in the paper can be applied to solve stiff stochastic differential systems with m-dimensional multiplicative noise.Compared to some other already reported split-step balanced methods,the drift increment function of the methods can be taken from any chosen one-step ordinary differential equations(ODEs)solver.The schemes is proved to be strong convergent with order one.For the mean-square stability analysis,the investigation is confined to two cases.Some numerical experiments are reported to testify the performance and the effectiveness of the methods.展开更多
文摘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.
基金Project supported by the National Natural Science Foundation of China(Nos.11972129 and12372008)the National Major Science and Technology Projects of China(No.2017-IV-0008-0045)+3 种基金the Natural Science Foundation of Heilongjiang Province of China(No.YQ2022A008)the Fundamental Research Funds for the Central Universities of China(No.HIT.OCEF.2023006)the Polish National Science Centre of Poland under the OPUS 18 grant(No.2019/35/B/ST8/00980)the Tianjin University Independent Innovation Foundation of China(No.2023XJS-0038)。
文摘The harmonic balance method(HBM)is one of the most widely used methods in solving nonlinear vibration problems,and its accuracy and computational efficiency largely depend on the number of the harmonics selected.The adaptive harmonic balance(AHB)method is an improved HBM method.This paper presents a modified AHB method with the asymptotic harmonic selection(AHS)procedure.This new harmonic selection procedure selects harmonics from the frequency spectra of nonlinear terms instead of estimating the contribution of each harmonic to the whole nonlinear response,by which the additional calculation is avoided.A modified continuation method is proposed to deal with the variable size of nonlinear algebraic equations at different values of path parameters,and then all solution branches of the amplitude-frequency response are obtained.Numerical experiments are carried out to verify the performance of the AHB-AHS method.Five typical nonlinear dynamic equations with different types of nonlinearities and excitations are chosen as the illustrative examples.Compared with the classical HBM and Runge-Kutta methods,the proposed AHB-AHS method is of higher accuracy and better convergence.The AHB-AHS method proposed in this paper has the potential to investigate the nonlinear vibrations of complex high-dimensional nonlinear systems.
基金supported by National Key Research and Development Program of China(No.2021YFA0717100)NationalNatural Science Foundation of China(Nos.12072270,U2013206).
文摘The harmonic balance(HB)method is one of the most commonly used methods for solving periodic solutions of both weakly and strongly nonlinear dynamical systems.However,it is confined to low-order approximations due to complex symbolic operations.Many variants have been developed to improve the HB method,among which the time domain HB-like methods are regarded as crucial improvements because of their fast computation and simple derivation.So far,there are two problems remaining to be addressed.i)A dozen of different versions of HB-like methods,in frequency domain or time domain or in hybrid,have been developed;unfortunately,misclassification pervades among them due to the unclear borderlines of different methods.ii)The time domain HB-like methods suffer from non-physical solutions,which have been shown to be caused by aliasing(mixture of the high-order into the low-order harmonics).Although a series of dealiasing techniques have been developed over the past two decades,the mechanism of aliasing and the final solution to dealiasing are still not well known to the academic community.This paper aims to provide a comprehensive review of the development of HB-like methods and enunciate their principal differences.In particular,the time domain methods are emphasized with the famous aliasing phenomenon clearly addressed.
基金support from the Youth Innovation Promotion Association of Chinese Academy of Sciences (Grant No.2019319)support from the Start-up Foundation of Suzhou Institute of Nano-Tech and Nano-Bionics,CAS,Suzhou (Grant No.Y9AAD110)。
文摘A Josephson traveling wave parametric amplifier(JTWPA),which is a quantum-limited amplifier with high gain and large bandwidth,is the core device of large-scale measurement and control systems for quantum computing.A typical JTWPA consists of thousands of Josephson junctions connected in series to form a transmission line and hundreds of shunt LC resonators periodically loaded along the line for phase matching.Because the variation of these capacitors and inductors can be detrimental to their high-frequency characteristics,the fabrication of a JTWPA typically necessitates precise processing equipment.To guide the fabrication process and further improve the design for manufacturability,it is necessary to understand how each electronic component affects the amplifier.In this paper,we use the harmonic balance method to conduct a comprehensive study on the impact of nonuniformity and fabrication yield of the electronic components on the performance of a JTWPA.The results provide insightful and scientific guidance for device design and fabrication processes.
文摘The dispersed phase in multiphase flows can be modeled by the population balance model(PBM). A typical population balance equation(PBE) contains terms for spatial transport, loss/growth and breakage/coalescence source terms. The equation is therefore quite complex and difficult to solve analytically or numerically. The quadrature-based moment methods(QBMMs) are a class of methods that solve the PBE by converting the transport equation of the number density function(NDF) into moment transport equations. The unknown source terms are closed by numerical quadrature. Over the years, many QBMMs have been developed for different problems, such as the quadrature method of moments(QMOM), direct quadrature method of moments(DQMOM),extended quadrature method of moments(EQMOM), conditional quadrature method of moments(CQMOM),extended conditional quadrature method of moments(ECQMOM) and hyperbolic quadrature method of moments(Hy QMOM). In this paper, we present a comprehensive algorithm review of these QBMMs. The mathematical equations for spatially homogeneous systems with first-order point processes and second-order point processes are derived in detail. The algorithms are further extended to the inhomogeneous system for multiphase flows, in which the computational fluid dynamics(CFD) can be coupled with the PBE. The physical limitations and the challenging numerical problems of these QBMMs are discussed. Possible solutions are also summarized.
基金supported by the National Natural Science Foundation of China(No.10632040)
文摘Comparisons of the common methods for obtaining the periodic responses show that the harmonic balance method with alternating frequency/time(HB-AFT) domain technique has some advantages in dealing with nonlinear problems of fractional exponential models. By the HB-AFT method, a rigid rotor supported by ball bearings with nonlinearity of Hertz contact and ball passage vibrations is considered. With the aid of the Floquet theory, the movement characteristics of interval stability are deeply studied. Besides, a simple strategy to determine the monodromy matrix is proposed for the stability analysis.
文摘A nonlinear transformation of the Whitham-Broer-Kaup (WBK) model equations in the shallow water small-amplitude regime is derived by using a simplified homogeneous balance method. The WBK model equations are linearized under the nonlinear transformation. Various exact solutions of the WBK model equations are obtained via the nonlinear transformation with the aid of solutions for the linear equation.
文摘In this study, forced nonlinear vibration of a circular micro-plate under two-sided electrostatic, two-sided Casimir and external harmonic forces is investigated analytically. For this purpose, at first, von Karman plate theory including geometrical nonlinearity is used to obtain the deflection of the micro-plate. Galerkin decomposition method is then employed, and nonlinear ordinary differential equations (ODEs) of motion are determined. A harmonic balance method (HBM) is applied to equations and analytical relation for nonlineaT frequency response (F-R) curves are derived for two categories (including and neglecting Casimir force) separately. The analytical results for three cases:(1) semi-linear vibration;(2) weakly nonlinear vibration;(3) highly non linear vibration, are validated by comparing with the numerical solutio ns. After validation, the effects of the voltage and Casimir force on the natural frequency of two-sided capacitor system are investigated. It is shown that by assuming Casimir force in small gap distances, reduction of the natural frequency is considerable. The influences of the applied voltage, damping, micro-plate thickness and Casimir force on the frequency response curves have been presented too. The results of this study can be useful for modeling circular parallel-plates in nano /microelectromechanical transducers such as microphones and pressure sensors.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11872254 and 11672191)
文摘The dynamical properties of fractional-order Duffing–van der Pol oscillator are studied, and the amplitude–frequency response equation of primary resonance is obtained by the harmonic balance method. The stability condition for steady-state solution is obtained based on Lyapunov theory. The comparison of the approximate analytical results with the numerical results is fulfilled, and the approximations obtained are in good agreement with the numerical solutions. The bifurcations of primary resonance for system parameters are analyzed. The results show that the harmonic balance method is effective and convenient for solving this problem, and it provides a reference for the dynamical analysis of similar nonlinear systems.
文摘The population balance modeling is regarded as a universally accepted mathematical framework for dynamic simulation of various particulate processes, such as crystallization, granulation and polymerization. This article is concerned with the application of the method of characteristics (MOC) for solving population balance models describing batch crystallization process. The growth and nucleation are considered as dominant phenomena, while the breakage and aggregation are neglected. The numerical solutions of such PBEs require high order accuracy due to the occurrence of steep moving fronts and narrow peaks in the solutions. The MOC has been found to be a very effective technique for resolving sharp discontinuities. Different case studies are carried out to analyze the accuracy of proposed algorithm. For validation, the results of MOC are compared with the available analytical solutions and the results of finite volume schemes. The results of MOC were found to be in good agreement with analytical solutions and superior than those obtained by finite volume schemes.
基金National Natural Science Foundations of China(Nos.11561028,11101101,11461032,11401267)Natural Science Foundations of Jiangxi Province,China(Nos.20151BAB201013,20151BAB201010,20151BAB201015)
文摘Positive results are proved here about the ability of balanced methods to reproduce the mean square stability of the impulsive stochastic differential equations. It is shown that the balanced methods with strong convergence can preserve the mean square stability with the sufficiently small stepsize. Weak variants and their mean square stability are also considered. Several numerical experiments are given for illustration and show that the fully implicit methods are superior to those of the explicit methods in terms of mean-square stabilities for relatively large stepsizes especially.
基金the NSF of China Grants 91937301 and 41675060,the National Key Scientific and Technological Infrastructure Project"EarthLab",and the ONR Grants N000141712375 and N000142012449 to the University of Oklahoma(OU)The numerical experiments were performed at the OU supercomputer SchoonerCIMMS by NOAA/Office of Oceanic and Atmospheric Research under NOAA-OU Cooperative Agreement#NA110AR4320072,U.S.Department of Commerce.
文摘Two types of existing iterative methods for solving the nonlinear balance equation(NBE)are revisited.In the first type,the NBE is rearranged into a linearized equation for a presumably small correction to the initial guess or the subsequent updated solution.In the second type,the NBE is rearranged into a quadratic form of the absolute vorticity with the positive root of this quadratic form used in the form of a Poisson equation to solve NBE iteratively.The two methods are rederived by expanding the solution asymptotically upon a small Rossby number,and a criterion for optimally truncating the asymptotic expansion is proposed to obtain the super-asymptotic approximation of the solution.For each rederived method,two iterative procedures are designed using the integral-form Poisson solver versus the over-relaxation scheme to solve the boundary value problem in each iteration.Upon testing with analytically formulated wavering jet flows on the synoptic,sub-synoptic and meso-αscales,the iterative procedure designed for the first method with the Poisson solver,named M1a,is found to be the most accurate and efficient.For the synoptic wavering jet flow in which the NBE is entirely elliptic,M1a is extremely accurate.For the sub-synoptic wavering jet flow in which the NBE is mostly elliptic,M1a is sufficiently accurate.For the meso-αwavering jet flow in which the NBE is partially hyperbolic so its boundary value problem becomes seriously ill-posed,M1a can effectively reduce the solution error for the cyclonically curved part of the wavering jet flow,but not for the anti-cyclonically curved part.
基金supported by the National Natural Science Foundation of China (10772202)Doctoral Program Foundation of Ministry of Education of China (20050558032)Guangdong Province Natural Science Foundation (07003680, 05003295)
文摘We have deduced an iteration scheme in theincremental harmonic balance (IHB) method using the harmonicbalance plus the Newton-Raphson method. Since theconvergence of the iteration is dependent upon the initial valuesin the iteration, the convergent region is greatly restrictedfor some cases. In this contribution, in order to enlarge theconvergent region of the IHB method, we constructed thezeroth-order deformation equation using the homotopy analysismethod, in which the IHB method is employed to solvethe deformation equation with an embedding parameter asthe active increment. Taking the Duffing and the van derPol equations as examples, we obtained the highly accuratesolutions. Importantly, the presented approach renders a convenientway to control and adjust the convergence.
基金the Project of the National Natural Science Foudation of China (Grant No.41372114,41340005,41172162,40972083)
文摘We applied the material balance principle of the denudation volume and sedimentary flux to study the denudation-accumulation system between the Longmen Mountains(Mts.) and the foreland basin. The amount of sediment in each sedimentation stage of the basin was estimated to obtain the denudation volume,erosion thickness and deposit thickness since the Late Triassic Epoch,to enable us to recover the paleoelevation of the provenance and the sedimentary area. The results show the following:(1) Since the Late Triassic Epoch,the elevation of the surface of the Longmen Mts. has uplifted from 0 m to 2751 m,and the crust of the Longmen Mts. has uplifted by 9.8 km. Approximately 72% of the materials introduced have been denuded from the mountains.(2) It is difficult to recover the paleoelevation of each stage of the Longmen Mts. foreland basin quantitatively by the present-day techniques and data.(3) The formation of the Longmen Mts. foreland basin consisted of three stages of thrust belt tectonic load and three stages of thrust belt erosional unload. During tectonic loading stages(Late Triassic Epoch,Late Jurassic–Early Cretaceous,Late Cretaceous–Miocene),the average elevation of Longmen Mts. was lower(approximately 700–1700 m). During erosional unloading stages(Early and Middle Jurassic,Middle Cretaceous and Jiaguan,Late Cenozoic),the average elevation of Longmen Mts. was high at approximately 2000–2800m.
基金supported by Science and Technology Major Project of Shanxi Province,China(No.20201101002)Science and Technology Major Project of China,China(No.2016ZX05043002)+1 种基金National Natural Science Foundation Project of China,China(No.51874319)Science Foundation of China University of Petroleum(Beijing),China(No.2462020QNXZ003)to support part of this work
文摘As the main unconventional natural gas reservoirs,shale gas reservoirs and coalbed methane(CBM)reservoirs belong to adsorptive gas reservoirs,i.e.,gas reservoirs containing adsorbed gas.Shale gas and CBM reservoirs usually have the characteristics of rich adsorbed gas and obvious dynamic changes of porosity and permeability.A generalized material balance equation and the corresponding reserve evaluation method considering all the mechanisms for both shale gas reservoirs and CBM reservoirs are necessary.In this work,a generalized material balance equation(GMBE)considering the effects of critical desorption pressure,stress sensitivity,matrix shrinkage,water production,water influx,and solubility of natural gas in water is established.Then,by converting the GMBE to a linear relationship between two parameter groups related with known formation/fluid properties and dynamic performance data,the straight-line reserve evaluation method is proposed.By using the slope and the y-intercept of this straight line,the original adsorbed gas in place(OAGIP),original free gas in place(OFGIP),original dissolved gas in place(ODGIP),and the original gas in place(OGIP)can be quickly calculated.Third,two validation cases for shale gas reservoir and CBM reservoir are conducted using commercial reservoir simulator and the coalbed methane dynamic performance analysis software,respectively.Finally,two field studies in the Fuling shale gas field and the Baode CBM field are presented.Results show that the GMBE and the corresponding straight-line reserve evaluation method are rational,accurate,and effective for both shale gas reservoirs and CBM reservoirs.More detailed information about reserves of shale gas and CBM reservoirs can be clarified,and only the straight-line fitting approach is used to determine all kinds of reserves without iteration,proving that the proposed method has great advantages compared with other current methods.
基金Project supported by the Ph. D. Programs Foundation of Ministry of Education of China (No.20050558032) the Natural Science Foundation of Guangdong Province of China (No.05003295) the Foundation of Sun Yat-sen University Advanced Research Center (No.06M8) the Young Teacher Scientific Research Foundation of Sun Sat-sen University (No.1131011)
文摘The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered in the flutter equations of two-dimensional airfoil. First, the equations were transferred into matrix form, then the vibration process was divided into the persistent incremental processes of vibration moments. And the expression of their solutions could be obtained by using a certain amplitude as control parameter in the harmonic balance process, and then the bifurcation, limit cycle flutter phenomena and the number of harmonic terms were analyzed. Finally, numerical results calculated by the Runge-Kutta method were given to verify the results obtained by the proposed procedure. It has been shown that the incremental harmonic method is effective and precise in the analysis of strongly nonlinear flutter with multiple structural nonlinearities.
基金Supported by the National Natural Science Foundation of China under Grant No. 11071209the Natural Science Foundation of the Higher Education Institutions of Jiangsu Province under Grant No. 10KJBll0011
文摘A modified homogeneous balance method is proposed by improving some key steps in the homogeneousbalance method.Bilinear equations of some nonlinear evolution equations are derived by using the modified homogeneousbalance method.Generalized Boussinesq equation,KP equation,and mKdV equation are chosen as examples to illustrateour method.This approach is also applicable to a large variety of nonlinear evolution equations.
基金Supported by the Natural Basic Research Program of China(No.2005CB422207)the Fund of Eco-enviromental Impacts and Protection in Devoloping and Utilizing of Oil-shale Resources(No.OSR-01-06)
文摘In order to identify the day and night pollution sources of PM10 in ambient air in Longyan City,the authors analyzed the elemental composition of respirable particulate matters in the day and night ambient air samples and various pollution sources which were collected in January 2010 in Longyan with inductivity coupled plasma-mass spectrometry(ICP-MS).Then chemical mass balance(CMB) model and factor analysis(FA) method were applied to comparatively study the inorganic components in the sources and receptor samples.The results of factor analysis show that the major sources were road dust,waste incineration and mixed sources which contained automobile exhaust,soil dust/secondary dust and coal dust during the daytime in Longyan City,China.There are two major sources of pollution which are soil dust and mixture sources of automobile exhaust and secondary dust during the night in Longyan.The results of CMB show that the major sources are secondary dust,automobile exhaust and road dust during the daytime in Longyan.The major sources are secondary dust,soil dust and automobile exhaust during the night in Longyan.The results of the two methods are similar to each other and the results will guide us to plan to control the PM10 pollution sources in Longyan.
基金National Natural Science Foundation of China(No.11171352)
文摘A class of general modified split-step balanced methods proposed in the paper can be applied to solve stiff stochastic differential systems with m-dimensional multiplicative noise.Compared to some other already reported split-step balanced methods,the drift increment function of the methods can be taken from any chosen one-step ordinary differential equations(ODEs)solver.The schemes is proved to be strong convergent with order one.For the mean-square stability analysis,the investigation is confined to two cases.Some numerical experiments are reported to testify the performance and the effectiveness of the methods.
