Phase-Sensitive Optical Time-Domain Reflectometer Based on a 120°-Phase-Difference Michelson Interferometer

A phase-sensitive optical time domain reflectometer （φ-OTDR） based on a 120°-phase-difference Michelson in- terferometer is proposed. The Michelson interferometer with arm difference of 4m is used to test the phase difference between the Rayleigh scattering from two sections of the fiber. A new demodulation method called the inverse transmission matrix demodulation scheme is utilized to demodulate the distributed phase from the backward scattering along the long fiber, The experimental results show that the 120°-phase-difference inter- ferometer φ-OTDR can detect the phase along the 3km fiber, and the acoustic signal within the whole human hearing range of 20 Hz-20 kHz is reproduced accurately and quickly.

Measurement of Turbulence Propagation Velocity Using Doppler Reflectometer on HL-2A Tokamak

O-mode Doppler reflectometer has been successfully developed as an important diagnostic system on HL-2A. It can be used to measure the turbulence propagation in both plasma edge and confinement zone. The Doppler reflectometer system consists of two fixed frequency homodyne receivers： 15 GHz （corresponding to cutoff density of 0.3×10^19 m^-3） and 33 GHz （corresponding to cutoff density of 1.35× 10^19 m^-3）. The Doppler reflectometry principle and the experimental arrangements on HL-2A are presented. Meanwhile, the experimental Doppler reflectometric spectra under different discharge conditions, with and without ECRH, were obtained. Furthermore, the turbulence propagation velocity change and the profile were also observed in different discharge conditions.

Reconstruction of the Density Profile for the EAST Tokamak Based on Polarimeter/Interferometer and Microwave Reflectometer Systems

A plasma density profile reconstruction procedure based on the Park matrix method has been developed for both circular and elongated plasma configuration on the Experimental Advanced Superconducting Tokamak（EAST）.This method incorporates the line integrated electron density measured by the HCN interferometer and polarimeter/interferometer（POINT） system,the equilibrium fit（EFIT） based on magnetic measurements and the edge electron density profile provided by the microwave reflectometer.It is shown that when the magnetic flux surfaces are slightly corrected,the fitting error is less than 5% in comparison with the measurement data.

A NOVEL METHOD FOR MEASURING COMPLEX REFLECTION COEFFICIENT USING A FOUR-PORT REFLECTOMETER

A novel method for precise measurement of complex reflection coeffcient using a four-port reflectometer is presented. First, three new complex system constants are introduced,which depend only on the scattering parameters of the four-port reflectometer. Therefore, the stability of the reflectometer is greatly improved. Then, these complex system constants are used to determine the complex reflection coeffcient F of the device under test by calibrating the reflectometer. Finally, a four-port reflectometer comprising a magic tee and a power detector is constructed and excellent experimental results are obtained.

Design and simulation of time-of-flight neutron reflectometer

In this paper,a time of flight neutron reflectometer model is proposed.According to the original objective inputting physical parameters and expectations,computer simulation was done by McStas package.With the open slit, neutron flux at sample position was about 1×10~5 n·cm^(-2)·S^(-1),and the optimal resolution could reach 0.5%.The results show that the model has a good performance.The design of this model can be an instruction for neutron reflectometer construction or project analysis in the future.

Preliminary Study of the Magnetic Perturbation Effects on the Edge Density Profiles and Fluctuations Using Reflectometers on EAST

The resonant magnetic perturbation （RMP） coils have been successfully designed and installed on the Experimental Advanced Superconducting Tokamak （EAST）. Using the reflectometer systems, the density profile and the density fluctuations during magnetic perturbations （MPs） phase have been investigated. During the experiments, two different cases are studied separately： steady MPs and rotating MPs. In both cases, a strongly density pump-out has been observed. In the steady MPs cases, an enhancement of the low frequency （〈60 kHz） density fluctuations in H-mode phase has been observed. The plasma density boundary out-shifts 5~ caused by the MPs. The pedestal density gradient is reduced by 50%, while the radial location nearly stays unchanged. In the rotating MPs, the line-averaged density, the D~ emission at the divertor region and the spectrum of the density fluctuations are modulated. The results suggest that the low frequency （〈60 kHz） density fluctuations may contribute to the strong density pump-out.

Plasma Boundary Position Measurements Using Frequency Sweep Microwave Reflectometer in LHD

A broadband frequency tunable microwave reflectometer system, which has the abil- ity of fast and stable frequency sweeping operation, is applied in the large helical device （LHD） to measure the boundary of the high-temperature plasma. When the microwave is launched, with extraordinary polarization, from the low field side of the magnetic field with its frequency swept from low to high, the microwave of minimum right-hand cut-off frequency is reflected most out- side. We can then estimate the plasma boundary position by measuring the change of the reflected power from the cut-off layer.

