Inversion of Seabed Geotechnical Properties in the Arctic Chukchi Deep Sea Basin Based on Time Domain Adaptive Search Matching Algorithm

The chirp sub-bottom profiler,for its high resolution,easy accessibility and cost-effectiveness,has been widely used in acoustic detection.In this paper,the acoustic impedance and grain size compositions were obtained based on the chirp sub-bottom profiler data collected in the Chukchi Plateau area during the 11th Arctic Expedition of China.The time-domain adaptive search matching algorithm was used and validated on our established theoretical model.The misfit between the inversion result and the theoretical model is less than 0.067%.The grain size was calculated according to the empirical relationship between the acoustic impedance and the grain size of the sediment.The average acoustic impedance of sub-seafloor strata is 2.5026×10^(6) kg(s m^(2))^(-1)and the average grain size(θvalue)of the seafloor surface sediment is 7.1498,indicating the predominant occurrence of very fine silt sediment in the study area.Comparison of the inversion results and the laboratory measurements of nearby borehole samples shows that they are in general agreement.

Effects of Standing Time during Pretreatment on the Nitrite Concentration Detected by Spectrophotometric Method

Food safety problems caused by excessive nitrite addition have been frequently reported and the detection of nitrite in food is particularly important. The standing time during the pretreatment of primary sample has a great influence on the concentration of nitrite tested by spectrophotometric method. In this context, three kinds of food samples are prepared, including canned mustard, canned fish and home-made pickled water. A series of standing times are placed during the sample pretreatments and the corresponding nitrite contents in these samples are detected by spectrophotometric method based on N-ethylenediamine dihydrochloride. This study aims to find out a reasonable standing time during the pretreatment of food sample, providing influence factor for precise detection of nitrite.

A Time-Domain Numerical Simulation for Free Motion Responses of Two Ships Advancing in Head Waves

The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems within a time domain framework,the free water surface needs to simultaneously satisfy both the kinematic and dynamic boundary conditions of the free water surface.This provides conditions for adding artificial damping layers.Using the Runge−Kutta method to solve equations related to time.An upwind differential scheme is used in the present method to deal with the convection terms on the free surface to prevent waves upstream.Through the comparison with the available experimental data and other numerical methods,the present method is proved to have good mesh convergence,and satisfactory results can be obtained.The constant panel method is applied to calculate the hydrodynamic interaction responses of two parallel ships advancing in head waves.Numerical simulations are conducted on the effects of forward speed,different longitudinal and lateral distances on the motion response of two modified Wigley ships in head waves.Then further investigations are conducted on the effects of different ship types on the motion response.

A Novel Method for Linear Systems of Fractional Ordinary Differential Equations with Applications to Time-Fractional PDEs

This paper presents an efficient numerical technique for solving multi-term linear systems of fractional ordinary differential equations(FODEs)which have been widely used in modeling various phenomena in engineering and science.An approximate solution of the system is sought in the formof the finite series over the Müntz polynomials.By using the collocation procedure in the time interval,one gets the linear algebraic system for the coefficient of the expansion which can be easily solved numerically by a standard procedure.This technique also serves as the basis for solving the time-fractional partial differential equations(PDEs).The modified radial basis functions are used for spatial approximation of the solution.The collocation in the solution domain transforms the equation into a system of fractional ordinary differential equations similar to the one mentioned above.Several examples have verified the performance of the proposed novel technique with high accuracy and efficiency.

Highly Accurate Golden Section Search Algorithms and Fictitious Time Integration Method for Solving Nonlinear Eigenvalue Problems

This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.

Research on fault time prediction method for high speed rail BTM unit based on multi method interactive validation

