Random Green's Function Method for Large-Scale Electronic Structure Calculation

We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.

Dynamic Characteristics of Functionally Graded Timoshenko Beams by Improved Differential Quadrature Method

This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node distribution.Firstly,based on the first-order shear deformation theory,the governing equation of free vibration of a functionally graded beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam axial displacement,transverse displacement,and cross-sectional rotation angle by considering the effects of shear deformation and rotational inertia of the beam cross-section.Then,ignoring the shear deformation of the beam section and only considering the effect of the rotational inertia of the section,the governing equation of the beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam transverse displacement.Based on the differential quadrature method theory,the eigenvalue problem of ordinary differential equations is transformed into the eigenvalue problem of standard generalized algebraic equations.Finally,the first several natural frequencies of the beam can be calculated.The feasibility and accuracy of the improved DQM are verified using the finite element method(FEM)and combined with the results of relevant literature.

A Hybrid Level Set Optimization Design Method of Functionally Graded Cellular Structures Considering Connectivity

With the continuous advancement in topology optimization and additive manufacturing(AM)technology,the capability to fabricate functionally graded materials and intricate cellular structures with spatially varying microstructures has grown significantly.However,a critical challenge is encountered in the design of these structures–the absence of robust interface connections between adjacent microstructures,potentially resulting in diminished efficiency or macroscopic failure.A Hybrid Level Set Method(HLSM)is proposed,specifically designed to enhance connectivity among non-uniform microstructures,contributing to the design of functionally graded cellular structures.The HLSM introduces a pioneering algorithm for effectively blending heterogeneous microstructure interfaces.Initially,an interpolation algorithm is presented to construct transition microstructures seamlessly connected on both sides.Subsequently,the algorithm enables the morphing of non-uniform unit cells to seamlessly adapt to interconnected adjacent microstructures.The method,seamlessly integrated into a multi-scale topology optimization framework using the level set method,exhibits its efficacy through numerical examples,showcasing its prowess in optimizing 2D and 3D functionally graded materials(FGM)and multi-scale topology optimization.In essence,the pressing issue of interface connections in complex structure design is not only addressed but also a robust methodology is introduced,substantiated by numerical evidence,advancing optimization capabilities in the realm of functionally graded materials and cellular structures.

Transfer matrix method for free and forced vibrations of multi-level functionally graded material stepped beams with different boundary conditions

Functionally graded materials(FGMs)are a novel class of composite materials that have attracted significant attention in the field of engineering due to their unique mechanical properties.This study aims to explore the dynamic behaviors of an FGM stepped beam with different boundary conditions based on an efficient solving method.Under the assumptions of the Euler-Bernoulli beam theory,the governing differential equations of an individual FGM beam are derived with Hamilton’s principle and decoupled via the separation-of-variable approach.Then,the free and forced vibrations of the FGM stepped beam are solved with the transfer matrix method(TMM).Two models,i.e.,a three-level FGM stepped beam and a five-level FGM stepped beam,are considered,and their natural frequencies and mode shapes are presented.To demonstrate the validity of the method in this paper,the simulation results by ABAQUS are also given.On this basis,the detailed parametric analyses on the frequencies and dynamic responses of the three-level FGM stepped beam are carried out.The results show the accuracy and efficiency of the TMM.

A Radial Basis Function Method with Improved Accuracy for Fourth Order Boundary Value Problems

Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.

Three-dimensional color particle image velocimetry based on a cross-correlation and optical flow method

