A REACTIVE DYNAMIC CONTINUUM USER EQUILIBRIUM MODEL FOR BI-DIRECTIONAL PEDESTRIAN FLOWS

In this paper, a reactive dynamic user equilibrium model is extended to simulate two groups of pedestrians traveling on crossing paths in a continuous walking facility. Each group makes path choices to minimize the travel cost to its destination in a reactive manner based on instantaneous information. The model consists of a conservation law equation coupled with an Eikonal-type equation for each group. The velocity-density relationship of pedestrian movement is obtained via an experimental method. The model is solved using a finite volume method for the conservation law equation and a fast-marching method for the Eikonal-type equation on unstructured grids. The numerical results verify the rationality of the model and the validity of the numerical method. Based on this continuum model, a number of results, e.g., the formation of strips or moving clusters composed of pedestrians walking to the same destination, are also observed.     

A GHOST FLUID BASED FRONT TRACKING METHOD FOR MULTIMEDIUM COMPRESSIBLE FLOWS

Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.     

A DIAGONAL SPLIT-CELL MODEL FOR THE OVERLAPPING YEE FDTD METHOD

In this paper, we present a nonorthogonal overlapping Yee method for solv- ing Maxwell＇s equations using the diagonal split-cell model. When material interface is presented, the diagonal split-cell model does not require permittivity averaging so that better accuracy can be achieved. Our numerical results on optical force computation show that the standard FDTD method converges linearly, while the proposed method achieves quadratic convergence and better accuracy.     

POSITIVE SOLUTIONS FOR WEAKLY COUPLED NONLINEAR ELLIPTIC SYSTEMS

In this article,we consider the existence of positive solutions for weakly cou-pled nonlinear elliptic systems {-△u+u (1+a(x))|u| p-1 u+μ|u| α-2 u|v|β+λv in R~N,-△v+v=(1+bx))|v|p-1v+μ|u|α|v|β-2v+λu in R N.(0.1) To find nontrivial solutions,we first investigate autonomous systems.In this case,results of bifurcation from semi-trivial solutions are obtained by the implicit function theorem.Next,the existence of positive solutions of problem(0.1) is obtained by variational methods.     

On the Spectrum of a Class of Strongly Coupled p-Laplacian Systems

Consider the nonlinear coupled elliptic system     

On sums of sparse prime subsets

For arbitrary c00, if A is a subset of the primes less than x with cardinality δx(logx)-1, δ≥(logx)-c0, then there exists a positive constant c such that the cardinality of A+A is larger than cδx(loglogx)-1.     

HOLDER ESTIMATES FOR A CLASS OF DEGENERATE ELLIPTIC EQUATIONS

In this paper we study the regularity theory of the solutions of a class of degenerate elliptic equations in divergence form. By introducing a proper distance and applying the compactness method we establish the HSlder type estimates for the weak solutions.     

A note on the existence of fractional f-factors in random graphs

Let G=Gn,p be a binomial random graph with n vertices and edge probability p=p(n),and f be a nonnegative integer-valued function defined on V(G) such that 0a≤f(x)≤bnp-2np ㏒n for every x ∈V(G). An fractional f-indicator function is an function h that assigns to each edge of a graph G a number h(e) in [0,1] so that for each vertex x,we have dh G(x)=f(x),where dh G(x) = x∈e h(e) is the fractional degree of x in G. Set Eh = {e:e ∈E(G) and h(e)=0}.If Gh is a spanning subgraph of G such that E(Gh)=Eh,then Gh is called an fractional f-factor of G. In this paper,we prove that for any binomial random graph Gn,p with p≥n-23,almost surely Gn,p contains an fractional f-factor.     

A CLASS OF MULTIWAVELETS AND PROJECTED FRAMES FROM TWO-DIRECTION WAVELETS

This article aims at studying two-direction refinable functions and two-direction wavelets in the setting Rs, s 1. We give a sufficient condition for a two-direction refinable function belonging to L2(Rs). Then, two theorems are given for constructing biorthogonal(orthogonal) two-direction refinable functions in L2(Rs) and their biorthogonal(orthogonal) two-direction wavelets, respectively. From the constructed biorthogonal(orthogonal)two-direction wavelets, symmetric biorthogonal(orthogonal) multiwaveles in L2(Rs) can be obtained easily. Applying the projection method to biorthogonal(orthogonal) two-direction wavelets in L2(Rs), we can get dual(tight) two-direction wavelet frames in L2(Rm), where m ≤ s. From the projected dual(tight) two-direction wavelet frames in L2(Rm), symmetric dual(tight) frames in L2(Rm) can be obtained easily. In the end, an example is given to illustrate theoretical results.     

SPECTRAL/HP ELEMENT METHOD WITH HIERARCHICAL RECONSTRUCTION FOR SOLVING NONLINEAR HYPERBOLIC CONSERVATION LAWS

The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectral volume schemes for solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element method for solving hyperbolic conservation laws. An orthogonal spectral basis written in terms of Jacobi polynomials is applied. High computational efficiency is obtained due to such matrix-free algorithm. The formulation is conservative, and essential nomoscillation is enforced by the HR limiter. We show that HR preserves the order of accuracy of the spectral/hp element method for smooth solution problems and generate essentially non-oscillatory solutions profiles for capturing discontinuous solutions without local characteristic decomposition. In addition, we introduce a postprocessing technique to improve HR for limiting high degree numerical solutions.     

A NEW REDUCED-ORDER FVE ALGORITHM BASED ON POD METHOD FOR VISCOELASTIC EQUATIONS

