High-order harmonic generation of ZnO crystals in chirped and static electric fields

High harmonic generation in ZnO crystals under chirped single-color field and static electric field are investigated by solving the semiconductor Bloch equation(SBE). It is found that when the chirp pulse is introduced, the interference structure becomes obvious while the harmonic cutoff is not extended. Furthermore, the harmonic efficiency is improved when the static electric field is included. These phenomena are demonstrated by the classical recollision model in real space affected by the waveform of laser field and inversion symmetry. Specifically, the electron motion in k-space shows that the change of waveform and the destruction of the symmetry of the laser field lead to the incomplete X-structure of the crystal-momentum-resolved(k-resolved) inter-band harmonic spectrum. Furthermore, a pre-acceleration process in the solid four-step model is confirmed.

High-order targeted essentially non-oscillatory scheme for two-fluid plasma model

The weakly ionized plasma flows in aerospace are commonly simulated by the single-fluid model,which cannot describe certain nonequilibrium phenomena by finite collisions of particles,decreasing the fidelity of the solution.Based on an alternative formulation of the targeted essentially non-oscillatory(TENO)scheme,a novel high-order numerical scheme is proposed to simulate the two-fluid plasmas problems.The numerical flux is constructed by the TENO interpolation of the solution and its derivatives,instead of being reconstructed from the physical flux.The present scheme is used to solve the two sets of Euler equations coupled with Maxwell's equations.The numerical methods are verified by several classical plasma problems.The results show that compared with the original TENO scheme,the present scheme can suppress the non-physical oscillations and reduce the numerical dissipation.

Time-varying gravity field model of Sichuan-Yunnan region based on the equivalent mass source model

High-precision time-varying gravity field is an effective way to study the internal mass movement and understanding the spatio-temporal evolution process of the geodynamic system.Compared to the satellite gravity measurement,the repeated terrestrial gravity observation can provide a more high-order signal related to the shallow crust and subsurface.However,the suitable and unified method for gravity model estimation is a key problem for further applications.In this study,we introduce the spherical hexahedron element to simulate the field source mass and forward model the change of gravity field located at the Sichuan-Yunnan region(99—104°E,23—29°N)in the four epochs from 2015 to 2017.Compared to the experimental results based on Slepian or spherical harmonics frequency domain method,this alternative approach is suitable for constructing the equivalent mass source model of regional-scale gravity data,by introducing the first-order smooth prior condition of gravity time-varying signal to suppress the high-frequency component of the signal.The results can provide a higher spatial resolution reference for regional gravity field modeling in the Sichuan-Yunnan region.

Calibration of quantitative rescattering model for simulating vortex high-order harmonic generation driven by Laguerre–Gaussian beam with nonzero orbital angular momentum

We calibrate the macroscopic vortex high-order harmonic generation(HHG)obtained by the quantitative rescattering(QRS)model to compute single-atom induced dipoles against that by solving the time-dependent Schr?dinger equation(TDSE).We show that the QRS perfectly agrees with the TDSE under the favorable phase-matching condition,and the QRS can accurately predict the main features in the spatial profiles of vortex HHG if the phase-matching condition is not good.We uncover that harmonic emissions from short and long trajectories are adjusted by the phase-matching condition through the time-frequency analysis and the QRS can simulate the vortex HHG accurately only when the interference between two trajectories is absent.This work confirms that it is an efficient way to employ the QRS model in the single-atom response for precisely simulating the macroscopic vortex HHG.

Research on Equivalent Modeling Method of AC-DC Power Networks Integrating with Renewable Energy Generation

