Average vector field methods for the coupled Schrdinger KdV equations

The energy preserving average vector field （AVF） method is applied to the coupled Schr6dinger-KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction discretization. In order to accelerate our simulation, the split-step technique is used. The numerical experiments show that the non-splitting scheme and splitting scheme are both effective, and have excellent long time numerical behavior. The comparisons show that the splitting scheme is faster than the non-splitting scheme, but it is not as good as the non-splitting scheme in preserving the invariants.

AN AVERAGING PRINCIPLE FOR STOCHASTIC DIFFERENTIAL DELAY EQUATIONS DRIVEN BY TIME-CHANGED LéVY NOISE

In this paper,we aim to derive an averaging principle for stochastic differential equations driven by time-changed Lévy noise with variable delays.Under certain assumptions,we show that the solutions of stochastic differential equations with time-changed Lévy noise can be approximated by solutions of the associated averaged stochastic differential equations in mean square convergence and in convergence in probability,respectively.The convergence order is also estimated in terms of noise intensity.Finally,an example with numerical simulation is given to illustrate the theoretical result.

EQUATION OF STATE CALCULATIONS FOR HOT, DENSE MATTER AT ARBITRARY DENSITIES AND TEMPERATURES

Within the approximations of spherical lattice cell, central-field, and relativistic Fermi statis- tics, an algorithm with average atom model is presented to calculate the electronic energy levels and equation of state for hot and dense matter at arbitrary densities and temperatures. Choosing Zink's analytical potential as initial potential, we have solved the Dirac-Slater equation which satisfies the Weigner-Seitz boundary condition. The electronic energy bands are not taken into account. Tak- ing energy level degeneracy as a continuous function of density, we have considered the pressure ionization effects for highly dense matter. Results for ^(13)Al atom are shown.

OSCILLATION CRITERIA FOR SECOND ORDER NONLINEAR DIFFERENTIAL EQUATION WITH DAMPING

By the generalized Riccati transformation and the integral averaging technique, some sufficient conditions of oscillation of the solutions for second order nonlinear differential equations with damping were discussed.Some sufficient oscillation criteria for previous equations were built up.Some oscillation criteria have been expanded and strengthened in some other known results.

Necessary Conditions for the Application of Moving Average Process of Order Three

Invertibility is one of the desirable properties of moving average processes. This study derives consequences of the invertibility condition on the parameters of a moving average process of order three. The study also establishes the intervals for the first three autocorrelation coefficients of the moving average process of order three for the purpose of distinguishing between the process and any other process (linear or nonlinear) with similar autocorrelation structure. For an invertible moving average process of order three, the intervals obtained are , -0.5ρ2ρ1<0.5.

Study of modeling unsteady blade row interaction in a transonic compressor stage part 2:influence of deterministic correlations on time-averaged flow prediction

The average-passage equation system （APES） provides a rigorous mathematical framework for account- ing for the unsteady blade row interaction through multistage compressors in steady state environment by introducing de- terministic correlations （DC） that need to be modeled to close the equation system. The primary purpose of this study was to provide insight into the DC characteristics and the in- fluence of DC on the time-averaged flow field of the APES. In Part 2 of this two-part paper, the influence of DC on the time-averaged flow field was systematically studied; Several time-averaging computations boundary conditions and DC were conducted with various for the downstream stator in a transonic compressor stage, by employing the CFD solver developed in Part 1 of this two-part paper. These results were compared with the time-averaged unsteady flow field and the steady one. The study indicat;d that the circumferential- averaged DC can take into account major part of the unsteady effects on spanwise redistribution of flow fields in compres- sors. Furthermore, it demonstrated that both deterministic stresses and deterministic enthalpy fluxes are necessary to reproduce the time-averaged flow field.

Dynamic flight stability of hovering model insects:theory versus simulation using equations of motion coupled with Navier-Stokes equations

