RKDG Methods with Multi-resolution WENO Limiters for Solving Steady-State Problems on Triangular Meshes

In this paper, we design high-order Runge-Kutta discontinuous Galerkin (RKDG) methods with multi-resolution weighted essentially non-oscillatory (multi-resolution WENO) limiters to compute compressible steady-state problems on triangular meshes. A troubled cell indicator extended from structured meshes to unstructured meshes is constructed to identify triangular cells in which the application of the limiting procedures is required. In such troubled cells, the multi-resolution WENO limiting methods are used to the hierarchical L^(2) projection polynomial sequence of the DG solution. Through using the RKDG methods with multi-resolution WENO limiters, the optimal high-order accuracy can be gradually reduced to first-order in the triangular troubled cells, so that the shock wave oscillations can be well suppressed. In steady-state simulations on triangular meshes, the numerical residual converges to near machine zero. The proposed spatial reconstruction methods enhance the robustness of classical DG methods on triangular meshes. The good results of these RKDG methods with multi-resolution WENO limiters are verified by a series of two-dimensional steady-state problems.

NUMERICAL SIMULATION OF UNSTEADY-STATE UNDEREXPANDED JET USING DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD

A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.

Interval finite difference method for steady-state temperature field prediction with interval parameters

A new numerical technique named interval finite difference method is proposed for the steady-state temperature field prediction with uncertainties in both physical parameters and boundary conditions. Interval variables are used to quantitatively describe the uncertain parameters with limited information. Based on different Taylor and Neumann series, two kinds of parameter perturbation methods are presented to approximately yield the ranges of the uncertain temperature field. By comparing the results with traditional Monte Carlo simulation, a numerical example is given to demonstrate the feasibility and effectiveness of the proposed method for solving steady-state heat conduction problem with uncertain-but-bounded parameters.

A Newton multigrid method for steady-state shallow water equations with topography and dry areas

A Newton multigrid method is developed for one-dimensional （1D） and two- dimensional （2D） steady-state shallow water equations （SWEs） with topography and dry areas. The nonlinear system arising from the well-balanced finite volume discretization of the steady-state SWEs is solved by the Newton method as the outer iteration and a geometric multigrid method with the block symmetric Gauss-Seidel smoother as the inner iteration. The proposed Newton multigrid method makes use of the local residual to regularize the Jacobian matrix of the Newton iteration, and can handle the steady- state problem with wet/dry transition. Several numerical experiments are conducted to demonstrate the efficiency, robustness, and well-balanced property of the proposed method. The relation between the convergence behavior of the Newton multigrid method and the distribution of the eigenvalues of the iteration matrix is detailedly discussed.

The Improved Element-Free Galerkin Method for Anisotropic Steady-State Heat Conduction Problems

In this paper,we considered the improved element-free Galerkin(IEFG)method for solving 2D anisotropic steadystate heat conduction problems.The improved moving least-squares(IMLS)approximation is used to establish the trial function,and the penalty method is applied to enforce the boundary conditions,thus the final discretized equations of the IEFG method for anisotropic steady-state heat conduction problems can be obtained by combining with the corresponding Galerkin weak form.The influences of node distribution,weight functions,scale parameters and penalty factors on the computational accuracy of the IEFG method are analyzed respectively,and these numerical solutions show that less computational resources are spent when using the IEFG method.

Insight Into the Separation-of-Variable Methods for the Closed-Form Solutions of Free Vibration of Rectangular Thin Plates

The separation-of-variable(SOV)methods,such as the improved SOV method,the variational SOV method,and the extended SOV method,have been proposed by the present authors and coworkers to obtain the closed-form analytical solutions for free vibration and eigenbuckling of rectangular plates and circular cylindrical shells.By taking the free vibration of rectangular thin plates as an example,this work presents the theoretical framework of the SOV methods in an instructive way,and the bisection–based solution procedures for a group of nonlinear eigenvalue equations.Besides,the explicit equations of nodal lines of the SOV methods are presented,and the relations of nodal line patterns and frequency orders are investigated.It is concluded that the highly accurate SOV methods have the same accuracy for all frequencies,the mode shapes about repeated frequencies can also be precisely captured,and the SOV methods do not have the problem of missing roots as well.