Along with the increasing integration of renewable energy generation in AC-DC power networks,investigating the dynamic behaviors of this complex system with a proper equivalent model is significant.This paper presents an equivalent modeling method for the AC-DC power networks with doubly-fed induction generator(DFIG)based wind farms to decrease the simulation scale and computational burden.For the AC-DC power networks,the equivalent modeling strategy in accordance with the physical structure simplification is stated.Regarding the DFIG-based wind farms,the equivalent modeling based on the sequential identification of multi-machine parameters using the improved chaotic cuckoo search algorithm(ICCSA)is conducted.In light of the MATLAB simulation platform,a two-zone four-DC interconnected power grid with wind farms is built to check the efficacy of the proposed equivalentmodelingmethod.Fromthe simulation analyses and comparative validation in different algorithms and cases,the proposed method can precisely reflect the steady and dynamic performance of the demonstrated system under N-1 and N-2 fault scenarios,and it can efficiently achieve the parameter identification of the wind farms and fulfill the equivalent modeling.Consequently,the proposed approach’s effectiveness and suitability are confirmed.

Retrieving chlorophyll content and equivalent water thickness of Moso bamboo(Phyllostachys pubescens) forests under Pantana phyllostachysae Chao-induced stress from Sentinel-2A/B images in a multiple LUTs-based PROSAIL framework

Biochemical components of Moso bamboo(Phyllostachys pubescens)are critical to physiological and ecological processes and play an important role in the material and energy cycles of the ecosystem.The coupled PROSPECT with SAIL(PROSAIL)radiative transfer model is widely used for vegetation biochemical component content inversion.However,the presence of leaf-eating pests,such as Pantana phyllostachysae Chao(PPC),weakens the performance of the model for estimating biochemical components of Moso bamboo and thus must be considered.Therefore,this study considered pest-induced stress signals associated with Sentinel-2A/B images and field data and established multiple sets of biochemical canopy reflectance look-up tables(LUTs)based on the PROSAIL framework by setting different parameter ranges according to infestation levels.Quantitative inversions of leaf area index(LAI),leaf chlorophyll content(LCC),and leaf equivalent water thickness(LEWT)were derived.The scale conversions from LCC to canopy chlorophyll content(CCC)and LEWT to canopy equivalent water thickness(CEWT)were calculated.The results showed that LAI,CCC,and CEWT were inversely related with PPC-induced stress.When applying multiple LUTs,the p-values were<0.01;the R2 values for LAI,CCC,and CEWT were 0.71,0.68,and 0.65 with root mean square error(RMSE)(normalized RMSE,NRMSE)values of 0.38(0.16),17.56μg cm-2(0.20),and 0.02 cm(0.51),respectively.Compared to the values obtained for the traditional PROSAIL model,for October,R2 values increased by 0.05 and 0.10 and NRMSE decreased by 0.09 and 0.02 for CCC and CEWT,respectively and RMSE decreased by 0.35μg cm-2 for CCC.The feasibility of the inverse strategy for integrating pest-induced stress factors into the PROSAIL model,while establishing multiple LUTs under different pest-induced damage levels,was successfully demonstrated and can potentially enhance future vegetation parameter inversion and monitoring of bamboo forest health and ecosystems.

Equivalent series system to model a multiple friction pendulum system with numerous sliding interfaces for seismic analyses

Current structural analysis software programs offer few if any applicable device-specifi c hysteresis rules or nonlinear elements to simulate the precise mechanical behavior of a multiple friction pendulum system(MFPS) with numerous sliding interfaces.Based on the concept of subsystems,an equivalent series system that adopts existing nonlinear elements with parameters systematically calculated and mathematically proven through rigorous derivations is proposed.The aim is to simulate the characteristics of sliding motions for an MFPS isolation system with numerous concave sliding interfaces without prior knowledge of detailed information on the mobilized forces at various sliding stages.An MFPS with numerous concave sliding interfaces and one articulated or rigid slider located between these interfaces is divided into two subsystems: the fi rst represents the concave sliding interfaces above the slider,and the second represents those below the slider.The equivalent series system for the entire system is then obtained by connecting those for each subsystem in series.The equivalent series system is validated by comparing numerical results for an MFPS with four sliding interfaces obtained from the proposed method with those from a previous study by Fenz and Constantinou.Furthermore,these numerical results demonstrate that an MFPS isolator with numerous concave sliding interfaces,which may have any number of sliding interfaces,is a good isolation device to protect structures from earthquake damage through appropriate designs with controllable mechanisms.

