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.

Compensated split-step balanced methods for nonlinear stiff SDEs with jump-diffusion and piecewise continuous arguments

This paper deals with numerical solutions of nonlinear stiff stochastic differential equations with jump-diffusion and piecewise continuous arguments.By combining compensated split-step methods and balanced methods,a class of compensated split-step balanced(CSSB)methods are suggested for solving the equations.Based on the one-sided Lipschitz condition and local Lipschitz condition,a strong convergence criterion of CSSB methods is derived.It is proved under some suitable conditions that the numerical solutions produced by CSSB methods can preserve the mean-square exponential stability of the corresponding analytical solutions.Several numerical examples are presented to illustrate the obtained theoretical results and the effectiveness of CSSB methods.Moreover,in order to show the computational advantage of CSSB methods,we also give a numerical comparison with the adapted split-step backward Euler methods with or without compensation and tamed explicit methods.

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.

A vertical track nonlinear energy sink

Eliminating the effects of gravity and designing nonlinear energy sinks(NESs)that suppress vibration in the vertical direction is a challenging task with numerous damping requirements.In this paper,the dynamic design of a vertical track nonlinear energy sink(VTNES)with zero linear stiffness in the vertical direction is proposed and realized for the first time.The motion differential equations of the VTNES coupled with a linear oscillator(LO)are established.With the strong nonlinearity considered of the VTNES,the steady-state response of the system is analyzed with the harmonic balance method(HBM),and the accuracy of the HBM is verified numerically.On this basis,the VTNES prototype is manufactured,and its nonlinear stiffness is identified.The damping effect and dynamic characteristics of the VTNES are studied theoretically and experimentally.The results show that the VTNES has better damping effects when strong modulation responses(SMRs)occur.Moreover,even for small-amplitude vibration,the VTNES also has a good vibration suppression effect.To sum up,in order to suppress the vertical vibration,an NES is designed and developed,which can suppress the vertical vibration within certain ranges of the resonance frequency and the vibration intensity.

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.

时域谐波平衡求解器中非物理解出现根源探析

The time domain harmonic balance method is an attractive reduced order method of analyzing unsteady flow for turbomachines. However, the method can admit non-physical solutions. Non-physical solutions were encountered from a three-blade-row compressor configuration in a time domain harmonic balance analysis. This paper aims to investigate the root cause of the non-physical solutions. The investigation involves several strategies, which include increasing the number of harmonics, increasing the number of time instants, including scattered modes,including the rotor-rotor interaction, and the use of a new method-the approximate time domain nonlinear harmonic method. Numerical analyses pertinent to each strategy are presented to reveal the root cause of the non-physical solution. It is found that the nonlinear interaction of unsteady flow components with different fundamental frequencies is the cause of the non-physical solution. The non-physical solution can be eliminated by incorporating extra scattered modes or using the approximate time domain nonlinear harmonic method.

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.

A new auxiliary equation method for finding travelling wave solutions to KdV equation

In this paper, a new auxiliary equation method is used to find exact travelling wave solutions to the （1＋1）-dimensional KdV equation. Some exact travelling wave solu- tions with parameters have been obtained, which cover the existing solutions. Compared to other methods, the presented method is more direct, more concise, more effective, and easier for calculations. In addition, it can be used to solve other nonlinear evolution equations in mathematical physics.

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.

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.

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.

Space decomposition based parallelization solutions for the combined finiteediscrete element method in 2D

The combined finiteediscrete element method （FDEM） belongs to a family of methods of computationalmechanics of discontinua. The method is suitable for problems of discontinua, where particles aredeformable and can fracture or fragment. The applications of FDEM have spread over a number of disciplinesincluding rock mechanics, where problems like mining, mineral processing or rock blasting canbe solved by employing FDEM. In this work, a novel approach for the parallelization of two-dimensional（2D） FDEM aiming at clusters and desktop computers is developed. Dynamic domain decompositionbased parallelization solvers covering all aspects of FDEM have been developed. These have beenimplemented into the open source Y2D software package and have been tested on a PC cluster. Theoverall performance and scalability of the parallel code have been studied using numerical examples. Theresults obtained confirm the suitability of the parallel implementation for solving large scale problems. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