A proper orthogonal decomposition （POD） technique is used to reduce the finite volume element （FVE） method for two-dimensional （2D） viscoelastic equations. A reduced-order fully discrete FVE algorithm with fewer degrees of freedom and sufficiently high accuracy based on POD method is established. The error estimates of the reduced- order fully discrete FVE solutions and the implementation for solving the reduced-order fully discrete FVE algorithm are provided. Some numerical examples are used to illus- trate that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced-order fully discrete FVE algorithm is one of the most effective numerical methods by comparing with corresponding numerical results of finite element formulation and finite difference scheme and that the reduced-order fully discrete FVE algorithm based on POD method is feasible and efficient for solving 2D viscoelastic equations.     

The Factorization Method to Solve a Class of Inverse Potential Scattering Problems for Schrdinger Equations

This paper is concerned with the inverse scattering problems for Schrdinger equations with compactly supported potentials.For purpose of reconstructing the support of the potential,we derive a factorization of the scattering amplitude operator A and prove that the ranges of （A＊ A） ^1/4 and G which maps more general incident fields than plane waves into the scattering amplitude coincide.As an application we characterize the support of the potential using only the spectral data of the operator A.     

FIXED POINTS AND STABILITY FOR QUARTIC MAPPINGS IN β-NORMED LEFT BANACH MODULES ON BANACH ALGEBRAS

The goal of the present paper is to investigate some new HUR-stability results by applying the alternative fixed point of generalized quartic functional equationin β-Banach modules on Banach algebras. The concept of Ulam-Hyers-Rassias stability （briefly, HUR-stability） originated from Th. M. Rassias stability theorem that appeared in his paper： On the stability of the linear mapping in Banach spaces, Proc. Amer. Math. Soc. 72 （1978）, 297-300.     

A REDUCED-ORDER MFE FORMULATION BASED ON POD METHOD FOR PARABOLIC EQUATIONS

In this paper, we extend the applications of proper orthogonal decomposition (POD) method, i.e., apply POD method to a mixed finite element (MFE) formulation naturally satisfied Brezz-Babu?ska for para...     

CHEBYSHEV METHODS WITH DISCRETE NOISE： THE τ-ROCK METHODS

Stabilized or Chebyshev explicit methods have been widely used in the past to solve stiff ordinary differential equations. Making use of special properties of Chebyshev-like polynomials, these methods have favorable stability properties compared to standard explicit methods while remaining explicit. A new class of such methods, called ROCK, introduced in [Numer. Math., 90, 1-18, 2001] has recently been extended to stiff stochastic differential equations under the name S-ROCK [C. R. Acad. Sci. Paris, 345（10）, 2007 and Commun. Math. Sci, 6（4）, 2008]. In this paper we discuss the extension of the S-ROCK methods to systems with discrete noise and propose a new class of methods for such problems, the τ-ROCK methods. One motivation for such methods is the simulation of multi-scale or stiff chemical kinetic systems and such systems are the focus of this paper, but our new methods could potentially be interesting for other stiff systems with discrete noise. Two versions of the τ-ROCK methods are discussed and their stability behavior is analyzed on a test problem. Compared to the τ-leaping method, a significant speed-up can be achieved for some stiff kinetic systems. The behavior of the proposed methods are tested on several numerical experiments.     

ON THE SINGULARITY OF LEAST SQUARES ESTIMATOR FOR MEAN-REVERTING α-STABLE MOTIONS

We study the problem of parameter estimation for mean-reverting α-stable motion, dXt = （a0 - θ0Xt）dt ＋ dZt, observed at discrete time instants. A least squares estimator is obtained and its asymptotics is discussed in the singular case （a0, θ0） = （0, 0）. If a0 = 0, then the mean-reverting α-stable motion becomes Ornstein-Uhlenbeck process and is studied in [7] in the ergodic case θ0 〉 0. For the Ornstein-Uhlenbeck process, asymptotics of the least squares estimators for the singular case （θ0 = 0） and for ergodic case （θ0 〉 0） are completely different.     

NUMERICAL SIMULATION OF THE XZ TAURI SUPERSONIC ASTROPHYSICAL JET

Computational gas dynamical simulations using the WENO-LF method are applied to modeling the high Mach number astrophysical jet XZ Tauri, including the effects of radiative cooling. Mach 55 simulations of the pulsed proto-jet are presented and analyzed in terms of interacting nonlinear waves： terminal Mach disks, bow shocks, and Meshkov-Richtmyer instabilities of the leading jet contact boundary.     

THE ACCELERATED SEARCH-EXTENSION METHOD FOR COMPUTING MULTIPLE SOLUTIONS OF SEMILINEAR PDEs

In this paper, we propose an accelerated search-extension method （ASEM） based on the interpolated coefficient finite element method, the search-extension method （SEM） and the two-grid method to obtain the multiple solutions for semilinear elliptic equations. This strategy is not only successfully implemented to obtain multiple solutions for a class of semilinear elliptic boundary value problems, but also reduces the expensive computation greatly. The numerical results in I-D and 2-D cases will show the efficiency of our approach.     

Nonradial Solutions of a Mixed Concave-Convex Elliptic Problem

We study the homogeneous Dirichlet problem in a ball for semi-linear elliptic problems derived from the Brezis-Nirenberg one with concave-convex nonlinearities. We are interested in determining non-radial solutions which are invariant with respect to some subgroup of the orthogonal group. We prove that unlike separated nonlinearities, there are two types of solutions, one converging to zero and one diverging. We conclude at the end on the classification of non radial solutions related to the nonlinearity used.     

The Tanh Method for Kink Solution of Some Modified Nonlinear Equation

The tanh method is a very powerful technique for computation of exact traveling wave, in this paper this method has been employed for special modified states of Burger, Klein-Gordon and Fisher-Burger equations and the solitary solution of these equations are derived.