Sensitivity Analysis of Electromagnetic Scattering from Dielectric Targets with Polynomial Chaos Expansion and Method of Moments

In this paper,an adaptive polynomial chaos expansion method(PCE)based on the method of moments(MoM)is proposed to construct surrogate models for electromagnetic scattering and further sensitivity analysis.The MoM is applied to accurately solve the electric field integral equation(EFIE)of electromagnetic scattering from homogeneous dielectric targets.Within the bistatic radar cross section(RCS)as the research object,the adaptive PCE algorithm is devoted to selecting the appropriate order to construct the multivariate surrogate model.The corresponding sensitivity results are given by the further derivative operation,which is compared with those of the finite difference method(FDM).Several examples are provided to demonstrate the effectiveness of the proposed algorithm for sensitivity analysis of electromagnetic scattering from homogeneous dielectric targets.

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.

Comparison of a Full Second-Order Moment Model and an Algebraic Stress Two-Phase Turbulence Model for Simulating Bubble-Liquid Flows in a Bubble Column

A full second-order moment (FSM) model and an algebraic stress (ASM) two-phase turbulence modelare proposed and applied to predict turbulent bubble-liquid flows in a 2D rectangular bubble column. Predictiongives the bubble and liquid velocities, bubble volume fraction, bubble and liquid Reynolds stresses and bubble-liquidvelocity correlation. For predicted two-phase velocities and bubble volume fraction there is only slight differencebetween these two models, and the simulation results using both two models are in good agreement with the particleimage velocimetry (PIV) measurements. Although the predicted two-phase Reynolds stresses using the FSM are insomewhat better agreement with the PIV measurements than those predicted using the ASM, the Reynolds stressespredicted using both two models are in general agreement with the experiments. Therefore, it is suggested to usethe ASM two-phase turbulence model in engineering application for saving the computation time.

A Second-Order Method for the Electromagnetic Scattering from a Large Cavity

In this paper,we study the electromagnetic scattering from a two dimen- sional large rectangular open cavity embedded in an infinite ground plane,which is modelled by Helmholtz equations.By introducing nonlocal transparent boundary con- ditions,the problem in the open cavity is reduced to a bounded domain problem.A hypersingular integral operator and a weakly singular integral operator are involved in the TM and TE cases,respectively.A new second-order Toeplitz type approximation and a second-order finite difference scheme are proposed for approximating the hyper- singular integral operator on the aperture and the Helmholtz in the cavity,respectively. The existence and uniqueness of the numerical solution in the TE case are established for arbitrary wavenumbers.A fast algorithm for the second-order approximation is pro- posed for solving the cavity model with layered media.Numerical results show the second-order accuracy and efficiency of the fast algorithm.More important is that the algorithm is easy to implement as a preconditioner for cavity models with more general media.

A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

Measurement of remanent magnetic moment using a torsion pendulum with single frequency modulation method

In Tian Qin spaceborne gravitational-wave detectors, the stringent requirements on the magnetic cleanliness of the test masses demand the high resolution ground-based characterization measurement of their magnetic properties. Here we present a single frequency modulation method based on a torsion pendulum to measure the remanent magnetic moment mr of 1.1 kg dummy copper test mass, and the measurement result is(6.45 ± 0.04(stat) ± 0.07(syst)) × 10^(-8)A · m^(2). The measurement precision of the mr is about 0.9 n A · m^(2), well below the present measurement requirement of Tian Qin. The method is particularly useful for measuring extremely low magnetic properties of the materials for use in the construction of space-borne gravitational wave detection and other precision scientific apparatus.

An Uncertainty Analysis and Reliability-Based Multidisciplinary Design Optimization Method Using Fourth-Moment Saddlepoint Approximation

