An adaptive finite-difference method for seismic traveltime modeling based on 3D eikonal equation

3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.

A wavelet finite-difference method for numerical simulation of wave propagation in fluid-saturated porous media

In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.

Investigation of three-pulse photon echo in thick crystal using finite-difference time-domain method

This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo＇s amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^（3＋）：YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp（α＇L）,which is previously employed to describe the relationship between echo＇s efficiency and thickness,should be modified as 1.3 · 0.09 exp（2.4 ·α＇L） with the propagation of echo considered.

Optical simulation of in-plane-switching blue phase liquid crystal display using the finite-difference time-domain method

The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell＇s equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.

USE OF FINITE-DIFFERENCE TIME-DOMAIN METHOD FOR CALCULATING EM ABSORPTION IN LOSSY DIELECTRIC SCATTERER

The problem for calculating EM energy absorption by lossy dielectric scatterer ir-radiated by plane wave are discussed.The factors affecting the accuracy of computation arediscussed.The calculated results of EM energy absorption and its distribution in homogeneousand layered homogenous lossy dielectric spheres are presented,and a comparison of these resultswith analytical solution is given.The calculation is carried out for dielectric cylinder on conduct-ing ground as well,and the results are compared with the image theory.All the computationsshew that the finite-difference time-domain method can give satisfactory results.

Numerical Solution of a Problem of Thermal Stresses of a Magnetothermoelastic Cylinder with Rotation by Finite-Difference Method

The present article deals with the investigation thermal stress of a magnetothermoelastic cylinder subjected to rotation,open or closed circuit,thermal and mechanical boundary conditions.The outer and inner surfaces of the cylinder are subjected to both mechanical and thermal boundary conditions.A The transient coupled thermoelasticity in an infinite cylinder with its base abruptly exposed to a heat flux of a decaying exponential function of time is devised solve by the finite-difference method.The fundamental equations’system is solved by utilizing an implicit finite-difference method.This current method is a second-order accurate in time and space;it is also unconditionally stable.To illustrate the present model’s efficiency,we consider a suitable material and acquire the numerical solution of temperature,displacement components,and the components of stresses with time t and through the radial of an infinite cylinder.The results indicate that the effect of coupled thermoelasticity,magnetic field,and rotation on the temperature,stresses,and displacement is quite pronounced.In order to illustrate and verify the analytical developments,the numerical solution of partial differential equations,stress components,displacement components and temperature is carried out and computer simulated results are presented graphically.This study is helpful in the development of piezoelectric devices.

P-and S-wavefield simulations using both the firstand second-order separated wave equations through a high-order staggered grid finite-difference method

In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.

An efficient locally one-dimensional finite-difference time-domain method based on the conformal scheme

An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.

IMPROVED LOCALLY CONFORMAL FINITE-DIFFERENCE TIME-DOMAIN METHOD FOR EDGE INCLINED SLOTS IN A FINITE WALL THICKNESS WAVEGUIDE

An Improved Locally Conformal Finite-Difference Time-Domain (ILC-FDTD) method is presented in this paper, which is used to analyze the edge inclined slots penetrating adjacent broadwalls of a finite wall thickness waveguide. ILC-FDTD not only removes tile instability of the original locally conformal FDTD algorithm, but also improves the computational accuracy by locally modifying magnetic field update equations and the virtual iterative electric fields accordlng to the complexity of tile slot fringe fields. The mutual coupling between two edge inclined slots can also be analyzed by ILC-FDTD effectively.

3D finite-difference modeling algorithm and anomaly features of ZTEM

The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.

Computer Processing of the Intersection of Two Cones with Their Axes Intersecting and a Method to Get the Special Points

The analytic expression of the special points on the intersection of two cones with their axes intersecting(ITCTAI) is given. It also presents a method to construct the special points graphically according to the analytic expression of them. Finally, with computer programming language, it gives a program to generate the intersection in several different cases.

Measurement Method of Graviton Velocity and Thought on Correcting Lorentz Transformation——The Superluminal Neutrino’s Influence on Special Relativity and Other Physical Theories

After discovery of the superluminal particle and consideration on development of contemporary physical theory research, also on the existing errors and omissions, the principle of constant light speed is found not a necessary condition in derivation of Lorentz Transformation;instead, this thesis proposes the velocity of graviton may feature superluminal, constant velocity in different directions, and independence of inertial reference frame speeds. This is an optional thought of correction. According serial hypothesis, an equation of graviton’s motion trace, i.e., the central curve of nebula density, is established for spiral galaxy. Thus we gain the method to measure velocity of graviton. If to totally avoid problem of limit speed, we have to search for independent of inertia frames, and relevant to space-time properties. Regarding current difficulties of singular points in the Theory of Limited Universe, this thesis points out that the document [1] is the best solution to these difficulties.

An Investigation of Buckled of Column on Liner Elastic foundation by Using Finite-Difference Mode

The use of columns on elastic foundation is very common in Civil Engineering, like bridge pier, the foundation of the buildings etc. So, it will be useful to find the critical load for the structure, the problem in this paper will be solved by Finite-Difference Mode, that＇ s simple and has an extensive use. The way it works is that by dividing the component into many units. Finite-difference methods （FDM） are numerical methods for anoroximating, the solutions to differential eauations usine finite difference equations to approximate derivatives.

Research on the theories and methods of the specialization of the college counselor＇s work

