Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations

The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.

Three-Dimensional Simulations of RESET Operation in Phase-Change Random Access Memory with Blade-Type Like Phase Change Layer by Finite Element Modeling

An optimized device structure for reducing the RESET current of phase-change random access memory （PCRAM） with blade-type like （BTL） phase change layer is proposed. The electrical thermal analysis of the BTL cell and the blade heater contactor structure by three-dimensional finite element modeling are compared with each other during RESET operation. The simulation results show that the programming region of the phase change layer in the BTL cell is much smaller, and thermal electrical distributions of the BTL cell are more concentrated on the TiN/GST interface. The results indicate that the BTL cell has the superiorities of increasing the heating efficiency, decreasing the power consumption and reducing the RESET current from 0.67mA to 0.32mA. Therefore, the BTL cell will be appropriate for high performance PCRAM device with lower power consumption and lower RESET current.

Similarities and differences in inherent mechanism and characteristic frequency between the one-dimensional poroelastic model and the layered White model

The similarities and differences in inherent mechanism and characteristic frequency between the onedimensional(1D)poroelastic model and the layered White model were investigated.This investigation was conducted under the assumption that the rock was homogenous and isotropic at the mesoscopic scale.For the inherent mechanism,both models resulted from quasi-static flow in a slow P-wave diffusion mode,and the differences between them originated from saturated fluids and boundary conditions.On the other hand,for the characteristic frequencies of the models,the characteristic frequency of the 1D poroelastic model was first modified because the elastic constant and formula for calculating it were misused and then compared to that of the layered White model.Both of them moved towards higher frequencies with increasing permeability and decreasing viscosity and diffusion length.The differences between them were due to the diffusion length.The diffusion length for the 1D poroelastic model was determined by the sample length,whereas that for the layered White model was determined by the length of the representative elementary volume(REV).Subsequently,a numerical example was presented to demonstrate the similarities and differences between the models.Finally,published experimental data were interpreted using the 1D poroelastic model combined with the Cole-Cole model.The prediction of the combined model was in good agreement with the experimental data,thereby validating the effectiveness of the 1D poroelastic model.Furthermore,the modified characteristic frequency in our study was much closer to the experimental data than the previous prediction,validating the effectiveness of our modification of the characteristic frequency of the 1D poroelastic model.The investigation provided insight into the internal relationship between wave-induced fluid flow(WIFF)models at macroscopic and mesoscopic scales and can aid in a better understanding of the elastic modulus dispersion and attenuation caused by the WIFF at different scales.

A three-layer structure model for the effect of a soft middle layer on Love waves propagating in layered piezoelectric systems

A three-layer structure model is proposed for investigating the effect of a soft elastic middle layer on the propagation behavior of Love waves in piezoelectric layered systems, with ＂soft＂ implying that the bulk-shear-wave velocity of the middle layer is smaller than that of the upper sensitive layer. Dispersion equations are obtained for unelectroded and traction-free upper surfaces which, in the limit, can be reduced to those for classical Love waves. Systematic parametric studies are subsequently carried out to quantify the effects of the soft middle layer upon Love wave propagation, including its thickness, mass density, dielectric constant and elastic coefficient. It is demonstrated that whilst the thickness and elastic coefficient of the middle layer affect significantly Love wave propagation, its mass density and dielectric constant have negligible influence. On condition that both the thickness and elastic coefficient of the middle layer are vanishingly small so that it degenerates into an imperfectly bonded interface, the three-layer model is also employed to investigate the influence of imperfect interfaces on Love waves propagating in piezoelectric layer/elastic sub- strate systems. Upon comparing with the predictions ob- tained by employing the traditional shear-lag model, the present three-layer structure model is found to be more ac- curate as it avoids the unrealistic displacement discontinuity across imperfectly bonded interfaces assumed by the shearlag model, especially for long waves when the piezoelectric layer is relatively thin.

STUDY ON STRATEGY OF DYNAMIC JOINT ROUTING AND RESOURCE ALLOCATION IN LAYERED OPTICAL TRANSPORT NETWORKS

A layered network model for optical transport networks is proposed in this paper,which involves Internet Protocol(IP) ,Synchronous Digital Hierarchy(SDH) and Wavelength Division Mul-tiplexing(WDM) layers. The strategy of Dynamic Joint Routing and Resource Allocation(DJRRA) and its algorithm description are also presented for the proposed layered network model. DJRRA op-timizes the bandwidth usage of interface links between different layers and the logic links inside all layers. The simulation results show that DJRRA can reduce the blocking probability and increase network throughput effectively,which is in contrast to the classical separate sequential routing and resource allocation solutions.

A Simplified Scheme of the Generalized Layered Radiative Transfer Model