In uncertainty analysis and reliability-based multidisciplinary design and optimization(RBMDO)of engineering structures,the saddlepoint approximation(SA)method can be utilized to enhance the accuracy and efficiency of reliability evaluation.However,the random variables involved in SA should be easy to handle.Additionally,the corresponding saddlepoint equation should not be complicated.Both of them limit the application of SA for engineering problems.The moment method can construct an approximate cumulative distribution function of the performance function based on the first few statistical moments.However,the traditional moment matching method is not very accurate generally.In order to take advantage of the SA method and the moment matching method to enhance the efficiency of design and optimization,a fourth-moment saddlepoint approximation(FMSA)method is introduced into RBMDO.In FMSA,the approximate cumulative generating functions are constructed based on the first four moments of the limit state function.The probability density function and cumulative distribution function are estimated based on this approximate cumulative generating function.Furthermore,the FMSA method is introduced and combined into RBMDO within the framework of sequence optimization and reliability assessment,which is based on the performance measure approach strategy.Two engineering examples are introduced to verify the effectiveness of proposed method.

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.

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.

Differential Quadrature Method for Steady Flow of an Incompressible Second-Order Viscoelastic Fluid and Heat Transfer Model

The two-dimensional steady flow of an incompressible second-order viscoelastic fluid between two parallel plates was studied in terms of vorticity, the stream function and temperature equations. The governing equations were expanded with respect to a snmll parameter to get the zeroth- and first-order approximate equations. By using the differenl2al quadrature method with only a few grid points, the high-accurate numerical results were obtained.

Two-Stream Approximation to the Radiative Transfer Equation:A New Improvement and Comparative Accuracy with Existing Methods

Mathematical modeling of the interaction between solar radiation and the Earth's atmosphere is formalized by the radiative transfer equation(RTE), whose resolution calls for two-stream approximations among other methods. This paper proposes a new two-stream approximation of the RTE with the development of the phase function and the intensity into a third-order series of Legendre polynomials. This new approach, which adds one more term in the expression of the intensity and the phase function, allows in the conditions of a plane parallel atmosphere a new mathematical formulation of γparameters. It is then compared to the Eddington, Hemispheric Constant, Quadrature, Combined Delta Function and Modified Eddington, and second-order approximation methods with reference to the Discrete Ordinate(Disort) method(δ –128 streams), considered as the most precise. This work also determines the conversion function of the proposed New Method using the fundamental definition of two-stream approximation(F-TSA) developed in a previous work. Notably,New Method has generally better precision compared to the second-order approximation and Hemispheric Constant methods. Compared to the Quadrature and Eddington methods, New Method shows very good precision for wide domains of the zenith angle μ 0, but tends to deviate from the Disort method with the zenith angle, especially for high values of optical thickness. In spite of this divergence in reflectance for high values of optical thickness, very strong correlation with the Disort method(R ≈ 1) was obtained for most cases of optical thickness in this study. An analysis of the Legendre polynomial series for simple functions shows that the high precision is due to the fact that the approximated functions ameliorate the accuracy when the order of approximation increases, although it has been proven that there is a limit order depending on the function from which the precision is lost. This observation indicates that increasing the order of approximation of the phase function of the RTE leads to a better precision in flux calculations. However, this approach may be limited to a certain order that has not been studied in this paper.

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.

Second-order two-scale computational method for ageing linear viscoelastic problem in composite materials with periodic structure

The correspondence principle is an important mathematical technique to compute the non-ageing linear viscoelastic problem as it allows to take advantage of the computational methods originally developed for the elastic case. However, the correspon- dence principle becomes invalid when the materials exhibit ageing. To deal with this problem, a second-order two-scale （SOTS） computational method in the time domain is presented to predict the ageing linear viscoelastic performance of composite materials with a periodic structure. First, in the time domain, the SOTS formulation for calcu- lating the effective relaxation modulus and displacement approximate solutions of the ageing viscoelastic problem is formally derived. Error estimates of the displacement ap- proximate solutions for SOTS method are then given. Numerical results obtained by the SOTS method are shown and compared with those by the finite element method in a very fine mesh. Both the analytical and numerical results show that the SOTS computational method is feasible and efficient to predict the ageing linear viscoelastic performance of composite materials with a periodic structure.

P-and S-wavefield simulations using both the firstand second-order separated wave equations through a high-order staggered grid finite-difference method

In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.

Stochastic Second-Order Two-Scale Method for Predicting the Mechanical Properties of Composite Materials with Random Interpenetrating Phase