In the present paper, the longitudinal dynamic flight stability properties of two model insects are predicted by an approximate theory and computed by numerical sim- ulation. The theory is based on the averaged model （which assumes that the frequency of wingbeat is sufficiently higher than that of the body motion, so that the flapping wings＇ degrees of freedom relative to the body can be dropped and the wings can be replaced by wingbeat-cycle-average forces and moments）; the simulation solves the complete equations of motion coupled with the Navier-Stokes equations. Comparison between the theory and the simulation provides a test to the validity of the assumptions in the theory. One of the insects is a model dronefly which has relatively high wingbeat frequency （164 Hz） and the other is a model hawkmoth which has relatively low wingbeat frequency （26 Hz）. The results show that the averaged model is valid for the hawkmoth as well as for the dronefly. Since the wingbeat frequency of the hawkmoth is relatively low （the characteristic times of the natural modes of motion of the body divided by wingbeat period are relatively large） compared with many other insects, that the theory based on the averaged model is valid for the hawkmoth means that it could be valid for many insects.

HOMOGENIZATION FOR NONLINEAR SCHRODINGER EQUATIONS WITH PERIODIC NONLINEARITY AND DISSIPATION IN FRACTIONAL ORDER SPACES

We study the nonlinear SchrSdinger equation with time-oscillating nonlinearity and dissipation originated from the recent studies of Bose-Einstein condensates and optical systems which reads iψt＋△ψ＋Ф（ωt）｜ψ｜αψ＋iξ （ωt）ψ= 0. Under some conditions, we show that as ω→∞ , the solution ψω will locally converge to the solution of the averaged equation iψt＋△ψ＋Ф（ωt）｜ψ｜αψ＋iξ （ωt）ψ= 0 with the same initial condition in Lq（（0, T）, B-S/T,2） for all admissible pairs （q, r）, where T∈ （0, Tmax）. We also show that if the dissipation coefficient ξ0 large enough, then, ψω is global if w is sufficiently large and ψω converges to ψ in Lq（（0, ∞）, B-S/T,2）, for all admissible pairs （q, r）. In particular, our results hold for both subcritical and critical nonlinearities.

A STUDY OF THE STATIC AND GLOBAL BIFURCATIONS FOR DUFFING EQUATION

In this paper, the static and global bifurcations of the forced Duffing equation have been studied by means of the averaged system. Bifurcation condition has been obtained in the whole parametric space. The change of the phase plane structure has been investigated.

Analytical solution of the Boltzmann-Poisson equation and its application to MIS tunneling junctions

In order to consider quantum transport under the influence of an electron-electron （e-e） interaction in a mesoscopic conductor,the Boltzmann equation and Poisson equation are investigated jointly.The analytical expressions of the distribution function for the Boltzmann equation and the self-consistent average potential concerned with e-e interaction are obtained,and the dielectric function appearing in the self-consistent average potential is naturally generalized beyond the Thomas-Fermi approximation.Then we apply these results to the tunneling junctions of a metal-insulator-semiconductor （MIS） in which the electrons are accumulated near the interface of the semiconductor,and we find that the e-e interaction plays an important role in the transport procedure of this system. The electronic density,electric current as well as screening Coulombic potential in this case are studied,and we reveal the time and position dependence of these physical quantities explicitly affected by the e-e interaction.

Oscillation theorems for second-order nonlinear perturbed differential equations with damping

A class of second-order nonlinear damped perturbed differential equations is considered and its oscillation theorems are studied.These theorems are more general and deal with the cases which are not covered by the known criteria.Particularly,these criteria extend and unify some existing results.An example is given to verify the results.

A TVD-WAF-based hybrid finite volume and finite difference scheme for nonlinearly dispersive wave equations

