A Method for Calculation of Low-Frequency Slow Drift Motions Based on NURBS for Floating Bodies

Through a higher-order boundary element method based on NURBS （Non-uniform Rational B-splines）, the calculation of second-order low-frequency forces and slow drift motions is conducted for floating bodies. In the floating body＇s inner domain, an auxiliary equation is obtained by applying a Green function which satisfies the solid surface condition. Then, the auxiliary equation and the velocity potential equation are combined in the fluid domain to remove the solid angle coefficient and the singularity of the double layer potentials in the integral equation. Thus, a new velocity potential integral equation is obtained. The new equation is extended to the inner domain to reheve the irregular frequency effects; on the basis of the order analysis, the comparison is made about the contribution of all integral terms with the result in the second-order tow-frequency problem; the higher-order boundary element method based on NURBS is apphed to calculate the geometric position and velocity potentials; the slow drift motions are calculated by the spectrum analysis method. Removing the solid angle coefficient can apply NURBS technology to the hydrodynamic calculation of floating bodies with complex surfaces, and the extended boundary integral method can reduce the irregular frequency effects. Order analysis shows that free surface integral can be neglected, and the numerical results can also prove the correctness of order analysis. The results of second-order low-frequency forces and slow drift motions and the comparison with the results from references show that the application of the NURBS technology to the second-order low-frequency problem is of high efficiency and credible results.

New methods to remove baseline drift in trapezoidal pulse shaping

Trapezoidal pulse shaping algorithm is widely applied to improve signal-to-noise ratio(SNR), throughput and energy resolution with the properties of noise suppression, pile-up pulse separation and ballistic deficit correction. The algorithm can be acquired by z transform method which is easier for derivation. However, the baseline drift of trapezoidal pulse appears because the noise superimposes on the input signal. In this paper,two new methods based on convergence analysis and noise suppression are proposed to remove the baseline drift resulting from trapezoidal pulse shaping. Simulations and experimental tests are carried out to verify the methods. The results demonstrate that the proposed methods can remove baseline drift in trapezoidal pulse shaping.

Study on the Drift Calibration Method and Device for the Tin Oxide pH Electrode

Drift phenomenon has been known as the drawback of sensors and causes inaccuracy on the long-term measurement. In general,there are two methods to reduce the drift problem.One is to tune the parameters of the fabrication process to improve the properties of the front-ended device.Another is to compensate the drift phenomenon by adding extra drift compensation circuit or software in the back-ended readout circuit.In this study,a drift calibration method for the potentiometric sensor was presented and the drift calibration method was performed by using the circuit.According to experimental results,the drift phenomenon of the SnO_2 pH electrode was reduced by the drift calibration device.

Effect of a two-phase wedge-sliding model on the ingredient drift of a stable mixed fluid and its computing method

A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-uniform stress field and to analyse its phase drift phenomenon. In the model, the drift-inhibition angle and the expansion-inhibition angle are also deduced and used as evaluating indexes to describe the drifting trend of different ingredients among the mixed fluids. For solving above two indexes of the model, a new calculation method is developed and used to compute the phase distributions of multiphase fluid at peak stress and gradient area stress, respectively. As an example, the flow process of grease in a pipe is analysed by simulation method and used to verify the validity of the model.

Dual-phase coexistence enables to alleviate resistance drift in phase-change films

The amorphous phase-change materials with spontaneous structural relaxation leads to the resistance drift with the time for phase-change neuron synaptic devices. Here, we modify the phase change properties of the conventional Ge_2Sb_2Te_5(GST) material by introducing an SnS phase. It is found that the resistance drift coefficient of SnS-doped GST was decreased from 0.06 to 0.01. It can be proposed that the origin originates from the precipitation of GST nanocrystals accompanied by the precipitation of SnS crystals compared to single-phase GST compound systems. We also found that the decrease in resistance drift can be attributed to the narrowed bandgap from 0.65 to 0.43 eV after SnS-doping. Thus, this study reveals the quantitative relationship between the resistance drift and the band gap and proposes a new idea for alleviating the resistance drift by composition optimization, which is of great significance for finding a promising phase change material.

Shift control method for the local time at descendi ngnode based on sun-synchronous orbit satellite

