The lack of soft magnetic composites with high power density in MHz frequency range has become an obstacle in the efficient operation of the electrical and electronic equipments.Here,a promising method to increase the...The lack of soft magnetic composites with high power density in MHz frequency range has become an obstacle in the efficient operation of the electrical and electronic equipments.Here,a promising method to increase the cut-off frequency of iron-based soft magnetic composites to hundreds of MHz is reported.The cut-off frequency is increased from 10 MHz to 1 GHz by modulating the height of the ring,the distribution of particles,and the particle size.The mechanism of cut-off frequency and permeability is the coherent rotation of domain modulated by inhomogeneous field due to the eddy current effect.An empirical formula for the cut-off frequency in a magnetic ring composed of iron-based particles is established from experimental data.This work provides an effective approach to fabricate soft magnetic composites with a cut-off frequency in hundreds of MHz.展开更多
BACKGROUND Glycated hemoglobin A1c(HbA1c)is considered the most suitable for diabetes mellitus diagnosis due to its accuracy and convenience.However,the effect of HbA1c on diabetic retinopathy(DR)in the Han and Korean...BACKGROUND Glycated hemoglobin A1c(HbA1c)is considered the most suitable for diabetes mellitus diagnosis due to its accuracy and convenience.However,the effect of HbA1c on diabetic retinopathy(DR)in the Han and Korean populations in Jilin,China,remains inconclusive.AIM To determine the best cut-off of HbA1c for diagnosing DR among the Chinese.METHODS This cross-sectional study included 1933 participants from the Yanbian area of Jilin Province,China.Trained investigators employed a questionnaire-based survey,physical examination,laboratory tests,and fundus photography for the investigation.The best cut-off value for HbA1c was established via the receiver operating characteristic curve.The factors associated with HbA1c-associated risk factors were determined via linear regression.RESULTS The analysis included 887 eligible Chinese Han and Korean participants,591 of whom were assigned randomly to the training set and 296 to the validation set.The prevalence of DR was 3.27% in the total population.HbA1c of 6.2% was the best cut-off value in the training set,while it was 5.9% in the validation set.In both Chinese Han and Korean populations,an HbA1c level of 6.2% was the best cut-off value.The optimal cut-off values of fasting blood glucose(FBG)≥7 mmol/L and<7 mmol/L were 8.1% and 6.2% respectively in Han populations,while those in Korean populations were 6.9%and 5.3%,respectively.Age,body mass index,and FBG were determined as the risk factors impacting HbA1c levels.CONCLUSION HbA1c may serve as a useful diagnostic indicator for DR.An HbA1c level of 6.2% may be an appropriate cut-off value for DR detection in the Chinese population.展开更多
Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumb...Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.展开更多
It is difficult to quantitatively detect defects by using the time domain or frequency domain features of Lamb wave signals due to their dispersion and multimodal characteristics.Therefore,it is important to discover ...It is difficult to quantitatively detect defects by using the time domain or frequency domain features of Lamb wave signals due to their dispersion and multimodal characteristics.Therefore,it is important to discover an intrinsical parameter of Lamb waves that could be used as a damage sensitive feature.In this paper,quantitative defect detection in aluminium plates is carried out by means of wavenumber analysis approach.The wavenumber of excited Lamb wave mode is a fixed value,given a frequency,a thickness and material properties of the target plate.When Lamb waves propagate to the structural discontinuity,new wavenumber components are created by abrupt wavefield change.The new wavenumber components can be identified in the frequency-wavenumber domain.To estimate spatially dependent wavenumber values,a short-space two-dimensional Fourier transform(FT)method is presented for processing wavefield data of Lamb waves.The results can be used to determine the location,size and depth of rectangular notch.The analysis techniques are demonstrated using simulation examples of an aluminium plate with a rectangular notch.Then,the wavenumber analysis method is applied to simulation data that are obtained through a range of notch depths and widths.The results are analyzed and rules of the technique with regards to estimating notch depth are determined.Based on simulation results,guidelines for using the technique are developed.Finally,experimental wavefield data are obtained in aluminium plates with rectangular notches by a full noncontact transceiving method,i.e.,laser-laser method.Band-pass filtering combined with continuous wavelet transform is used to extract a certain frequency component from the full laser-induced wavefield with wide band.Shortspace two-dimensional FT method is used for further processing full wavefield data at a certain frequency to estimate spatially dependent wavenumber values.The consistency of simulation and experimental results shows the effectiveness of proposed wavenumber method for quantitative rectangular notch detection.展开更多
An exact solution based on the wavenumber integration method is proposed and implemented in a numerical model for the acoustic field in a Pekeris waveguide excited by either a point source in cylindrical geometry or a...An exact solution based on the wavenumber integration method is proposed and implemented in a numerical model for the acoustic field in a Pekeris waveguide excited by either a point source in cylindrical geometry or a line source in plane geometry. Besides, an unconditionally stable numerical solution is also presented, which entirely resolves the stability problem in previous methods. Generally the branch line integral contributes to the total field only at short ranges, and hence is usually ignored in traditional normal mode models. However, for the special case where a mode lies near the branch cut, the branch line integral can contribute to the total field significantly at all ranges. The wavenumber integration method is well-suited for such problems. Numerical results are also provided, which show that the present model can serve as a benchmark for sound propagation in a Pekeris waveguide.展开更多
The current local wavenumber methods for the interpretation of magnetic anomalies compute the locations of geological bodies by solving complex matrices. Presently, such methods require to know the structural index, w...The current local wavenumber methods for the interpretation of magnetic anomalies compute the locations of geological bodies by solving complex matrices. Presently, such methods require to know the structural index, which is a parameter that represents the source type. The structural index is hard to know in real data; consequently, the precision of current methods is low. We present the fast local wavenumber (FLW) method, and define the squared sum of the horizontal and vertical local wavenumbers as the cumulative local wavenumber. The FLW method is the linear combination of the umulative local wavenumberand other wavenumbers, and is used to compute the locations and structural index of the source without a priori information and matrix solution. We apply the FLW method to synthetic magnetic anomalies, and the results suggest that the FLW method is insensitive to background and oblique magnetization. Next, we apply the FLW method to real magnetic data to obtain the location and structural index of the source.展开更多
The wavenumber spectral components WN4 at the mesosphere and low thermosphere(MLT)altitudes(70–10 km)and in the latitude range between±45°are obtained from temperature data(T)observed by the Sounding of the...The wavenumber spectral components WN4 at the mesosphere and low thermosphere(MLT)altitudes(70–10 km)and in the latitude range between±45°are obtained from temperature data(T)observed by the Sounding of the Atmosphere using Broadband Emission Radiometry(SABER)instruments on board the National Aeronautics and Space Administration(NASA)’s Thermosphere–Ionosphere–Mesosphere Energetics and Dynamics(TIMED)spacecraft during the 11-year solar period from 2002 to 2012.We analyze in detail these spectral components WNk and obtain the main properties of their vertical profiles and global structures.We report that all of the wavenumber spectral components WNk occur mainly around 100 km altitude,and that the most prominent component is the wavenumber spectral component WN4 structure.Comparing these long duration temperature data with results of previous investigations,we have found that the yearly variation of spectral component WN4 is similar to that of the eastward propagating non-migrating diurnal tide with zonal wavenumber 3(DE3)at the low latitudes,and to that of the semi-diurnal tide with zonal wavenumber 2(SE2)at the mid-latitudes:the amplitudes of the A4 are larger during boreal summer and autumn at the low-latitudes;at the mid-latitudes the amplitudes have a weak peak in March.In addition,the amplitudes of component WN4 undergo a remarkable short period variation:significant day-to-day variation of the spectral amplitudes A4 occurs primarily in July and September at the low-latitudes.In summary,we conclude that the non-migrating tides DE3 and SE2 are likely to be the origins,at the low-latitudes and the mid-latitudes in the MLT region,respectively,of the observed wavenumber spectral component WN4.展开更多
In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over ...In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over large-scale and complex terrain requires additional methods. To this end, we have proposed a new topography-capable By performing a two-dimensional Fourier transform in the horizontal directions, threedimensional partial differential equations in the spatial domain were transformed into a group of independent, one-dimensional differential equations engaged with different wave numbers. These independent differential equations are highly parallel across different wave numbers. differential equations with different wave numbers, and the efficiency of solving fixedbandwidth linear equations was further improved by a chasing method. In a synthetic test, a prism model was used to verify the accuracy and reliability of the proposed algorithm by comparing the numerical solution with the analytical solution. We studied the computational precision and efficiency with and without topography using different Fourier transform methods. The results showed that the Guass-FFT method has higher numerical precision, while the standard FFT method is superior, in terms of computation time, for inversion and quantitative interpretation under complicated terrain.展开更多
Capillary and capillary-gravity waves possess a random character, and the slope wavenumber spectra of them can be used to represent mean distributions of wave energy with respect to spatial scale of variability. But s...Capillary and capillary-gravity waves possess a random character, and the slope wavenumber spectra of them can be used to represent mean distributions of wave energy with respect to spatial scale of variability. But simple and practical models of the slope wavenumber spectra have not been put forward so far. In this article, we address the accurate definition of the slope wavenumber spectra of water surface capillary and capillary-gravity waves. By combining the existing slope wavenumber models and using the dispersion relation of water surface waves, we derive the slope wavenumber spectrum models of capillary and capillary-gravity waves. Simultaneously, by using the slope wavenumber models, the dependence of the slope wavenumber spectrum on wind speed is analyzed using data obtained in an experiment which was performed in a laboratory wind wave tank. Generally speaking, the slope wavenumber spectra are influenced profoundly by the wind speed above water surface. The slope wavenumber spectrum increases with wind speed obviously and do not cross each other for different wind speeds. But, for the same wind speed, the slope wavenumber spectra are essentially identical, even though the capillary and capillary-gravity waves are excited at different times and locations. Furthermore, the slope wavenumber spectra obtained from the models agree quite well with experimental results as regards both the values and the shape of the curve.展开更多
It is complicated to model the acoustic field in stratified ocean for airborne aircraft,due to high speed of the source and air-to-water sound transmission.To our knowledge,there are very few papers in the open litera...It is complicated to model the acoustic field in stratified ocean for airborne aircraft,due to high speed of the source and air-to-water sound transmission.To our knowledge,there are very few papers in the open literature dealing with this complicated problem;but,in our opinion,they all require great amount of computation.We now propose a different method that requires much less computation.We improve the wavenumber integration method to model the received temporal signal for a moving source in stratified ocean and sum up in a concise form the core of our paper as follows:(A) Eq.(11) can be calculated by means of fast Chirp Z transform and the signals at all N time points are generated simultaneously;(B) direct numerical evaluation of the wavenumber integral in Eq.(4) produces large numerical errors;so it is necessary to shift the integration slightly below the real axis;(C) we compare the computation cost of direct calculation method with that of our fast calculation method;from the results presented in table 1,we can see that the fast calculation method consumes much less computation time,particularly for long duration signals;(D) for an airborne rapidly moving source,we compute the Doppler-shifted signals in shallow water and analyze their short-time Fourier transform;from Fig.1b,we can see that the received signals have multiple frequency components for a tonal source due to source motion and that each component corresponds to an arrival path.展开更多
Northeast Asian cut-off lows are crucial cyclonic systems that can bring temperature and precipitation extremes over large areas.Skillful subseasonal forecasting of Northeast Asian cut-off lows is of great importance....Northeast Asian cut-off lows are crucial cyclonic systems that can bring temperature and precipitation extremes over large areas.Skillful subseasonal forecasting of Northeast Asian cut-off lows is of great importance.Using two dynamical forecasting systems,one from the Beijing Climate Center(BCC-CSM2-HR)and the other from the Met Office(GloSea5),this study assesses simulation ability and subseasonal prediction skill for early-summer Northeast Asian cut-off lows.Both models are shown to have good ability in representing the spatial structure of cut-off lows,but they underestimate the intensity.The skillful prediction time scales for cut-off low intensity are about 10.2 days for BCC-CSM2-HR and 11.4 days for GloSea5 in advance.Further examination shows that both models can essentially capture the initial Rossby wave train,rapid growth and decay processes responsible for the evolution of cut-off lows,but the models show weaker amplitudes for the three-stage processes.The underestimated simulated strength of both the Eurasian midlatitude and East Asian subtropical jets may lead to the weaker local eddy-mean flow interaction responsible for the cut-off low evolution.展开更多
The prediction on small disturbance propagation in complex three-dimensional(3D) boundary layers is of great significance in transition prediction methodology, especially in the aircraft design. In this paper, the lin...The prediction on small disturbance propagation in complex three-dimensional(3D) boundary layers is of great significance in transition prediction methodology, especially in the aircraft design. In this paper, the linear stability theory(LST) with the equivalent spanwise wavenumber correction(ESWC) is proposed in order to accurately predict the linear evolution of a disturbance in a kind of boundary layer flow with a vital variation in the spanwise direction. The LST with the ESWC takes not only the scale of the mean flow with the significant variation but also the wavenumber evolution of the disturbance itself. Compared with the conventional LST, the results obtained by the new method are in excellent agreement with those of the numerical simulations. The LST with the ESWC is an effective method on the prediction of the disturbance evolution in 3D boundary layers, which improves the prediction of the LST in the applications to complex 3D boundary layers greatly.展开更多
Large-scale,fine,and efficient numerical simulation of a geothermal field plays an important role in geothermal energy development.Confronted with the problem of large computation and high storage requirements for com...Large-scale,fine,and efficient numerical simulation of a geothermal field plays an important role in geothermal energy development.Confronted with the problem of large computation and high storage requirements for complex underground models in a three-dimensional(3-D)numerical simulation of a geothermal fi eld,a mixed space-wavenumber domain 3-D numerical simulation algorithm is proposed in this paper.According to the superposition principle of temperature field,the geothermal field is decomposed into background and abnormal temperature fi elds for calculation.The uniform layered model is used to solve the background field.When the abnormal field is solved,the horizontal two-dimensional(2-D)Fourier transform is used to transform the 3-D diff erential equation satisfi ed by an abnormal field into a series of one-dimensional ordinary differential equations with diff erent wavenumbers,which greatly reduces the calculation and storage.The unit division of an ordinary diff erential equation is fl exible,and the calculation amount is small.The algorithm fully takes advantage of the effi ciency of the Fourier transform and the quickness of the catch-up method to solve linear equations with a fixed bandwidth,which effectively improves the computational efficiency.Compared with the COMSOL Multiphysics professional simulation finite element software,the time consumption and memory requirements of the algorithm proposed in this paper are reduced by multiple orders of magnitude in terms of ensuring accuracy and the same mesh division.The more the number of calculated nodes is,the more obvious is the advantage.We design models to study the thermal conductivity,heat fl ux boundary,regional tectonic morphology,and topographic relief of the geothermal fi eld distribution.A 3-D geophysical model is developed based on topographic elevation data,geothermal geology,and geophysical exploration data in the Qiabuqia area of Gonghe Basin,Qinghai Province,China.Numerical simulation of the geothermal fi eld in this area is realized,which shows that the algorithm is suitable for precise and effi cient simulation of an arbitrary complex terrain and geological conditions.展开更多
Marangoni Benard convection, which is mainly driven by the thermocapillary (Marangoni) effect, occurs in a thin liquid layer heated uniformly from the bottom. The wavenumber of supercritical convection is studied ex...Marangoni Benard convection, which is mainly driven by the thermocapillary (Marangoni) effect, occurs in a thin liquid layer heated uniformly from the bottom. The wavenumber of supercritical convection is studied experimentally in a 160×160-mm^2 cavity that & heated from the bottom block. The convection pattern & visualized by an infrared thermography camera. It is shown that the onset of the Benard cell is consistent with theoretical analysis. The wavenumber decreases obviously with increasing temperature, except for a slight increase near the onset. The wavenumber gradually approaches the minimum when the supercritical number e is larger than 10. Finally, a formula is devised to describe the wavenumber selection in supercritical convection.展开更多
Laplace's tidal equations are of great importance in various fields of geophysics. Here, the special case of zonal symmetry (zonal wavenumber m = 0) is investigated, where degenerate sets of eigensolutions appear....Laplace's tidal equations are of great importance in various fields of geophysics. Here, the special case of zonal symmetry (zonal wavenumber m = 0) is investigated, where degenerate sets of eigensolutions appear. New results are presented for the inclusion of dissipative processes and the case of unstable conditions. In both instances the (nonzero) eigenfrequencies are complex. In the latter case, additional stable (i.e. real) eigenfrequencies appear in the numerical results for the absolute value of the Lambparameter ε being larger than a critical value εc. Further, it is shown that any degeneracies are removed through the inclusion of dissipation. Moreover, asymptotic relations are derived employing the relation of Laplace's tidal equations for m = 0 to the spheroidal differential equation. The implications of these findings to numerical techniques are demonstrated and results of computations are presented.展开更多
The silicon on glasses process is a common preparation method of micro-electro-mechanical system inertial devices,which can realize the processing of thick silicon structures.This paper proposes that indium tin oxides...The silicon on glasses process is a common preparation method of micro-electro-mechanical system inertial devices,which can realize the processing of thick silicon structures.This paper proposes that indium tin oxides(ITO)film can serve as a deep silicon etching cut-off layer because ITO is less damaged under the attack of fluoride ions.ITO has good electrical conductivity and can absorb fluoride ions for silicon etching and reduce the reflection of fluoride ions,thus reducing the foot effect.The removal and release of ITO use an acidic solution,which does not damage the silicon structure.Therefore,the selection of the sacrificial layer has an excellent effect in maintaining the shape of the MEMS structure.This method is used in the preparation of MEMS accelerometers with a structure thickness of 100μm and a feature size of 4μm.The over-etching of the bottom of the silicon structure caused by the foot effect is negligible.The difference between the simulated value and the designed value of the device characteristic frequency is less than 5%.This indicates that ITO is an excellent deep silicon etch stopper material.展开更多
The transverse relaxation time (T_(2)) cut-off value plays a crucial role in nuclear magnetic resonance for identifying movable and immovable boundaries, evaluating permeability, and determining fluid saturation in pe...The transverse relaxation time (T_(2)) cut-off value plays a crucial role in nuclear magnetic resonance for identifying movable and immovable boundaries, evaluating permeability, and determining fluid saturation in petrophysical characterization of petroleum reservoirs. This study focuses on the systematic analysis of T_(2) spectra and T_(2) cut-off values in low-permeability reservoir rocks. Analysis of 36 low-permeability cores revealed a wide distribution of T_(2) cut-off values, ranging from 7 to 50 ms. Additionally, the T_(2) spectra exhibited multimodal characteristics, predominantly displaying unimodal and bimodal morphologies, with a few trimodal morphologies, which are inherently influenced by different pore types. Fractal characteristics of pore structure in fully water-saturated cores were captured through the T_(2) spectra, which were calculated using generalized fractal and multifractal theories. To augment the limited dataset of 36 cores, the synthetic minority oversampling technique was employed. Models for evaluating the T_(2) cut-off value were separately developed based on the classified T_(2) spectra, considering the number of peaks, and utilizing generalized fractal dimensions at the weight <0 and the singular intensity range. The underlying mechanism is that the singular intensity and generalized fractal dimensions at the weight <0 can detect the T_(2) spectral shift. However, the T_(2) spectral shift has negligible effects on multifractal spectrum function difference and generalized fractal dimensions at the weight >0. The primary objective of this work is to gain insights into the relationship between the kurtosis of the T_(2) spectrum and pore types, as well as to predict the T_(2) cut-off value of low-permeability rocks using machine learning and data augmentation techniques.展开更多
Patients living with chronic kidney disease (CKD) are at high risk of cardiovascular events. Our aim in this study was to assess the cut-off value for lipoprotein (a) (Lp(a)) in CKD patients with a history of cardiova...Patients living with chronic kidney disease (CKD) are at high risk of cardiovascular events. Our aim in this study was to assess the cut-off value for lipoprotein (a) (Lp(a)) in CKD patients with a history of cardiovascular disease (CVD). This was a cross-sectional study. Variables including age, sex, history of CVD, body mass index and CKD stage, were collected during CKD patient’s first admission in the nephrology dialysis department. Blood samples were collected for quantitative determination of Lp(a) by immunoturbidimetric method. They were divided into two groups: CKD patients without history of CVD and CKD patients with history of CVD. Fisher’s exact test was used to assess associations with a significance level of 0.05%. Area under the curve (AUC) and new cut-off value for Lp(a) were identified by drawing Receiver Operating Characteristic (ROC) curve. A total of seventy CKD patients with median age of 43 years [minimum-maximum = 15 - 78 years] were included. Patients with history of CVD were 65.71% (46/70). New Lp(a) cut-off point in CKD patients with history of CVD was 66.50 nmol/L [sensitivity, 87.00%;specificity, 58.30%;AUC = 0.727;p = 0.000]. ROC curve demonstrated good performance of Lp(a) to screen CKD patients with history of CVD. Further research is needed to determine an LPA gene polymorphism’s contribution to increasing risk for CVD at each kidney disease stage.展开更多
基金the National Natural Science Foun-dation of China(Grant Nos.91963201 and 12174163)the 111 Project(Grant No.B20063).
文摘The lack of soft magnetic composites with high power density in MHz frequency range has become an obstacle in the efficient operation of the electrical and electronic equipments.Here,a promising method to increase the cut-off frequency of iron-based soft magnetic composites to hundreds of MHz is reported.The cut-off frequency is increased from 10 MHz to 1 GHz by modulating the height of the ring,the distribution of particles,and the particle size.The mechanism of cut-off frequency and permeability is the coherent rotation of domain modulated by inhomogeneous field due to the eddy current effect.An empirical formula for the cut-off frequency in a magnetic ring composed of iron-based particles is established from experimental data.This work provides an effective approach to fabricate soft magnetic composites with a cut-off frequency in hundreds of MHz.
基金Supported by National Key R&D Program of China,No.2016YFC1305700.
文摘BACKGROUND Glycated hemoglobin A1c(HbA1c)is considered the most suitable for diabetes mellitus diagnosis due to its accuracy and convenience.However,the effect of HbA1c on diabetic retinopathy(DR)in the Han and Korean populations in Jilin,China,remains inconclusive.AIM To determine the best cut-off of HbA1c for diagnosing DR among the Chinese.METHODS This cross-sectional study included 1933 participants from the Yanbian area of Jilin Province,China.Trained investigators employed a questionnaire-based survey,physical examination,laboratory tests,and fundus photography for the investigation.The best cut-off value for HbA1c was established via the receiver operating characteristic curve.The factors associated with HbA1c-associated risk factors were determined via linear regression.RESULTS The analysis included 887 eligible Chinese Han and Korean participants,591 of whom were assigned randomly to the training set and 296 to the validation set.The prevalence of DR was 3.27% in the total population.HbA1c of 6.2% was the best cut-off value in the training set,while it was 5.9% in the validation set.In both Chinese Han and Korean populations,an HbA1c level of 6.2% was the best cut-off value.The optimal cut-off values of fasting blood glucose(FBG)≥7 mmol/L and<7 mmol/L were 8.1% and 6.2% respectively in Han populations,while those in Korean populations were 6.9%and 5.3%,respectively.Age,body mass index,and FBG were determined as the risk factors impacting HbA1c levels.CONCLUSION HbA1c may serve as a useful diagnostic indicator for DR.An HbA1c level of 6.2% may be an appropriate cut-off value for DR detection in the Chinese population.
基金supported by the National Science Foundation of China(No.41174100)the Large-scale Oil and Gas Field and Coalbed Methane Development Major Projects(No.2011ZX05019-008-08)the China National Petroleum Corporation(No.2014A-3609)
文摘Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.
基金supported by the National Natural Science Foundation of China(Nos.51475012,11772014,and 11272021)
文摘It is difficult to quantitatively detect defects by using the time domain or frequency domain features of Lamb wave signals due to their dispersion and multimodal characteristics.Therefore,it is important to discover an intrinsical parameter of Lamb waves that could be used as a damage sensitive feature.In this paper,quantitative defect detection in aluminium plates is carried out by means of wavenumber analysis approach.The wavenumber of excited Lamb wave mode is a fixed value,given a frequency,a thickness and material properties of the target plate.When Lamb waves propagate to the structural discontinuity,new wavenumber components are created by abrupt wavefield change.The new wavenumber components can be identified in the frequency-wavenumber domain.To estimate spatially dependent wavenumber values,a short-space two-dimensional Fourier transform(FT)method is presented for processing wavefield data of Lamb waves.The results can be used to determine the location,size and depth of rectangular notch.The analysis techniques are demonstrated using simulation examples of an aluminium plate with a rectangular notch.Then,the wavenumber analysis method is applied to simulation data that are obtained through a range of notch depths and widths.The results are analyzed and rules of the technique with regards to estimating notch depth are determined.Based on simulation results,guidelines for using the technique are developed.Finally,experimental wavefield data are obtained in aluminium plates with rectangular notches by a full noncontact transceiving method,i.e.,laser-laser method.Band-pass filtering combined with continuous wavelet transform is used to extract a certain frequency component from the full laser-induced wavefield with wide band.Shortspace two-dimensional FT method is used for further processing full wavefield data at a certain frequency to estimate spatially dependent wavenumber values.The consistency of simulation and experimental results shows the effectiveness of proposed wavenumber method for quantitative rectangular notch detection.
基金supported by the National Natural Science Foundation of China(Grant No.11125420)the Knowledge Innovation Program of the Chinese Academy of Sciences+2 种基金the China Postdoctoral Science Foundation(Grant No.2014M561882)the Doctoral Fund of Shandong ProvinceChina(Grant No.BS2012HZ015)
文摘An exact solution based on the wavenumber integration method is proposed and implemented in a numerical model for the acoustic field in a Pekeris waveguide excited by either a point source in cylindrical geometry or a line source in plane geometry. Besides, an unconditionally stable numerical solution is also presented, which entirely resolves the stability problem in previous methods. Generally the branch line integral contributes to the total field only at short ranges, and hence is usually ignored in traditional normal mode models. However, for the special case where a mode lies near the branch cut, the branch line integral can contribute to the total field significantly at all ranges. The wavenumber integration method is well-suited for such problems. Numerical results are also provided, which show that the present model can serve as a benchmark for sound propagation in a Pekeris waveguide.
基金This work was supported by the National Key Research and Development Program of China (Nos. 2017YFC0601305, 2017YFC0602203, and 2017YFC0601606), National Science and Technology Major Project task (No. 2016ZX05027-002-03), National Natural Science Foundation of China (No. 41604098), and State Key Program of National Natural Science of China (No. 41430322).
文摘The current local wavenumber methods for the interpretation of magnetic anomalies compute the locations of geological bodies by solving complex matrices. Presently, such methods require to know the structural index, which is a parameter that represents the source type. The structural index is hard to know in real data; consequently, the precision of current methods is low. We present the fast local wavenumber (FLW) method, and define the squared sum of the horizontal and vertical local wavenumbers as the cumulative local wavenumber. The FLW method is the linear combination of the umulative local wavenumberand other wavenumbers, and is used to compute the locations and structural index of the source without a priori information and matrix solution. We apply the FLW method to synthetic magnetic anomalies, and the results suggest that the FLW method is insensitive to background and oblique magnetization. Next, we apply the FLW method to real magnetic data to obtain the location and structural index of the source.
基金The present work is supported by National Science Foundation of China(41604138,41427901,41621063,41474133,41674158,41874179,41322030).
文摘The wavenumber spectral components WN4 at the mesosphere and low thermosphere(MLT)altitudes(70–10 km)and in the latitude range between±45°are obtained from temperature data(T)observed by the Sounding of the Atmosphere using Broadband Emission Radiometry(SABER)instruments on board the National Aeronautics and Space Administration(NASA)’s Thermosphere–Ionosphere–Mesosphere Energetics and Dynamics(TIMED)spacecraft during the 11-year solar period from 2002 to 2012.We analyze in detail these spectral components WNk and obtain the main properties of their vertical profiles and global structures.We report that all of the wavenumber spectral components WNk occur mainly around 100 km altitude,and that the most prominent component is the wavenumber spectral component WN4 structure.Comparing these long duration temperature data with results of previous investigations,we have found that the yearly variation of spectral component WN4 is similar to that of the eastward propagating non-migrating diurnal tide with zonal wavenumber 3(DE3)at the low latitudes,and to that of the semi-diurnal tide with zonal wavenumber 2(SE2)at the mid-latitudes:the amplitudes of the A4 are larger during boreal summer and autumn at the low-latitudes;at the mid-latitudes the amplitudes have a weak peak in March.In addition,the amplitudes of component WN4 undergo a remarkable short period variation:significant day-to-day variation of the spectral amplitudes A4 occurs primarily in July and September at the low-latitudes.In summary,we conclude that the non-migrating tides DE3 and SE2 are likely to be the origins,at the low-latitudes and the mid-latitudes in the MLT region,respectively,of the observed wavenumber spectral component WN4.
基金supported by the Natural Science Foundation of China(No.41574127)the China Postdoctoral Science Foundation(No.2017M622608)the project for the independent exploration of graduate students at Central South University(No.2017zzts008)
文摘In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over large-scale and complex terrain requires additional methods. To this end, we have proposed a new topography-capable By performing a two-dimensional Fourier transform in the horizontal directions, threedimensional partial differential equations in the spatial domain were transformed into a group of independent, one-dimensional differential equations engaged with different wave numbers. These independent differential equations are highly parallel across different wave numbers. differential equations with different wave numbers, and the efficiency of solving fixedbandwidth linear equations was further improved by a chasing method. In a synthetic test, a prism model was used to verify the accuracy and reliability of the proposed algorithm by comparing the numerical solution with the analytical solution. We studied the computational precision and efficiency with and without topography using different Fourier transform methods. The results showed that the Guass-FFT method has higher numerical precision, while the standard FFT method is superior, in terms of computation time, for inversion and quantitative interpretation under complicated terrain.
基金Supported by the National Natural Science Foundation of China (No. 60372077)
文摘Capillary and capillary-gravity waves possess a random character, and the slope wavenumber spectra of them can be used to represent mean distributions of wave energy with respect to spatial scale of variability. But simple and practical models of the slope wavenumber spectra have not been put forward so far. In this article, we address the accurate definition of the slope wavenumber spectra of water surface capillary and capillary-gravity waves. By combining the existing slope wavenumber models and using the dispersion relation of water surface waves, we derive the slope wavenumber spectrum models of capillary and capillary-gravity waves. Simultaneously, by using the slope wavenumber models, the dependence of the slope wavenumber spectrum on wind speed is analyzed using data obtained in an experiment which was performed in a laboratory wind wave tank. Generally speaking, the slope wavenumber spectra are influenced profoundly by the wind speed above water surface. The slope wavenumber spectrum increases with wind speed obviously and do not cross each other for different wind speeds. But, for the same wind speed, the slope wavenumber spectra are essentially identical, even though the capillary and capillary-gravity waves are excited at different times and locations. Furthermore, the slope wavenumber spectra obtained from the models agree quite well with experimental results as regards both the values and the shape of the curve.
文摘It is complicated to model the acoustic field in stratified ocean for airborne aircraft,due to high speed of the source and air-to-water sound transmission.To our knowledge,there are very few papers in the open literature dealing with this complicated problem;but,in our opinion,they all require great amount of computation.We now propose a different method that requires much less computation.We improve the wavenumber integration method to model the received temporal signal for a moving source in stratified ocean and sum up in a concise form the core of our paper as follows:(A) Eq.(11) can be calculated by means of fast Chirp Z transform and the signals at all N time points are generated simultaneously;(B) direct numerical evaluation of the wavenumber integral in Eq.(4) produces large numerical errors;so it is necessary to shift the integration slightly below the real axis;(C) we compare the computation cost of direct calculation method with that of our fast calculation method;from the results presented in table 1,we can see that the fast calculation method consumes much less computation time,particularly for long duration signals;(D) for an airborne rapidly moving source,we compute the Doppler-shifted signals in shallow water and analyze their short-time Fourier transform;from Fig.1b,we can see that the received signals have multiple frequency components for a tonal source due to source motion and that each component corresponds to an arrival path.
基金supported by the National Key Research and Development Program of China(2021YFA0718000)NSF of China under Grant No.42175075the UK-China Research&Innovation Partnership Fund through the Met Office Climate Science for Service Partnership(CSSP)China as part of the Newton Fund.
文摘Northeast Asian cut-off lows are crucial cyclonic systems that can bring temperature and precipitation extremes over large areas.Skillful subseasonal forecasting of Northeast Asian cut-off lows is of great importance.Using two dynamical forecasting systems,one from the Beijing Climate Center(BCC-CSM2-HR)and the other from the Met Office(GloSea5),this study assesses simulation ability and subseasonal prediction skill for early-summer Northeast Asian cut-off lows.Both models are shown to have good ability in representing the spatial structure of cut-off lows,but they underestimate the intensity.The skillful prediction time scales for cut-off low intensity are about 10.2 days for BCC-CSM2-HR and 11.4 days for GloSea5 in advance.Further examination shows that both models can essentially capture the initial Rossby wave train,rapid growth and decay processes responsible for the evolution of cut-off lows,but the models show weaker amplitudes for the three-stage processes.The underestimated simulated strength of both the Eurasian midlatitude and East Asian subtropical jets may lead to the weaker local eddy-mean flow interaction responsible for the cut-off low evolution.
基金Project supported by the National Key Research and Development(R&D)Program of China(No.2016YFA0401200)the National Natural Science Foundation of China(Nos.11402167,11332007,11672204,11672205,and 11732011)
文摘The prediction on small disturbance propagation in complex three-dimensional(3D) boundary layers is of great significance in transition prediction methodology, especially in the aircraft design. In this paper, the linear stability theory(LST) with the equivalent spanwise wavenumber correction(ESWC) is proposed in order to accurately predict the linear evolution of a disturbance in a kind of boundary layer flow with a vital variation in the spanwise direction. The LST with the ESWC takes not only the scale of the mean flow with the significant variation but also the wavenumber evolution of the disturbance itself. Compared with the conventional LST, the results obtained by the new method are in excellent agreement with those of the numerical simulations. The LST with the ESWC is an effective method on the prediction of the disturbance evolution in 3D boundary layers, which improves the prediction of the LST in the applications to complex 3D boundary layers greatly.
基金supported by National Natural Science Foundation of China (No. 41574127, 42174080)Innovation research team project of Guangxi Natural Science Foundation (No. GXNSFGA380004)Central South University independent exploration and innovation project for Postgraduates (Nos. 2021zzts0831, 2021zzts0271)
文摘Large-scale,fine,and efficient numerical simulation of a geothermal field plays an important role in geothermal energy development.Confronted with the problem of large computation and high storage requirements for complex underground models in a three-dimensional(3-D)numerical simulation of a geothermal fi eld,a mixed space-wavenumber domain 3-D numerical simulation algorithm is proposed in this paper.According to the superposition principle of temperature field,the geothermal field is decomposed into background and abnormal temperature fi elds for calculation.The uniform layered model is used to solve the background field.When the abnormal field is solved,the horizontal two-dimensional(2-D)Fourier transform is used to transform the 3-D diff erential equation satisfi ed by an abnormal field into a series of one-dimensional ordinary differential equations with diff erent wavenumbers,which greatly reduces the calculation and storage.The unit division of an ordinary diff erential equation is fl exible,and the calculation amount is small.The algorithm fully takes advantage of the effi ciency of the Fourier transform and the quickness of the catch-up method to solve linear equations with a fixed bandwidth,which effectively improves the computational efficiency.Compared with the COMSOL Multiphysics professional simulation finite element software,the time consumption and memory requirements of the algorithm proposed in this paper are reduced by multiple orders of magnitude in terms of ensuring accuracy and the same mesh division.The more the number of calculated nodes is,the more obvious is the advantage.We design models to study the thermal conductivity,heat fl ux boundary,regional tectonic morphology,and topographic relief of the geothermal fi eld distribution.A 3-D geophysical model is developed based on topographic elevation data,geothermal geology,and geophysical exploration data in the Qiabuqia area of Gonghe Basin,Qinghai Province,China.Numerical simulation of the geothermal fi eld in this area is realized,which shows that the algorithm is suitable for precise and effi cient simulation of an arbitrary complex terrain and geological conditions.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11502271 and 11372328the Strategic Priority Research Program on Space Science of the Chinese Academy of Sciences under Grant Nos XDA04020405 and XDA04020202-05the China Manned Space Engineering Program
文摘Marangoni Benard convection, which is mainly driven by the thermocapillary (Marangoni) effect, occurs in a thin liquid layer heated uniformly from the bottom. The wavenumber of supercritical convection is studied experimentally in a 160×160-mm^2 cavity that & heated from the bottom block. The convection pattern & visualized by an infrared thermography camera. It is shown that the onset of the Benard cell is consistent with theoretical analysis. The wavenumber decreases obviously with increasing temperature, except for a slight increase near the onset. The wavenumber gradually approaches the minimum when the supercritical number e is larger than 10. Finally, a formula is devised to describe the wavenumber selection in supercritical convection.
文摘Laplace's tidal equations are of great importance in various fields of geophysics. Here, the special case of zonal symmetry (zonal wavenumber m = 0) is investigated, where degenerate sets of eigensolutions appear. New results are presented for the inclusion of dissipative processes and the case of unstable conditions. In both instances the (nonzero) eigenfrequencies are complex. In the latter case, additional stable (i.e. real) eigenfrequencies appear in the numerical results for the absolute value of the Lambparameter ε being larger than a critical value εc. Further, it is shown that any degeneracies are removed through the inclusion of dissipation. Moreover, asymptotic relations are derived employing the relation of Laplace's tidal equations for m = 0 to the spheroidal differential equation. The implications of these findings to numerical techniques are demonstrated and results of computations are presented.
基金the Laboratory Open Fund of Beijing Smart-chip Microelectronics Technology Co.,Ltd and Chinese National Science Foundation(Contract No.52075519 and 61974136).
文摘The silicon on glasses process is a common preparation method of micro-electro-mechanical system inertial devices,which can realize the processing of thick silicon structures.This paper proposes that indium tin oxides(ITO)film can serve as a deep silicon etching cut-off layer because ITO is less damaged under the attack of fluoride ions.ITO has good electrical conductivity and can absorb fluoride ions for silicon etching and reduce the reflection of fluoride ions,thus reducing the foot effect.The removal and release of ITO use an acidic solution,which does not damage the silicon structure.Therefore,the selection of the sacrificial layer has an excellent effect in maintaining the shape of the MEMS structure.This method is used in the preparation of MEMS accelerometers with a structure thickness of 100μm and a feature size of 4μm.The over-etching of the bottom of the silicon structure caused by the foot effect is negligible.The difference between the simulated value and the designed value of the device characteristic frequency is less than 5%.This indicates that ITO is an excellent deep silicon etch stopper material.
基金supported by National Natural Science Foundation of China(Nos.42002171,42172159)China Postdoctoral Science Foundation(Nos.2020TQ0299,2020M682520)Postdoctoral Innovation Science Foundation of Hubei Province of China.
文摘The transverse relaxation time (T_(2)) cut-off value plays a crucial role in nuclear magnetic resonance for identifying movable and immovable boundaries, evaluating permeability, and determining fluid saturation in petrophysical characterization of petroleum reservoirs. This study focuses on the systematic analysis of T_(2) spectra and T_(2) cut-off values in low-permeability reservoir rocks. Analysis of 36 low-permeability cores revealed a wide distribution of T_(2) cut-off values, ranging from 7 to 50 ms. Additionally, the T_(2) spectra exhibited multimodal characteristics, predominantly displaying unimodal and bimodal morphologies, with a few trimodal morphologies, which are inherently influenced by different pore types. Fractal characteristics of pore structure in fully water-saturated cores were captured through the T_(2) spectra, which were calculated using generalized fractal and multifractal theories. To augment the limited dataset of 36 cores, the synthetic minority oversampling technique was employed. Models for evaluating the T_(2) cut-off value were separately developed based on the classified T_(2) spectra, considering the number of peaks, and utilizing generalized fractal dimensions at the weight <0 and the singular intensity range. The underlying mechanism is that the singular intensity and generalized fractal dimensions at the weight <0 can detect the T_(2) spectral shift. However, the T_(2) spectral shift has negligible effects on multifractal spectrum function difference and generalized fractal dimensions at the weight >0. The primary objective of this work is to gain insights into the relationship between the kurtosis of the T_(2) spectrum and pore types, as well as to predict the T_(2) cut-off value of low-permeability rocks using machine learning and data augmentation techniques.
文摘Patients living with chronic kidney disease (CKD) are at high risk of cardiovascular events. Our aim in this study was to assess the cut-off value for lipoprotein (a) (Lp(a)) in CKD patients with a history of cardiovascular disease (CVD). This was a cross-sectional study. Variables including age, sex, history of CVD, body mass index and CKD stage, were collected during CKD patient’s first admission in the nephrology dialysis department. Blood samples were collected for quantitative determination of Lp(a) by immunoturbidimetric method. They were divided into two groups: CKD patients without history of CVD and CKD patients with history of CVD. Fisher’s exact test was used to assess associations with a significance level of 0.05%. Area under the curve (AUC) and new cut-off value for Lp(a) were identified by drawing Receiver Operating Characteristic (ROC) curve. A total of seventy CKD patients with median age of 43 years [minimum-maximum = 15 - 78 years] were included. Patients with history of CVD were 65.71% (46/70). New Lp(a) cut-off point in CKD patients with history of CVD was 66.50 nmol/L [sensitivity, 87.00%;specificity, 58.30%;AUC = 0.727;p = 0.000]. ROC curve demonstrated good performance of Lp(a) to screen CKD patients with history of CVD. Further research is needed to determine an LPA gene polymorphism’s contribution to increasing risk for CVD at each kidney disease stage.