In this paper,a stochastic second-order two-scale(SSOTS)method is proposed for predicting the non-deterministic mechanical properties of composites with random interpenetrating phase.Firstly,based on random morphology description functions(RMDF),the randomness of the material properties of the constituents as well as the correlation among these random properties are fully characterized through the topologies of the constituents.Then,by virtue of multiscale asymptotic analysis,the random effective quantities such as stiffness parameters and strength parameters along with their numerical computation formulae are derived by a SSOTS strategy combined with the Monte-Carlo method.Finally,the SSOTS method developed in this paper shows an excellent computational accuracy,and therefore present an important advance towards computationally efficient multiscale modeling frameworks considering microstructure uncertainties.

SECOND-ORDER ACCURATE DIFFERENCE METHOD FOR THE SINGULARLY PERTURBED PROBLEM OF FOURTH-ORDER ORDINARY DIFFERENTIAL EQUATIONS

In this paper, we construct a uniform second-order difference scheme for a class of boundary value problems of fourth-order ordinary differential equations. Finally, a numerical example is given.

Anti-plane deformations around arbitrary-shaped canyons on a wedge-shape half-space:moment method solutions

The wave propagation behavior in an elastic wedge-shaped medium with an arbitrary shaped cylindrical canyon at its vertex has been studied.Numerical computation of the wave displacement field is carried out on and near the canyon surfaces using weighted-residuals(moment method).The wave displacement fields are computed by the residual method for the cases of elliptic,circular,rounded-rectangular and flat-elliptic canyons,The analysis demonstrates that the resulting surface displacement depends,as in similar previous analyses,on several factors including,but not limited,to the angle of the wedge,the geometry of the vertex,the frequencies of the incident waves,the angles of incidence,and the material properties of the media.The analysis provides intriguing results that help to explain geophysical observations regarding the amplification of seismic energy as a function of site conditions.

SEMI-IMPLICIT SPECTRAL DEFERRED CORRECTION METHODS BASED ON SECOND-ORDER TIME INTEGRATION SCHEMES FOR NONLINEAR PDES

In[20],a semi-implicit spectral deferred correction(SDC)method was proposed,which is efficient for highly nonlinear partial differential equations(PDEs).The semi-implicit SDC method in[20]is based on first-order time integration methods,which are corrected iteratively,with the order of accuracy increased by one for each additional iteration.In this paper,we will develop a class of semi-implicit SDC methods,which are based on second-order time integration methods and the order of accuracy are increased by two for each additional iteration.For spatial discretization,we employ the local discontinuous Galerkin(LDG)method to arrive at fully-discrete schemes,which are high-order accurate in both space and time.Numerical experiments are presented to demonstrate the accuracy,efficiency and robustness of the proposed semi-implicit SDC methods for solving complex nonlinear PDEs.

Quadrature-based moment methods for the population balance equation: An algorithm review

The dispersed phase in multiphase flows can be modeled by the population balance model(PBM). A typical population balance equation(PBE) contains terms for spatial transport, loss/growth and breakage/coalescence source terms. The equation is therefore quite complex and difficult to solve analytically or numerically. The quadrature-based moment methods(QBMMs) are a class of methods that solve the PBE by converting the transport equation of the number density function(NDF) into moment transport equations. The unknown source terms are closed by numerical quadrature. Over the years, many QBMMs have been developed for different problems, such as the quadrature method of moments(QMOM), direct quadrature method of moments(DQMOM),extended quadrature method of moments(EQMOM), conditional quadrature method of moments(CQMOM),extended conditional quadrature method of moments(ECQMOM) and hyperbolic quadrature method of moments(Hy QMOM). In this paper, we present a comprehensive algorithm review of these QBMMs. The mathematical equations for spatially homogeneous systems with first-order point processes and second-order point processes are derived in detail. The algorithms are further extended to the inhomogeneous system for multiphase flows, in which the computational fluid dynamics(CFD) can be coupled with the PBE. The physical limitations and the challenging numerical problems of these QBMMs are discussed. Possible solutions are also summarized.