A coupled Legendre-Laguerre polynomial method with analytical integration for the Rayleigh waves in a quasicrystal layered half-space with an imperfect interface

The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.

Effect of microalloying on wettability and interface characteristics of Zr-based bulk metallic glasses with W substrate

The infiltration casting method is widely employed for the preparation of ex-situ composite materials.However,the production of composite materials using this method must necessitates a comprehensive understanding of the wettability and interface characteristics between the reinforcing phase and the bulk metallic glasses(BMGs).This work optimized the composition of Zr-based BMGs through microalloying methods,resulting in a new set of Zr-based BMGs with excellent glass-forming ability.Wetting experiments between the Zr-based BMGs melts and W substrates were conducted using the traditional sessile drop method,and the interfaces were characterized utilizing a scanning electron microscope(SEM)equipped with energy dispersive X-ray spectroscopy(EDS).The work demonstrates that the microalloying method substantially enhances the wettability of the Zr-based BMGs melt.Additionally,the incorporation of Nb element impedes the formation of W-Zr phases,but the introduction of Nb element does not alter the extent of interdiffusion between the constituent elements of the amorphous matrix and W element,indicating that the influence of Nb element on the diffusion of individual elements is minute.

Comparison of the interface reaction behaviors of CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) solid-state systems based on the diffusion couple method

The formation mechanism of calcium vanadate and manganese vanadate and the difference between calcium and manganese in the reaction with vanadium are basic issues in the calcification roasting and manganese roasting process with vanadium slag.In this work,CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) diffusion couples were prepared and roasted for different time periods to illustrate and compare the diffusion reaction mechanisms.Then,the changes in the diffusion product and diffusion coefficient were investigated and calculated based on scanning electron microscopy (SEM) with energy dispersive X-ray spectroscopy (EDS) analysis.Results show that with the extension of the roasting time,the diffusion reaction gradually proceeds among the CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) diffusion couples.The regional boundaries of calcium and vanadium are easily identifiable for the CaO–V_(2)O_(5) diffusion couple.Meanwhile,for the MnO_(2)–V_(2)O_(5) diffusion couple,MnO_(2) gradually decomposes to form Mn_(2)O_(3),and vanadium diffuses into the interior of Mn_(2)O_(3).Only a part of vanadium combines with manganese to form the diffusion production layer.CaV_(2)O_(6) and MnV_(2)O_(6) are the interfacial reaction products of the CaO–V_(2)O_(5) and MnO_(2)–V_(2)O_(5) diffusion couples,respectively,whose thicknesses are 39.85 and 32.13μm when roasted for 16 h.After 16 h,both diffusion couples reach the reaction equilibrium due to the limitation of diffusion.The diffusion coefficient of the CaO–V_(2)O_(5) diffusion couple is higher than that of the MnO_(2)–V_(2)O_(5) diffusion couple for the same roasting time,and the diffusion reaction between vanadium and calcium is easier than that between vanadium and manganese.

Growth interface of CdMnTe crystal by traveling heater method

The growth interfaces of CdMnTe（CMT） crystals grown by traveling heater method（THM） were studied. Two types of polycrystalline CMT feed ingots synthesized in a traditional rocking furnace and vertical Bridgman（VB） furnace were adopted in THM growth, and the effects of the polycrystalline feed on the growth interface were revealed. The morphology of the growth interface of CMT crystal（CMT2） grown from the feed by vertical Bridgman was smoother with lower curvature compared with that of CMT crystal（CMT1） from the feed by rocking furnace. The radial Mn composition and Te inclusion distribution of the CMT wafers were analyzed and correlated to the growth interface. The Mn segregation along the radial direction and Te inclusion density of CMT2 were lower than those of CMT1. The VB method synthesized polycrystalline feed could improve the growth interface morphology, which is beneficial for decreasing the Te inclusions and Mn segregation in CMT wafers.

Two-grid methods for semi-linear elliptic interface problems by immersed finite element methods

In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of finite element (FE) formulation and the simple structure of Cartesian grids,the IFE discretization is used in this paper.Two-grid schemes are formulated to linearize the FE equations.It is theoretically and numerically illustrated that the coarse space can be selected as coarse as H =O(h^1/4)(or H =O(h^1/8)),and the asymptotically optimal approximation can be achieved as the nonlinear schemes.As a result,we can settle a great majority of nonlinear equations as easy as linearized problems.In order to estimate the present two-grid algorithms,we derive the optimal error estimates of the IFE solution in the L^p norm.Numerical experiments are given to verify the theorems and indicate that the present two-grid algorithms can greatly improve the computing efficiency.

