Contribution to the Full 3D Finite Element Modelling of a Hybrid Stepping Motor with and without Current in the Coils

The paper presents our contribution to the full 3D finite element modelling of a hybrid stepping motor using COMSOL Multiphysics software. This type of four-phase motor has a permanent magnet interposed between the two identical and coaxial half stators. The calculation of the field with or without current in the windings (respectively with or without permanent magnet) is done using a mixed formulation with strong coupling. In addition, the local high saturation of the ferromagnetic material and the radial and axial components of the magnetic flux are taken into account. The results obtained make it possible to clearly observe, as a function of the intensity of the bus current or the remanent induction, the saturation zones, the lines, the orientations and the magnetic flux densities. 3D finite element modelling provide more accurate numerical data on the magnetic field through multiphysics analysis. This analysis considers the actual operating conditions and leads to the design of an optimized machine structure, with or without current in the windings and/or permanent magnet.

Second-Order MaxEnt Predictive Modelling Methodology. I: Deterministically Incorporated Computational Model (2nd-BERRU-PMD)

This work presents a comprehensive second-order predictive modeling (PM) methodology designated by the acronym 2nd-BERRU-PMD. The attribute “2nd” indicates that this methodology incorporates second-order uncertainties (means and covariances) and second-order sensitivities of computed model responses to model parameters. The acronym BERRU stands for “Best- Estimate Results with Reduced Uncertainties” and the last letter (“D”) in the acronym indicates “deterministic,” referring to the deterministic inclusion of the computational model responses. The 2nd-BERRU-PMD methodology is fundamentally based on the maximum entropy (MaxEnt) principle. This principle is in contradistinction to the fundamental principle that underlies the extant data assimilation and/or adjustment procedures which minimize in a least-square sense a subjective user-defined functional which is meant to represent the discrepancies between measured and computed model responses. It is shown that the 2nd-BERRU-PMD methodology generalizes and extends current data assimilation and/or data adjustment procedures while overcoming the fundamental limitations of these procedures. In the accompanying work (Part II), the alternative framework for developing the “second- order MaxEnt predictive modelling methodology” is presented by incorporating probabilistically (as opposed to “deterministically”) the computed model responses.

Second-Order MaxEnt Predictive Modelling Methodology. II: Probabilistically Incorporated Computational Model (2nd-BERRU-PMP)

This work presents a comprehensive second-order predictive modeling (PM) methodology based on the maximum entropy (MaxEnt) principle for obtaining best-estimate mean values and correlations for model responses and parameters. This methodology is designated by the acronym 2nd-BERRU-PMP, where the attribute “2nd” indicates that this methodology incorporates second- order uncertainties (means and covariances) and second (and higher) order sensitivities of computed model responses to model parameters. The acronym BERRU stands for “Best-Estimate Results with Reduced Uncertainties” and the last letter (“P”) in the acronym indicates “probabilistic,” referring to the MaxEnt probabilistic inclusion of the computational model responses. This is in contradistinction to the 2nd-BERRU-PMD methodology, which deterministically combines the computed model responses with the experimental information, as presented in the accompanying work (Part I). Although both the 2nd-BERRU-PMP and the 2nd-BERRU-PMD methodologies yield expressions that include second (and higher) order sensitivities of responses to model parameters, the respective expressions for the predicted responses, for the calibrated predicted parameters and for their predicted uncertainties (covariances), are not identical to each other. Nevertheless, the results predicted by both the 2nd-BERRU-PMP and the 2nd-BERRU-PMD methodologies encompass, as particular cases, the results produced by the extant data assimilation and data adjustment procedures, which rely on the minimization, in a least-square sense, of a user-defined functional meant to represent the discrepancies between measured and computed model responses.

Second-Order MaxEnt Predictive Modelling Methodology. III: Illustrative Application to a Reactor Physics Benchmark

