Numerical method of the Riemann problem for two-dimensional multi-fluid flows with general equation of state

Based on the classical Roe method, we develop an interface capture method according to the general equation of state, and extend the single-fluid Roe method to the two-dimensional （2D） multi-fluid flows, as well as construct the continuous Roe matrix for the whole flow field. The interface capture equations and fluid dynamic conservative equations are coupled together and solved by using any high-resolution schemes that usually suit for the single-fluid flows. Some numerical examples are given to illustrate the solution of 1D and 2D multi-fluid Riemann problems.

Some New Systems of Exponentially General Equations

Some new systems of exponentially general equations are introduced and investigated, which can be used to study the odd-order, non-positive and nonsymmetric exponentially boundary value problems. Some important and interesting results such as Riesz-Frechet representation theorem, Lax-Milgram lemma and system of absolute values equations can be obtained as special cases. It is shown that the system of exponentially general equations is equivalent to nonlinear optimization problem. The auxiliary principle technique is used to prove the existence of a solution to the system of exponentially general equations. This technique is also used to suggest some new iterative methods for solving the system of the exponentially general equations. The convergence analysis of the proposed methods is analyzed. Ideas and techniques of this paper may stimulate further research.

An Interface-Capturing Method for Resolving Compressible Two-Fluid Flows with General Equation of State

In this study,a stable and robust interface-capturing method is developed to resolve inviscid,compressible two-fluid flows with general equation of state(EOS).The governing equations consist of mass conservation equation for each fluid,momentum and energy equations for mixture and an advection equation for volume fraction of one fluid component.Assumption of pressure equilibrium across an interface is used to close the model system.MUSCL-Hancock scheme is extended to construct input states for Riemann problems,whose solutions are calculated using generalized HLLC approximate Riemann solver.Adaptive mesh refinement(AMR)capability is built into hydrodynamic code.The resulting method has some advantages.First,it is very stable and robust,as the advection equation is handled properly.Second,general equation of state can model more materials than simple EOSs such as ideal and stiffened gas EOSs for example.In addition,AMR enables us to properly resolve flow features at disparate scales.Finally,this method is quite simple,time-efficient and easy to implement.

Solution of general dynamic equation for nanoparticles in turbulent flow considering fluctuating coagulation

A new averaged general dynamic equation （GDE） for nanoparticles in the turbulent flow is derived by considering the combined effect of convection, Brownian diffusion, turbulent diffusion, turbulent coagulation, and fluctuating coagulation. The equation is solved with the Taylor-series expansion moment method in a turbulent pipe flow. The experiments are performed. The numerical results of particle size distribu- tion correlate well with the experimental data. The results show that, for a turbulent nanoparticulate flow, a fluctuating coagulation term should be included in the averaged particle GDE. The larger the Schmidt number is and the lower the Reynolds number is, the smaller the value of ratio of particle diameter at the outlet to that at the inlet is. At the outlet, the particle number concentration increases from the near-wall region to the near-center region. The larger the Schmidt number is and the higher the Reynolds num- ber is, the larger the difference in particle number concentration between the near-wall region and near-center region is. Particle polydispersity increases from the near-center region to the near-wall region. The particles with a smaller Schmidt number and the flow with a higher Reynolds number show a higher polydispersity. The degree of particle polydispersity is higher considering fluctuating coagulation than that without considering fluctuating coagulation.

EXACT SOLUTIONS FOR GENERAL VARIABLE-COEFFICIENT KdV EQUATION

By asing the nonclassical method of symmetry reductions, the exact solutions for general variable coefficient KdV equation with dissipative loss and nonuniformity terms are obtained. When the dissipative loss and nonuniformity terms don't exist, the multisoliton solutions are found and the corresponding Painleve II type equation for the variable coefficient KdV equation is given.

The extended symmetry approach for studying the general Korteweg-de Vries-type equation

The extended symmetry approach is used to study the general Korteweg-de Vries-type （KdV-type） equation. Several variable-coefficient equations are obtained. The solutions of these resulting equations can be constructed by the solutions of original models if their solutions are well known, such as the standard constant coefficient KdV equation and the standard compound KdV--Burgers equation, and so on. Then any one of these variable-coefficient equations can be considered as an original model to obtain new variable-coefficient equations whose solutions can also be known by means of transformation relations between solutions of the resulting new variable-coefficient equations and the original equation.

THE GENERALIZED HYPERSTABILITY OF GENERAL LINEAR EQUATION IN QUASI-2-BANACH SPACE

In this paper,we study the hyperstability for the general linear equation f(ax+by)=Af(x)+Bf(y)in the setting of complete quasi-2-Banach spaces.We first extend the main fixed point result of Brzdek and Ciepliński(Acta Mathematica Scientia,2018,38 B(2):377-390)to quasi-2-Banach spaces by defining an equivalent quasi-2-Banach space.Then we use this result to generalize the main results on the hyperstability for the general linear equation in quasi-2-Banach spaces.Our results improve and generalize many results of literature.

