Distinguishing the main components of submarine groundwater and estimating the corresponding fluxes based on radium tracing method—taking the Maowei Sea for example

Submarine groundwater discharge(SGD)is an important part in the land-sea interactions,which mainly contains three components:submarine fresh groundwater discharge(SFGD),tidal flat recirculated saline groundwater discharge(tidal flat RSGD)and subtidal recirculated saline groundwater discharge(subtidal RSGD).In order to make a more accurate assessment of the impact of SGD on coastal ecological environment,it is necessary to distinguish the main components of SGD.In this study,the Maowei Sea,located in the northern part of the Beibu Gulf,was selected as the study area.Based on the radium(Ra)tracing method,we present a new analytical method for distinguishing the three main components of SGD in this area combined with field data.The average daily flow along the coastline of the Maowei Sea of tidal flat RSGD was slightly higher than that of SFGD,and both two were on the magnitude of 1×10^(5)m^(3)/d.The average daily flow for the subtidal RSGD of the entire subtidal zone of the Maowei Sea reached to the magnitude of 1×10^(6)-1×10^(7)m^(3)/d.The long-term variation trend of terrestrial SGD is a valuable information for the study of the influence of terrigenous material on the coastal ecological environment.Based on the results of four sampling periods,it is found that the fluxes of SFGD and tidal flat RSGD in the Maowei Sea had good linear correlation with the net precipitation.As an example,January 2015 to August 2022 were selected as the study periods,and the variation trends of SFGD and tidal flat RSGD were calculated by linear function with net precipitation as the independent variable.The results showed that the flux of tidal flat RSGD was slightly higher than that of SFGD,and the difference between the two is larger in flood season while smaller in dry season.In general,in the coastal range of China,the total SGD flux in the Maowei Sea area is at a high level,and the SFGD flux is at a medium level.

Deviation of the Lagrangian particle tracing method in the evaluation of the Southern Hemisphere annual subduction rate

The classical Lagrangian particle tracing method is widely used in the evaluation of the ocean annual subduction rate.However,our analysis indicates that in addition to neglecting the effect of mixing,there are two possible deviations in the method:one is an overestimation due to not considering that the amount of subducted water at the source location may be inadequate during the late winter of the first year when the mixed layer becomes shallow;the other one is an underestimation due to the neglect of the effective subduction caused by strong vertical pumping.Quantitative analysis shows that these two deviations mainly exist in the low-latitude subduction areas of the South Pacific and South Atlantic.The two deviations have very similar distribution areas and can partially off set each other.However,the overall deviation is still large,and the maximum relative deviation ratio can reach 50%;therefore,it cannot be ignored.

An improved algorithm based on equivalent sound velocity profile method at large incident angle

With the development of ultra-wide coverage technology,multibeam echo-sounder(MBES)system has put forward higher requirements for localization accuracy and computational efficiency of ray tracing method.The classical equivalent sound speed profile(ESSP)method replaces the measured sound velocity profile(SVP)with a simple constant gradient SVP,reducing the computational workload of beam positioning.However,in deep-sea environment,the depth measurement error of this method rapidly increases from the central beam to the edge beam.By analyzing the positioning error of the ESSP method at edge beam,it is discovered that the positioning error increases monotonically with the incident angle,and the relationship between them could be expressed by polynomial function.Therefore,an error correction algorithm based on polynomial fitting is obtained.The simulation experiment conducted on an inclined seafloor shows that the proposed algorithm exhibits comparable efficiency to the original ESSP method,while significantly improving bathymetry accuracy by nearly eight times in the edge beam.

Low-complexity signal detection for massive MIMO systems via trace iterative method

Linear minimum mean square error(MMSE)detection has been shown to achieve near-optimal performance for massive multiple-input multiple-output(MIMO)systems but inevitably involves complicated matrix inversion,which entails high complexity.To avoid the exact matrix inversion,a considerable number of implicit and explicit approximate matrix inversion based detection methods is proposed.By combining the advantages of both the explicit and the implicit matrix inversion,this paper introduces a new low-complexity signal detection algorithm.Firstly,the relationship between implicit and explicit techniques is analyzed.Then,an enhanced Newton iteration method is introduced to realize an approximate MMSE detection for massive MIMO uplink systems.The proposed improved Newton iteration significantly reduces the complexity of conventional Newton iteration.However,its complexity is still high for higher iterations.Thus,it is applied only for first two iterations.For subsequent iterations,we propose a novel trace iterative method(TIM)based low-complexity algorithm,which has significantly lower complexity than higher Newton iterations.Convergence guarantees of the proposed detector are also provided.Numerical simulations verify that the proposed detector exhibits significant performance enhancement over recently reported iterative detectors and achieves close-to-MMSE performance while retaining the low-complexity advantage for systems with hundreds of antennas.

