Full-vectorial finite-difference beam propagation method based on the modified alternating direction implicit method

A modified alternating direction implicit algorithm is proposed to solve the full-vectorial finite-difference beam propagation method formulation based on H fields. The cross-coupling terms are neglected in the first sub-step, but evaluated and doubly used in the second sub-step. The order of two sub-steps is reversed for each transverse magnetic field component so that the cross-coupling terms are always expressed in implicit form, thus the calculation is very efficient and stable. Moreover, an improved six-point finite-difference scheme with high accuracy independent of specific structures of waveguide is also constructed to approximate the cross-coupling terms along the transverse directions. The imaginary-distance procedure is used to assess the validity and utility of the present method. The field patterns and the normalized propagation constants of the fundamental mode for a buried rectangular waveguide and a rib waveguide are presented. Solutions are in excellent agreement with the benchmark results from the modal transverse resonance method.

The efficiency of direct integration methods in elastic contact-impact problems

A comparison of direct integration methods is madeand their efficiency is investigated for impact problems.New-mark,Wilson-θ,Central Difference and Houbolt Methodsare used as direct integration methods.Impact analysisincludes that of elastic and large deformation based uponupdated Lagrangian including buckling check.The resultsshow that the direct integration methods give differentresults in different contact-impact cases.

Structural Reliability Analysis Based on Support Vector Machine and Dual Neural Network Direct Integration Method

Aiming at the reliability analysis of small sample data or implicit structural function,a novel structural reliability analysis model based on support vector machine(SVM)and neural network direct integration method(DNN)is proposed.Firstly,SVM with good small sample learning ability is used to train small sample data,fit structural performance functions and establish regular integration regions.Secondly,DNN is approximated the integral function to achieve multiple integration in the integration region.Finally,structural reliability was obtained by DNN.Numerical examples are investigated to demonstrate the effectiveness of the present method,which provides a feasible way for the structural reliability analysis.

3D elastic wave equation forward modeling based on the precise integration method

The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.

OPTIMIZATION OF ADAPTIVE DIRECT METHOD FOR APPROXIMATE SOLUTION OF INTEGRAL EQUATIONS OF SEVERAL VARIABLES

This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic Sobolev classes, and also gives an optimal algorithm.

Synthesis of Hierarchically Porous FAU/γ-Al2O3 Composites with Different Morphologies via Directing Agent Induced Method

Zeolite FAU composites with a macro/meso-microporous hierarchical structure were hydrothermally synthesized using macro-mesoporous γ-Al_2O_3 monolith as the substrate by means of the liquid crystallization directing agent(LCDA) induced method. No template was needed throughout the synthesis processes. The structure and porosity of zeolite composites were analyzed by means of X-ray powder diffraction(XRD), scanning electron microscopy(SEM) and N_2adsorption-desorption isotherms. The results showed that the supported zeolite composites with varied zeolitic crystalline phases and different morphologies can be obtained by adjusting the crystallization parameters, such as the crystallization temperature, the composition and the alkalinity of the precursor solution. The presence of LCDA was defined as a determinant for synthesizing the zeolite composites. The mechanisms for formation of the hierarchically porous FAU zeolite composites in the LCDA induced synthesis process were discussed. The resulting monolithic zeolite with a trimodal-porous hierarchical structure shows potential applicability where facile diffusion is required.

Direct method of finding first integral of two-dimensional autonomous systems in polar coordinates

A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous Hamiltonian system is multiplied by a set of multiplicative functions, then the general expression of the first integral can be obtained, An example is given to illustrate the application of the results.

Compact implicit integration factor methods for some complex-valued nonlinear equations

The compact implicit integration factor （cIIF） method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger （NLS） equation and the complex Ginzburg-Landau （GL） equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.

RETHINKING TO FINITE DIFFERENCE TIME-STEP INTEGRATIONS

The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integration schemes derived by the singlepoint precise integration method given in this paper are proved unconditionally stable.Comparisons between the schemes derived by the finite difference method and theschemes by the method employed in the present paper are made for diffusion andconvective-diffusion equations. Nunierical examples show the superiority of the singlepoint integration method.

ON DIRECT METHOD OF SOLUTION FOR A CLASS OF SINGULAR INTEGRAL EQUATIONS