In this paper, firstly, a simplified version （SGRTM） of the generalized layered radiative transfer model （GRTM） within the canopy, developed by us, is presented. It reduces the information requirement of inputted sky diffuse radiation, as well as of canopy morphology, and in turn saves computer resources. Results from the SGRTM agree perfectly with those of the GRTM. Secondly, by applying the linear superposition principle of the optics and by using the basic solutions of the GRTM for radiative transfer within the canopy under the condition of assumed zero soil reflectance, two sets of explicit analytical solutions of radiative transfer within the canopy with any soil reflectance magnitude are derived： one for incident diffuse, and the other for direct beam radiation. The explicit analytical solutions need two sets of basic solutions of canopy reflectance and transmittance under zero soil reflectance, run by the model for both diffuse and direct beam radiation. One set of basic solutions is the canopy reflectance αf （written as α1 for direct beam radiation） and transmittance βf （written as β1 for direction beam radiation） with zero soil reflectance for the downward radiation from above the canopy （i.e. sky）, and the other set is the canopy reflectance （αb） and transmittance βb for the upward radiation from below the canopy （i.e., ground）. Under the condition of the same plant architecture in the vertical layers, and the same leaf adaxial and abaxial optical properties in the canopies for the uniform diffuse radiation, the explicit solutions need only one set of basic solutions, because under this condition the two basic solutions are equal, i.e., αf = αb and βf = βb. Using the explicit analytical solutions, the fractions of any kind of incident solar radiation reflected from （defined as surface albedo, or canopy reflectance）, transmitted through （defined as canopy transmittance）, and absorbed by （defined as canopy absorptance） the canopy and other properties pertinent to the radiative transfer within the canopy can be estimated easily on the ground surface below the canopy （soil or snow surface） with any reflectance magnitudes. The simplified transfer model is proven to have a similar accuracy compared to the detailed model, as well as very efficient computing.

Characterization of the adsorption behavior of aqueous cadmium on nanozero-valent iron based on orthogonal experiment and surface complexation modeling

Polyvinylpyrrolidone K-30(PVP) was introduced into the preparation of nanozero-valent iron(n ZVI) and the traditional liquid-phase reduction was improved. The introduction of PVP simplified the traditional method.The n ZVI prepared with this new approach showed excellent surface characters and high performance on the removal of cadmium. TEM results showed that the aggregates of n ZVI can reach to several micrometers in length but less than 100 nm in diameter. The iron particles that were enclosed by a layer of oxide film that is less than10 nm, demonstrated that the n ZVI possesses a core–shell structure. BET results indicate that the specific surface area of the n ZVI was 20.3159 m^2g^(-1). A three factor and three level orthogonal experiment was employed to find out the dominant factor that affects the removal rate of cadmium by n ZVI. Based on the range values, the prominence order of each factor was: initial p H of the solution N initial concentration of cadmium N dosage of n ZVI, the range was 96.453, 3.294 and 1.747, respectively. A simulation was performed under the same condition and a same conclusion was derived, this consistence confirmed the validity of the conclusion that p H is the most significant factor that affects the adsorption efficiency.

A numerical study for boundary layer current and sheet flow transport induced by a skewed asymmetric wave

An analytical model with essential parameters given by a two-phase numerical model is utilized to study the net boundary layer current and sediment transport under skewed asymmetric oscillatory sheet flows. The analytical model is the first instantaneous type model that can consider phase-lag and asymmetric boundary layer development. The two-phase model supplies the essential phase-lead, instantaneous erosion depth and boundary layer development for the analytical model to enhance the understanding of velocity skewness and acceleration skewness in sediment flux and transport rate. The sediment transport difference between onshore and offshore stages caused by velocity skewness or acceleration skewness is shown to illustrate the determination of net sediment transport by the analytical model. In previous studies about sediment transport in skewed asymmetric sheet flows, the generation of net sediment transport is mainly concluded to the phase-lag effect.However, the phase-lag effect is shown important but not enough for the net sediment transport, while the skewed asymmetric boundary layer development generated net boundary layer current and mobile bed effect are key important in the transport process.

Wind Structure in an Intermediate Boundary Layer Model Based on Ekman Momentum Approximation

A quasi three–dimensional, intermediate planetary boundary layer (PBL) model is developed by including inertial acceleration with the Ekman momentum approximation, but a nonlinear eddy viscosity based on Blackadar’s scheme was included to improve the theoretical model proposed by Tan and Wu (1993). The model could keep the same complexity as the classical Ekman model in numerical, but extends the conventional Ekman model to include the horizontal accelerated flow with the Ekman momentum approximation. A comparison between this modified Ekman model and other simplified accelerating PBL models is made. Results show that the Ekman model overestimates (underestimates) the wind speed and pumping velocity in the cyclonic (anticyclonic) shear flow due to the neglect of the acceleration flow, however, the semi–geostrophic Ekman model overestimates the acceleration effects resulting from the underestimating (overestimating) of the wind speed and pumping velocity in the cyclonic (anticyclonic) shear flow. The Ekman momentum approximation boundary layer model could be applied to the baroclinic atmosphere. The baroclinic Ekman momentum approximation boundary layer solution has both features of classical baroclinic Ekman layer and the Ekman momentum approximate boundary lager.