Rainbow particle image velocimetry(PIV)can restore the three-dimensional velocity field of particles with a single camera;however,it requires a relatively long time to complete the reconstruction.This paper proposes a hybrid algorithm that combines the fast Fourier transform(FFT)based co-correlation algorithm and the Horn–Schunck(HS)optical flow pyramid iterative algorithm to increase the reconstruction speed.The Rankine vortex simulation experiment was performed,in which the particle velocity field was reconstructed using the proposed algorithm and the rainbow PIV method.The average endpoint error and average angular error of the proposed algorithm were roughly the same as those of the rainbow PIV algorithm;nevertheless,the reconstruction time was 20%shorter.Furthermore,the effect of velocity magnitude and particle density on the reconstruction results was analyzed.In the end,the performance of the proposed algorithm was verified using real experimental single-vortex and double-vortex datasets,from which a similar particle velocity field was obtained compared with the rainbow PIV algorithm.The results show that the reconstruction speed of the proposed hybrid algorithm is approximately 25%faster than that of the rainbow PIV algorithm.

The influences of large earthquake signals on the recovery of surface waves from ambient noise cross-correlation functions

Ambient noise tomography(ANT)has been widely used to image crust and upmost mantle structures.ANT assumes that sources of ambient noise are diffuse and evenly distributed in space and the energy of different modes is equipartitioned.At present,the sources of the primary and the secondary microseisms are well studied,but there are only a few on the studies of long-period ambient noise sources.In this study,we study the effects of large earthquake signals on the recovery of surface waves from seismic ambient noise data recorded by seismic stations from the US permanent networks and Global Seismographic Network(GSN).Our results show that large earthquake signals play an important role on the recovery of long-period surface waves from ambient noise cross-correlation functions.Our results are consistent with previous studies that suggest the contribution of earthquake signals to the recovery of surface waves from cross-correlations of ambient noise is dominant at periods larger than 20–40 s.

Intensity Correlation Function of a Single-Mode Laser Driven by Colored Cross-Correlation Noises in Modulation of a Bias Signal

By using the linear approximation method, the intensity correlation function is calculated for a single-mode laser modulated by a bias signal and driven by colored pump and quantum noises with colored cross-correlation. We found that, when the correlation time between the two noises is very short, the behavior of the intensity correlation function versus the time, in addition to decreasing monotonously, also exhibits several cases, such as one maximum, one minimum, and two extrema. When the correlation time between the two noises is very long, the behavior of the intensity correlation function exhibits oscillation and the envelope is similar to the case of short cross-correlation time.

Study on Characteristics of 3-D Translating-Pulsating Source Green Function of Deep-Water Havelock Form and Its Fast Integration Method

The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.

SOME PROBLEMS WITH THE METHOD OF FUNDAMENTAL SOLUTION USING RADIAL BASIS FUNCTIONS

The present work describes the application of the method of fundamental solutions （MFS） along with the analog equation method （AEM） and radial basis function （RBF） approximation for solving the 2D isotropic and anisotropic Helmholtz problems with different wave numbers. The AEM is used to convert the original governing equation into the classical Poisson＇s equation, and the MFS and RBF approximations are used to derive the homogeneous and particular solutions, respectively. Finally, the satisfaction of the solution consisting of the homogeneous and particular parts to the related governing equation and boundary conditions can produce a system of linear equations, which can be solved with the singular value decomposition （SVD） technique. In the computation, such crucial factors related to the MFS-RBF as the location of the virtual boundary, the differential and integrating strategies, and the variation of shape parameters in multi-quadric （MQ） are fully analyzed to provide useful reference.

Function Combined Method for Design Innovation of Children's Bike

As children mature, bike products for children in the market develop at the same time, and the conditions are frequently updated. Certain problems occur when using a bike, such as cycle overlapping, repeating function, and short life cycle, which go against the principles of energy conservation and the environmental protection intensive design concept. In this paper, a rational multi-function method of design through functional superposition, transformation, and technical implementation is proposed. An organic combination of frog-style scooter and children’s tricycle is developed using the multi-function method. From the ergonomic perspective, the paper elaborates on the body size of children aged 5 to 12 and effectively extracts data for a multi-function children’s bike, which can be used for gliding and riding. By inverting the body, parts can be interchanged between the handles and the pedals of the bike. Finally, the paper provides a detailed analysis of the components and structural design, body material, and processing technology of the bike. The study of Industrial Product Innovation Design provides an effective design method to solve the bicycle problems, extends the function problems, improves the product market situation, and enhances the energy saving feature while implementing intensive product development effectively at the same time.