This work illustrates the innovative results obtained by applying the recently developed the 2nd-order predictive modeling methodology called “2nd- BERRU-PM”, where the acronym BERRU denotes “best-estimate results with reduced uncertainties” and “PM” denotes “predictive modeling.” The physical system selected for this illustrative application is a polyethylene-reflected plutonium (acronym: PERP) OECD/NEA reactor physics benchmark. This benchmark is modeled using the neutron transport Boltzmann equation (involving 21,976 uncertain parameters), the solution of which is representative of “large-scale computations.” The results obtained in this work confirm the fact that the 2nd-BERRU-PM methodology predicts best-estimate results that fall in between the corresponding computed and measured values, while reducing the predicted standard deviations of the predicted results to values smaller than either the experimentally measured or the computed values of the respective standard deviations. The obtained results also indicate that 2nd-order response sensitivities must always be included to quantify the need for including (or not) the 3rd- and/or 4th-order sensitivities. When the parameters are known with high precision, the contributions of the higher-order sensitivities diminish with increasing order, so that the inclusion of the 1st- and 2nd-order sensitivities may suffice for obtaining accurate predicted best- estimate response values and best-estimate standard deviations. On the other hand, when the parameters’ standard deviations are sufficiently large to approach (or be outside of) the radius of convergence of the multivariate Taylor-series which represents the response in the phase-space of model parameters, the contributions stemming from the 3rd- and even 4th-order sensitivities are necessary to ensure consistency between the computed and measured response. In such cases, the use of only the 1st-order sensitivities erroneously indicates that the computed results are inconsistent with the respective measured response. Ongoing research aims at extending the 2nd-BERRU-PM methodology to fourth-order, thus enabling the computation of third-order response correlations (skewness) and fourth-order response correlations (kurtosis).

NUMERICAL SIMULATION OF METHANE-AIR TURBULENT JET FLAME USING A NEW SECOND-ORDER MOMENT MODEL

A new second-order moment model for turbulent combustion is applied in the simulation of methane-air turbulent jet flame. The predicted results are compared with the experimental results and with those predicted using the well-known EBU-Arrhenius model and the original second-order moment model. The comparison shows the advantage of the new model that it requires almost the same computational storage and time as that of the original second-order moment model, but its modeling results are in better agreement with experiments than those using other models. Hence, the new second-order moment model is promising in modeling turbulent combustion with NOx formation with finite reaction rate for engineering application.

Closed-form steady-state solutions for forced vibration of second-order axially moving systems

Second-order axially moving systems are common models in the field of dynamics, such as axially moving strings, cables, and belts. In the traditional research work, it is difficult to obtain closed-form solutions for the forced vibration when the damping effect and the coupling effect of multiple second-order models are considered.In this paper, Green's function method based on the Laplace transform is used to obtain closed-form solutions for the forced vibration of second-order axially moving systems. By taking the axially moving damping string system and multi-string system connected by springs as examples, the detailed solution methods and the analytical Green's functions of these second-order systems are given. The mode functions and frequency equations are also obtained by the obtained Green's functions. The reliability and convenience of the results are verified by several examples. This paper provides a systematic analytical method for the dynamic analysis of second-order axially moving systems, and the obtained Green's functions are applicable to different second-order systems rather than just string systems. In addition, the work of this paper also has positive significance for the study on the forced vibration of high-order systems.

A two-scale second-order moment two-phase turbulence model for simulating dense gas-particle flows

A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concepts of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision and through a unified treatment of these two kinds of fluctuations. The proposed model is used to simulate gas-particle flows in a channel and in a downer. Simulation results are in agreement with the experimental results reported in references and are near the results obtained using the sin- gle-scale second-order moment two-phase turbulence model superposed with a particle collision model （USM-θ model） in most regions.

Comparisons of Wave Force Model Effects on the Structural Responses and Fatigue Loads of a Semi-Submersible Floating Wind Turbine

The selection of wave force models will significantly impact the structural responses of floating wind turbines.In this study,comparisons of wave force model effects on the structural responses and fatigue loads of a semi-submersible floating wind turbine(SFWT)were conducted.Simulations were performed by employing the Morison equation(ME)with linear or second-order wave kinematics and potential flow theory(PFT)with first-or second-order wave forces.A comparison of regular waves,irregular waves,and coupled wind/waves analyses with the experimental data showed that many of the simulation results and experimental data are relatively consistent.However,notable discrepancies are found in the response amplitude operators for platform heave,tower base bending moment,and tension in mooring lines.PFT models give more satisfactory results of heave but more significant discrepan-cies in tower base bending moment than the ME models.In irregular wave analyses,low-frequency resonances were captured by PFT models with second-order difference-frequency terms,and high-frequency resonances were captured by the ME models or PFT models with second-order sum-frequency terms.These force models capture the response frequencies but do not reasonably predict the response amplitudes.The coupled wind/waves analyses showed more satisfactory results than the wave-only analyses.However,an important detail to note is that this satisfactory result is based on the overprediction of wind-induced responses.

Derivation of a second-order model for Reynolds stress using renormalization group analysis and the two-scale expansion technique