The (G'/G, 1/G)-expansion method and its application to travelling wave solutions of the Zakharov equations

The （G＇/G, 1/G）-expansion method for finding exact travelling wave solutions of nonlinear evolution equations, which can be thought of as an extension of the （G＇/G）-expansion method proposed recently, is presented. By using this method abundant travelling wave so- lutions with arbitrary parameters of the Zakharov equations are successfully obtained. When the parameters are replaced by special values, the well-known solitary wave solutions of the equations are rediscovered from the travelling waves.

The Spatial Sensitivity Analysis of Evapotranspiration using Penman-Monteith Method at Grid Scale

The need to allocate the existing water in a sustainable manner, even with the projected population growth, has made to assess the consumptive use or evapotranspiration (ET), which determines the irrigation demand. As underscored in the literature, Penman-Monteith method which is a combination of aerodynamic and energy balance method is widely used and accepted as the method of estimation of ET. However, the application of Penman-Monteith relies on many climate parameters such as relative humidity, solar radiation, temperature, and wind speed. Therefore, there exists a need to determine the parameters that are most sensitive and correlated with dependent variable (i.e., ET), to strengthen the knowledge base. However, the sensitivity of ET using Penman-Monteith is oftentimes estimated using meteorological data from climate stations. Such estimation of sensitivity may vary spatially and thus there exists a need to estimate sensitivity of ET spatially. Thus, in this paper, based on One-AT-A-Time (OAT) method, a spatial sensitivity tool that can geographically encompass all the best available climate datasets to produce ET and its sensitivity at different spatial scales is developed. The spatial tool is developed as a Python toolbox in ArcGIS using Python, an open source programming language, and the ArcPy site-package of ArcGIS. The developed spatial tool is demonstrated using the meteorological data from Automated Weather Data Network in Nebraska in 2010. To summarize the outcome of the sensitivity analysis using OAT method, sensitivity indices are developed for each raster cell. The demonstration of the tool shows that, among the considered parameters, the computed ET using Penman-Monteith is highly sensitive to solar radiation followed by temperature for the state of Nebraska, as depicted by the sensitivity index. The computed sensitivity index of wind speed and the relative humidity are not that significant compared to the sensitivity index of solar radiation and temperature.

A Simplified Hirota Method and Its Application

By means of the homogeneous balance principle, a nonlinear transformation to the well known breaking soliton equation with physical interest was given. The original equation was turned into a homogeneity differential equation with this nonlinear transformation. By solving the homogeneity equation via the simplified Hirota method and applying the nonlinear transformation, one soliton, two soliton and three soliton solutions as well as some other types of explicit solutions to the breaking soliton equation were obtained with the assistance of Maple.

Auto-Backlund Transformation and Extended Tanh-Function Methods to Solve the Time-Dependent Coefficients Calogero-Degasperis Equation

In this paper, the Auto-B?cklund transformation connected with the homogeneous balance method (HB) and the extended tanh-function method are used to construct new exact solutions for the time-dependent coefficients Calogero-Degasperis (VCCD) equation. New soliton and periodic solutions of many types are obtained. Furthermore, the soliton propagation is discussed under the effect of the variable coefficients.

Data correction of the force balanced accelerometer

ＩｎｔｒｏｄｕｃｔｉｏｎＤｉｇｉｔａｌｓｔｒｏｎｇｇｒｏｕｎｄａｃｅｌｅｒａｔｉｏｎｏｂｓｅｒｖａｔｉｏｎｉｎｓｔｒｕｍｅｎｔ，ｓｕｃｈａｓＰＤＲ１，ＳＳＡ１ａｎｄＳＳＲ１ｐｒｏｄｕｃｅｄｂｙＫｉｎｅｍｅｔｒｉｃｓＩｎｃ．，ＵＳＡａｎｄＳＣＱ?..

Chief Physician Defen Wang’s Clinical Experience in Dialectical Treatment with Mediation Method of Perimenopausal Sweating

The objective of this paper is to outline Chief Physician Defen Wang’s clinical experience in perimenopausal sweating syndrome differentiation and treatment,to inherit Mr.Gao’s academic concept of mediation technique,and to integrate multiple prescriptions under Mr.Gao’s thinking.The deficiency is based on the deficiency of liver and kidney Yin,the weakness of spleen and stomach,and the standard is the imbalance of Ying and Wei,the imbalance of Qi,blood and body fluid,and the main method of treatment is nourishing liver and kidney,soothing liver and invigorating spleen,regulating Ying and Wei,as well as tonifying Qi and absorbing fluid.

Adjacent mode resonance of a hydraulic pipe system consisting of parallel pipes coupled at middle points

The coupling vibration of a hydraulic pipe system consisting of two pipes is studied.The pipes are installed in parallel and fixed at their ends,and are restrained by clips to one bracket at their middle points.The pipe subjected to the basement excitation at the left end is named as the active pipe,while the pipe without excitation is called the passive pipe.The clips between the two pipes are the bridge for the vibration energy.The adjacent natural frequencies will enhance the vibration coupling.The governing equation of the coupled system is deduced by the generalized Hamilton principle,and is discretized to the modal space.The modal correction is used during the discretization.The investigation on the natural characters indicates that the adjacent natural frequencies can be adjusted by the stiffness of the two clips and bracket.The harmonic balance method(HBM)is used to study the responses in the adjacent natural frequency region.The results show that the vibration energy transmits from the active pipe to the passive pipe swimmingly via the clips together with a flexible bracket,while the locations of them are not node points.The adjacent natural frequencies may arouse wide resonance curves with two peaks for both pipes.The stiffness of the clip and bracket can release the vibration coupling.It is suggested that the stiffness of the clip on the passive pipe should be weak and the bracket should be strong enough.In this way,the vibration energy is reflected by the almost rigid bracket,and is hard to transfer to the passive pipe via a soft clip.The best choice is to set the clips at the pipe node points.The current work gives some suggestions for weakening the coupled vibration during the dynamic design of a coupled hydraulic pipe system.