This article analyzes the shift factors of the descending node local time for sun-synchronous satellites and proposes a shift control method to keep the local time shift within an allowance range. It is found that the satellite orbit design and the orbit injection deviation are the causes for the initial shift velocity, whereas the atmospheric drag and the sun gravitational perturbation produce the shift acceleration. To deal with these shift factors, a shift control method is put forward, through such methods as orbit variation design, orbit altitude, and inclination keeping control. The simulation experiment and practical application have proved the effectiveness of this control method.

Using a Time-domain Higher-order Boundary Element Method to Simulate Wave and Current Diffraction from a 3-D Body

To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.

Wave Drift Forces on Two Floating Bodies Arranged Side by Side

An innovative hydrodynamic theory and numerical model were developed to help improve the efficiency, accuracy, and convergence of the numerical prediction of wave drift forces on two side-by-side deepwater floating bodies. The wave drift forces were expressed by the double integration of source strength and the corresponding Green function on the body surface, which is consistent with the far field formula based on momentum conservation and sharing the advantage of near field calculations providing the drift force on each body. Numerical results were validated through comparing the general far field model and pressure integral model, as well as the middle field model developed usin^z the software HydroStar.

Investigation Responses of the Diagrid Structural System of High-rise Buildings Equipped with Tuned Mass Damper Using New Dynamic Method

Due to the shortage of land in cities and population growth,the significance of high rise buildings has risen.Controlling lateral displacement of structures under different loading such as an earthquake is an important issue for designers.One of the best systems is the diagrid method which is built with diagonal elements with no columns for manufacturing tall buildings.In this study,the effect of the distribution of the tuned mass damper(TMD)on the structural responses of diagrid tall buildings was investigated using a new dynamic method.So,a diagrid structural systems with variable height with TMDs was solved as an example of structure.The reason for the selection of the diagrid system was the formation of a stiffness matrix for the diagonal and angular elements.Therefore,the effect of TMDs distribution on the story drift,base shear and structural behaviour were studied.The obtained outcomes showed that the TMDs distribution does not significantly affect on improving the behaviour of the diagrid structural system during an earthquake.Furthermore,the new dynamic scheme represented in this study has good performance for analyzing different systems.Abbreviation:TMD-tuned mass damper;SATMD-semiactive-tuned mass dampers;MDOF-multiple degrees of freedom;m_(i)-mass of ith story of the building;c_(i)-damping coefficient of the ith story of the building;k_(i)-stiffness of ith story of the building;x_(i)-displacement of the ith story of the building;md-mass of damper;c_(d)-damping coefficient of the damper;k_(d)-stiffness of damper;x_(d)-displacement of TMD;M_(i)-generalized mass of the ith normal mode;C_(i)-generalized damping of the ith normal mode;K_(i)-generalized stiffness of the ith normal mode;K_(i)(t)-generalized load of the ith normal mode;Y_(i)(t)-generalized displacement of the ith normal mode;[M]-matrices of mass;[C]-matrices of damping;{P(t)}-consequence external forces;N_(i)(τ)-interpolation functions;[Ai]-mechanical properties of the structure.

Combined Active Islanding Detection Method for Grid-connected Photovoltaic System

Unintentional islanding phenomenon has been one of the most important problems of gridconnected photovoltaic inverters. To prevent this phenomenon, all kinds of anti-islanding methods have been discussed. This paper presents a combined active islanding detection method, which consists of active frequency drift method and automatic phase-shift method. The traditional active anti-islanding methods of grid-connected PV inverters bear nondetection zone possibilities for certain paralleled RLC loads. The combined method shows islanding detection ability effectively, and it can eliminate nondetection zones even in the worst case conditions. Simulation in different load conditions is performed for verification.

DRIFT MOTION OF FREE-ROTOR GYROSCOPE WITH RADIAL MASS-UNBALANCE

The motion of a rigid body about fixed point with small radial mass-unbalance in homogeneous gravitational field was discussed. The dynamical equations described by state variables of the body were established, and approximate analytical solutions for a spinning body with high speed were obtained by use of the average method. The influence of the radial mass-unbalance of the rotor to the precession character of a free-rotor gyroscope was analyzed. And a physical explanation of the drift phenomenon of the gyro was given. An applicable formula of gyro's constant drift in analytical form was obtained, which is perfectly coincident with the numerical calculation.

Dimension-reduction of FPK equation via equivalent drift coefficient