With the two-scale expansion technique proposed by Yoshizawa,the turbulent fluctuating field is expanded around the isotropic field.At a low-order two-scale expansion,applying the mode coupling approximation in the Yakhot-Orszag renormalization group method to analyze the fluctuating field,the Reynolds-average terms in the Reynolds stress transport equation,such as the convective term,the pressure-gradient-velocity correlation term and the dissipation term,are modeled.Two numerical examples：turbulent flow past a backward-facing step and the fully developed flow in a rotating channel,are presented for testing the efficiency of the proposed second-order model.For these two numerical examples,the proposed model performs as well as the Gibson-Launder （GL） model,giving better prediction than the standard k-ε model,especially in the abilities to calculate the secondary flow in the backward-facing step flow and to capture the asymmetric turbulent structure caused by frame rotation.

EFFECT OF EMPIRICAL COEFFICIENTS ON SIMULATION IN TWO-SCALE SECOND-ORDER MOMENT PARTICLE-PHASE TURBULENCE MODEL

A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concept of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision. The proposed model is used to simulate gas-particle downer reactor flows. The computational results of both particle volume fraction and mean velocity are in agreement with the experimental results. After analyzing effects of empirical coefficient on prediction results, we can come to a conclusion that, inside the limit range of empirical coefficient, the predictions do not reveal a large sensitivity to the empirical coefficient in the downer reactor, but a relatively great change of the constants has important effect on the prediction.

SECOND-ORDER MOMENT MODEL FOR DENSE TWO-PHASE TURBULENT FLOW OF BINGHAM FLUID WITH PARTICLES

The USM-θ model of Bingham fluid for dense two-phase turbulent flow was developed, which combines the second-order moment model for two-phase turbulence with the particle kinetic theory for the inter-particle collision. In this model, phases interaction and the extra term of Bingham fluid yield stress are taken into account. An algorithm for USM-θ model in dense two-phase flow was proposed, in which the influence of particle volume fraction is accounted for. This model was used to simulate turbulent flow of Bingham fluid single-phase and dense liquid-particle two-phase in pipe. It is shown USM-θ model has better prediction result than the five-equation model, in which the particle-particle collision is modeled by the particle kinetic theory, while the turbulence of both phase is simulated by the two-equation turbulence model. The USM-θ model was then used to simulate the dense two-phase turbulent up flow of Bingham fluid with particles. With the increasing of the yield stress, the velocities of Bingham and particle decrease near the pipe centre. Comparing the two-phase flow of Bingham-particle with that of liquid-particle, it is found the source term of yield stress has significant effect on flow.

Static output feedback stabilization for second-order singular systems using model reduction methods

In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.

Dynamics of Rabi model under second-order Born Oppenheimer approximation

We apply the second-order Born-Oppenheimer （BO） approximation to investigate the dynamics of the Rabi model, which describes the interaction between a two-level system and a single bosonic mode beyond the rotating wave approxi- mation. By comparing with the numerical results, we find that our approach works well when the frequency of the two-level system is much smaller than that of the bosonic mode.

Effect of Time Step Size and Turbulence Model on the Open Water Hydrodynamic Performance Prediction of Contra-Rotating Propellers

A growing interest has been devoted to the contra-rotating propellers （CRPs） due to their high propulsive efficiency, torque balance, low fuel consumption, low cavitations, low noise performance and low hull vibration. Compared with the single-screw system, it is more difficult for the open water performance prediction because forward and aft propellers interact with each other and generate a more complicated flow field around the CRPs system. The current work focuses on the open water performance prediction of contra-rotating propellers by RANS and sliding mesh method considering the effect of computational time step size and turbulence model. The validation study has been performed on two sets of contra-rotating propellers developed by David W Taylor Naval Ship R ＆ D center. Compared with the experimental data, it shows that RANS with sliding mesh method and SST k-ω turbulence model has a good precision in the open water performance prediction of contra-rotating propellers, and small time step size can improve the level of accuracy for CRPs with the same blade number of forward and aft propellers, while a relatively large time step size is a better choice for CRPs with different blade numbers.

Breakdown voltage model and structure realization of a thin silicon layer with linear variable doping on a silicon on insulator high voltage device with multiple step field plates

Based on the theoretical and experimental investigation of a thin silicon layer（TSL） with linear variable doping（LVD） and further research on the TSL LVD with a multiple step field plate（MSFP）,a breakdown voltage（BV） model is proposed and experimentally verified in this paper.With the two-dimensional Poisson equation of the silicon on insulator（SOI） device,the lateral electric field in drift region of the thin silicon layer is assumed to be constant.For the SOI device with LVD in the thin silicon layer,the dependence of the BV on impurity concentration under the drain is investigated by an enhanced dielectric layer field（ENDIF）,from which the reduced surface field（RESURF） condition is deduced.The drain in the centre of the device has a good self-isolation effect,but the problem of the high voltage interconnection（HVI） line will become serious.The two step field plates including the source field plate and gate field plate can be adopted to shield the HVI adverse effect on the device.Based on this model,the TSL LVD SOI n-channel lateral double-diffused MOSFET（nLDMOS） with MSFP is realized.The experimental breakdown voltage（BV） and specific on-resistance（R on,sp） of the TSL LVD SOI device are 694 V and 21.3 ·mm 2 with a drift region length of 60 μm,buried oxide layer of 3 μm,and silicon layer of 0.15 μm,respectively.