Sub-harmonic Resonance Solutions of Generalized Strongly Nonlinear Van der Pol Equation with Parametric and External Excitations

In this paper a modified L-P method and multiple scale method are used to solve sub-harmonic resonance solutions of strong and nonlinear resonance of general Van der Pol equation with parametric and external excitations by parametric transformation. Bifurcation response equation and transition sets of sub-harmonic resonance with strong nonlinearity of general Van der Pol equation with parametric and external excitation are worked out.Besides, transition sets and bifurcation graphs are drawn to help to analysis the problems theoretically. Conclusions show that the transition sets of general and nonlinear Van der Pol equation with parametric and external excitations are more complex than those of general and nonlinear Van der Pol equation only with parametric excitation, which is helpful for the qualitative and quantitative reference for engineering and science applications.

GENERAL DYNAMIC EQUATION AND DYNAMICAL CHARACTERISTICS OF VISCOELASTIC TIMOSHENKO BEAMS

In this paper, a governing differential equation of viscoelastic Timoshenko beam including both extension and shear viscosity is developed in the time domain by direct method. To measure the complex moduli and three parameters of standard linear solid, the forced vibration technique of beam is successfully used for PCL and PMMA specimens. The dynamical characteristics of viscoelastic Timoshenko beams, especially the damping properties, are derived from a considerable number of numerical computations. The analyses show that the viscosity of materials has great influence on dynamical characteristics of structures, especially on damping, and the standard linear solid model is the better one for describing the dynamic behavior of high viscous materials.

DEFORMATION OF STRUCTURE AND SPECTRUM OF EVOLUTION EQUATIONS

In this paper, we study the general structure of evolution equations of the AKNS eigenvalue problem q(x,t), r(x,t) with the spectrum varying asand AV BV CV are all positive or negative power polynomials of where q, r are not limited with any additional conditions at infinity.

L^P SOLUTION OF GENERAL MEAN-FIELD BSDES WITH CONTINUOUS COEFFICIENTS

In this paper we consider one dimensional mean-field backward stochastic differential equations(BSDEs)under weak assumptions on the coefficient.Unlike[3],the generator of our mean-field BSDEs depends not only on the solution(Y,Z)but also on the law PY of Y.The first part of the paper is devoted to the existence and uniqueness of solutions in Lp,1<p≤2,where the monotonicity conditions are satisfied.Next,we show that if the generator/is uniformly continuous in(μ,y,z),uniformly with respect to(t,ω) and if the terminal valueξbelongs to Lp(Ω,F,P)with 1<p≤2,the mean-field BSDE has a unique Lp solution.

The Conservation Law and the Radiation of General Relativity

This paper presents the conservation law of Einstein's General Relativity, in whichthe continuity equation of matter and gravitational field is implicitly derived from theequation of motion of matter. Although the obtained energy-momentum tensor is thesame as the Landau- Lifshitz pseudotensor. the physical and conceptual foundation aredifferent. Two alternative methods to obtain the gravitational radiation are proposed inthis paper as well The radiation will be derived from the equation of motion of matterand geodesic equation separately.

A New Framework of Convergence Analysis for Solving the General Nonlinear Schrodinger Equation using the Fourier Pseudo-Spectral Method in Two Dimensions

This paper aims to build a new framework of convergence analysis of conservative Fourier pseudo-spectral method for the general nonlinear Schr¨odinger equation in two dimensions,which is not restricted that the nonlinear term is mere cubic.The new framework of convergence analysis consists of two steps.In the first step,by truncating the nonlinear term into a global Lipschitz function,an alternative numerical method is proposed and proved in a rigorous way to be convergent in the discrete L2 norm;followed in the second step,the maximum bound of the numerical solution of the alternative numerical method is obtained by using a lifting technique,as implies that the two numerical methods are the same one.Under our framework of convergence analysis,with neither any restriction on the grid ratio nor any requirement of the small initial value,we establish the error estimate of the proposed conservative Fourier pseudo-spectral method,while previous work requires the certain restriction for the focusing case.The error bound is proved to be of O(h^(r)+t^(2))with grid size h and time step t.In fact,the framework can be used to prove the unconditional convergence of many other Fourier pseudo-spectral methods for solving the nonlinear Schr¨odinger-type equations.Numerical results are conducted to indicate the accuracy and efficiency of the proposed method,and investigate the effect of the nonlinear term and initial data on the blow-up solution.

Oscillatory Behavior of Solutions of General Third Order Mixed Neutral Difference Equations