Studying the mechanical properties of the soil-root interface using the pullout test method

It is important to quantify the effect of the root diameter, the embedment length of the root and load speed on the soil-root interface mechanical properties for studying the root anchorage. The soilroot interface mechanical properties can be obtained through the pullout force and root slippage curve(F-S curve). About 120 Pinus tabulaeformis single roots whose diameters ranged from 1 mm to 10 mm divided into 6 groups based on different root embedment length(50 mm, 100 mm and 150 mm) and different load velocity(10 mm·min^(-1), 50 mm·min^(-1), 100 mm·min^(-1) and 300 mm·min^(-1)) were investigated using the pullout method. This study aims to explore the mechanical properties of the soil-root interface in the real conditions using the pullout test method. The results showed two kinds of pullout test failure modes during the experimental process: breakage failure and pullout failure. The results showed that the roots were easier to be broken when the root diameter was smaller or the loading speed was larger. The relationship between the maximum anchorage force and root diameter was linear and the linearly dependent coefficient(R^2) was larger than 0.85. The anchorage force increased with the root embedment length. An increase of 10%^(-1)5% for the maximumanchorage force was found when load speed increased from 10 to 300 mm.min^(-1). The mean peak slippage of the root was from 13.81 to 35.79 mm when the load velocity varied from 10 to 300 mm.min^(-1). The study will be helpful for the design of slopes reinforced by vegetation and in predicting risk of uprooting of trees, and will have practical benefits for understanding the mechanism of landslide.

Experimental investigation of interface curing stresses between PMMA and composite using digital speckle correlation method

This paper studies the interface curing stresses between polymethyl methacrylate (PMMA) and composite by means of digital speckle correlation method (DSCM).A new method by combining DSCM with the marker points is developed to measure the interface curing stresses,and the measurement principle is introduced.The interface curing stresses between PMMA and composite with different curing bonding conditions are measured and analyzed,this indicates that the residual stress for furnace heating and furnace cooling is the smallest.Finally,the measurement error is discussed by means of finite element method,the influences of glass microsphere between adhesive and PMMA can be ignored.

LOCAL AND PARALLEL FINITE ELEMENT METHOD FOR THE MIXED NAVIER-STOKES/DARCY MODEL WITH BEAVERS-JOSEPH INTERFACE CONDITIONS

In this paper, we consider the mixed Navier-Stokes/Darcy model with BeaversJoseph interface conditions. Based on two-grid discretizations, a local and parallel finite element algorithm for this mixed model is proposed and analyzed. Optimal errors are obtained and numerical experiments are presented to show the efficiency and effectiveness of the local and parallel finite element algorithm.

Immersed Interface Finite Element Methods for Elasticity Interface Problems with Non-Homogeneous Jump Conditions

In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.

P_1-nonconforming triangular finite element method for elliptic and parabolic interface problems

The lowest order Pl-nonconforming triangular finite element method （FEM） for elliptic and parabolic interface problems is investigated. Under some reasonable regularity assumptions on the exact solutions, the optimal order error estimates are obtained in the broken energy norm. Finally, some numerical results are provided to verify the theoretical analysis.

A NEW TREATMENT FOR INTERFACE OF METAL AND DIELECTRIC IN FDTD METHOD

A new treatment for interface of metal and dielectric in finite-difference time domain (FDTD) method was presented, further how and why this treatment should be as it is was explained. The FDTD method was developed to compute the low-frequency electromagnetic fields in the metal instead of the traditional high-frequency ones as it was originally proposed. The computational results agreed well with the experimental ones.

Modeling Methods and Test Verification of Root Insert Contact Interface for Wind Turbine Blade

Two modeling methods of the root insert for wind turbine blade are presented,i.e.,the local mesh optimization method(LMOM)and the global modeling method(GMM).Based on the optimized mesh of the local model for the metal contact interface,LMOM is proposed to analyze the load path and stress distribution characteristics,while GMM is used to calculate and analyze the stress distribution characteristics of the resin layer established between the bushing and composite layers of root insert.To validate the GMM,a tension test is carried out.The result successfully shows that the shear strain expresses a similar strain distribution tendency with the GMM′s results.