Seepage simulation of high concrete-faced rockfill dams based on generalized equivalent continuum model

This research focused on the three-dimensional(3 D) seepage field simulation of a high concrete-faced rockfill dam(CFRD) under complex hydraulic conditions. A generalized equivalent continuum model of fractured rock mass was used for equivalent continuous seepage field analysis based on the improved node virtual flow method. Using a high CFRD as an example, the generalized equivalent continuum range was determined, and a finite element model was established based on the terrain and geological conditions, as well as structural face characteristics of the dam area. The equivalent seepage coefficients of different material zones or positions in the dam foundation were calculated with the Snow model or inverse analysis. Then, the 3 D seepage field in the dam area was calculated under the normal water storage conditions, and the corresponding water head distribution, seepage flow, seepage gradient, and seepage characteristics in the dam area were analyzed. The results show that the generalized equivalent continuum model can effectively simulate overall seepage patterns of the CFRD under complex hydraulic conditions and provide a reference for seepage analysis of similar CFRDs.

Hybrid Model Testing Technique for Deep-Sea Platforms Based on Equivalent Water Depth Truncation

In this paper, an inner turret moored FPSO which works in the water of 320 m depth, is selected to study the socalled ＂passively-truncated ＋ numerical-simulation＂ type of hybrid model testing technique while the tnmcated water depth is 160 m and the model scale ）, = 80. During the investigation, the optimization design of the equivalent-depth truncated system is performed by using the similarity of the static characteristics between the truncated system and the full depth one as the objective function. According to the truncated system, the corresponding physical test model is made. By adopting the coupling time domain simulation method, the tnmcated system model test is numerically reconstructed to carefully verify the computer simulation software and to adjust the corresponding hydrodynamic parameters. Based on the above work, the numerical extrapolation to the full depth system is performed by using the verified computer software and the adjusted hydrodyrmmic parameters. The full depth system model test is then performed in the basin and the results are compared with those from the numerical extrapolation. At last, the implementation procedure and the key technique of the hybrid model testing of the deep-sea platforms are summarized and printed. Through the above investigations, some beneficial conclusions are presented.

EQUIVALENT NORMAL CURVATURE APPROACH MILLING MODEL OF MACHINING FREEFORM SURFACES

A new milling methodology with the equivalent normal curvature milling model machining freeform surfaces is proposed based on the normal curvature theorems on differential geometry. Moreover, a specialized whirlwind milling tool and a 5-axis CNC horizontal milling machine are introduced. This new milling model can efficiently enlarge the material removal volume at the tip of the whirlwind milling tool and improve the producing capacity. The machining strategy of this model is to regulate the orientation of the whirlwind milling tool relatively to the principal directions of the workpiece surface at the point of contact, so as to create a full match with collision avoidance between the workpiece surface and the symmetric rotational surface of the milling tool. The practical results show that this new milling model is an effective method in machining complex three- dimensional surfaces. This model has a good improvement on finishing machining time and scallop height in machining the freeform surfaces over other milling processes. Some actual examples for manufacturing the freeform surfaces with this new model are given.

On iterative learning control design for tracking iteration-varying trajectories with high-order internal model

In this paper, iterative learning control （ILC） design is studied for an iteration-varying tracking problem in which reference trajectories are generated by high-order internal models （HOLM）. An HOlM formulated as a polynomial operator between consecutive iterations describes the changes of desired trajectories in the iteration domain and makes the iterative learning problem become iteration varying. The classical ILC for tracking iteration-invariant reference trajectories, on the other hand, is a special case of HOlM where the polynomial renders to a unity coefficient or a special first-order internal model. By inserting the HOlM into P-type ILC, the tracking performance along the iteration axis is investigated for a class of continuous-time nonlinear systems. Time-weighted norm method is utilized to guarantee validity of proposed algorithm in a sense of data-driven control.