In this article, by introducing characteristic singular integral operator and associate singular integral equations （SIEs）, the authors discuss the direct method of solution for a class of singular integral equations with certain analytic inputs. They obtain both the conditions of solvability and the solutions in closed form. It is noteworthy that the method is different from the classical one that is due to Lu.

High accuracy compact finite difference methods and their applications

Numerical simulation of complex flow fields with multi-scale structures is one of the most important and challenging branches of computational fluid dynamics. From linear analysis and numerical experiments it has been discovered that the higher-order accurate method can give reliable and efficient computational results, as well as better resolution of the complex flow fields with multi-scale structures. Compact finite difference schemes, which feature higher-order accuracy and spectral-like resolution with smaller stencils and easier application of boundary conditions, has attracted more and more interest and attention.

A Straightforward Direct Traction Boundary Integral Method for Two-Dimensional Crack Problems Simulation of Linear Elastic Materials

This paper presents a direct traction boundary integral equation method(DTBIEM)for two-dimensional crack problems of materials.The traction boundary integral equation was collocated on both the external boundary and either side of the crack surfaces.The displacements and tractions were used as unknowns on the external boundary,while the relative crack opening displacement(RCOD)was chosen as unknowns on either side of crack surfaces to keep the single-domain merit.Only one side of the crack surfaces was concerned and needed to be discretized,thus the proposed method resulted in a smaller system of algebraic equations compared with the dual boundary element method(DBEM).A new set of crack-tip shape functions was constructed to represent the strain field singularity exactly,and the SIFs were evaluated by the extrapolation of the RCOD.Numerical examples for both straight and curved cracks are given to validate the accuracy and efficiency of the presented method.

Solving Stiff Reaction-Diffusion Equations Using Exponential Time Differences Methods

Reaction-diffusion equations modeling Predator-Prey interaction are of current interest. Standard approaches such as first-order (in time) finite difference schemes for approximating the solution are widely spread. Though, this paper shows that recent advance methods can be more favored. In this work, we have incorporated, throughout numerical comparison experiments, spectral methods, for the space discretization, in conjunction with second and fourth-order time integrating methods for approximating the solution of the reaction-diffusion differential equations. The results have revealed that these methods have advantages over the conventional methods, some of which to mention are: the ease of implementation, accuracy and CPU time.

Multidomain Hybrid Direct DG and Central Difference Methods for Viscous Terms in Hyperbolic-Parabolic Equations

A class of multidomain hybrid methods of direct discontinuous Galerkin(DDG)methods and central difference(CD)schemes for the viscous terms is pro-posed in this paper.Both conservative and nonconservative coupling modes are dis-cussed.To treat the shock wave,the nonconservative coupling mode automatically switch to conservative coupling mode to preserve the conservative property when discontinuities pass through the artificial interface.To maintain the accuracy of the hybrid methods,the Lagrange interpolation polynomials and their derivatives are reconstructed to handle the coupling cells in the DDG subdomain,while the values of ghost points for the CD subdomain are calculated by the approximate polynomials from the DDG methods.The linear stabilities of these methods are demonstrated in detail through von-Neumann analysis.The multidomain hybrid DDG and CD meth-ods are then extended to one-and two-dimensional hyperbolic-parabolic equations.Numerical results validate that the multidomain hybrid methods are high-order ac-curate in the smooth regions,robust for viscous shock simulations and capable to save computational cost.

COMPARISON OF TDOA LOCATION ALGORITHMS WITH DIRECT SOLUTION METHOD

For Time Difference Of Arrival(TDOA) location based on multi-ground stations scene,two direct solution methods are proposed to solve the target position in TDOA location.Therein,the solving methods are realized in the rectangular and polar coordinates.On the condition of rectangular coordinates,first of all,it solves the radial range between the target and reference station,then cal-culates the location of the target.In the case of polar coordinates,the azimuth between the target and reference station is solved first,then the radial range between the target and reference station is figured out,finally the location of the target is obtained.Simultaneously,the simulation and comparison analysis are given in detail,and show that the polar solving method has the better fuzzy performance than that of rectangular coordinate.

Underdetermined direction of arrival estimation with nonuniform linear motion sampling based on a small unmanned aerial vehicle platform