Reconstruction of sub cross-correlation cancellation technique for unambiguous acquisition of BOC(kn, n) signals

In order to solve the problem of ambiguous acquisition of BOC signals caused by its property of multiple peaks,an unambiguous acquisition algorithm named reconstruction of sub cross-correlation cancellation technique(RSCCT)for BOC(kn,n)signals is proposed.In this paper,the principle of signal decomposition is combined with the traditional acquisition algorithm structure,and then based on the method of reconstructing the correlation function.The method firstly gets the sub-pseudorandom noise(PRN)code by decomposing the local PRN code,then uses BOC(kn,n)and the sub-PRN code cross-correlation to get the sub cross-correlation function.Finally,the correlation peak with a single peak is obtained by reconstructing the sub cross-correlation function so that the ambiguities of BOC acquisition are removed.The simulation shows that RSCCT can completely eliminate the side peaks of BOC(kn,n)group signals while maintaining the narrow correlation of BOC,and its computational complexity is equivalent to sub carrier phase cancellation(SCPC)and autocorrelation side-peak cancellation technique(ASPeCT),and it reduces the computational complexity relative to BPSK-like.For BOC(n,n),the acquisition sensitivity of RSCCT is 3.25 dB,0.81 dB and 0.25 dB higher than binary phase shift keying(BPSK)-like,SCPC and ASPeCT at the acquisition probability of 90%,respectively.The peak to average power ratio is 1.91,3.0 and 3.7 times higher than ASPeCT,SCPC and BPSK-like at SNR=–20 dB,respectively.For BOC(2n,n),the acquisition sensitivity of RSCCT is 5.5 dB,1.25 dB and 2.69 dB higher than BPSK-like,SCPC and ASPeCT at the acquisition probability of 90%,respectively.The peak to average power ratio is 1.02,1.68 and 2.12 times higher than ASPeCT,SCPC and BPSK-like at SNR=–20 dB,respectively.

Construction of doubly-periodic solutions to nonlinear partial differential equations using improved Jacobi elliptic function expansion method and symbolic computation

Some doubly-periodic solutions of the Zakharov-Kuznetsov equation are presented. Our approach is to introduce an auxiliary ordinary differential equation and use its Jacobi elliptic function solutions to construct doubly-periodic solutions of the Zakharov-Kuznetsov equation, which has been derived by Gottwald as a two-dimensional model for nonlinear Rossby waves. When the modulus k →1, these solutions reduce to the solitary wave solutions of the equation.

Free vibration of non-uniform axially functionally graded beams using the asymptotic development method

The asymptotic development method is applied to analyze the free vibration of non-uniform axially functionally graded(AFG) beams, of which the governing equations are differential equations with variable coefficients. By decomposing the variable flexural stiffness and mass per unit length into reference invariant and variant parts, the perturbation theory is introduced to obtain an approximate analytical formula of the natural frequencies of the non-uniform AFG beams with different boundary conditions.Furthermore, assuming polynomial distributions of Young's modulus and the mass density, the numerical results of the AFG beams with various taper ratios are obtained and compared with the published literature results. The discussion results illustrate that the proposed method yields an effective estimate of the first three order natural frequencies for the AFG tapered beams. However, the errors increase with the increase in the mode orders especially for the cases with variable heights. In brief, the asymptotic development method is verified to be simple and efficient to analytically study the free vibration of non-uniform AFG beams, and it could be used to analyze any tapered beams with an arbitrary varying cross width.

Temporal-spatial cross-correlation analysis of non-stationary near-surface wind speed time series

