A Fixed-Point Fast Sweeping WENO Method with Inverse Lax-Wendroff Boundary Treatment for Steady State of Hyperbolic Conservation Laws

Fixed-point fast sweeping WENO methods are a class of efficient high-order numerical methods to solve steady-state solutions of hyperbolic partial differential equations(PDEs).The Gauss-Seidel iterations and alternating sweeping strategy are used to cover characteristics of hyperbolic PDEs in each sweeping order to achieve fast convergence rate to steady-state solutions.A nice property of fixed-point fast sweeping WENO methods which distinguishes them from other fast sweeping methods is that they are explicit and do not require inverse operation of nonlinear local systems.Hence,they are easy to be applied to a general hyperbolic system.To deal with the difficulties associated with numerical boundary treatment when high-order finite difference methods on a Cartesian mesh are used to solve hyperbolic PDEs on complex domains,inverse Lax-Wendroff(ILW)procedures were developed as a very effective approach in the literature.In this paper,we combine a fifthorder fixed-point fast sweeping WENO method with an ILW procedure to solve steadystate solution of hyperbolic conservation laws on complex computing regions.Numerical experiments are performed to test the method in solving various problems including the cases with the physical boundary not aligned with the grids.Numerical results show highorder accuracy and good performance of the method.Furthermore,the method is compared with the popular third-order total variation diminishing Runge-Kutta(TVD-RK3)time-marching method for steady-state computations.Numerical examples show that for most of examples,the fixed-point fast sweeping method saves more than half CPU time costs than TVD-RK3 to converge to steady-state solutions.

Normalized Ground State Solutions of Nonlinear Choquard Equations with Nonconstant Potential

