Numerical Investigations on the Resonance Errors of Multiscale Discontinuous Galerkin Methods for One-Dimensional Stationary Schrödinger Equation

In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.

Mobility and dynamic erosion process of granular flow:insights from numerical investigation using material point method

In order to understand the dynamics of granular flow on an erodible base soil,in this paper,a series of material point method-based granular column collapse tests were conducted to investigate numerically the mobility and dynamic erosion process of granular flow subjected to the complex settings,i.e.,the aspect ratio,granular mass,friction and dilatancy resistance,gravity and presence of water.A set of power scaling laws were proposed to describe the final deposit characteristics of granular flow by the relations of the normalized run-out distance and the normalized final height of granular flow against the aspect ratio,being greatly affected by the complex geological settings,e.g.,granular mass,the friction and dilatancy resistance of granular soil,and presence of water in granular flow.An index of the coefficient of friction of granular soil was defined as a ratio of the target coefficient of friction over the initial coefficient of friction to quantify the scaling extent of friction change(i.e.,friction strengthening or weakening).There is a characteristic aspect ratio of granular column corresponding to the maximum mobility of granular flow with the minimum index of the apparent coefficient of friction.The index of the repose coefficient of friction of granular flow decreased gradually with the increase in aspect ratio because higher potential energy of granular column at a larger aspect ratio causes a larger kinetic energy of granular soil to weaken the friction of granular soil as a kind of velocity-related friction weakening.An increase in granular mass reduces gradually the indexes of the apparent and repose coefficients of friction of granular soil to enhance the mobility of granular flow.The mobility of granular flow increases gradually with the decrease in friction angle or increase in dilatancy angle of granular soil.However,the increase of gravity accelerates granular flow but showing the same final deposit profile without any dependence on gravity.The mobility of granular flow increases gradually by lowering the indexes of the apparent and repose coefficients of friction of granular flow while changing the surroundings,in turn,the dry soil,submerged soil and saturated soil,implying a gradually increased excessive mobility of granular flow with the friction weakening of granular soil.Presence of water in granular flow may be a potential catalyzer to yield a long run-out granular flow,as revealed in comparison of water-absent and water-present granular flows.In addition,the dynamic erosion and entrainment of based soil induced by granular flow subjected to the complex geological settings,i.e.,the aspect ratio,granular mass,gravity,friction and dilatancy resistance,and presence of water,were comprehensively investigated as well.

Experimental investigation of closed-loop active control to modulate coherent structures by mu-level method

The experimental research on zero-net-mass-flux jet closed-loop active control was conducted in the wind tunnel.The mu-level method successfully detected burst events of the coherent structures. The streamwise velocity signals in the turbulent boundary layer were measured by HWA. The drag reduction rate of 16.7% is obtained comparable to that of the open-loop control and saves 75% of the input energy at the asynchronous 100 V/160 Hz control case, which reflects the advantages of the closed-loop control. The experimental findings indicate that the intensity increases in the near-wall region.The perturbation of the PZT vibrators on the skewness factor is concentrated in the region y+< 60. The generation of highspeed fluids is depressed and the downward effect of high-speed fluids weakens. The alteration of energy distribution and the discernible impact of modulation between structures of varying scales are observed. The correlation coefficient exhibits a strong positive correlation, which indicates that the large-scale structures produce modulation effect on small-scale ones.The occurrence of burst events is effectively suppressed. The disturbance has the characteristics of stable periodicity,positive and negative symmetry, low intermittency, and high pulsation strength. The conditional phase waveform shows that the fluctuation amplitude increases, indicating amplitude modulation effects on coherent structures.

NUMERICAL INVESTIGATION OF TOROIDAL SHOCK WAVES FOCUSING USING DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD

A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.

INVESTIGATION OF FLIGHT DYNAMICS OF THRUST VECTORING AIRCRAFT USING EXTENDED CONTINUATION METHODS

This paper presents the flight dynamical behavior of the thrust vectoring aircraft with extended bifurcation and continuation methods. In contrast to the standard bifurcation and continuation methods, the extended methods are capable of calculating the continuation curves of the equilibrium points for the particular type of trimming flight. Therefore, these methods can not only give the performance measures of aircraft, but also determine the stability of trimming points. In this paper, the methods are used to verify the effectiveness of the thrust vectoring control law, to define the flight envelope boundary, to analyze the stability and controllability of trimming flight, and to predict the departures of the instable flight. The result shows that the extended methods provide more flight dynamic information and are useful in preliminary design of the thrust vectoring aircraft.

Offshore fouling: investigation methods

On the hasis of the reason that the offshore oil development areas in the northern South China Sea are at greater depth,more hazardous conditions and distance from shore, the methods and technologies used in coastal waters are notfeasible because of the higher risk of losing investigation equipment or facilities. A series of methods and technologiesfor offshore fouling investigation that have been successfully applied is given in detail and it is hoped that theexperience can be helpful for further studies.