Design of a variable frequency comb reflectometer system for the ASDEX Upgrade tokamak

Comb reflectometers offer the advantage of measuring several radial positions in plasma simultaneously.This allows for the investigation of fast timescales during L-H transitions,I-phases,I-mode bursts,transients during heat wave propagation,etc.A drawback of many present-day systems is that they use a fixed frequency difference between the probing frequencies.Hence,although the central probing frequency can be varied,the probing frequency difference is usually fixed.The new design presented in this work uses an advanced microwave generation and detection scheme,which allows for arbitrary probing frequencies and probing frequency separations.

Comparison of methods for turbulence Doppler frequency shift calculation in Doppler reflectometer

The Doppler reflectometer(DR),a powerful diagnostic for the plasma perpendicular velocity(u⊥)and turbulence measurement,has been widely used in various fusion devices.Many efforts have been put into extracting the Doppler shift from the DR signal.There are several commonly used methods for Doppler shift extraction,such as the phase derivative,the center of gravity,and symmetric fitting(SFIT).However,the strong zero-order reflection component around 0 kHz may interfere with the calculation of the Doppler shift.To avoid the influence of the zerofrequency peak,the asymmetric fitting(AFIT)method was designed to calculate the Doppler shift.Nevertheless,the AFIT method may lead to an unacceptable error when the Doppler shift is relatively small compared to the half width at half maximum(HWHM).Therefore,an improved method,which can remove the zero-frequency peak and fit the remaining Doppler peak with a Gaussian function,is devised to extract the Doppler shift.This method can still work reliably whether the HWHM is larger than the Doppler shift or not.

Time-domain dynamic constitutive model suitable for mucky soil site seismic response

Soil nonlinear behavior displays noticeable effects on the site seismic response.This study proposes a new functional expression of the skeleton curve to replace the hyperbolic skeleton curve.By integrating shear modulus and combining the dynamic skeleton curve and the damping degradation coefficient,the constitutive equation of the logarithmic dynamic skeleton can be obtained,which considers the damping effect in a soil dynamics problem.Based on the finite difference method and the multi-transmitting boundary condition,a 1D site seismic response analysis program called Soilresp1D has been developed herein and used to analyze the time-domain seismic response in three types of sites.At the same time,this study also provides numerical simulation results based on the hyperbolic constitutive model and the equivalent linear method.The results verify the rationality of the new soil dynamic constitutive model.It can analyze the mucky soil site nonlinear seismic response,reflecting the deformation characteristics and damping effect of the silty soil.The hysteresis loop area is more extensive,and the residual strain is evident.

Time-Domain Analysis of Body Freedom Flutter Based on 6DOF Equation

The reduced weight and improved efficiency of modern aeronautical structures result in a decreasing separation of frequency ranges of rigid and elastic modes.Particularly,a high-aspect-ratio flexible flying wing is prone to body freedomflutter(BFF),which is a result of coupling of the rigid body short-periodmodewith 1st wing bendingmode.Accurate prediction of the BFF characteristics is helpful to reflect the attitude changes of the vehicle intuitively and design the active flutter suppression control law.Instead of using the rigid body mode,this work simulates the rigid bodymotion of the model by using the six-degree-of-freedom(6DOF)equation.A dynamicmesh generation strategy particularly suitable for BFF simulation of free flying aircraft is developed.An accurate Computational Fluid Dynamics/Computational Structural Dynamics/six-degree-of-freedom equation(CFD/CSD/6DOF)-based BFF prediction method is proposed.Firstly,the time-domain CFD/CSD method is used to calculate the static equilibrium state of the model.Based on this state,the CFD/CSD/6DOF equation is solved in time domain to evaluate the structural response of themodel.Then combinedwith the variable stiffnessmethod,the critical flutter point of the model is obtained.This method is applied to the BFF calculation of a flyingwing model.The calculation results of the BFF characteristics of the model agree well with those fromthe modalmethod andNastran software.Finally,the method is used to analyze the influence factors of BFF.The analysis results show that the flutter speed can be improved by either releasing plunge constraint or moving the center ofmass forward or increasing the pitch inertia.

PC-based artif icial neural network inversion for airborne time-domain electromagnetic data

Traditionally, airborne time-domain electromagnetic （ATEM） data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components （PCs） to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network （ANN） to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody＇s methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.

Solution and Analysis of Chatter Stability for End Milling in the Time-domain