A total variation diminishing-weighted average flux （TVD-WAF）-based hybrid numerical scheme for the enhanced version of nonlinearly dispersive Boussinesq-type equations was developed. The one-dimensional governing equations were rewritten in the conservative form and then discretized on a uniform grid. The finite volume method was used to discretize the flux term while the remaining terms were approximated with the finite difference method. The second-order TVD-WAF method was employed in conjunction with the Harten-Lax-van Leer （HLL） Riemann solver to calculate the numerical flux, and the variables at the cell interface for the local Riemann problem were reconstructed via the fourth- order monotone upstream-centered scheme for conservation laws （MUSCL）. The time marching scheme based on the third-order TVD Runge- Kutta method was used to obtain numerical solutions. The model was validated through a series of numerical tests, in which wave breaking and a moving shoreline were treated. The good agreement between the computed results, documented analytical solutions, and experimental data demonstrates the correct discretization of the governing equations and high accuracy of the proposed scheme, and also conforms the advantages of the proposed shock-capturing scheme for the enhanced version of the Boussinesq model, including the convenience in the treatment of wave breaking and moving shorelines and without the need for a numerical filter.

Numerical research on flow and thermal transport in cooling pool of electrical power station using three depth-averaged turbulence models

This paper describes a numerical simulation of thermal discharge in the cooling pool of an electrical power station, aiming to develop general-purpose computational programs for grid generation and flow/pollutant transport in the complex domains of natural and artificial waterways. Three depth-averaged two-equation closure turbulence models, k-ε, k- w, and k- w, were used to close the quasi three-dimensional hydrodynamic model. The k- w model was recently established by the authors and is still in the testing process. The general-purpose computational programs and turbulence models will be involved in a software that is under development. The SIMPLE （Semi-Implicit Method for Pressure-Linked Equation） algorithm and multi-grid iterative method are used to solve the hydrodynamic fundamental governing equations, which are discretized on non-orthogonal boundary-fitted grids with a variable collocated arrangement. The results calculated with the three turbulence models were compared with one another. In addition to the steady flow and thermal transport simulation, the unsteady process of waste heat inpouring and development in the cooling pool was also investigated.

High Accuracy Arithmetic Average Discretization for Non-Linear Two Point Boundary Value Problems with a Source Function in Integral Form

In this article, we report the derivation of high accuracy finite difference method based on arithmetic average discretization for the solution of Un=F(x,u,u′)+∫K(x,s)ds , 0 x s < 1 subject to natural boundary conditions on a non-uniform mesh. The proposed variable mesh approximation is directly applicable to the integro-differential equation with singular coefficients. We need not require any special discretization to obtain the solution near the singular point. The convergence analysis of a difference scheme for the diffusion convection equation is briefly discussed. The presented variable mesh strategy is applicable when the internal grid points of the solution space are both even and odd in number as compared to the method discussed by authors in their previous work in which the internal grid points are strictly odd in number. The advantage of using this new variable mesh strategy is highlighted computationally.

Average Rainfall Estimation: Methods Performance Comparison in the Brazilian Semi-Arid

Considering the rainfall’s importance in hydrological modeling, the objective of this study was the performance comparison, in convergence terms, of techniques often used to estimate the average rainfall over an area: Thiessen Polygon (TP) Method;Reciprocal Distance Squared (RDS) Method;Kriging Method (KM) and Multiquadric Equations (ME) Method. The comparison was done indirectly, using GORE and BALANCE index to assess the convergence results from each method by increasing the rain gauges density in a region, through six scenarios. The Coremas/Mae D’água Watershed employed as study area, with an area of 8385 km2, is situated on Brazilian semi-arid. The results showed the TP, as RDS and ME techniques to be employed successfully to obtain the average rainfall over an area, highlighting the MEM. On the other hand, KM, using two variograms models, had an unstable behavior, pointing the prior study of data and variogram’s choice as a need to practical applying.

Modeling continuous traffic flow with the average velocity effect of multiple vehicles ahead on gyroidal roads