Uniform linear array(ULA)radars are widely used in the collision-avoidance radar systems of small unmanned aerial vehicles(UAVs).In practice,a ULA's multi-target direction of arrival(DOA)estimation performance suffers from significant performance degradation owing to the limited number of physical elements.To improve the underdetermined DOA estimation performance of a ULA radar mounted on a small UAV platform,we propose a nonuniform linear motion sampling underdetermined DOA estimation method.Using the motion of the UAV platform,the echo signal is sampled at different positions.Then,according to the concept of difference co-array,a virtual ULA with multiple array elements and a large aperture is synthesized to increase the degrees of freedom(DOFs).Through position analysis of the original and motion arrays,we propose a nonuniform linear motion sampling method based on ULA for determining the optimal DOFs.Under the condition of no increase in the aperture of the physical array,the proposed method obtains a high DOF with fewer sampling runs and greatly improves the underdetermined DOA estimation performance of ULA.The results of numerical simulations conducted herein verify the superior performance of the proposed method.

A PREDICT-CORRECT NUMERICAL INTEGRATION SCHEME FOR SOLVING NONLINEAR DYNAMIC EQUATIONS

A new numerical integration scheme incorporating a predict-correct algorithm for solving the nonlinear dynamic systems was proposed in this paper. A nonlinear dynamic system governed by the equation v = F（v, t） was transformed into the form as v = Hv ＋ f（v, t）. The nonlinear part f（v, t） was then expanded by Taylor series and only the first-order term retained in the polynomial. Utilizing the theory of linear differential equation and the precise time-integration method, an exact solution for linearizing equation was obtained. In order to find the solution of the original system, a third-order interpolation polynomial of v was used and an equivalent nonlinear ordinary differential equation was regenerated. With a predicted solution as an initial value and an iteration scheme, a corrected result was achieved. Since the error caused by linearization could be eliminated in the correction process, the accuracy of calculation was improved greatly. Three engineering scenarios were used to assess the accuracy and reliability of the proposed method and the results were satisfactory.

Adaptive Reduced Basis Methods Applied to Structural Dynamic Analysis

The reduced basis methods (RBM) have been demonstrated as a promising numerical technique for statics problems and are extended to structural dynamic problems in this paper. Direct step-by-step integration and mode superposition are the most widely used methods in the field of the finite element analysis of structural dynamic response and solid mechanics. Herein these two methods are both transformed into reduced forms according to the proposed reduced basis methods. To generate a reduced surrogate model with small size, a greedy algorithm is suggested to construct sample set and reduced basis space adaptively in a prescribed training parameter space. For mode superposition method, the reduced basis space comprises the truncated eigenvectors from generalized eigenvalue problem associated with selected sample parameters. The reduced generalized eigenvalue problem is obtained by the projection of original generalized eigenvalue problem onto the reduced basis space. In the situation of direct integration, the solutions of the original increment formulation corresponding to the sample set are extracted to construct the reduced basis space. The reduced increment formulation is formed by the same method as mode superposition method. Numerical example is given in Section 5 to validate the efficiency of the presented reduced basis methods for structural dynamic problems.

A two-step unconditionally stable explicit method with controllable numerical dissipations

A family of unconditionally stable direct integration algorithm with controllable numerical dissipations is proposed. The numerical properties of the new algorithms are controlled by three parameters α, β and γ. By the consistent and stability analysis, the proposed algorithms achieve the second-order accuracy and are unconditionally stable under the condition that α≥-0.5, β≤ 0.5 and γ≥-(1+α)/2. Compared with other unconditionally stable algorithms, such as Chang's algorithms and CR algorithm, the proposed algorithms are found to be superior in terms of the controllable numerical damping ratios. The unconditional stability and numerical damping ratios of the proposed algorithms are examined by three numerical examples. The results demonstrate that the proposed algorithms have a superior performance and can be used expediently in solving linear elastic dynamics problems.

Removing Random-Valued Impulse Noises by a Two-Staged Nonlinear Filtering Method

Digital images are frequently contaminated by impulse noise(IN)during acquisition and transmission.The removal of this noise from images is essential for their further processing.In this paper,a two-staged nonlinear filtering algorithm is proposed for removing random-valued impulse noise(RVIN)from digital images.Noisy pixels are identified and corrected in two cascaded stages.The statistics of two subsets of nearest neighbors are employed as the criterion for detecting noisy pixels in the first stage,while directional differences are adopted as the detector criterion in the second stage.The respective adaptive median values are taken as the replacement values for noisy pixels in each stage.The performance of the proposed method was compared with that of several existing methods.The experimental results show that the performance of the suggested algorithm is superior to those of the compared methods in terms of noise removal,edge preservation,and processing time.