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.

Adaptive hp finite element method for fluorescence molecular tomography with simplified spherical harmonics approximation

Recently,the simplified spherical harmonics equations(SP)model has at tracted much att entionin modeling the light propagation in small tissue ggeometriesat visible and near-infrared wave-leng ths.In this paper,we report an eficient numerical method for fluorescence moleeular tom-ography(FMT)that combines the advantage of SP model and adaptive hp finite elementmethod(hp-FEM).For purposes of comparison,hp-FEM and h-FEM are,respectively applied tothe reconstruction pro cess with diffusion approximation and SPs model.Simulation experiments on a 3D digital mouse atlas and physical experiments on a phantom are designed to evaluate thereconstruction methods in terms of the location and the reconstructed fluorescent yield.Theexperimental results demonstrate that hp-FEM with SPy model,yield more accurate results thanh-FEM with difusion approximation model does.The phantom experiments show the potentialand feasibility of the proposed approach in FMT applications.

Spherical harmonics method for neutron transport equation based on unstructured-meshes

Based on a new second-order neutron transport equation, self-adjoint angular flux (SAAF) equation, the spherical harmonics (PN) method for neutron transport equation on unstructured-meshes is derived. The spherical harmonics function is used to expand the angular flux. A set of differential equations about the spatial variable, which are coupled with each other, can be obtained. They are solved iteratively by using the finite element method on un- structured-meshes. A two-dimension transport calculation program is coded according to the model. The numerical results of some benchmark problems demonstrate that this method can give high precision results and avoid the ray effect very well.

Utility and Application of a Versatile Analytical Method for MMF Calculation in AC Machines

A versatile analytical method(VAM) for calculating the harmonic components of the magnetomotive force(MMF) generated by diverse armature windings in AC machines has been proposed, and the versatility of this method has been established in early literature. However, its practical applications and significance in advancing the analysis of AC machines need further elaboration. This paper aims to complement VAM by augmenting its theory, offering additional insights into its conclusions, as well as demonstrating its utility in assessing armature windings and its application of calculating torque for permanent magnet synchronous machines(PMSM). This work contributes to advancing the analysis of AC machines and underscores the potential for improved design and performance optimization.

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.

Harmonics analysis of the ITER poloidal field converter based on a piecewise method

Poloidal field(PF) converters provide controlled DC voltage and current to PF coils. The many harmonics generated by the PF converter flow into the power grid and seriously affect power systems and electric equipment. Due to the complexity of the system, the traditional integral operation in Fourier analysis is complicated and inaccurate. This paper presents a piecewise method to calculate the harmonics of the ITER PF converter. The relationship between the grid input current and the DC output current of the ITER PF converter is deduced. The grid current is decomposed into the sum of some simple functions. By calculating simple function harmonics based on the piecewise method, the harmonics of the PF converter under different operation modes are obtained.In order to examine the validity of the method, a simulation model is established based on Matlab/Simulink and a relevant experiment is implemented in the ITER PF integration test platform.Comparative results are given. The calculated results are found to be consistent with simulation and experiment. The piecewise method is proved correct and valid for calculating the system harmonics.

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

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.

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.

Wave reflection, transmission and harmonics due to different submerged structures

Comparisons of wave reflection, transmission and harmonics due to different types of sub merged structures are investigated by a numerical method, the boundary-fitted coordinate (BFC) method. The types of submerged structures include a submerged horizontal plate, submerged breakwa ters (rectangular and trapezoidal) and a step-type structure (topography). First, the BFC method is ver ified by comparing the computed results with the experimental data, including wave surface elevations, reflected and transmitted wave heights, and amplitudes of higher harmonics, showing that the method is a reasonable one to predict wave deformations due to the submerged structures. Secondly, the wave sur face elevations and the higher harmonics over different submerged structures are compared. Thirdly, re flected and transmitted waves due to different submerged structures are investigated.

Experimental and Numerical Study of Structural Identification Using Non-Linear Resonant Decay Method

One of the practical approaches in identifying structures is the non-linear resonant decay method which identifies a non-linear dynamic system utilizing a model based on linear modal space containing the underlying linear system and a small number of extra terms that exhibit the non-linear effects.In this paper,the method is illustrated in a simulated system and an experimental structure.The main objective of the non-linear resonant decay method is to identify the non-linear dynamic systems based on the use of a multi-shaker excitation using appropriated excitation which is obtained from the force appropriation approach.The experimental application of the method is indicated to provide suitable estimates of modal parameters for the identification of non-linear models of structures.