Topology Optimization for Steady-State Navier-Stokes Flow Based on Parameterized Level Set Based Method

In this paper,we consider solving the topology optimization for steady-state incompressibleNavier-Stokes problems via a new topology optimization method called parameterized level set method,which can maintain a relatively smooth level set function with a local optimality condition.The objective of topology optimization is to􀀀nd an optimal con􀀀guration of theuid and solid materials that minimizes power dissipation under a prescribeduid volume fraction constraint.An arti􀀀cial friction force is added to the Navier-Stokes equations to apply the no-slip boundary condition.Although a great deal of work has been carried out for topology optimization ofuidow in recent years,there are few researches on the topology optimization ofuidow with physical body forces.To simulate theuidow in reality,the constant body force(e.g.,gravity)is considered in this paper.Several 2D numerical examples are presented to discuss the relationships between the proposed method with Reynolds number and initial design,and demonstrate the feasibility and superiority of the proposed method in dealing with unstructuredmesh problems.Three 3D numerical examples demonstrate the proposedmethod is feasible in three-dimensional.

EFFICIENT STEADY-STATE ANALYSIS METHOD FOR CLOSED-LOOP PWM SWITCHING CONVERTERS

An analysis technique of steady state and stability for closed-loop PWM DC/DC switching converters is presented. Using this method, the closed-loop switching converter is transformed into an open-loop system. By means of the fact that in steady state, the two boundary values are equal in one switching period. The exponential matrix is evaluated by precise time-domain-integration method, and then the related curve between feedback duty cycle and the input one is obtained. Not only can the steady-state duty cycle be found from the curve, but also the stability and stable domain of the system. Compared with other methods, it features with simplicity and less calculation, and fit for numerical simulation and analysis for closed-loop switching converters. The simulation results of examples indicate the correctness of the presented method.

Wave Number Method for Three-Dimensional Steady-State Acoustic Problems

Based on the indirect Trefftz approach, a wave number method （WNM） is proposed to deal with three-dimensional steady-state acoustic problems. In the WNM, the dynamic pressure response variable is approximated by a set of wave functions, which exactly satisfy the Helmholtz equation. The set of wave functions comprise the exact solutions of the homogeneous part of the governing equations and some particular solution functions. The unknown coefficients of the wave functions can be obtained by enforcing the pressure approximation to satisfy the boundary conditions. Compared with the boundary element method （BEM）, the WNM have a smaller system matrix, and is applicable to the radiation problems since the wave functions are independent of the domain size. A 3D acoustic cavity is exemplified to show the properties of the method. The results show that the wave number method is more efficient than the BEM, and it is fairly accurate.

Application of silver nitrate colorimetric method to non-steady-state diffusion test

NT build 443, or profile fitting method, is often used to measure the diffusion coefficient of chloride in concrete. However, this method is quite laborious and needs special equipment. Colorimetric method is a quick and simple method to measure the penetration depth of chloride by spraying 0.1 mol/L silver nitrate solution. The objective of this work is to study the possibility of the use of colorimetric method in the calculation of non-steady-state diffusion coefficient. Twelve concrete mixtures with different supplementary cementitious materials and water-to-cement ratios of 0.35, 0.48 and 0.6 were used for study. According to NT build 443, the concrete specimens were immersed in 165 g/L NaC1 （2.8 mol/L） solution for 42 d. Both water-soluble （convert to free chloride） chloride and acid-soluble chloride at different layers of specimens were measured. The results show that the mean value of free chloride concentration at the color change boundary Cd was 0.306 mol/L. The surface free chloride concentration cs was obtained by profile fitting method, which was 40% lower than the chloride concentration of exposure solution after an immersion period of 42 d. Chloride diffusion coefficients obtained by the colorimetric method was not well correlated with those obtained by profile fitting method.

Isogeometric Boundary Element Method for Two-Dimensional Steady-State Non-Homogeneous Heat Conduction Problem