In this paper,we mainly focus on the following Choquard equation-{△u-V(x)(I_(a*)|u|^(p))|u|^(p-2)u=λu,x∈R^(N),u∈H^(1)(R^(N))where N≥1,λ∈R will arise as a Lagrange multiplier,0<a<N and N+a/N<p<N+a+2/N Under appropriate hypotheses on V(x),we prove that the above Choquard equation has a normalized ground state solution by utilizing variational methods.

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.

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.

Real-time embedded software testing method based on extended finite state machine

The reliability of real-time embedded software directly determines the reliability of the whole real-time embedded sys- tem, and the effective software testing is an important way to ensure software quality and reliability. Based on the analysis of the characteristics of real-time embedded software, the formal method is introduced into the real-time embedded software testing field and the real-time extended finite state machine （RT-EFSM） model is studied firstly. Then, the time zone division method of real-time embedded system is presented and the definition and description methods of time-constrained transition equivalence class （timeCTEC） are presented. Furthermore, the approaches of the testing sequence and test case generation are put forward. Finally, the proposed method is applied to a typical avionics real- time embedded software testing practice and the examples of the timeCTEC, testing sequences and test cases are given. With the analysis of the testing result, the application verification shows that the proposed method can effectively describe the real-time embedded software state transition characteristics and real-time requirements and play the advantages of the formal methods in accuracy, effectiveness and the automation supporting. Combined with the testing platform, the real-time, closed loop and automated simulation testing for real-time embedded software can be realized effectively.

A circular zone counting method of identifying a Duffing oscillator state transition and determining the critical value in weak signal detection

Identifying state transition and determining the critical value of the Duffing oscillator are crucial to indicating external signal existence and have a great influence on detection accuracy in weak signal detection. A circular zone counting （CZC） method is proposed in this paper, by combining the Duffing oscillator＇s phase trajectory feature and numerical calculation for quickly and accurately identifying state transition and determining the critical value, to realize a high- efficiency weak signal detection. Detailed model analysis and method construction of the CZC method are introduced. Numerical experiments into the reliability of the proposed CZC method compared with the maximum Lyapunov exponent （MLE） method are carried out. The CZC method is demonstrated to have better detecting ability than the MLE method, and furthermore it is simpler and clearer in calculation to extend to engineering application.

ONE-CELL-STATE METHOD FOR DETERMINATION OF ELECTRONIC STRUCTURE OF INTERMETALLIC COMPOUNDS

ＩＮＴＲＯＤＵＣＴＩＯＮＩｎｒｅｃｅｎｔｙｅａｒｓ，ｔｈｅＰａｕｌｉｎｇ’ｓｖａｌｅｎｃｅｂｏｎｄ（ＶＢ）ｔｈｅｏｒｙｈａｓｂｅｎｍａｄｅａｇｒｅａｔｐｒｏｇｒｅｓｉｎｄｅｓｉｇｎｉｎｇａｎｄｄｅｔｅｒｍｉｎｉｎｇｔｈｅａｔｏｍｉｃｓｔａｔｅ，ａｎｄ...

Application of the state space model to the system of wave energy conversion Analytical method for wave forces on singleoscillating rectangular buoy

A new analytical method is proposed to analyze the force acting on a rectangular oscillating buoy due to linear waves.In the method a new analytical expression for the diffraction velocity potential is obtained first by use of theeigenfunction expansion method and then the wave excitation force is calculated by use of the known incident wavepotential and the diffraction potential. Compared with the classical analytical method, it can be seen that the presentmethod is simpler for a two-dimensional problem due to the comparable effort needed for the computation ofdiffraction potential and for that of radiated potential. To verify the correctness of the method, a classical example inthe reference is recomputed and the obtained results are in good accordance with those by use of other methods,which shows that the present method is correct.

Mixed H_2/H_∞ State Feedback Attitude Control of Microsatellite Based on Extended LMI Method

For the appearance of the additive perturbation of controller gain when the controller parameter has minute adjustment at the initial running stage of system,to avoid the adverse effects,this paper investigates the mixed H_2/H_∞ state feedback attitude control problem of microsatellite based on extended LMI method.Firstly,the microsatellite attitude control system is established and transformed into corresponding state space form.Then,without the equivalence restriction of the two Lyapunov variables of H_2 and H∞performance,this paper introduces additional variables to design the mixed H_2/H_∞ control method based on LMI which can also reduce the conservatives.Finally,numerical simulations are analyzed to show that the proposed method can make the satellite stable within 20 s whether there is additive perturbation of the controller gain or not.The comparative analysis of the simulation results between extended LMI method and traditional LMI method also demonstrates the effectiveness and feasibility of the proposed method in this paper.

Present state,basic theories,methods and progresses of investigation and assessment on marine hazardous geology in China

The achievement progresses of investigation and studies on marine hazardous geology are summarized and presentsd in the late 20 century in China. The importance, research value and present-day studies of marine hazardous geology, a newly developing branch of geoscience, are well expatiated. Several often confused concepts and theories are explained and redefined here. The comment on the means of investigations, assessment of marine hazardous geology, as well as its evolution, innovation, existing questions and future tasks are also introduced and presented. The concepts of 'hazard geology', geohazard', 'map of marine hazard geology', 'integrated evaluaton on seafloor stablity' are respectively discussed, including their definition, research objects, methods and contents. The types and classification of marine hazardous geology, principles and methods of marine hazardous geology map compilation, the assessment methods and models of marine hazardous geology environment and seafloor stability and so on are also discussed.

BCS Ground State and XXZ Antiferromagnetic Model as SU（2）,SU（1,1） Coherent States：AN Algebraic Diagonalization Method

An algebraic diagonalization method is proposed. As two examples, the Hamiltonians of BCS ground stateunder mean-field approximation and XXZ antiferromagnetic model in linear spin-wave frame have been diagonalized byusing SU(2), SU(1,1) Lie algebraic method, respectively. Meanwhile, the eigenstates of the above two models are revealedto be SU(2), SU(1,1) coherent states, respectively. The relation between the usual Bogoliubov Valatin transformationand the algebraic method in a special case is also discussed.

An improved data space inversion method to predict reservoir state fields via observed production data

A data-space inversion(DSI)method has been recently proposed and successfully applied to the history matching and production prediction of reservoirs.Based on Bayesian theory,DSI can directly and effectively obtain good posterior flow predictions without inversion of geological parameters of reservoir model.This paper presents an improved DSI method to fast predict reservoir state fields(e.g.saturation and pressure profiles)via observed production data.Firstly,a large number of production curves and state data are generated by reservoir model simulation to expand the data space of original DSI.Then,efficient history matching only on the observed production data is carried out via the original DSI to obtain related parameters which reflects the weight of the real reservoir model relative to prior reservoir models.Finally,those parameters are used to predict the oil saturation and pressure profiles of the real reservoir model by combining large amounts of state data of prior reservoir models.Two examples including conventional heterogeneous and unconventional fractured reservoir are implemented to test the performances of predicting saturation and pressure profiles of this improved DSI method.Besides,this method is also tested in a real field and the obtained results show the high computational efficiency and high accuracy of the practical application of this method.

Probabilistic limit state design methodof the axial resistance of post grouting pile

The reliability of post grouting pile axial resistance was studied by proposing a design method for its probabilistic limit state,which is represented by the partial coefficients of load,end,and side resistance.The hyperbolic,modified hyperbolic,and polynomial models were employed to predict the ultimate bearing capacity of test piles that were not loaded to damage in field tests.The results were used for the calculation and calibration of the reliability index.The reliability of the probabilistic limit state design method was verified by an engineering case.The results show that the prediction results obtained from the modified hyperbolic model are closest to those obtained through the static load test.The proposed corresponding values of total,side,and end resistance partial coefficients are 1.84,1.66,and 2.73 when the dead and live load partial coefficients are taken as 1.1 and 1.4,respectively.Meanwhile,the corresponding partial coefficients of total,side,and end resistance are 1.70,1.56,and 2.34 when the dead and live load partial coefficients are taken as 1.2 and 1.4,respectively.

Power System State Estimation Solution With Zero Injection Constraints Using Modified Newton Method and Fast Decoupled Method in Polar Coordinate

如何保证零注入节点的注入功率在状态估计结果中严格为0是电力系统状态估计研究中的重要问题。在直角坐标下,由于零注入约束为线性约束,可使用修正牛顿法来有效地解决这一问题。因此,借鉴直角坐标下修正牛顿法的思路,提出了极坐标下的修正牛顿法和修正快速解耦估计。这些方法的计算流程与传统的极坐标下的牛顿法和快速解耦估计非常相似,计算速度与大权重法相当,同时能够保证零注入约束严格满足。仿真结果验证了所得结论。

DYNAMIC RESPONSE OF ELASTIC RECTANGULAR PLATES BY SPLINE STATE VARIABLE METHOD

Application of spline element and state space method for analysis of dynamic response of elastic rectangular plates is presented. The spline element method is used for space domain and the state space method in control theory of system is used for time domain. A state variable recursive scheme is developed, then the dynamic response of structure can he calculated directly. Several numerical examples are given. The results which are presented to demonstrate the accuracy and efficiency of the present method are quite satisfactory.

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.

Barut-Girardello and Gilmore-Perelomov coherent states for pseudoharmonic oscillators and their nonclassical properties: Factorization method

In this paper we try to introduce the ladder operators associated with the pseudoharmonic oscillator, after solving the corresponding Schrrdinger equation by using the factorization method. The obtained generalized raising and lowering operators naturally lead us to the Dirac representation space of the system which is much easier to work with, in comparison to the functional Hilbert space. The SU（1, 1） dynamical symmetry group associated with the considered system is exactly established through investigating the fact that the deduced operators satisfy appropriate commutation relations. This result enables us to construct two important and distinct classes of Barut-Girardello and Gilmore-Perelomov coherent states associated with the system. Finally, their identities as the most important task are exactly resolved and some of their nonclassical properties are illustrated, numerically.

STATE SPACE METHODS FOR DECENTRALIZED H_∞ CONTROL

This paper presents state space methods for decentralized Hoe control, which contain two respects: a parametrization approach and an iterative algorithm. For large scale systems with N subsystems, decentralized Hoe con trollers can be derived by a parametrization result for centralized Her: controllers and designed by an iterative algorithm with structured constraint to the controllers.

Generalized Lanczos method for systematic optimization of tensor network states

We propose a generalized Lanczos method to generate the many-body basis states of quantum lattice models using tensor-network states （TNS）. The ground-state wave function is represented as a linear superposition composed from a set of TNS generated by Lanczos iteration. This method improves significantly the accuracy of the tensor-network algorithm and provides an effective way to enlarge the maximal bond dimension of TNS. The ground state such obtained contains significantly more entanglement than each individual TNS, reproducing correctly the logarithmic size dependence of the entanglement entropy in a critical system. The method can be generalized to non-Hamiltonian systems and to the calculation of low-lying excited states, dynamical correlation functions, and other physical properties of strongly correlated systems.