Numerical modeling of wave equation by a truncated high-order finite-difference method

Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with increased order of accuracy. Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note that there exist some very small coefficients. With the order increasing, the number of these small coefficients increases, however, the values decrease sharply. An error analysis demonstrates that omitting these small coefficients not only maintain approximately the same level of accuracy of finite difference but also reduce computational cost significantly. Moreover, it is easier to truncate for the high-order finite-difference formulas than for the pseudospectral for- mulas. Thus this study proposes a truncated high-order finite-difference method, and then demonstrates the efficiency and applicability of the method with some numerical examples.

Equivalent continuum modeling of beam-like truss structures with flexible joints

The paper investigated the equivalent continuum modeling of beam-like repetitive truss structures considering the flexibility of joints,which models the contact between the truss member and joint by spring-damper with six directional stiffnesses and dampings.Firstly,a two-node hybrid joint-beam element was derived for modeling the truss member with flexible end joints,and a condensed model for the repeating element with flexible joints was obtained.Then,the energy equivalence method was adopted to equivalently model the truss structure with flexible joints and material damping as a spatial viscoelastic anisotropic beam model.Afterwards,the equations of motion for the equivalent beam model were derived and solved analytically in the frequency domain.In the numerical studies,the correctness of the presented method was verified by comparisons of the natural frequencies and frequency responses evaluated by the equivalent beam model with the results of the finite element method model.

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.

Optimum design of equivalent accelerated life testing plans based on proportional hazards-proportional odds model

The optimum design of equivalent accelerated life testing plan based on proportional hazards-proportional odds model using D-optimality is presented. The defined equivalent test plan is the test plan that has the same value of the determinant of Fisher information matrix. The equivalent test plan of step stress accelerated life testing （SSALT） to a baseline optimum constant stress accelerated life testing （CSALT） plan is obtained by adjusting the censoring time of SSALT and solving the optimization problem for each case to achieve the same value of the determinant of Fisher information matrix as in the baseline optimum CSALT plan. Numer- ical examples are given finally which demonstrate the equivalent SSALT plan to the baseline optimum CSALT plan reduces almost half of the test time while achieving approximately the same estimation errors of model parameters.

Modelling of Quasi 3-D Groundwater Flow and Studying of Equivalent Drainage Boundary in Beihai Peninsula,Guangxi

This paper describes a quasi 3-D finite element model of the groundwater flow in two -aquifer system which is constructed from a sequence of aquifer flow equations coupled by leakage terms representing flow through the aquitard . It is applied to evaluate the maximum rate of groundwater resources exploited from the coastal aquifer without seawater intrusion . The main task in this model is to determine the drainage boundary of the aquifer extending under the sea . The information of the boundary can be obtained from the fluctuations of the groundwater level caused by sea-tide fluctuations . A new idea, Equivalent Drainage Boundary (EDB), is proposed and the corresponding methods , determining the EDB, are developed with tidal fluctuations data observed in boreholes . The quasi 3-D model and the methods determining EDB have been applied to the aquifer system of Beihai peninsula , Guangxi Autonomous Region of China for calculating the available groundwater resources .

A High-Order Conservative Semi-Lagrangian Solver for 3D Free Surface Flows with Sediment Transport on Voronoi Meshes

In this paper,we present a conservative semi-Lagrangian scheme designed for the numeri-cal solution of 3D hydrostatic free surface flows involving sediment transport on unstruc-tured Voronoi meshes.A high-order reconstruction procedure is employed for obtaining a piecewise polynomial representation of the velocity field and sediment concentration within each control volume.This is subsequently exploited for the numerical integration of the Lagrangian trajectories needed for the discretization of the nonlinear convective and viscous terms.The presented method is fully conservative by construction,since the transported quantity or the vector field is integrated for each cell over the deformed vol-ume obtained at the foot of the characteristics that arises from all the vertexes defining the computational element.The semi-Lagrangian approach allows the numerical scheme to be unconditionally stable for what concerns the advection part of the governing equations.Furthermore,a semi-implicit discretization permits to relax the time step restriction due to the acoustic impedance,hence yielding a stability condition which depends only on the explicit discretization of the viscous terms.A decoupled approach is then employed for the hydrostatic fluid solver and the transport of suspended sediment,which is assumed to be passive.The accuracy and the robustness of the resulting conservative semi-Lagrangian scheme are assessed through a suite of test cases and compared against the analytical solu-tion whenever is known.The new numerical scheme can reach up to fourth order of accu-racy on general orthogonal meshes composed by Voronoi polygons.