Application of the High-resolution EM Method in the Investigation of B Tunnel

A successful case history of exploring for concealed structure using the high-resolution EM method in the investigation of the West-East Gas Pipeline Project＇s B Tunnel is presented in this paper. The high frequency electromagnetic image system named STRATAGEM EH4, operating at frequencies ranging from 90KHz to 1Hz, was used for data acquisition. The orthogonal components of the electromagnetic field were measured during the field acquisition and the relevant electromagnetic attributes of the object body were extracted from the electromagnetic data. Hybrid sources, consisting of natural and full tensor-controlled sources, were utilized to produce high-quality electromagnetic field data. B Tunnel lies in the western part of Hubei province, at depths of less than 200m. The geologic setting of B tunnel is very complex. Following an initial geologic investigation, an outcrop considered to be a bedrock interface by investigators, collapsed during tunneling operations. A second investigation applied high-resolution EM and seismic refraction methods to reveal a more complex geologic structure along the tunnel route. The predicted rock classes and fault were encountered during the subsequent tunneling operations.

Hybrid method for investigation of electromagnetic scattering from conducting target above the randomly rough surface

A current based hybrid method （HM） is proposed which combines the method of moment （MOM） with the Kirchhoff approximation （KA） for the analysis of scattering interaction between a two-dimensional （2D） infinitely long conducting target with arbitrary cross section and a one-dimensional （1D） Gaussian rough surface. The electromagnetic scattering region in the HM is split into KA region and MOM region. The electric field integral equation （EFIE） in MOM region （target） is derived, the computational time of the HM depends mainly on the number of unknowns of the target. The bistatic scattering coefficient for the infinitely long cylinder above the rough surface with Gaussian roughness spectrum is calculated, and the numerical results are compared and verified with those obtained by the conventional MOM, which shows the high efficiency of the HM. Finally, the influence of the size, location of the target, the rms height and correlation length of the rough surface on the bistatic scattering coefficient with different polarizations is discussed in detail.

Performance Investigation of Cone Dielectric Elastomer Actuator Using Taguchi Method

Dielectric elastomer actuator （DEA） show promise for mechatronic applications due to the advantages of dielectric elastomer, such as lightweight, flexible, low cost, high strain, etc, and many configurations of DEAs have been demonstrated. As a kind of linear actuator, cone DEAs are studied in some laboratory prototypes due to easy manufacturing, however, their performance have not been exploited fully. Based on the working principle of DEA, a four-bar linkage mechanism is designed to provide negative stiffness preload, which can increase displacement output of actuator （outer diameter 100 mm） to 17 mm. Three cone actuating units are assembled in parallel to enhance the maximum force output to 5.07 N. Loading experiments of actuator in forward and backward strokes are performed, the experimental results show that backward stroke has stronger actuating capability than forward stroke, accordingly application of actuator is recommended. Four factors rather than applied voltage, i.e., number of actuating units, pre-stretch ratio, inner diameter, and outer diameter, are determined as influencing factors for Taguchi method. Then the performance objectives of actuator, i.e., displacement output, maximum force output, and maximum work in backward stroke, are investigated based on L9（34） Taguchi orthogonal design. The mean signal-to-noise （S/N） ratio based on the larger-the-better criterion is calculated according to the acquired displacement and force output. Analytical results show that outer diameter has the most significant influence on displacement output, and maximum force out and work output are influenced most by number of actuating units. Inner diameter also has an important effect on the performance objectives of actuator, while pre-stretch ratio has the least influence. The proposed performance investigation is helpful for the design and application of cone actuator in mechatronic system.

Present state,basic theories,methods and progresses of investigation and assessment on marine hazardous geology in China

The achievement progresses of investigation and studies on marine hazardous geology are summarized and presentsd in the late 20 century in China. The importance, research value and present-day studies of marine hazardous geology, a newly developing branch of geoscience, are well expatiated. Several often confused concepts and theories are explained and redefined here. The comment on the means of investigations, assessment of marine hazardous geology, as well as its evolution, innovation, existing questions and future tasks are also introduced and presented. The concepts of 'hazard geology', geohazard', 'map of marine hazard geology', 'integrated evaluaton on seafloor stablity' are respectively discussed, including their definition, research objects, methods and contents. The types and classification of marine hazardous geology, principles and methods of marine hazardous geology map compilation, the assessment methods and models of marine hazardous geology environment and seafloor stability and so on are also discussed.

Potato farmers’ preference for agricultural insurance in China: An investigation using the choice experimental method