In this paper, the instantaneous undeformed chip thickness is modeled to include the dynamic modulation caused by the tool vibration while the dynamic regenerative effects are taken into account. The numerical method is used to solve the differential equations goveming the dynamics of the milling system. Several chatter detection criteria are applied synthetically to the simulated signals and the stability diagram is obtained in time-domain. The simulation results in time-domain show a good agreement with the analytical prediction, which is validated by the cutting experiments. By simulating the chatter stability lobes in the time-domain and analyzing the influences of different spindle speeds on the vibration amplitudes of the tool under a Fixed chip-load condition, conclusions could be drawn as follows： In rough milling, higher machining efficiency can be achieved by selecting a spindle speed corresponding to the axial depth of cut in accordance with the simulated chatter stability lobes, and in Fmish milling, lower surface roughness can be achieved by selecting a spindle speed well beyond the resonant frequency of machining system.

Comparison of the time-domain electromagnetic field from an infinitesimal point charge and dipole source

An electromagnetic field is generated through the accelerating movement of two equal but opposite charges of a single dipole. An electromagnetic field can also be generated by a time-varying infinitesimal point charge. In this study, a comparison between the electromagnetic fields of an infinitesimal point charge and a dipole has been presented. First, the time-domain potential function of a point source in a 3D conductive medium is derived. Then the electric and magnetic fields in a 3D homogeneous lossless space are derived via the relation between the potential and field. The field differences between the infinitesimal point charge and the dipole in the step-off time, far-source, and near-source zones are analyzed, and the accuracy of the solutions from these sources is investigated. It is also shown that the field of the infinitesimal point charge in the near-source zone is different from that of the dipole, whereas the far-source zone fields of these two sources are identical. The comparison of real and simulated data shows that the infinitesimal point charge represents the real source better than the divole source.

Numerical modeling of the 2D time-domain transient electromagnetic secondary field of the line source of the current excitation

To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force （Emf）, and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.

Three-dimensional arbitrarily anisotropic modeling for time-domain airborne electromagnetic surveys

Electrically anisotropic strata are abundant in nature, so their study can help our data interpretation and our understanding of the processes of geodynamics. However, current data processing generally assumes isotropic conditions when surveying anisotropic structures, which may cause discrepancies between reality and electromagnetic data interpretation. Moreover, the anisotropic interpretation of the time-domain airborne electromagnetic （TDAEM） method is still confined to one dimensional （1D） cases, and the corresponding three-dimensional （3D） numerical simulations are still in development. In this study, we expanded the 3D TDAEM modeling of arbitrarily anisotropic media. First, through coordinate rotation of isotropic conductivity, we obtained the conductivity tensor of an arbitrary anisotropic rock. Next, we incorporated this into Maxwell＇s equations, using a regular hexahedral grid of vector finite elements to subdivide the solution area. A direct solver software package provided the solution for the sparse linear equations that resulted. Analytical solutions were used to verify the accuracy and feasibility of the algorithm. The proven model was then applied to analyze the effects of arbitrary anisotropy in 3D TDAEM via the distribution of responses and amplitude changes, which revealed that different anisotropy situations strongly affected the responses of TDAEM.

A 1D time-domain method for in-plane wave motions in a layered half-space

A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the layered half-space is discretized on the basis of the propagation characteristic of elastic wave according to the Snell law. Then, the finite element method with lumped mass and the central difference method are incorporated to establish 2D wave motion equations, which can be transformed into 1D equations by discretization principle and explicit finite element method. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained, and the wave motions in layered half-space are finally determined based on the characteristic of traveling wave. Both the theoretical analysis and the numerical results demonstrate that the proposed method has high accuracy and good stability.

Investigation of copper sulfate pentahydrate dehydration by terahertz time-domain spectroscopy

Copper sulfate pentahydrate is investigated by terahertz time-domain spectroscopy. It is shown that the terahertz absorption coefficients are correlated with the particle size of the samples, as well as the heating rates of the ambient temperature. Furthermore, the water molecules of copper sulfate pentahydrate can be quantitatively characterized due to the high sensitivity of the terahertz wave to water molecules. Based on such results, the status of water incorporated in mineral opal is also characterized using terahertz time-domain spectroscopy. It indicates that terahertz technology can be considered as an efficient method to detect the dehydration of minerals.

Second-Order Wave Diffraction Around 3-D Bodies by A Time-Domain Method

A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.

Investigating path loss characteristics of UWB signals in vacancy indoor environment based on time-domain technique

The path loss analysis model based on 5 rays in vacancy indoor environment is proposed. The relationship between multipath overlapping and the path loss is analyzed mathematically. Time-domain technique is introduced to compute reflection coefficient in a very short time interval. A 5 rays path loss calculation method, which is satisfactory accurate, is developed. 5 typical environments are involved to analyze and generalize the common path loss characteristics in vacancy indoor environment. The simulation result shows that the path loss can be characterized as 3 zones with different path loss exponent as distance between transmitter and receiver increasing.