Indirect-approach method for specified-end points seismic ray tracing in three dimen-sional inhomogeneous media

Generally speaking, the factors of both medium and tectonic that give rise to heterogeneity of the earth crust and mantle structures should be taken into account simultaneously in three dimensional seismic ray tracing. In this paper, the three dimensional structure models are constructed with the model similar to generation system in computer aid design and manufacturing (CAD/CAM). Based on the algorithm proposed by Cerveny et al . for complete seismic ray tracing in complex three dimensional structures, a new technique called the indirect approach method for two point seismic ray tracing in three dimensional laterally heterogeneous media has been put forward, and its numerical computing examples were given.

Measurement and verification of concentration-dependent diffusion coefficient:Ray tracing imagery of diffusion process

Ray tracing method is used to study the propagation of collimated beams in a liquid-core cylindrical lens(LCL),which has dual functions of diffusion cell and image formation.The diffusion images on the focal plane of the used LCL are simulated by establishing and solving both linear and nonlinear ray equations,the calculated results indicate that the complex imaging results of LCL in inhomogeneous media can be treated by the law of ray propagation in homogeneous media under the condition of small refractive index gradient of diffusion solution.Guided by the calculation conditions,the diffusion process of triethylene glycol aqueous solution is experimentally studied at room temperature by using the LCL in this paper.The spatial and temporal concentration profile Ce(z,t)of diffusion solution is obtained by analyzing diffusion image appearing on the focal plane of the LCL;Then,the concentration-dependent diffusion coefficient is assumed to be a polynomial D(C)=D0×(1+α1C+α2C2+α3C3+…).The finite difference method is used to solve the Fick diffusion equation for calculating numerically the concentration profiles Cn(z,t).The D(C)of triethylene glycol aqueous solution is obtained by comparing the Cn(z,t)with Ce(z,t).Finally,the obtained polynomial D(C)is used to calculate the refractive index profiles nn(z,t)s of diffusion solution in the used LCL.Based on the ray propagation law in inhomogeneous media and the calculated n(z,t),the ray tracing method is used again to simulate the dynamic images of the whole experimental diffusion process to varify the correctness of the calculated D(C).The method presented in this work opens up a new way for both measuring and verifying the concentration-dependent liquid diffusion coefficients.

A FLUORESCENCE QUENCHING IMNUNOASSAY METHOD FOR DETERMINATION OF TRACE ALBUMIN USING UNLABELED EUROPIUM CHELATE

A fluoroimmunoassay method using unlabeled europium chelate is described.The principle is similar to that of fluoroimmunoassay method using lanthanide chelate as labels.The procedure is simple because labeling process is omitted.The detection limit is about 10^(10) mol/L antigen.The relative standard deviation of immunoassay is less than 10%.The recoveries of human serum albumin and estradiol protein conjugate are 96-105% and 111% respectively.

A QUENCHING FLUORESCENCE IMMUNOASSAY METHOD FOR DETERMINATION OF TRACE ALBUMIN USING UNLABELED TERBIUM CHELATE

A fluoroimmunoassay method using unlabeled Terbium chelate is described.The principle is similar to that of fluoroimmunoassay method using lanthanide chelate as labels.The procedure is simpte because labeling process is unnecessary.The recovery of HSA and albumin in urine is 107% and 95% respectively.The standard deviation is tess than 10%.

A simple and effective model for Jones matrix to evaluate the transmission performance of Savart polariscope

Starting from the wave normal tracing treatment in the Savart polariscope that relates the Jones matrix to its transmission performance, this paper establishes a simple and effective model for the Jones matrix at an arbitrary incidence in the spatial domain. Analytical expressions of all the matrix components are determined with the consideration of all the main impact factors. This model needs only a few parameters hence it is convenient to be employed to evaluate the propagation performance of any birefringent optical system. The simulated results obtained with it demonstrate that this model gives a precise representation of the characteristic of light propagation in the Savart polariscope. This would provide useful suggestions for the design, calibration, and performance improvement of any other birefringent polarisation element and optical system.

Relation between viscoplastic flow, periodic feature and Joule heat validity for electrically assisted friction stir welding

Electrically assisted friction stir welding is a promising way to improve liquidity of viscoplastic metal. However, this local heat source only can take into effect in a certain contact interface status. Semi-analytical model and FEM are adopted to analyze materials flow in electrically assisted FSW. Semi-analysis results show that Joule heating is validity with rotation speed 1 400 r/min and travel speed 400 mm/min. The maximum flow velocity is 0.033 m/s when travel speed is 400 mm/min, and it increases to 0.043 m/s with current assistant. Visco-plasticizing efficiency is a key factor to decide whether Joule heat is significant, and viscoplastic flow can periodically make circuit short. Periodic torque, axial force and onion rings surface show dynamic flow process, which can be used to illustrate the dynamic contact resistance feature and correct the calculation equation of Joule heat. This paper is aim to point out the potential relation between viscoplastic flow, periodic torque and Joule heating validity.