Temporal-spatial cross-correlation analysis of non-stationary wind speed time series plays a crucial role in wind field reconstruction as well as in wind pattern recognition.Firstly,the near-surface wind speed time series recorded at different locations are studied using the detrended fluctuation analysis(DFA),and the corresponding scaling exponents are larger than 1.This indicates that all these wind speed time series have non-stationary characteristics.Secondly,concerning this special feature( i.e.,non-stationarity)of wind signals,a cross-correlation analysis method,namely detrended cross-correlation analysis(DCCA) coefficient,is employed to evaluate the temporal-spatial cross-correlations between non-stationary time series of different anemometer pairs.Finally,experiments on ten wind speed data synchronously collected by the ten anemometers with equidistant arrangement illustrate that the method of DCCA cross-correlation coefficient can accurately analyze full-scale temporal-spatial cross-correlation between non-stationary time series and also can easily identify the seasonal component,while three traditional cross-correlation techniques(i.e.,Pearson coefficient,cross-correlation function,and DCCA method) cannot give us these information directly.

Detrended cross-correlation analysis of electroencephalogram

In the paper we use detrended cross-correlation analysis （DCCA） to study the electroencephalograms of healthy young subjects and healthy old subjects. It is found that the cross-correlation between different leads of a healthy young subject is larger than that of a healthy old subject. It was shown that the cross-correlation relationship decreases with the aging process and the phenomenon can help to diagnose whether the subject＇s brain function is healthy or not.

A STUDY ON THE WEIGHT FUNCTION OF THE MOVING LEAST SQUARE APPROXIMATION IN THE LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.

An improved influence function method for predicting subsidence caused by longwall mining operations in inclined coal seams

Prediction of surface subsidence caused by longwall mining operation in inclined coal seams is often very challenging. The existing empirical prediction methods are inflexible for varying geological and mining conditions. An improved influence function method has been developed to take the advantage of its fundamentally sound nature and flexibility. In developing this method, the original Knothe function has been transformed to produce a continuous and asymmetrical subsidence influence function. The empirical equations for final subsidence parameters derived from col- lected longwall subsidence data have been incorporated into the mathematical models to improve the prediction accuracy. A number of demonstration cases for longwall mining operations in coal seams with varying inclination angles, depths and panel widths have been used to verify the applicability of the new subsidence prediction model.

Uncertainty analysis of strain modal parameters by Bayesian method using frequency response function

Structural strain modes are able to detect changes in local structural performance, but errors are inevitably intermixed in the measured data. In this paper, strain modal parameters are considered as random variables, and their uncertainty is analyzed by a Bayesian method based on the structural frequency response function （FRF）. The estimates of strain modal parameters with maximal posterior probability are determined. Several independent measurements of the FRF of a four-story reinforced concrete flame structural model were performed in the laboratory. The ability to identify the stiffness change in a concrete column using the strain mode was verified. It is shown that the uncertainty of the natural frequency is very small. Compared with the displacement mode shape, the variations of strain mode shapes at each point are quite different. The damping ratios are more affected by the types of test systems. Except for the case where a high order strain mode does not identify local damage, the first order strain mode can provide an exact indication of the damage location.

An influence function method based subsidence prediction program for longwall mining operations in inclined coal seams

The distribution of the final surface subsidence basin induced by longwall operations in inclined coal seam could be significantly different from that in flat coal seam and demands special prediction methods. Though many empirical prediction methods have been developed, these methods are inflexible for varying geological and mining conditions. An influence function method has been developed to take the advantage of its fundamentally sound nature and flexibility. In developing this method, significant modifications have been made to the original Knothe function to produce an asymmetrical influence function. The empirical equations for final subsidence parameters derived from US subsidence data and Chinese empirical values have been incorpo- rated into the mathematical models to improve the prediction accuracy. A corresponding computer program is developed. A number of subsidence cases for longwall mining operations in coal seams with varying inclination angles have been used to demonstrate the applicability of the developed subsidence prediction model.