Critical properties of XY model dimensional layered magnetic on two- films

Using Monte Carlo simulations, we have investigated the classical XY model on triangular lattices of ultra-thin film structures with middle ferromagnetic layers sandwiched between two antiferromagnetic layers. The internal energy, the specifc heat, the chirality and the chiral susceptibility are calculated in order to clarify phase transitions and critical phenomena. Prom the finite-size scaling analyses, the values of critical exponents are determined. In a range of interaction parameters, we find that the chirality steeply goes up as temperature increases in a temperature range; correspondingly the value of a critical exponent for this change is estimated.

A 4-layer method of developing integrated sensor systems with LabVIEW

System integrity is important for fast and accurate measurement and control.LabVIEW is widely used in education and industry.Many LabVIEW codes are hard to be read and shown because of their 2D topology.In order to simplify the programming,a 4-layer model of developing sensor or measurement systems with LabVIEW is proposed in this paper.The purpose of this paper is to show the readers how to design simple,clear and strong automated systems with LabVIEW.Using a Sensirion SHT75 humidity sensor and an NI USB6008 DAQ board as an example,this paper describes the steps of developing a sensor system from the physical layer to application layer in detail.In layer 2,port selection and signal regulation are demonstrated.In layer 3,timing waveform analysis and synthesis,state diagram analysis,instruction set design,micro operation and fault tolerance designs are demonstrated.In layer 4,data visualization is covered with a vivid example.Programmers found it was hard to show readers their LabVIEW codes because many LabVIEW codes occupy several screens.A layer model simplifies the programming,so one is able to show a medium size LabVIEW code easily.

A rapid and accurate two-point ray tracing method in horizontally layered velocity model

A rapid and accurate method for two-point ray tracing in horizontally layered velocity model is presented in this paper. Numerical experiments show that this method provides stable and rapid convergence with high accuracies, regardless of various 1-D velocity structures, takeoff angles and epicentral distances. This two-point ray tracing method is compared with the pseudobending technique and the method advanced by Kim and Baag (2002). It turns out that the method in this paper is much more efficient and accurate than the pseudobending technique, but is only applicable to 1-D velocity model. Kim's method is equivalent to ours for cases without large takeoff angles, but it fails to work when the takeoff angle is close to 90°. On the other hand, the method presented in this paper is applicable to cases with any takeoff angles with rapid and accurate convergence. Therefore, this method is a good choice for two-point ray tracing problems in horizontally layered velocity model and is efficient enough to be applied to a wide range of seismic problems.

Contouring error modeling and simulation of a four-axis motion control system

A layered modeling method is proposed to resolve the problems resulting from the complexity of the error model of a multi-axis motion control system. In this model, a low level layer can be used as a virtual axis by the high level layer. The first advantage of this model is that the complex error model of a four-axis motion control system can be divided into several simple layers and each layer has different coupling strength to match the real control system. The second advantage lies in the fact that the controller in each layer can be designed specifically for a certain purpose. In this research, a three-layered cross coupling scheme in a four-axis motion control system is proposed to compensate the contouring error of the motion control system. Simulation results show that the maximum contouring error is reduced from 0.208 mm to 0.022 mm and the integration of absolute error is reduced from 0.108 mm to 0.015 mm, which are respectively better than 0.027 mm and 0.037 mm by the traditional method. And in the bottom layer the proposed method also has remarkable ability to achieve high contouring accuracy.

RECONSTRUCTION OF LAYER DATA WITH DEFORMABLE B-SPLINES

A new B-spline surface reconstruction method from layer data based on deformable model is presented. An initial deformable surface, which is represented as a closed cylinder, is firstly given. The surface is subject to internal forces describing its implicit smoothness property and external forces attracting it toward the layer data points. And then finite element method is adopted to solve its energy minimization problem, which results a bicubic closed B-spline surface with C^2 continuity. The proposed method can provide a smoothness and accurate surface model directly from the layer data, without the need to fit cross-sectional curves and make them compatible. The feasibility of the proposed method is verified by the experimental results.

Physical layer impairments tolerance based lightpath provision in software defined optical network