The isogeometric boundary element technique(IGABEM)is presented in this study for steady-state inhomogeneous heat conduction analysis.The physical unknowns in the boundary integral formulations of the governing equations are discretized using non-uniform rational B-spline(NURBS)basis functions,which are utilized to build the geometry of the structures.To speed up the assessment of NURBS basis functions,the Bezier extraction´approach is used.To solve the extra domain integrals,we use a radial integration approach.The numerical examples show the potential of IGABEM for dimension reduction and smooth integration of CAD and numerical analysis.

Goal-oriented error estimation applied to direct solution of steady-state analysis with frequency-domain finite element method

Based on the concept of the constitutive relation error along with the residuals of both the origin and the dual problems, a goal-oriented error estimation method with extended degrees of freedom is developed. It leads to the high quality locM error bounds in the problem of the direct-solution steady-state dynamic analysis with a frequency-domain finite element, which involves the enrichments with plural variable basis functions. The solution of the steady-state dynamic procedure calculates the harmonic response directly in terms of the physical degrees of freedom in the model, which uses the mass, damping, and stiffness matrices of the system. A three-dimensional finite element example is carried out to illustrate the computational procedures.

Analytical Expressions for Steady-State Concentrations of Substrate and Product in an Amperometric Biosensor with the Substrate Inhibition—The Adomian Decomposition Method

A mathematical model of an amperometric biosensor with the substrate inhibition for steady-state condition is discussed. The model is based on the system of non-stationary diffusion equation containing a non-linear term related to non-Michaelis–Menten kinetics of the enzymatic reaction. This paper presents the analytical expression of concentrations and current for all values of parameters φ2s φ2s α and β . Here the Adomian decomposition method (ADM) is used to find the analytical expressions for substrate, product concentration and current. A comparison of the analytical approximation and numerical simulation is also presented. A good agreement between theoretical predictions and numerical results is observed.

Reliable Method for Steady-State Concentrations and Current over the Diagnostic Biosensor Transducers

A mathematical modelling of diagnostic biosensors system at three basic types of enzyme kinetics is discussed in the presence of diffusion. Enzyme kinetics is adopted to be first order, Michaelis-Menten and ping-pong mechanism. In this paper, approximate analytical solutions are obtained for the non-linear equations under steady-state conditions by using the new Homotopy perturbation method. Simple and closed forms of analytical expressions for concentrations of substrate, product and co-substrate and corresponding current response have been derived for all possible values of parameters. Furthermore, the numerical simulation of the problem is also reported here by using Matlab program. Good agreement between analytical and numerical results is noted.

Structural Modal Parameter Recognition and Related Damage Identification Methods under Environmental Excitations: A Review

Modal parameters can accurately characterize the structural dynamic properties and assess the physical state of the structure.Therefore,it is particularly significant to identify the structural modal parameters according to the monitoring data information in the structural health monitoring(SHM)system,so as to provide a scientific basis for structural damage identification and dynamic model modification.In view of this,this paper reviews methods for identifying structural modal parameters under environmental excitation and briefly describes how to identify structural damages based on the derived modal parameters.The paper primarily introduces data-driven modal parameter recognition methods(e.g.,time-domain,frequency-domain,and time-frequency-domain methods,etc.),briefly describes damage identification methods based on the variations of modal parameters(e.g.,natural frequency,modal shapes,and curvature modal shapes,etc.)and modal validation methods(e.g.,Stability Diagram and Modal Assurance Criterion,etc.).The current status of the application of artificial intelligence(AI)methods in the direction of modal parameter recognition and damage identification is further discussed.Based on the pre-vious analysis,the main development trends of structural modal parameter recognition and damage identification methods are given to provide scientific references for the optimized design and functional upgrading of SHM systems.

A Disturbance Source Location Method on the Low Frequency Oscillation With Time-varying Steady-state Points