Mass‑Spring‑Damping Theory Based Equivalent Mechanical Model for Cylindrical Lithium‑ion Batteries under Mechanical Abuse

An equivalent mechanical model with the equivalent physical meaning of mass-spring-damping is proposed for cylindrical lithium-ion batteries through experiments and theory.The equivalent mechanical model of a cylindrical lithium-ion battery consists of a spring-damping parallel unit.Therefore,a spring-damping parallel unit connecting a damping unit in series is selected to construct the constitutive characteristics of the battery under mechanical abuse.Comparison results show that the equivalent mechanical model can more effectively describe the mechanical properties of the batteries than most cubic fitting models,of which the average relative error of the equivalent mechanical model under different states-of-charge is less than 6.75%.Combined with the proposed equivalent mechanical model,the failure process of the batteries was simulated and analyzed using LS-Dyna and HyperWorks.Under rigid rod tests,failure occurred at the core and bottom of the batteries;under hemispherical punch tests,failure occurred at the core and top,consistent with the experimental results.The average prediction error for the failure displacement under different abuse conditions is less than 4%in the simulations.The equivalent mechanical model requires only a few parameters and can be recognized easily.In the future,the model can be used in safety warning devices based on mechanical penetration.

Multi-layer Quasi Three-dimensional Equivalent Model of Axial-Flux Permanent Magnet Synchronous Machine

Axial-flux permanent magnet synchronous machine(AFPMSM)enjoys the merits of high torque density and high efficiency,which make it one good candidate in the direct-drive application.The AFPMSM is usually analyzed based on the three-dimensional finite element method(3D FEM)due to its three-dimensional magnetic field distribution.However,the 3D FEM suffers large amount of calculation,time-consuming and is not suitable for the optimization of AFPMSM.Addressing this issue,a multi-layer quasi three-dimensional equivalent model of the AFPMSM is investigated in this paper,which could take the end leakage into consideration.Firstly,the multi-layer quasi three-dimensional equivalent model of the AFPMSM with single stator and single rotor is derived in details,including the equivalent processes and conversions of structure dimensions,motion conditions and electromagnetic parameters.Then,to consider the influence of end leakage on the performance,a correction factor is introduced in the multi-layer quasi three-dimensional equivalent model.Finally,the proposed multi-layer quasi three-dimensional equivalent model is verified by the 3D FEM based on an AFPMSM under different structure parameters.It demonstrates that the errors of flux linkage and average torque obtained by the multi-layer quasi three-dimensional equivalent model and 3D FEM are only around 2%although the structure parameters of the AFPMSM are varied.Besides,the computation time of one case based on the multi-layer quasi three-dimensional equivalent model is only 6 min,which is much less than that of the 3D FEM,1.8 h,under the same conditions.Thus,the proposed multi-layer quasi three-dimensional equivalent model could be used to optimize the AFPMSM and much time could be saved by this method compared with the 3D FEM.

An equivalent model of discharge instability in the discharge chamber of Kaufman ion thruster

The industrial application of the Kaufman ion thruster in its arc stage is limited owing to the instability of the discharge pulse.Presently,a complete prediction model that can predict the discharge pulse in the high-current stage does not exist.In this study,a complete prediction model for the pulse in the ion thruster is established using the zero-dimensional plasma discharge model and equivalent circuit model.The zero-dimensional plasma discharge model is used to obtain the corresponding plasma parameters by calculating the beam current,discharge current,voltage,and gas flow under actual working conditions.The input parameters of the equivalent circuit model are calculated using empirical formulae to acquire the estimated discharge waveforms.The pulse waveforms obtained using the model are found to be consistent with the experimental results.The model is used to evaluate the process of rapid changes in plasma density.Additionally,this model is employed to predict changes in the pulse waveforms when the volume of the discharge chamber and grid plate transmittance are changed.