AN EMBEDDED BOUNDARY METHOD FOR ELLIPTIC AND PARABOLIC PROBLEMS WITH INTERFACES AND APPLICATION TO MULTI-MATERIAL SYSTEMS WITH PHASE TRANSITIONS

The embedded boundary method for solving elliptic and parabolic problems in geometrically complex domains using Cartesian meshes by Johansen and Colella （1998, J. Comput. Phys. 147, 60） has been extended for elliptic and parabolic problems with interior boundaries or interfaces of discontinuities of material properties or solutions. Second order accuracy is achieved in space and time for both stationary and moving interface problems. The method is conservative for elliptic and parabolic problems with fixed interfaces. Based on this method, a front tracking algorithm for the Stefan problem has been developed. The accuracy of the method is measured through comparison with exact solution to a two-dimensional Stefan problem. The algorithm has been used for the study of melting and solidification problems.

Effect of SDBS on interfacial electron transfer at the liquid/liquid interface by thin layer method

The effect of an adsorbed anionic surfactant sodium dodecyl benzene sulfonate （SDBS） on electron transfer （ET） reaction between ferricyanide aqueous solution and decamethylferrocene （DMFc） located on the adjacent organic phase was investigated for the first time by thin layer method. The adsorption of SDBS at the interface resulted in a decay in the cathodic plateau current of bimolecular reaction with increasing concentrations of SDBS in aqueous phase. However, the rate constant of electron transfer （ket） increased monotonically as the SDBS concentrations increased from 0 to 200 p, moFL. The experimental results showed that SDBS formed patches on the interface and influenced the structure of electrical double layer. 2009 Xiao Quan Lu. Published by Elsevier B.V. on behalf of Chinese Chemical Society. All rights reserved.

Distributed Lagrange Multiplier/Fictitious Domain Finite Element Method for a Transient Stokes Interface Problem with Jump Coefficients

The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface,where the arbitrary Lagrangian-Eulerian(ALE)technique is employed to deal with the moving and immersed subdomain.Stability and optimal convergence properties are obtained for both schemes.Numerical experiments are carried out for different scenarios of jump coefficients,and all theoretical results are validated.

Immersed Interface Method for Fokker-Planck Equation with Discontinuous Drift

The Immersed Interface Method (IIM) is derived to solve the corresponding Fokker-Planck equation of Brownian motion with pure dry friction, which is one of the simplest models of piecewise-smooth stochastic systems. The IIM is capable of treating a discontinuity in the drift of Fokker-Planck equation and it is readily extended to the dry and viscous friction model. Analytic results of the considered model are used to confirm the effectiveness and design accuracy of the scheme.

Quantitative analysis on influencing factors for interface propagation-based thermal conductivity measurement method during solid-liquid transition

The recently proposed interface propagation-based method has shown its advantages in obtaining the thermal conductivity of phase change materials during solid-liquid transition over conventional techniques. However, in previous investigation, the analysis on the measurement error was qualitative and only focused on the total effects on the measurement without decoupling the influencing factors. This paper discusses the effects of influencing factors on the measurement results for the interface propagation-based method. Numerical simulations were performed to explore the influencing factors, namely model simplification, subcooling and natural convection, along with their impact on the measurement process and corresponding measurement results. The numerical solutions were provided in terms of moving curves of the solid-liquid interface and the predicted values of thermal conductivity. Results indicated that the impact of simplified model was strongly dependent on Stefan number of the melting process. The degree of subcooling would lead to underestimated values for thermal conductivity prediction. The natural convection would intensify the heat transfer rate in the liquid region, thereby overestimating the obtained results of thermal conductivity. Correlations and experimental guidelines are provided. The relative errors are limited in ±1.5%,±3%and ±2% corresponding to the impact of simplified model, subcooling and natural convection, respectively.

Wall-modeling for large-eddy simulation of flows around an axisymmetric body using the diffuse-interface immersed boundary method

A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.

TRACKING AND RECONSTRUCTION METHODS FOR MOVING-INTERFACES

A simple,efficient and accurate high resolution method to tracking moving-interfaces-the characteristic integral-averaging finite volume method on unstructured meshes is proposed. And some numerical tests and evaluation of six main efficient methods for interface reconstruction are made. Through strict numerical simulation,their characters,advantages and shortcomings are compared,analyzed and commended in particular.

Superconvergence of Direct Discontinuous Galerkin Methods:Eigen-structure Analysis Based on Fourier Approach

This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.