Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

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.

Harmonic balance simulation of the influence of component uniformity and reliability on the performance of a Josephson traveling wave parametric amplifier

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.

A novel adaptive harmonic balance method with an asymptotic harmonic selection

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.

Harmonic Balance Methods:A Review and Recent Developments

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.

Quadrature-based moment methods for the population balance equation: An algorithm review

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.

Harmonic balance method with alternating frequency/time domain technique for nonlinear dynamical system with fractional exponential

Comparisons of the common methods for obtaining the periodic responses show that the harmonic balance method with alternating frequency/time （HB-AFT） do- main 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.

Improving convergence of incremental harmonic balance method using homotopy analysis method

We have deduced incremental harmonic balance an iteration scheme in the （IHB） method using the harmonic balance plus the Newton-Raphson method. Since the convergence of the iteration is dependent upon the initial values in the iteration, the convergent region is greatly restricted for some cases. In this contribution, in order to enlarge the convergent region of the IHB method, we constructed the zeroth-order deformation equation using the homotopy analysis method, in which the IHB method is employed to solve the deformation equation with an embedding parameter as the active increment. Taking the Duffing and the van der Pol equations as examples, we obtained the highly accurate solutions. Importantly, the presented approach renders a convenient way to control and adjust the convergence.

Determination of initial cable force of cantilever casting concrete arch bridge using stress balance and influence matrix methods

Cantilever casting concrete arch bridge using form traveller has a broad application prospect.However,it is difficult to obtain reasonable initial cable force in construction stage.In this study,stress balance and influence matrix methods were developed to determine the initial cable force of cantilever casting concrete arch bridge.The stress balance equation and influence matrix of arch rib critical section were established,and the buckle cable force range was determined by the allowable stress of arch rib critical section.Then a group of buckle cable forces were selected and substituted into the stress balance equation,and the reasonable initial buckle cable force was determined through iteration.Based on the principle of force balance,the initial anchor cable force was determined.In an engineering application example,it is shown that the stress balance and influence matrix methods for the determination of initial cable force are feasible and reliable.The initial cable forces of arch rib segments only need to be adjusted once in the corresponding construction process,which improves the working efficiency and reduces the construction risk.It is found that the methods have great advantages for determining initial cable force in cantilever casting construction process of concrete arch bridge.

Nonlinear vibration analysis of a circular micro-plate in two-sided NEMS/MEMS capacitive system by using harmonic balance method

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.

Simplified Homogeneous Balance Method and Its Applications to the Whitham-Broer-Kaup Model Equations

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.

Application of the Material Balance Method in Paleoelevation Recovery: A Case Study of the Longmen Mountains Foreland Basin on the Eastern Margin of the Tibetan Plateau

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.

Application of the Method of Characteristics to Population Balance Models Considering Growth and Nucleation Phenomena

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.

Iterative Methods for Solving the Nonlinear Balance Equation with Optimal Truncation

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.

Primary resonance of fractional-order Duffing–van der Pol oscillator by harmonic balance method

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.

INCREMENTAL HARMONIC BALANCE METHOD FOR AIRFOIL FLUTTER WITH MULTIPLE STRONG NONLINEARITIES

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.

Numerical Analysis of Balanced Methods for the Impulsive Stochastic Differential Equations

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.

Source Apportionment of Ambient PM_(10) in the Urban Area of Longyan City,China:a Comparative Study Based on Chemical Mass Balance Model and Factor Analysis Method

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 Modified Homogeneous Balance Method and Its Applications

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.

General Modified Split-Step Balanced Methods for Stiff Stochastic Differential Equations

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 ane-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.

Homogenous Balance Method and Exact Analytical Solutions for Whitham-Broer-Kaup Equations in Shallow Water

Based on the homogenous balance method and with the help of mathematica, the Backlund transformation and the transfer heat equation are derived. Analyzing the heat-transfer equation, the multiple soliton solutions and other exact analytical solution for Whitham-Broer-Kaup equations(WBK) are derived. These solutions contain Fan's, Xie's and Yan's results and other new types of analytical solutions, such as rational function solutions and periodic solutions. The method can also be applied to solve more nonlinear differential equations.