Potato insurance plays an important role in transferring agricultural risks to promote the potato staple strategy. Understanding farmers’ real preferences for potato insurance attributes is important to improve the potato insurance. In this study, a choice experiment was designed with attributes such as peril, minimum compensation ratio for production loss, types of crops covered by insurance, annual premium per mu after subsidy, and complexity of claims procedures. We constructed a mixed logit model based on a questionnaire survey of 362 potato farmers’ choices of insurance attributes from 24 villages in Dingxi City of Gansu Province using random sampling. The results show that:(1) farmers prefer agricultural insurance with widely perils including output price drop and input cost risk;(2) farmers who have suffered plant disease and insect and pest damage are willing to pay a high premium for the insurance with a high compensation ratio and several types of crops covered;full-time farmers and large-scale farmers preferred insurance products with low complexity of claims procedures. Therefore, new potato insurance design could be considered specifically for different farmer groups such as those who had previous disaster experiences, and large-or small-scale farm operations.

Analytical investigation of Jeffery-Hamel flow with high magnetic field and nanoparticle by Adomian decomposition method

In this study, the effects of magnetic field and nanoparticle on the Jeffery- Hamel flow are studied using a powerful analytical method called the Adomian decomposition method （ADM）. The traditional Navier-Stokes equation of fluid mechanics and Maxwell＇s electromagnetism governing equations are reduced to nonlinear ordinary differential equations to model the problem. The obtained results are well agreed with that of the Runge-Kutta method. The present plots confirm that the method has high accuracy for different a, Ha, and Re numbers. The flow field inside the divergent channel is studied for various values of Hartmann ：number and angle of channel. The effect of nanoparticle volume fraction in the absence of magnetic field is investigated.

Investigation on mechanism of magnetization reversal for nanocrystalline Pr-Fe-B permanent magnets by micromagnetic finite element methods

Magnetization configurations were calculated under various magnetic fields for nanocrystalline Pr-Fe-B permanent magnets by micromagnetic finite element method.According to the configurations during demagnetization process, the mechanism of magnetization reversal was analyzed.For the Pr2Fe14B with 10 nm grains or its composite with 10vol.% α-Fe, the coercivity was determined by nucleation of reversed domain that took place at grain boundaries.However, for Pr2Fe14B with 30 nm grains, coercivity was controlled by pinning of the nucle-ated domain.For Pr2Fe14B/α-Fe with 30vol.% α-Fe, the demagnetization behavior was characterized by continuous reversal of α-Fe moment.

Experimental investigation of interface curing stresses between PMMA and composite using digital speckle correlation method

This paper studies the interface curing stresses between polymethyl methacrylate (PMMA) and composite by means of digital speckle correlation method (DSCM).A new method by combining DSCM with the marker points is developed to measure the interface curing stresses,and the measurement principle is introduced.The interface curing stresses between PMMA and composite with different curing bonding conditions are measured and analyzed,this indicates that the residual stress for furnace heating and furnace cooling is the smallest.Finally,the measurement error is discussed by means of finite element method,the influences of glass microsphere between adhesive and PMMA can be ignored.

Investigation on charge parameters of underwater contact explosion based on axisymmetric SPH method

This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to improve the pressure fluctuation and numerical accuracy near the symmetric axis. An in-depth study is carried out over the influence of charge shapes and detonation modes on the near-field loads in terms of the peak pressure and impulse of shock waves. For different charge shapes, the cylindrical charge with different length-diameter ratios may cause strong directivity of peak pressure and impulse in the near field. Compared with spherical charge, the peak pressure of cylindrical charge may be either weakened or enhanced in different directions. Within a certain range, the greater the length-diameter ratio is, the more obvious the effect will be. The weakened ratio near the detonation end may reach 25% approximately, while the enhanced ratio may reach around 20% in the opposite direction. However, the impulse in different directions seems to be uniform. For different detonation modes, compared with point-source explosion, the peak pressure of plane-source explosion is enhanced by about 5%. Besides, the impulse of plane-source explosion is enhanced by around 5% near the detonation end, but close to those of the point-source explosion in other directions. Based on the material constitutive relation in the axisymmetric coordinates, a simple case of underwater contact explosion is simulated to verify the above conclusions, showing that the charge parameters of underwater contact explosion should not be ignored.

Numerical Investigation on Vortex-Induced Rotations of A Triangular Cylinder Using An Immersed Boundary Method