As all-optical networks grow with ever increasing ultra-high speed,the communication quality suffers seriously from physical layer impairments( PLIs). The same problem still exists in software defined optical network( SDON) controlled by OpenFlow. Aimed to solve this problem,a PLIs tolerance based lightpath provision scheme is proposed for OpenFlow controlled optical networks. This proposed approach not only takes the OSNR model to represent those linear PLIs factors,but also introduces those nonlinear factors into the OSNR model. Thus,the proposed scheme is able to cover most PLIs factors of each optical link and conduct optical lightpath provison with better communication quality. Moreover,PLIs tolerance model is also set up and considered in this work with some necessary extension to OpenFlow protocols to achieve better compatibility between physical layer impairments factors and various services connections. Simulation results show that the proposed scheme is able to get better performance in terms of packet loss rate and connection setup time.

A robust behavior of Feed Forward Back propagation algorithm of Artificial Neural Networks in the application of vertical electrical sounding data inversion

The applications of intelligent techniques have increased exponentially in recent days to study most of the non-linear parameters. In particular, the behavior of earth resembles the non- linearity applications. An efficient tool is needed for the interpretation of geophysical parameters to study the subsurface of the earth. Artificial Neural Networks （ANN） perform certain tasks if the structure of the network is modified accordingly for the purpose it has been used. The three most robust networks were taken and comparatively analyzed for their performance to choose the appropriate network. The single- layer feed-forward neural network with the back propagation algorithm is chosen as one of the well- suited networks after comparing the results. Initially, certain synthetic data sets of all three-layer curves have been taken tk^r training the network, and the network is validated by the field datasets collected from Tuticorin Coastal Region （78°7＇30＂E and 8°48＇45＂N）, Tamil Nadu, India. The interpretation has been done successfully using the corresponding learning algorithm in the present study. With proper training of back propagation networks, it tends to give the resistivity and thickness of the subsurface layer model of the field resistivity data concerning the synthetic data trained earlier in the appropriate network. The network is trained with more Vertical Electrical Sounding （VES） data, and this trained network is demon- strated by the field data. Groundwater table depth also has been modeled.

Building Morphological Characteristics and Their Effect on the Wind in Beijing

An urban boundary layer model （UBLM） is improved by incorporating the effect of buildings with a sectional drag coefficient and a height-distributed canopy drag length scale. The improved UBLM is applied to simulate the wind fields over three typical urban blocks over the Beijing area with different height-towidth ratios. For comparisons, the wind fields over the same blocks are simulated by an urban sub-domain scale model resolving the buildings explicitly. The wind fields simulated from the two different methods are in good agreement. Then, two-dimensional building morphological characteristics and urban canopy parameters for Beijing are derived from detailed building height data. Finally, experiements are conducted to investigate the effect of buildings on the wind field in Beijing using the improved UBLM.

Science Letters:Mask synthesis and verification based on geometric model for surface micro-machined MEMS

Traditional MEMS （microelectromechanical system） design methodology is not a structured method and has become an obstacle for MEMS creative design. In this paper, a novel method of mask synthesis and verification for surface micro-machined MEMS is proposed, which is based on the geometric model of a MEMS device. The emphasis is focused on synthesizing the masks at the basis of the layer model generated from the geometric model of the MEMS device. The method is comprised of several steps： the correction of the layer model, the generation of initial masks and final masks including multi-layer etch masks, and mask simulation. Finally some test resuhs are given.

Simulation study of a magnetocardiogram based on a virtual heart model:effect of a cardiac equivalent source and a volume conductor

In this paper, we present a magnetocardiogram （MCG） simulation study using the boundary element method （BEM） and based on the virtual heart model and the realistic human volume conductor model. The different contributions of cardiac equivalent source models and volume conductor models to the MCG are deeply and comprehensively investigated. The single dipole source model, the multiple dipoles source model and the equivalent double layer （EDL） source model are analysed and compared with the cardiac equivalent source models. Meanwhile, the effect of the volume conductor model on the MCG combined with these cardiac equivalent sources is investigated. The simulation results demonstrate that the cardiac electrophysiological information will be partly missed when only the single dipole source is taken, while the EDL source is a good option for MCG simulation and the effect of the volume conductor is smallest for the EDL source. Therefore, the EDL source is suitable for the study of MCG forward and inverse problems, and more attention should be paid to it in future MCG studies.

A virtual lock-in amplifier,spectrum analyzer,impedance meter and semiconductor analyzer implemented on an SR7265 hardware target

Lock-in amplifiers are used to detect and measure very small alternating current(AC)signals down to the range of nVs.Accurate measurements can be made even when the small signals are buried by noise thousands of times larger.With the digital signal processing(DSP)technology involved in modern instrumentation,a lock-in amplifier is more versatile in sensing and recovering small signals.Combining the virtual instrumentation technology,we reorganize the functional blocks of a programmable lock-in amplifier and build it as a virtual spectrum analyzer,virtual impedance meter,virtual network analyzer,virtual semiconductor parameter analyzer,signal generator,etc.A 4 layer model is used to implement these virtual instruments.The same virtual instrument can also be implemented on a general purpose FPGA developing board.