The Fokker–Planck–Kolmogorov（FPK） equation plays an essential role in nonlinear stochastic dynamics. However, neither analytical nor numerical solution is available as yet to FPK equations for high-dimensional systems. In the present paper, the dimension reduction of FPK equation for systems excited by additive white noise is studied. In the proposed method, probability density evolution method（PDEM）, in which a decoupled generalized density evolution equation is solved, is employed to reproduce the equivalent flux of probability for the marginalized FPK equation. A further step of constructing an equivalent coefficient finally completes the dimension-reduction of FPK equation. Examples are illustrated to verify the proposed method.

A numerical method based on probability theory

By using the connections between Brownian family with drift and elliptic differential equations, an efficient probabilistic computing method is given. This method is applied to a wide-range Dirichlet problem. Detail ana- lysis and deduction of solving the problem are offered. The stochastic representation of the solution to the problem makes a 3-dimensional problem turned into a 2-dimensional problem. And an auxiliary ball is constructed. The strong Markov property and the joint distributions of the time and place of hitting spheres for Brownian family with drift are employed. Finally, good convergence of the numerical solution to the problem over domain with arbitrary boundary is obtained.

An improved optical flow method to estimate Arctic sea ice velocity(winter 2014-2016)

Sea ice velocity impacts the distribution of sea ice,and the flux of exported sea ice through the Fram Strait increases with increasing ice velocity.Therefore,improving the accuracy of estimates of the sea ice velocity is important.We introduce a pyramid algorithm into the Horn-Schunck optical flow(HS-OF)method(to develop the PHS-OF method).Before calculating the sea ice velocity,we generate multilayer pyramid images from an original brightness temperature image.Then,the sea ice velocity of the pyramid layer is calculated,and the ice velocity in the original image is calculated by layer iteration.Winter Arctic sea ice velocities from 2014 to 2016 are obtained and used to discuss the accuracy of the HS-OF method and PHS-OF(specifically the 2-layer PHS-OF(2 LPHS-OF)and 4-layer PHS-OF(4 LPHS-OF))methods.The results prove that the PHS-OF method indeed improves the accuracy of sea ice velocity estimates,and the 2 LPHS-OF scheme is more appropriate for estimating ice velocity.The error is smaller for the 2 LPHS-OF velocity estimates than values from the Ocean and Sea Ice Satellite Application Facility and the Copernicus Marine Environment Monitoring Service,and estimates of changes in velocity by the 2 LPHS-OF method are consistent with those from the National Snow and Ice Data Center.Sea ice undergoes two main motion patterns,i.e.,transpolar drift and the Beaufort Gyre.In addition,cyclonic and anticyclonic ice drift occurred during winter 2016.Variations in sea ice velocity are related to the open water area,sea ice retreat time and length of the open water season.

Monte Carlo simulation based thinning and calculating method for noninvasive blood glucose sensing by near-infrared spectroscopy

Previous results show that the floating reference theory(FRT)is an effective tool to reduce the infuence of interference factors on noninvasive blood glucose sensing by near infrared spectros-copy(NIRS).It is the key to measure the floating reference point(FRP)precisely for the application of FRT.Monte Carlo(MC)simulation has been introduced to quantitatively in-vestigate the effects of positioning errors and light source drifts on measuring FRP.In this article,thinning and calculating method(TCM)is proposed to quantify the positioning error.Mean-while,the normalization process(NP)is developed to significantly reduce the error induced by light source drift.The results according to TCM show that 7 purm deviations in positioning can generate about 10.63%relative error in FRP.It is more noticeable that 1%fluctuation in light source intensity may lead to 12.21%relative errors.Gratifyingly,the proposed NP model can effectively reduce the error caused by light source drift.Therefore,the measurement system for FRPs must meet that the positioning error is less than 7 purm,and the light source drift is kept within 1%.Furthermore,an improvement for measurement system is proposed in order to take advantage of the NP model.

The effect of forced oscillations on the kinetics of wave drift in an inhomogeneous plasma

The kinetics is analyzed of the drift of non-potential plasma waves in spatial positions and wavevectors due to plasma's spatial inhomogeneity. The analysis is based on highly informative kinetic scenarios of the drift of electromagnetic waves in a cold ionized plasma in the absence of a magnetic field(Erofeev 2015 Phys. Plasmas 22 092302) and the drift of long Langmuir waves in a cold magnetized plasma(Erofeev 2019 J. Plasma Phys. 85 905850104). It is shown that the traditional concept of the wave kinetic equation does not account for the effects of the forced plasma oscillations that are excited when the waves propagate in an inhomogeneous plasma.Terms are highlighted that account for these oscillations in the kinetic equations of the abovementioned highly informative wave drift scenarios.

