Superconvergence of Direct Discontinuous Galerkin Methods:Eigen-structure Analysis Based on Fourier Approach

This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.

Evaluation of bronchoscopic direct vision glottis anesthesia method in bronchoscopy

BACKGROUND Fibrobronchoscopy is a common adjunct tool that requires anesthesia and is widely used in the diagnosis and treatment of various respiratory diseases.However,current anesthesia methods,such as spray,nebulized inhalation,and cricothyroid membrane puncture,have their own advantages and disadvantages.Recently,studies have shown that bronchoscopic direct-view glottis anesthesia is a simple and inexpensive method that shortens the examination time and provides excellent anesthetic results.AIM To evaluate the effectiveness of bronchoscopic direct vision glottis anesthesia for bronchoscopy.METHODS The study included 100 patients who underwent bronchoscopy during thoracic surgery.A random number table method was used to divide the patients into control and observation groups(50 patients each).The control and observation groups were anesthetized using the nebulized inhalation and bronchoscopic direct vision glottis method,respectively.Hemodynamic indices[systolic blood pressure(SBP),diastolic blood pressure(DBP),heart rate(HR),and oxygen saturation(SpO_(2))before(T1),5 min after anesthesia(T2),and at the end of the operation(T3)]serum stress hormone indices[norepinephrine(NE),epinephrine(E),adrenocorticotropic hormone(ACTH),and cortisol(Cor)before and after treatment]were compared between the 2 groups.Adverse effects were also RESULTS At T2 and T3,SBP,DBP,and HR were lower in the observation group than the control group,whereas SpO_(2) was higher than the control group[(119.05±8.01)mmHg vs(127.05±7.83)mmHg,(119.35±6.66)mmHg vs(128.39±6.56)mmHg,(84.68±6.04)mmHg vs(92.42±5.57)mmHg,(84.53±4.97)mmHg compared to(92.57±6.02)mmHg,(74.25±5.18)beats/min compared to(88.32±5.72)beats/min,(74.38±5.31)beats/min compared to(88.42±5.69)beats/min,(97.36±2.21)%vs(94.35±2.16)%,(97.42±2.36)%vs(94.38±2.69%],with statistically significant differences(all P<0.05).After treatment,NE,E,ACTH,and Cor were significantly higher in both groups than before treatment,but were lower in the observation group than in the control group[(68.25±8.87)ng/mL vs(93.35±14.00)ng/mL,(53.59±5.89)ng/mL vs(82.32±10.70)ng/mL,(14.32±1.58)pg/mL vs(20.35±3.05)pg/mL,(227.35±25.01)nmol/L vs(322.28±45.12)nmol/L],with statistically significant differences(all P<0.05).The incidence of adverse reactions was higher in the control group than in the observation group[12.00%(12/50)vs 6.00%(3/50)](P<0.05).CONCLUSION The use of bronchoscopic direct vision glottis anesthesia method for bronchoscopy patients is beneficial for stabilizing hemodynamic indices during bronchoscopy and reducing the level of patient stress,with good safety and practicality.

Greedy Randomized Gauss-Seidel Method with Oblique Direction

For the linear least squares problem with coefficient matrix columns being highly correlated, we develop a greedy randomized Gauss-Seidel method with oblique direction. Then the corresponding convergence result is deduced. Numerical examples demonstrate that our proposed method is superior to the greedy randomized Gauss-Seidel method and the randomized Gauss-Seidel method with oblique direction.

Optimization of Random Feature Method in the High-Precision Regime

Machine learning has been widely used for solving partial differential equations(PDEs)in recent years,among which the random feature method(RFM)exhibits spectral accuracy and can compete with traditional solvers in terms of both accuracy and efficiency.Potentially,the optimization problem in the RFM is more difficult to solve than those that arise in traditional methods.Unlike the broader machine-learning research,which frequently targets tasks within the low-precision regime,our study focuses on the high-precision regime crucial for solving PDEs.In this work,we study this problem from the following aspects:(i)we analyze the coeffcient matrix that arises in the RFM by studying the distribution of singular values;(ii)we investigate whether the continuous training causes the overfitting issue;(ii)we test direct and iterative methods as well as randomized methods for solving the optimization problem.Based on these results,we find that direct methods are superior to other methods if memory is not an issue,while iterative methods typically have low accuracy and can be improved by preconditioning to some extent.

Impact Force Localization and Reconstruction via ADMM-based Sparse Regularization Method

In practice,simultaneous impact localization and time history reconstruction can hardly be achieved,due to the illposed and under-determined problems induced by the constrained and harsh measuring conditions.Although l_(1) regularization can be used to obtain sparse solutions,it tends to underestimate solution amplitudes as a biased estimator.To address this issue,a novel impact force identification method with l_(p) regularization is proposed in this paper,using the alternating direction method of multipliers(ADMM).By decomposing the complex primal problem into sub-problems solvable in parallel via proximal operators,ADMM can address the challenge effectively.To mitigate the sensitivity to regularization parameters,an adaptive regularization parameter is derived based on the K-sparsity strategy.Then,an ADMM-based sparse regularization method is developed,which is capable of handling l_(p) regularization with arbitrary p values using adaptively-updated parameters.The effectiveness and performance of the proposed method are validated on an aircraft skin-like composite structure.Additionally,an investigation into the optimal p value for achieving high-accuracy solutions via l_(p) regularization is conducted.It turns out that l_(0.6)regularization consistently yields sparser and more accurate solutions for impact force identification compared to the classic l_(1) regularization method.The impact force identification method proposed in this paper can simultaneously reconstruct impact time history with high accuracy and accurately localize the impact using an under-determined sensor configuration.

A modified stochastic finite-fault method for estimating strong ground motion:Validation and application

We developed a modified stochastic finite-fault method for estimating strong ground motions.An adjustment to the dynamic corner frequency was introduced,which accounted for the effect of the location of the subfault relative to the hypocenter and rupture propagation direction,to account for the influence of the rupture propagation direction on the subfault dynamic corner frequency.By comparing the peak ground acceleration(PGA),pseudo-absolute response spectra acceleration(PSA,damping ratio of 5%),and duration,the results of the modified and existing methods were compared,demonstrating that our proposed adjustment to the dynamic corner frequency can accurately reflect the rupture directivity effect.We applied our modified method to simulate near-field strong motions within 150 km of the 2008 MW7.9 Wenchuan earthquake rupture.Our modified method performed well over a broad period range,particularly at 0.04-4 s.The total deviations of the stochastic finite-fault method(EXSIM)and the modified EXSIM were 0.1676 and 0.1494,respectively.The modified method can effectively account for the influence of the rupture propagation direction and provide more realistic ground motion estimations for earthquake disaster mitigation.

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.

Direction finding of bistatic MIMO radar in strong impulse noise

For bistatic multiple-input multiple-output(MIMO)radar,this paper presents a robust and direction finding method in strong impulse noise environment.By means of a new lower order covariance,the method is effective in suppressing impulse noise and achieving superior direction finding performance using the maximum likelihood(ML)estimation method.A quantum equilibrium optimizer algorithm(QEOA)is devised to resolve the corresponding objective function for efficient and accurate direc-tion finding.The results of simulation reveal the capability of the presented method in success rate and root mean square error over existing direction-finding methods in different application situations,e.g.,locating coherent signal sources with very few snapshots in strong impulse noise.Other than that,the Cramér-Rao bound(CRB)under impulse noise environment has been drawn to test the capability of the presented method.

Finite element numerical simulation of 2.5D direct current method based on mesh refinement and recoarsement

To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direct current method, we propose a new mesh refinement and recoarsement method for a two-dimensional point source. We introduce the mesh refinement and mesh recoarsement into the traditional structured mesh subdivision. By refining the horizontal grids, the singularity owing to the point source is minimized and the topography is simulated. By recoarsening the horizontal grids, the number of grid cells is reduced significantly and computational efficiency is improved. Model tests show that the proposed method solves the singularity problem and reduces the number of grid cells by 80% compared to the uniform grid refinement.

THREE-DIMENSIONAL TRAJECTORY OPTIMIZATION WITH DIRECT METHOD

The principle of direct method used in optimal control problem is introduced. Details of applying this method to flight trajectory generation are presented including calculation of velocity and controls histories. And capabilities of flight and propulsion systems are considered also. Combined with digital terrain map technique, the direct method is applied to the three dimensional trajectory optimization for low altitude penetration, and simplex algorithm is used to solve the parameters in optimization. For the small number of parameters, the trajectory can be optimized in real time on board.

Fast alternating direction method of multipliers for total-variation-based image restoration

A novel algorithm, i.e. the fast alternating direction method of multipliers （ADMM）, is applied to solve the classical total-variation （ TV ）-based model for image reconstruction. First, the TV-based model is reformulated as a linear equality constrained problem where the objective function is separable. Then, by introducing the augmented Lagrangian function, the two variables are alternatively minimized by the Gauss-Seidel idea. Finally, the dual variable is updated. Because the approach makes full use of the special structure of the problem and decomposes the original problem into several low-dimensional sub-problems, the per iteration computational complexity of the approach is dominated by two fast Fourier transforms. Elementary experimental results indicate that the proposed approach is more stable and efficient compared with some state-of-the-art algorithms.

Direct computing methods for turn flows in traffic assignment

Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks. The first one is designed on the path-turn incidencerelationship, and it is similar to the computational procedure of link flows. It applies to thetraffic assignment algorithms that can provide detailed path structures. The second utilizes thelink-turn incidence relationship and the conservation of flow on links, a law deriving from thisrelationship. It is actually an improved version of Dial's logit assignment algorithm. The proposedapproaches can avoid the shortcomings both of the estimation methods, e. g. Furness's model andFrator's model, and of the network-expanding method in precision, stability and computation scale.Finally, they are validated by numerical examples.

DIRECTED HYPERGRAPH THEORY AND DECOMPOSITION CONTRACTION METHOD

A new branch of hypergraph theory-directed hyperaph theory and a kind of new methods-dicomposition contraction(DCP, PDCP and GDC) methods are presented for solving hypernetwork problems.lts computing time is lower than that of ECP method in several order of magnitude.

Preparation of Al-Al_3Ti in situ composites by direct reaction method

Al-Al3Ti composites were prepared by a direct reaction method, in which Al3Ti was formed by the reaction of Ti and Al in aluminum alloy melt. The morphology of Al3Ti changes apparently from the fine particle, needle-like to large block with the increase of Al3Ti content. The addition of magnesium can markedly change the morphology of Al3Ti and reduce their size. Short rod-like Al3Ti was formed and homogeneous distribution was obtained with the addition of 3 wt.% Mg. The effect of Al3Ti and Mg on the microstructure of Al-Al3Ti composites and the mechanism were also discussed.

Image reconstruction based on total-variation minimization and alternating direction method in linear scan computed tomography

Linear scan computed tomography （LCT） is of great benefit to online industrial scanning and security inspection due to its characteristics of straight-line source trajectory and high scanning speed. However, in practical applications of LCT, there are challenges to image reconstruction due to limited-angle and insufficient data. In this paper, a new reconstruction algorithm based on total-variation （TV） minimization is developed to reconstruct images from limited-angle and insufficient data in LCT. The main idea of our approach is to reformulate a TV problem as a linear equality constrained problem where the objective function is separable, and then minimize its augmented Lagrangian function by using alternating direction method （ADM） to solve subproblems. The proposed method is robust and efficient in the task of reconstruction by showing the convergence of ADM. The numerical simulations and real data reconstructions show that the proposed reconstruction method brings reasonable performance and outperforms some previous ones when applied to an LCT imaging problem.

Study on Measurement Methods of Multi-Directional Waves—Wave Gauge Array Method

Measuring multi-directional waves with the wave gauge array is one of the fundamental and easily realised methods. In this paper, the wave gauge array is described and the effects of the gauge spacing, the array orientations, etc. of the three array arrangements, i. e., linear array, T-type array and pentagon array, on the resolution of the directional spreading of waves, are investigated experimentally. This study can be used as a reference in the experimental study and the field measurement of directional waves.

Direct discontinuous Galerkin method for the generalized Burgers-Fisher equation

In this study, we use the direct discontinuous Galerkin method to solve the generalized Burgers-Fisher equation. The method is based on the direct weak formulation of the Burgers-Fisher equation. The two adjacent cells are jointed by a numerical flux that includes the convection numerical flux and the diffusion numerical flux. We solve the ordinary differential equations arising in the direct Galerkin method by using the strong stability preserving Runge^Kutta method. Numerical results are compared with the exact solution and the other results to show the accuracy and reliability of the method.

A laboratory study of directional spectra with maximum likelihood method─I Developing wind wave

From measurements by a circular array consisting of 18 wave gauges in a large wave tank, directional spectra of wind-generated waves in deep water are systematically determined by using maximum likehood method.The investigations reveal that the angular spreading of the wave energy is consistent with cos2s(θ/2) proposed by Longuet-Higgins et al. (1963, Ocean Wad Spectra,11~136), if the bimodal distributions of wave energy are not taken into account. Bimodality occurring on higher frequency than peak frequency is too rare to affect our whole results. Surprisingly, a much broader directional spreading than that of the field, which is interpreted by the strongly nonlinear energy transfer because of the very young waves in laboratory, is found. The parameter s depends on frequency in the same way as observed by Mitsuyasu et al. (1975, Journal of Physical Oceanography, 5, 750~760)and Hasselmann et al. (1980, Journal of physical Oceanography, 10, 1264~1280) in the field, and the relationship between the four nondimensional parameters sm, fo, b1 and b2, determining the directional width, and (corresponding to the inverse of wave age) are given respectively. The observed distributions are found to agree well with the suggestion of Donelan et al. (1985, Philosophical Transaction of Royal Society of London, A315, 509~562) when applied to field waves.

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.

Advances in traction methods for endoscopic submucosal dissection:What is the best traction method and traction direction?

Endoscopic submucosal dissection(ESD)has been developed as a treatment for superficial gastrointestinal neoplasms,which can achieve en bloc resection regardless of the lesion size.However,ESD is technically difficult because endoscopists cannot bring their hand into the gastrointestinal tract,unlike surgeons in regular surgery.It is difficult to obtain sufficient tension in the dissection plane and a good field of vision.Therefore,ESD is associated with a long procedure time and a high risk of adverse events in comparison with endoscopic mucosal resection.Traction methods have been developed to provide sufficient tension for the dissection plane and a good field of vision during the ESD procedure.However,traction direction is limited in most traction methods,resulting in insufficient effect in some cases.Although traction direction is considered important,there have been few investigations of its effect.In the first half of this review,important traction methods are discussed,including traction direction.In second half,appropriate traction methods for each organ are considered.Other important considerations for traction method,such as ability to adjust traction strength,interference between traction device and endoscope,and the need for specialized devices are also discussed.