Incremental harmonic balance method for periodic forced oscillation of a dielectric elastomer balloon

Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices call for an efficient and reliable method to analyze the dynamic response of DE. This remains to be a challenge since the resultant vibration equation of DE, for example, the vibration of a DE balloon considered here is highly nonlinear with higher-order power terms and time-dependent coefficients. Previous efforts toward this goal use largely the numerical integration method with the simple harmonic balance method as a supplement. The numerical integration and the simple harmonic balance method are inefficient for large parametric analysis or with difficulty in improving the solution accuracy. To overcome the weakness of these two methods,we describe formulations of the incremental harmonic balance(IHB) method for periodic forced solutions of such a unique system. Combined with an arc-length continuation technique, the proposed strategy can capture the whole solution branches, both stable and unstable, automatically with any desired accuracy.

Improving convergence of incremental harmonic balance method using homotopy analysis method

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.

Quantitative Measurement of the Proportions of High-Order Harmonics for the 4B7B Soft-X-Ray Source at Beijing Synchrotron Radiation Facility

A transmission grating coupled with an X-ray charge coupled device(CCD)is used to quantitatively measure the proportion of high-order harmonics of the soft-X-ray source of beam line 4B7B.The results show that the monochromatic X-ray has third-order and second-order harmonics.The proportion of second-order harmonic of 4B7B is less than 9.0%and the thirdorder harmonic is below 0.7%when no suppressing method is applied.When suppression methods are used,the proportion of second-order harmonic is less than 1.7%and the third-order harmonic is ignorable.

Spectral Method for Three-Dimensional Nonlinear Klein-Gordon Equation by Using Generalized Laguerre and Spherical Harmonic Functions

In this paper,a generalized Laguerre-spherical harmonic spectral method is proposed for the Cauchy problem of three-dimensional nonlinear Klein-Gordon equation. The goal is to make the numerical solutions to preserve the same conservation as that for the exact solution.The stability and convergence of the proposed scheme are proved.Numerical results demonstrate the efficiency of this approach.We also establish some basic results on the generalized Laguerre-spherical harmonic orthogonal approximation,which play an important role in spectral methods for various problems defined on the whole space and unbounded domains with spherical geometry.

Computer Methodologies for the Comparison of Some Efficient Derivative FreeSimultaneous Iterative Methods for Finding Roots of Non-Linear Equations

In this article,we construct the most powerful family of simultaneous iterative method with global convergence behavior among all the existing methods in literature for finding all roots of non-linear equations.Convergence analysis proved that the order of convergence of the family of derivative free simultaneous iterative method is nine.Our main aim is to check out the most regularly used simultaneous iterative methods for finding all roots of non-linear equations by studying their dynamical planes,numerical experiments and CPU time-methodology.Dynamical planes of iterative methods are drawn by using MATLAB for the comparison of global convergence properties of simultaneous iterative methods.Convergence behavior of the higher order simultaneous iterative methods are also illustrated by residual graph obtained from some numerical test examples.Numerical test examples,dynamical behavior and computational efficiency are provided to present the performance and dominant efficiency of the newly constructed derivative free family of simultaneous iterative method over existing higher order simultaneous methods in literature.

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.

Fibonacci Harmonics: A New Mathematical Model of Synchronicity

This article aims to provide a brief overview of the relevance of new findings about the Fibonacci Life Chart Method (FLCM) for understanding synchronicity. The FLCM is reviewed first, including an exposition of the golden section model, and elaboration of a new harmonic model. The two models are then compared to illuminate several strengths and weaknesses in connection with the following four major criteria regarding synchronicity: explanatory adequacy;predictability of future synchronicities;simplicity of the model;and generalizability to other branches of knowledge. The review indicates that both models appear capable of simulating nonlinear and fractal dynamics. Hybrid approaches that combine both models are feasible and necessary for projects that aim to experimentally address synchronicity.

A SCR method for uncertainty estimation in geodesy non-linear error propagation: Comparisons and applications

We review three derivative-free methods developed for uncertainty estimation of non-linear error propagation, namely, MC(Monte Carlo), SUT(scaled unscented transformation), and SI(sterling interpolation). In order to avoid preset parameters like as these three methods need, we introduce a new method to uncertainty estimation for the first time, namely, SCR(spherical cubature rule), which is no need for setting parameters. By theoretical derivation, we prove that the precision of uncertainty obtained by SCR can reach second-order. We conduct four synthetic experiments, for the first two experiments, the results obtained by SCR are consistent with the other three methods with optimal setting parameters, but SCR is easier to operate than other three methods, which verifies the superiority of SCR in calculating the uncertainty. For the third experiment, real-time calculation is required, so the MC is hardly feasible. For the forth experiment, the SCR is applied to the inversion of seismic fault parameter which is a common problem in geophysics, and we study the sensitivity of surface displacements to fault parameters with errors. Our results show that the uncertainty of the surface displacements is the magnitude of ±10 mm when the fault length contains a variance of 0.01 km^(2).