The Balise Transmission Module(BTM)unit of the on-board train control system is a crucial component.Due to its unique installation position and complex environment,this unit has a higher fault rate within the on-board train control system.To conduct fault prediction for the BTM unit based on actual fault data,this study proposes a prediction method combining reliability statistics and machine learning,and achieves the fusion of prediction results from different dimensions through multi-method interactive validation.Firstly,a method for predicting equipment fault time targeting batch equipment is introduced.This method utilizes reliability statistics to construct a model of the remaining faultless operating time distribution considering uncertainty,thereby predicting the remaining faultless operating probability of the BTM unit.Secondly,considering the complexity of the BTM unit’s fault mechanism,the small sample size of fault cases,and the potential presence of multiple fault features in fault text records,an individual-oriented fault prediction method based on Bayesian-optimized Gradient Boosting Regression Tree(Bayes-GBRT)is proposed.This method achieves better prediction results compared to linear regression algorithms and random forest regression algorithms,with an average absolute error of only 0.224 years for predicting the fault time of this type of equipment.Finally,a multi-method interactive validation approach is proposed,enabling the fusion and validation of multi-dimensional results.The results indicate that the predicted fault time and the actual fault time conform to a log-normal distribution,and the parameter estimation results are basically consistent,verifying the accuracy and effectiveness of the prediction results.The above research findings can provide technical support for the maintenance and modification of BTM units,effectively reducing maintenance costs and ensuring the safe operation of high-speed railway,thus having practical engineering value for preventive maintenance.

Viscoacoustic prestack reverse time migration based onthe optimal time-space domain high-order finite-difference method

Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.

A Time Neighborhood Method for the Verification of Landfalling Typhoon Track Forecast

Landfalling typhoons can cause disasters over large regions.The government and emergency responders need to take measures to mitigate disasters according to the forecast of landfall position,while slight timing error can be ignored.The reliability of operational model forecasts of typhoon landfall position needs to be evaluated beforehand,according to the forecasts and observation of historical cases.In the evaluation of landfalling typhoon track,the traditional method based on point-to-point matching methods could be influenced by the predicted typhoon translation speed.Consequently,the traditional track evaluation method may result in a large track error even if the predicted landfall position is close to observation.The purpose of this paper is to address the above issue using a simple evaluation method of landfalling typhoon track forecast based on the time neighborhood approach.In this new method,the timing error was lessened to highlight the importance of the position error during the landfall of typhoon.The properties of the time neighborhood method are compared with the traditional method based on numerical forecast results of 12 landfalling typhoon cases.Results demonstrated that the new method is not sensitive to the sampling frequency,and that the difference between the time neighborhood and traditional method will be more obvious when the moving speed of typhoon is moderate(between 15−30 km h^(−1)).The time neighborhood concept can be easily extended to a broader context when one attempts to examine the position error more than the timing error.

Efficient Finite Difference/Spectral Method for the Time Fractional Ito Equation Using Fast Fourier Transform Technic

A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.

Dynamic interaction numerical models in the time domain based on the high performance scaled boundary finite element method

Consideration of structure-foundation-soil dynamic interaction is a basic requirement in the evaluation of the seismic safety of nuclear power facilities. An efficient and accurate dynamic interaction numerical model in the time domain has become an important topic of current research. In this study, the scaled boundary finite element method （SBFEM） is improved for use as an effective numerical approach with good application prospects. This method has several advantages, including dimensionality reduction, accuracy of the radial analytical solution, and unlike other boundary element methods, it does not require a fundamental solution. This study focuses on establishing a high performance scaled boundary finite element interaction analysis model in the time domain based on the acceleration unit-impulse response matrix, in which several new solution techniques, such as a dimensionless method to solve the interaction force, are applied to improve the numerical stability of the actual soil parameters and reduce the amount of calculation. Finally, the feasibility of the time domain methods are illustrated by the response of the nuclear power structure and the accuracy of the algorithms are dynamically verified by comparison with the refinement of a large-scale viscoelastic soil model.

Stochastic transient analysis of thermal stresses in solids by explicit time-domain method

Stochastic heat conduction and thermal stress analysis of structures has received considerable attention in recent years.The propagation of uncertain thermal environments will lead to stochastic variations in temperature fields and thermal stresses.Therefore,it is reasonable to consider the variability of thermal environments while conducting thermal analysis.However,for ambient thermal excitations,only stationary random processes have been investigated thus far.In this study,the highly efficient explicit time-domain method(ETDM)is proposed for the analysis of non-stationary stochastic transient heat conduction and thermal stress problems.The explicit time-domain expressions of thermal responses are first constructed for a thermoelastic body.Then the statistical moments of thermal displacements and stresses can be directly obtained based on the explicit expressions of thermal responses.A numerical example involving non-stationary stochastic internal heat generation rate is investigated.The accuracy and efficiency of the proposed method are validated by comparison with the Monte-Carlo simulation.