The research proposition of this article is to do a brief analysis of the specialized theories and the systems of the methods of the college counselors＇ work under the perspective of the professional theories of the college counselors, and discuss the way of the realization of the specialization of the college counseling work, to achieve the professional orientation of the work of the college counselors. With the use of the specialized ways of work, we can make the college counselors develop themselves in the strict professional training and the continuous autonomous learning. They should use their own professional knowledge to equip themselves, to be the teachers and management personnel able to correctly guide the healthy growth of the university students.

Construction of Conservative Numerical Fluxes for the Entropy Split Method

The entropy split method is based on the physical entropies of the thermally perfect gas Euler equations.The Euler flux derivatives are approximated as a sum of a conservative portion and a non-conservative portion in conjunction with summation-by-parts(SBP)difference boundary closure of(Gerritsen and Olsson in J Comput Phys 129:245-262,1996;Olsson and Oliger in RIACS Tech Rep 94.01,1994;Yee et al.in J Comp Phys 162:33-81,2000).Sj?green and Yee(J Sci Comput)recently proved that the entropy split method is entropy conservative and stable.Stand-ard high-order spatial central differencing as well as high order central spatial dispersion relation preserving(DRP)spatial differencing is part of the entropy stable split methodol-ogy framework.The current work is our first attempt to derive a high order conservative numerical flux for the non-conservative portion of the entropy splitting of the Euler flux derivatives.Due to the construction,this conservative numerical flux requires higher oper-ations count and is less stable than the original semi-conservative split method.However,the Tadmor entropy conservative(EC)method(Tadmor in Acta Numerica 12:451-512,2003)of the same order requires more operations count than the new construction.Since the entropy split method is a semi-conservative skew-symmetric splitting of the Euler flux derivative,a modified nonlinear filter approach of(Yee et al.in J Comput Phys 150:199-238,1999,J Comp Phys 162:3381,2000;Yee and Sj?green in J Comput Phys 225:910934,2007,High Order Filter Methods for Wide Range of Compressible flow Speeds.Proceedings of the ICOSAHOM09,June 22-26,Trondheim,Norway,2009)is proposed in conjunction with the entropy split method as the base method for problems containing shock waves.Long-time integration of 2D and 3D test cases is included to show the com-parison of these new approaches.

Food Safety Detection Methods Applied to National Special Rectification of Product Quality and Food Safety

Afour-month period of national special rectification for product quality and food safety officially started on August 25, and was focused on eight fields, including those of agricultural products and processed foods.……

新工科视域下环境工程课程思政建设策略研究

课程思政和新工科建设人才培养的新要求相切合,二者的融合协同作用将对我国高校的环境工程专业课程建设和未来人才培养发挥重要的指导作用。当前我国环境工程专业的课程思政仍处在初步建设阶段,尚缺乏清晰的课程思政建设原则和与专业紧密结合的指导方针。该文以新工科建设为背景,依据2020年《高等学校课程思政建设指导纲要》和《环境科学与工程类教学质量国家标准(2018年)》,通过对环境工程专业课程思政建设理念和可行性的梳理,提出环境工程专业课程思政的构建原则和具体建设策略,为提升高校环境工程专业教学效果提供新的理论指导。

临床视角下医务社会工作专业化发展路径探析

目的:从多学科医学人文团队中临床医务人员的视角梳理医务社会工作的服务内容、临床服务流程和专业化发展路径,为医务社会工作标准化研究提供依据。方法:以江苏省人民医院为个案,运用质性研究法对16名参与多学科医学人文团队的临床医务人员进行结构式访谈。结果:运用Colaizzi七步分析法对访谈资料归纳整理,就医务社会工作的专业化发展提炼出5个主题,即需求、内容、优势、问题和期待。结论:医务社工作为临床工作中的专业角色已得到多学科医学人文团队中医护人员的认可,但在专业技能强化、服务流程优化、社会资源链接、科研水平提升等方面还需要进一步加强;可参照医务社会工作的通用过程模式,从“点-线-面”3个阶段逐步向规范化、科学化、专科化发展。

激光粒度仪在特殊土颗粒分析试验中的应用

采用激光粒度仪对膨胀土、黄土和红土进行了颗粒分析试验,研究了取样质量和超声波分散时间对试验结果的影响,发现遮光率与取样质量呈现很好的线性正相关,给出了特殊土在激光粒度仪法试验中适宜的遮光率范围(15%~25%)、取样质量(膨胀土取0.1~0.2 g,黄土取0.4~0.8 g,红土取0.2~0.4 g)和超声波分散时间(膨胀土取6 min,黄土和红土取2 min);激光粒度仪法试验结果的离散性,对于黄土和膨胀土较小,但对红土相对较大,红土平行试验中粒组含量极差的最大值为10.1%;与密度计法相比,激光粒度仪法测得的胶粒(0,0.002]mm含量偏小,而黏粒(0.002,0.005]mm或更大粒组的含量则偏大,可能是因为两者试验原理的不同而导致的;2种方法对红土的测试结果差异较大,与红土的特殊性质有一定关系。研究成果为激光粒度仪在特殊土颗粒分析试验中的应用提供了借鉴。

介入专科护士岗位胜任力模型的构建

目的 构建介入室专科护士岗位胜任力模型。方法2023年3月至2024年3月,组建课题小组,查阅国内外文献,运用德尔菲法对7名护理专家、9名医疗专家进行2轮问卷函询。结果 2轮函询问卷的有效回收率均为100%。经过16名专家2轮咨询,构建了4个一级指标及28个二级指标。结论 探索专科护士岗位胜任力模型构建,构建结果可信,具有一定的科学性和适用性,但指标体系仍需在实践中不断验证、修改、完善。