Combined Electromagnetic and Drift Diffusion Models for Microwave Semiconductor Device

In this work, we present a numerical model to solve the drift diffusion equations coupled with electromagnetic model, where all simulations codes are implemented using MATLAB code software. As first, we present a one-dimensional (1-D) PIN diode structure simulation achieved by solving the drift diffusion model (DDM). Backward Euler algorithm is used for the discretization of the proposed model. The aim is to accomplish time-domain integration. Also, finite different method (FDM) is considered to achieve space-Domain mesh. We introduced an iterative scheme to solve the obtained matrix systems, which combines the Gummel’s iteration with an efficient direct numerical UMFPACK method. The obtained solutions of the proposed algorithm provide the time and space distribution of the unknown functions like electrostatic potential and carrier’s concentration for the PIN diode. As second case, the finite-difference time-domain (FDTD) technique is adopted to analyze the entire 3-D structure of the stripline circuit including the lumped element PIN diode. The microwave circuit is located in an unbounded medium, requiring absorbing boundaries to avoid nonphysical reflections. Active device results were presented and show a good agreement with other reference. Electromagnetic results are qualitatively in agreement with other results obtained using SILVACO-TCAD.

粉煤灰堆积体大断面连拱隧道失稳破坏模式及施工力学研究

国内目前并无在粉煤灰堆积体中采用暗挖法修建大断面连拱隧道的先例,给隧道的设计和施工带来了巨大挑战。依托盐坪坝隧道工程,开展了室内相似模型试验,探讨了该地层下大断面连拱隧道塌落拱演化规律及开挖施工力学行为。研究结果表明:连拱隧道在粉煤灰堆积体中开挖极易发生失稳,将失稳过程划分为4个演化阶段,破坏模式具有突发性、剧烈性、对称性、同步性特征;塌落拱高度约为隧道开挖跨度的40%。采用单侧壁导坑法开挖时粉煤灰堆积体承载力严重不足,根据相似换算关系左、右洞拱顶沉降分别为16.86,14.91cm,需对地层实施注浆加固;地层沉降呈对称“双峰状”,先行洞相较于后行洞受开挖影响较大,中隔墙的承载作用能够有效降低地层位移,单洞压力拱向双洞压力拱转换过程中,中隔墙受到偏压作用容易向先行洞侧偏转,建议在先行洞开挖时中隔墙添加支撑物。

厦门海域落水集装箱漂移路径的数值模拟

对落水集装箱的受力和进水淹没过程进行分析,细化落水集装箱风致漂移系数与淹没率的线性关系。并以“SITC PYEONGTAEK”轮落水集装箱为例,结合区域海洋模式模拟的海流数据,用拉格朗日法和蒙特卡洛随机扩散相结合的粒子追踪法建立集装箱追踪模型模拟集装箱的漂移轨迹,再利用网格化统计方案对零散分布的粒子进行概率统计,分析掉落集装箱的漂移位置。通过实际案例中集装箱落水后的漂移轨迹对模型进行验证,表明该模型模拟得到的轨迹点与海上发现的漂浮集装箱点位基本吻合,说明该模型有助于快速锁定漂移区域,减少由落水集装箱带来的风险,增加航行单位对落水集装箱漂移轨迹的预测准确性,为集装箱落水港区及周边水域的通航安全提供保障。

浅埋无中导洞连拱隧道先行洞衬砌承担围岩荷载确定方法

不同于常规隧道衬砌结构,无中导洞连拱隧道先行洞衬砌施作稳定后还要受到邻近后行洞开挖的二次扰动,后行洞开挖对先行洞衬砌受力的影响不可忽视。依托杨家寨公路隧道工程,基于考虑了后行洞开挖影响的隧道围岩压力计算方法,结合现场监测位移数据采用荷载-结构法进行位移反分析,提出了一种确定浅埋无中导洞连拱隧道先行洞衬砌承担围岩荷载的方法,最后和衬砌背后围岩荷载现场监测值进行对比验证。研究结果表明:后行洞开挖对先行洞衬砌承担围岩荷载的影响较为显著,监测点增长幅度均值高达142%;先行洞衬砌受力呈现“拱腰较小,拱顶及拱肩较大”的分布规律,且靠近后行洞一侧拱肩位置受力最大;通过现场监测位移反分析得到的先行洞衬砌受力结果和现场监测结果较为一致,论证了该方法的合理性和实用性,研究成果可为无中导洞连拱隧道结构设计提供参考。