A Novel Absorbing Boundary Condition for the Frequency-Dependent Finite-Difference Time-Domain Method

A new absorbing boundary condition (ABC) for frequency dependent finite difference time domain algorithm for the arbitrary dispersive media is presented. The concepts of the digital systems are introduced to the (FD) 2TD method. On the basis of digital filter designing and vector algebra, the absorbing boundary condition under arbitrary angle of incidence are derived. The transient electromagnetic problems in two dimensions and three dimensions are calculated and the validity of the ABC is verified.

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.

Multi-scale seismic full waveform inversion in the frequency-domain with a multi-grid method

Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based on the optimal steep descent methods, we present an algorithm which combines the preconditioned bi-conjugated gradient stable method and the multi-grid method to compute the wave propagation and the gradient space. The multiple scale prosperity of the waveform inversion and the multi-grid method can overcome the inverse problems local minima defect and accelerate convergence. The local inhomogeneous three-hole model simulated results and the Marmousi model certify the algorithm effectiveness.

High-Order Leap-Frog Based Discontinuous Galerkin Method for the Time-Domain Maxwell Equations on Non-Conforming Simplicial Meshes

A high-order leap-frog based non-dissipative discontinuous Galerkin time- domain method for solving Maxwell＇s equations is introduced and analyzed. The pro- posed method combines a centered approximation for the evaluation of fluxes at the in- terface between neighboring elements, with a Nth-order leap-frog time scheme. More- over, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwelrs equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with high- order elements show the potential of the method.

Numerical Solutions of Finite Well in Two Dimensions Using the Finite Difference Time Domain Method

The higher excited states for two dimensional finite rectangular well potential are calculated numerically,by solving the Schrödinger equation using the finite difference time domain method.Although,this method is suitable to calculate the ground state of the quantum systems,it has been improved to calculate the higher excited states directly.The improvement is based on modifying the iterative process involved in this method to include two procedures.The first is known as cooling steps and the second is known as a heating step.By determining the required length of the cooling iteration steps using suitable excitation energy estimate,and repeating these two procedures using suitable initial guess function for sufficient times.This modified iteration will lead automatically to the desired excited state.In the two dimensional finite rectangular well potential problem both of the suitable excitation energy and the suitable initial guess wave function are calculated analytically using the separation of variables technique.

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.

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.

Distributed Dynamic Load in Structural Dynamics by the Impulse-Based Force Estimation Algorithm

This paper proposes a novel approach for identifying distributed dynamic loads in the time domain.Using polynomial andmodal analysis,the load is transformed intomodal space for coefficient identification.This allows the distributed dynamic load with a two-dimensional form in terms of time and space to be simultaneously identified in the form of modal force,thereby achieving dimensionality reduction.The Impulse-based Force Estimation Algorithm is proposed to identify dynamic loads in the time domain.Firstly,the algorithm establishes a recursion scheme based on convolution integral,enabling it to identify loads with a long history and rapidly changing forms over time.Secondly,the algorithm introduces moving mean and polynomial fitting to detrend,enhancing its applicability in load estimation.The aforementioned methodology successfully accomplishes the reconstruction of distributed,instead of centralized,dynamic loads on the continuum in the time domain by utilizing acceleration response.To validate the effectiveness of the method,computational and experimental verification were conducted.

Monte Carlo method for evaluation of surface emission rate measurement uncertainty

The aim of this study is to evaluate the uncertainty of 2πα and 2πβ surface emission rates using the windowless multiwire proportional counter method.This study used the Monte Carlo method (MCM) to validate the conventional Guide to the Expression of Uncertainty in Measurement (GUM) method.A dead time measurement model for the two-source method was established based on the characteristics of a single-channel measurement system,and the voltage threshold correction factor measurement function was indirectly obtained by fitting the threshold correction curve.The uncertainty in the surface emission rate was calculated using the GUM method and the law of propagation of uncertainty.The MCM provided clear definitions for each input quantity and its uncertainty distribution,and the simulation training was realized with a complete and complex mathematical model.The results of the surface emission rate uncertainty evaluation for four radioactive plane sources using both methods showed the uncertainty’s consistency E_(n)<0.070 for the comparison of each source,and the uncertainty results of the GUM were all lower than those of the MCM.However,the MCM has a more objective evaluation process and can serve as a validation tool for GUM results.