In the future connected vehicle environment,the information of multiple vehicles ahead can be readily collected in real-time,such as the velocity or headway,which provides more opportunities for information exchange and cooperative control.Meanwhile,gyroidal roads are one of the fundamental road patterns prevalent in mountainous areas.To effectively control the system,it is therefore significant to explore the evolution mechanism of traffic flow on gyroidal roads under a connected vehicle environment.In this paper,we present a new continuum model with the average velocity of multiple vehicles ahead on gyroidal roads.The stability criterion and KdV-Burger equation are deduced via linear and nonlinear stability analysis,respectively.Solving the above KdV-Burger equation yields the density wave solution,which explores the formation and propagation property of traffic jams near the neutral stability curve.Simulation examples verify that the model can reproduce complex phenomena,such as shock waves and rarefaction waves.The analysis of the local cluster effect shows that the number of vehicles ahead and the radius information,and the slope information of gyroidal roads can exert a great influence on traffic jams.The effect of the first and second terms are positive,while the last term is negative.

A Scalar Acoustic Equation for Gases, Liquids, and Solids, Including Viscoelastic Media

The work deals with a mathematical model for real-time acoustic monitoring of material parameters of media in multi-state viscoelastic engineering systems continuously operating in irregular external environments (e.g., wind turbines in cold climate areas, aircrafts, etc.). This monitoring is a high-reliability time-critical task. The work consistently derives a scalar wave PDE of the Stokes type for the non-equilibrium part (NEP) of the average normal stress in a medium. The explicit expression for the NEP of the corresponding pressure and the solution-adequateness condition are also obtained. The derived Stokes-type wave equation includes the stress relaxation time and is applicable to gases, liquids, and solids.

Macroscopic Equation and Its Application for Free Ways

Effective transportation systems lead to the efficient movement of goods and people, which significantly contribute to the quality of life in every society. In the heart of every economic and social development, there is always a transportation system. Mathematically the problem of modeling vehicle traffic flow can be solved at two main observation scales: The microscopic and the macroscopic levels. In the microscopic level, every vehicle is considered individually, and therefore, for every vehicle, we have an equation that is usually an ordinary differential equation (ODE). At a macroscopic level, we use from the dynamics models, where we have a system of partial differential equation, which involves variables such as density, speed, and flow rate of traffic stream with respect to time and space. Therefore, considering above content, this study has tried to compare solution of equation of macroscopic flow considering linear form (speed-density) and applying boundary condition that resulting to form solved is non-linear one-order partial differential equation (sharpy method) with non-linear assuming (speed and density) and consequently homographic nonlinear relation (speed-density). The recent case clearly gives more significant speeds than linear case of speed and density that can be a good scientific basis. In terms of safety for accidents and traffic signal, just as a reminder, but it is resulted of the reality that generally solutions of partial differential equations can have different forms. Therefore, the solution of partial differential equation (macroscopic flow) can have different answers and solutions so that all of these solutions apply in PDE (equation of macroscopic flow). Thus, under this condition, we can have solution of linear equation similar to greenberg or greenshield & android that are explained in logarithm and exponential function, but this article is based mostly on nonlinear solution of macroscopic equation, provided that existing nonlinear relationship between speed and density (homographic the second degree function). As mentioned above, as it gives more reliable and reasonable speeds than greenshield case, it will have more safety. This article has been provided in this field.

Novel energy dissipative method on the adaptive spatial discretization for the Allen–Cahn equation

We propose a novel energy dissipative method for the Allen–Cahn equation on nonuniform grids.For spatial discretization,the classical central difference method is utilized,while the average vector field method is applied for time discretization.Compared with the average vector field method on the uniform mesh,the proposed method can involve fewer grid points and achieve better numerical performance over long time simulation.This is due to the moving mesh method,which can concentrate the grid points more densely where the solution changes drastically.Numerical experiments are provided to illustrate the advantages of the proposed concrete adaptive energy dissipative scheme under large time and space steps over a long time.

Periodic Solutions of a Class of Second-Order Differential Equation

We study the periodic solutions of the second-order differential equations of the form where the functions, , and are periodic of period in the variable t.