A numerical study of vortex-induced rotations(VIRs) of an equivalent triangular cylinder, which is free to rotate in the azimuthal direction in a uniform flow, is presented. Based on an immersed boundary method, the numerical model is established, and is verified through the benchmark problem of flow past a freely rotating rectangular body.The computation is performed for a fixed reduced mass of m~*=2.0 and the structural stiffness and damping ratio are set to zero. The effects of Reynolds number(Re=25-180) on the characteristics of VIR are studied. It is found that the dynamic response of the triangular cylinder exhibits four distinct modes with increasing Re: a rest position,periodic rotational oscillation, random rotation and autorotation. For the rotational oscillation mode, the cylinder undergoes a periodic vibration around an equilibrium position with one side facing the incoming flow. Since the rotation effect, the outset of vortex shedding from cylinder shifts to a much lower Reynolds number. Further increase in Re leads to 2 P and P+S vortex shedding modes besides the typical 2 S pattern. Our simulation results also elucidate that the free rotation significantly changes the drag and lift forces. Inspired by these facts, the effect of free rotation on flow-induced vibration of a triangular cylinder in the in-line and transverse directions is investigated. The results show that when the translational vibration is coupled with rotation, the triangular cylinder presents a galloping response instead of vortex-induced vibration(VIV).

Numerical investigations on mechanical characteristics and failure mechanism of outwash deposits based on random meso-structures using discrete element method

Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this work, random meso-structure of outwash deposits was constructed by the technique of computer random simulation based on characteristics of its meso-structure in the statistical sense and some simplifications, and a series of large direct shear tests on numerical samples of outwash deposits with stone contents of 15%, 30%, 45% and 60% were conducted using the discrete element method to further investigate its mechanical characteristics and failure mechanism under external load. The results show that the deformation characteristics and shear strength of outwash deposits are to some extent improved with the increase of stone content, and the shear stress–shear displacement curves of outwash deposits show great differences at the post-peak stage due to the random spatial distribution and content of stones. From the mesoscopic view, normal directions of contacts between "soil" and "stone" particles undergo apparent deflection as the shear displacement continues during the shearing process, accompanying redistribution of the magnitude of contact forces during the shearing process. For outwash deposits, the shear zone formed after shear failure is an irregular stripe due to the movements of stones near the shear zone, and it expands gradually with the increase of stone content. In addition, there is an approximately linear relation between the mean increment of internal friction angle and the stone content lying between 30% and 60%, and a concave nonlinear relation between the mean increment of cohesion and stone content, which are in good agreement with the existing research results.

Numerical and experimental investigation of rock breaking method under free surface by TBM disc cutter

To study the rock breaking method under the free surface induced by disc cutter,the rock breaking simulations were first conducted based on the discrete element method,and the dynamic process of rock breaking under the free surface was studied including stressed zone,crush zone,crack initiation and propagation.Then the crack propagation conditions,specific energy,etc.under different free surface distance(S)were also investigated combined with linear cutting experiments.The results show that the rock breaking process under the free surface induced by disc cutter is dominated by tension failure mode.There exists a critical S to promote crack propagation to free surface effectively.And this rock breaking method can improve the rock breaking force and breaking efficiency significantly when proper.

Thermoelastic Stress Field Investigation of GaN Material for Laser Lift-off Technique based on Finite Element Method

The transient thermoelastic stress fields of GaN films is analyzed by the finite element method for the laser lift-off （LLO） technique. Stress distributions in GaN films irradiated by pulse laser with different energy densities as functions of time and depth are simulated. The results show that the high thermoelastic stress distributions in GaN films localize within about 1 μm below the GaN/Al2O3 interface using proper laser parameters. It is also found that GaN films can avoid the thermal deformation because the maximum thermoelastic stress 4.28 GPa is much smaller than the yield strength of GaN 15GPa. The effects of laser beam dimension and the thickness of GaN films on stress distribution are also analyzed. The variation range of laser beam dimension as a function of the thickness of GaN films is simulated to keep the GaN films free of thermal deformation. LLO experiments are also carried out. GaN-based light-emitting diodes （LEDs） are separated from sapphire substrates using the parameters obtained from the simulation. Compared with devices before LLO, P-I-V measurements of GaN- based LEDs after LLO show that the electrical and optical characteristics improve greatly, indicating that no stress damage is brought to GaN films using proper parameters obtained by calculation during LLO.

Investigation on Hard Magnetic Properties of Nanocomposite Permanent Magnets with Precipitate-Typed Microstructure by Micromagnetic Finite-Element Methods

The demagnetization curves were calculated using micromagnetic finite-element method for nanocomposite Pr 2Fe 14B/α-Fe permanent magnets with precipitate-typed microstructure. Due to intergrain exchange coupling, both remanence enhancement and a single magnetic phase behavior in demagnetization curve were found. For the samples with the hard phase as the precipitate and the soft one as the matrix, a coercivity μ 0H c of 0.78 T, a remanence J r of 1.18 T and a large energy product (BH) max of 200 kJ·m -3 are obtained for the sample with hard grain size being 23 nm. While, for the sample with the soft phase as the precipitate, μ 0H c of 0.95 T, J r of 1.24 T and (BH) max of 240 kJ·m -3 are obtained for the sample with soft grain size being 10 nm. The calculating results were compared with the experimental Pr 8Fe 87B 5 ribbons. The dependence of remanence and coercivity on the microstructure was discussed extensively.