Multi-point Contact of the High-speed Vehicle-turnout System Dynamics

The wheel-rail contact problems, such as the number, location and the track of contact patches, are very important for optimizing the spatial structure of the rails and lowering the vehicle-turnout system dynamics. However, the above problems are not well solved currently because of having the difficulties in how to determine the multi-contact, to preciously present the changeable profiles of the rails and to establish an accurate spatial turnout system dynamics model. Based on a high-speed vehicle-turnout coupled model in which the track is modeled as flexible with rails and sleepers represented by beams, the line tracing extreme point method is introduced to investigate the wheel-rail multiple contact conditions and the key sections of the blade rail, longer nose rail, shorter rail in the switch and nose rail area are discretized to represent the varying profiles of rails in the turnout. The dynamic interaction between the vehicle and turnout is simulated for cases of the vehicle divergently passing the turnout and the multi-point contact is obtained. The tracks of the contact patches on the top of the rails are presented and the wheel-rail impact forces are offered in comparison with the contact patches transference on the rails. The numerical simulation results indicate that the length of two-point contact occurrence of a worn wheel profile and rails is longer than that of the new wheel profile and rails; The two-point contact definitely occurs in the switch and crossing area. Generally, three-point contact doesn’t occur for the new rail profile, which is testified by the wheel-rails interpolation distance and the first order derivative function of the tracing line extreme points. The presented research is not only helpful to optimize the structure of the turnout, but also useful to lower the dynamics of the high speed vehicle-turnout system.

An efficient adaptive space partitioning algorithm for electromagnetic scattering calculation of complex 3D models

The space partitioning algorithm based on the rounding and addressing operations has been proved to be an efficient space partitioning algorithm with the potential for real-time calculation.An improvement on this kind of space partitioning algorithms for solving complex 3D models is presented.Numerical examples show that the efficiency of the improved algorithm is better than that of the original method.When the size of most target elements is smaller than the size of spatial grids,the efficiency of the improved method can be more than four times of that of the original method.An adaptive method of space partitioning based on the improved algorithm is developed by taking the surface element density or the curvature as the threshold for deep partitioning and conducting the deep partitioning using the octree method.A computer program implementation for applying the method in some typical applications is discussed,and the performance in terms of the efficiency,reliability,and resource use is evaluated.Application testing shows that the results of the adaptive spacing partitioning are more convenient for the follow-up use than that of the basic uniform space partitioning.Furthermore,when it is used to calculate the electromagnetic scattering of complex targets by the ray tracing(RT)method,the adaptive space partitioning algorithm can reduce the calculation time of the RT process by more than 40%compared with the uniform space segmentation algorithm.

Thermal vacuum state for the two-coupled-oscillator model at finite temperature: Derivation and application

Following the spirit of thermo field dynamics initiated by Takahashi and Umezawa, we employ the technique of integration within an ordered product of operators to derive the thermal vacuum state （TVS） for the Hamiltonian H of the two-coupled-oscillator model. The ensemble averages of the system are derived conveniently by using the TVS. In addition, the entropy for this system is discussed based on the relation between the generalized Hellmann-Feynman theorem and the entroy variation in the context of the TVS.

A Tracing Algorithm for Surface—Surface Intersections on Surface Boundaries

In this paper we present an algorithm with a new trace-terminating conditionfor tracing along surface-surface intersection curves on surface boundaries, while several tracingmethods and embedding methods that include tracing scheme may cause false termination with atraditional trace-terminating condition: tracing stops when the surface-domain's boundary isreached. And we also suggest another iterative method to calculate intersection points on surfaceboundaries with parallel surface normal. Some numerical examples with these two ideas andcomparisons to 'DESIGNBASE', 'ACIS', and 'Parasolid' are included to demonstrate the effectivenessof our algorithm.

TRACE AVERAGING DOMAIN DECOMPOSITION METHOD WITH NONCONFORMING FINITE ELEMENTS

We consider, in this paper, the trace averaging domain decomposition method for the second order self-adjoint elliptic problems discretized by a class of nonconforming finite elements, which is only continuous at the nodes of the quasi-uniform mesh. We show its geometric convergence and present the dependence of the convergence factor on the relaxation factor, the subdomain diameter H and the mesh parameter h. In essence;, this method is equivalent to the simple iterative method for the preconditioned capacitance equation. The preconditioner implied in this iteration is easily invertible and can be applied to preconditioning the capacitance matrix with the condition number no more than O((1 + In H/h)max(1 + H-2, 1 + In H/h)).