Second-Order Adjoint Sensitivity Analysis Methodology for Computing Exactly Response Sensitivities to Uncertain Parameters and Boundaries of Linear Systems: Mathematical Framework

This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2nd-CASAM)” for the efficient and exact computation of 1st- and 2nd-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (i.e., model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2nd-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1st-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2nd-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (N2/2 + 3 N/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2nd-LASS requires very little additional effort beyond the construction of the 1st-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1st-LASS and 2nd-LASS requires the same computational solvers as needed for solving (i.e., “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1st-LASS and the 2nd-LASS. Since neither the 1st-LASS nor the 2nd-LASS involves any differentials of the operators underlying the original system, the 1st-LASS is designated as a “first-level” (as opposed to a “first-order”) adjoint sensitivity system, while the 2nd-LASS is designated as a “second-level” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2nd-LASS that involve the imprecisely known boundary parameters. Notably, the 2nd-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.

Model Test and Numerical Analysis on Long-Term Mechanical Properties of Stepped Reinforced Retaining Wall

Model tests and numerical analyses of stepped reinforced retaining wall were performed to investigate the effects of rheology of backfill and creep of geogrids on the long-term performance of the structure.The geogrid tensions,soil pressures,wall deformations and foundation pressure were measured during model construction and loading.A visco-elasto-plastic model and an empirical nonlinear visco-elastic model were utilized to simulate the stresses and deformations of geogrid-reinforced earth-retaining wall under long-term loads.By comparing test data with numerical results,it is shown that the foundation pressure distribution is nonlinear,and the lateral constraint of geogrids for backfill can cause a redistribution of foundation pressure.The curve of soil pressure is outside convex at each step initially,and it is close to the distribution for the case of vertical wall subsequently.The variation trend of geogrid tensions at different heights is obtained.Moreover,the failure mechanism and development mode of potential slip surface in retaining wall are proposed.

Development and Experimental Evaluation of a Steady-state Model for the Step-feed Biological Nitrogen Removal Process

In this article,a steady-state mathematical model was developed and experimentally evaluated to inves- tigate the effect of influent flow distribution and volume ratios of anoxic and aerobic zones in each stage on the to- tal nitrogen concentration of the effluent in the step-feed biological nitrogen removal process.Unlike the previous modeling methods,this model can be used to calculate the removal rates of ammonia and nitrate in each stage and thereby predict the concentrations of ammonia,nitrate,and total nitrogen in the effluent.To verify the simulation results,pilot-scale experimental studies were carried out in a four-stage step feed process.Good correlations were achieved between the measured data and the simulation results,which proved the validity of the developed model. The sensitivity of the model predictions was analyzed.After verification of the validity,the step feed process was optimally operated for five months using the model and the criteria developed for the design and operation.During the pilot-scale experimental period,the effluent total nitrogen concentrations were all below 5mg·L -1 ,with more than 90%removal efficiency.

Physical modelling and scale effects of air-water flows on stepped spillways

During the last three decades, the introduction of new construction materials (e.g. RCC (Roller Compacted Concrete), strengthened gabions) has increased the interest for stepped channels and spillways. However stepped chute hydraulics is not simple, because of different flow regimes and importantly because of very-strong interactions between entrained air and turbu- lence. In this study, new air-water flow measurements were conducted in two large-size stepped chute facilities with two step heights in each facility to study experimental distortion caused by scale effects and the soundness of result extrapolation to pro- totypes. Experimental data included distributions of air concentration, air-water flow velocity, bubble frequency, bubble chord length and air-water flow turbulence intensity. For a Froude similitude, the results implied that scale effects were observed in both facilities, although the geometric scaling ratio was only Lr=2 in each case. The selection of the criterion for scale effects is a critical issue. For example, major differences (i.e. scale effects) were observed in terms of bubble chord sizes and turbulence levels al- though little scale effects were seen in terms of void fraction and velocity distributions. Overall the findings emphasize that physical modelling of stepped chutes based upon a Froude similitude is more sensitive to scale effects than classical smooth-invert chute studies, and this is consistent with basic dimensional analysis developed herein.