The objective of this paper is to study the oscillatory and asymptotic properties of the general mixed type third order neutral difference equation of the form △（aαn△2（xn＋ bnxn-τ1＋ cnxn＋τ2）） ＋ qnx^a n＋1-σ＋ pnx^βn＋1＋=σ2 0,where {an}, {bn}, {cn}, {qn} and {pn} are positive real sequences, both α and β are ratios of odd positive integers, τ1, τ2, σ1 and σ2 are positive integers. We establish some sufficient conditions which ensure all solutions are either oscillatory or converge to zero.

Suboptimal Robust Stabilization of Discrete-time Mismatched Nonlinear System

This paper proposes a discrete-time robust control technique for an uncertain nonlinear system. The uncertainty mainly affects the system dynamics due to mismatched parameter variation which is bounded by a predefined known function. In order to compensate the effect of uncertainty, a robust control input is derived by formulating an equivalent optimal control problem for a virtual nominal system with a modified costfunctional. To derive the stabilizing control law for a mismatched system, this paper introduces another control input named as virtual input. This virtual input is not applied directly to stabilize the uncertain system, rather it is used to define a sufficient condition. To solve the nonlinear optimal control problem, a discretetime general Hamilton-Jacobi-Bellman(DT-GHJB) equation is considered and it is approximated numerically through a neural network(NN) implementation. The approximated solution of DTGHJB is used to compute the suboptimal control input for the virtual system. The suboptimal inputs for the virtual system ensure the asymptotic stability of the closed-loop uncertain system. A numerical example is illustrated with simulation results to prove the efficacy of the proposed control algorithm.

Bivariate association analysis for quantitative traits using generalized estimation equation

Quantitative traits often underlie risk for complex diseases. Many studies collect multiple correlated quantitative phenotypes and perform univariate analyses on each of them respectively. However, this strategy may not be powerful and has limitations to detect plei- otropic genes that may underlie correlated quantitative traits. In addition, testing multiple traits individually will exacerbate perplexing problem of multiple testing. In this study, generalized estimating equation 2 （GEE2） is applied to association mapping of two correlated quantitative traits. We suppose that a quantitative trait locus is located in a chromosome region that exerts pleiotropic effects on multiple quantitative traits. In that region, multiple SNPs are genotyped. Genotypes of these SNPs and the two quantitative traits affected by a causal SNP were simulated under various parameter values： residual correlation coefficient between two traits, causal SNP heritability, minor allele frequency of the causal SNP, extent of linkage disequilibrium with the causal SNP, and the test sample size. By power ana- lytical analyses, it is showed that the bivariate method is generally more powerful than the univariate method. This method is robust and yields false-positive rates close to the pre-set nominal significance level. Our real data analyses attested to the usefulness of the method.

Error Estimates ofMixedMethods forOptimal Control Problems Governed by General Elliptic Equations

In this paper,we investigate the error estimates of mixed finite element methods for optimal control problems governed by general elliptic equations.The state and co-state are approximated by the lowest order Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions.We derive L2 and H−1-error estimates both for the control variable and the state variables.Finally,a numerical example is given to demonstrate the theoretical results.

Study on General Governing Equations of Computational Heat Transfer and Fluid Flow

The governing equations for heat transfer and fluid flow are often formulated in a general formfor the simplification of discretization and programming,which has achieved great success in thermal science and engineering.Based on the analysis of the popular general form of governing equations,we found that energy conservation cannot be guaranteed when specific heat capacity is not constant,which may lead to unreliable results.A new concept of generalized density is put forward,based on which a new general form of governing equations is proposed to guarantee energy conservation.A number of calculation examples have been employed to verify validation and feasibility.

Efcient Iterative Solutions to General Coupled Matrix Equations

Linear matrix equations are encountered in many systems and control applications.In this paper,we consider the general coupled matrix equations（including the generalized coupled Sylvester matrix equations as a special case）l t=1EstYtFst = Gs,s = 1,2,···,l over the generalized reflexive matrix group（Y1,Y2,···,Yl）.We derive an efcient gradient-iterative（GI） algorithm for fnding the generalized reflexive solution group of the general coupled matrix equations.Convergence analysis indicates that the algorithm always converges to the generalized reflexive solution group for any initial generalized reflexive matrix group（Y1（1）,Y2（1）,···,Yl（1））.Finally,numerical results are presented to test and illustrate the performance of the algorithm in terms of convergence,accuracy as well as the efciency.

Efficient Chebyshev Spectral Method for Solving Linear Elliptic PDEs Using Quasi-Inverse Technique

We present a systematic and efficient Chebyshev spectral method using quasiinverse technique to directly solve the second order equation with the homogeneous Robin boundary conditions and the fourth order equation with the first and second boundary conditions.The key to the efficiency of the method is to multiply quasiinverse matrix on both sides of discrete systems,which leads to band structure systems.We can obtain high order accuracy with less computational cost.For multi-dimensional and more complicated linear elliptic PDEs,the advantage of this methodology is obvious.Numerical results indicate that the spectral accuracy is achieved and the proposed method is very efficient for 2-D high order problems.