The low frequency oscillation is a serious threat to security and stability of a power grid.How to locate the disturbance source accurately is an important issue to low frequency oscillation disposal.Existing methods have poor adaptability to the low frequency oscillation with time-varying steady-state points because of the limitations in the location criterion derivation.A disturbance source location method on a low frequency oscillation with good generality is presented in the paper.Firstly,the reasons why the steady-state points are time-varying on a low frequency oscillation are analyzed.Then,based on the energy function construction form,the branch transmission energy is decomposed into state energy,reciprocating energy and dissipation energy by mathematical derivation.The flow direction of the dissipation energy shows the source and destination of the disturbance energy,and the specific location of a disturbance source can be identified according to its flow direction.Meanwhile,to meet the needs of energy calculation,a recognition method on the electrical quantities steady-state points is also presented by using the cubic spline interpolation.Simulation results show the correctness of the derivation and analysis on energy structure in the paper,and the disturbance source can be located accurately according to the dissipation energy.

A Steady-State Errors Correction Method for Eliminating Measurement Errors

This paper proposes a steady-state errors correction(SSEC)method for eliminating measurement errors.This method is based on the detections of error signal E(s)and output C(s)which generate an expected output R(s).In comparison with the conventional solutions which are based on detecting the expected output R(s)and output C(s)to obtain error signal E(s),the measurement errors are eliminated even the error might be at a significant level.Moreover,it is possible that the individual debugging by regulating the coefficient K for every member of the multiple objectives achieves the optimization of the open loop gain.Therefore,this simple method can be applied to the weak coupling and multiple objectives system,which is usually controlled by complex controller.The principle of eliminating measurement errors is derived analytically,and the advantages comparing with the conventional solutions are depicted.Based on the SSEC method analysis,an application of this method for an active power filter(APF)is investigated and the effectiveness and viability of the scheme are demonstrated through the simulation and experimental verifications.

Meshless analysis of three-dimensional steady-state heat conduction problems

Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson＇s equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method （RKPM）. A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.

Observing the steady-state visual evoked potentials with a compact quad-channel spin exchange relaxation-free magnetometer

We observed the steady-state visually evoked potential(SSVEP) from a healthy subject using a compact quad-channel potassium spin exchange relaxation-free(SERF) optically pumped magnetometer(OPM). To this end, 30 s of data were collected, and SSVEP-related magnetic responses with signal intensity ranging from 150 fT to 300 f T were observed for all four channels. The corresponding signal to noise ratio(SNR) was in the range of 3.5–5.5. We then used different channels to operate the sensor as a gradiometer. In the specific case of detecting SSVEP, we noticed that the short channel separation distance led to a strongly diminished gradiometer signal. Although not optimal for the case of SSVEP detection, this set-up can prove to be highly useful for other magnetoencephalography(MEG) paradigms that require good noise cancellation.Considering its compactness, low cost, and good performance, the K-SERF sensor has great potential for biomagnetic field measurements and brain-computer interfaces(BCI).

Steady-state responses of a belt-drive dynamical system under dual excitations

The stable steady-state periodic responses of a belt-drive system with a one-way clutch are studied. For the first time, the dynamical system is investigated under dual excitations. The system is simultaneously excited by the firing pulsations of the engine and the harmonic motion of the foundation. Nonlinear discrete-continuous equations are derived for coupling the transverse vibration of the belt spans and the rotations of the driving and driven pulleys and the accessory pulley. The nonlinear dynamics is studied under equal and multiple relations between the frequency of the fir- ing pulsations and the frequency of the foundation motion. Furthermore, translating belt spans are modeled as axially moving strings. A set of nonlinear piecewise ordinary differ- ential equations is achieved by using the Galerkin truncation. Under various relations between the excitation frequencies, the time histories of the dynamical system are numerically simulated based on the time discretization method. Further- more, the stable steady-state periodic response curves are calculated based on the frequency sweep. Moreover, the convergence of the Galerkin truncation is examined. Numer- ical results demonstrate that the one-way clutch reduces the resonance amplitude of the rotations of the driven pul- ley and the accessory pulley. On the other hand, numerical examples prove that the resonance areas of the belt spans are decreased by eliminating the torque-transmitting in the opposite direction. With the increasing amplitude of the foun- dation excitation, the damping effect of the one-way clutch will be reduced. Furthermore, as the amplitude of the firing pulsations of the engine increases, the jumping phenomena in steady-state response curves of the belt-drive system with or without a one-way clutch both occur.