The modified characteristic finite difference fractional steps method for the coupled system of fluid dynamics in porous media and its analysis

For the coupled system of multilayer fluid dynamics in porous media, the modified characteristic finite difference fractional steps method applicable to parallel arithmetic is put forward and two‐dimensional and three‐dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, energy method, piecewise biquadratic interpolation, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in L² norm are derived to determine the error in the approximate solution. This method has already been applied to the numerical simulation of multilayer fluid dynamics in porous media. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 665–681, 2003.     

Time-Reversed Dynamical Entropy and Irreversibility in Markovian Random Processes

A concept of time-reversed entropy per unit time is introduced in analogy with the entropy per unit time by Shannon, Kolmogorov, and Sinai. This time-reversed entropy per unit time characterizes the dynamical randomness of a stochastic process backward in time, while the standard entropy per unit time characterizes the dynamical randomness forward in time. The difference between the time-reversed and standard entropies per unit time is shown to give the entropy production of Markovian processes in nonequilibrium steady states.     

Chaotic resonance: A simulation

Stochastic resonance is a statistical phenomenon that has been observed in periodically modulated, noise-driven, bistable systems. The characteristic signatures of the effect include an increase in the signal-to-noise of the output as noise is added to the system, and exponentially decreasing peaks in the probability density as a function of residence times in one state. Presented are the results of a numerical simulation where these same signatures were observed by adding achaotic driving term instead of a white noise term. Although the probability distributions of the noise and chaos inputs were significantly different, the stochastic and chaotic resonances were equal within the experimental error.     

Analysis of nonlinear sliding structures by modified stochastic linearization methods

Probabilistic characteristics of a sliding structure is investigated by using new versions of stochastic linearization technique. The structure is composed of base part and upper part, which are connected to each other in a spring-damping system. Coulomb friction between the base structure and earth ground is considered. Two alternative versions of stochastic linearization approach, suggested by X. Zhang and I. Elishakoff, respectively, are applied to such a sliding structure to evaluate its statistical properties. Compared with the results of Monte Carlo simulation, the two new approaches are performing much better than the conventional one in their applications to the sliding structure. Moreover, numerical results indicate that the criterion proposed by Elishakoff turns out to be superior to all other versions in the problem under study. Numerical results also suggest that the entire structure may be replaced by the rigid body in the sliding problem as long as the difference of velocity responses are considered less important than those of displacement responses.     

The dynamics of a two-component fluid in the presence of capillary forces

In the present paper we study the qualitative behavior ast→∞ of the solution of the Cauchy problem for a system of equations describing a dynamics of a two-component viscous fluid. The model under consideration takes into account the mutual diffusion of the fluid components as well as their capillary interaction. We describe the ω-limit set of trajectories of the dynamical system generated by the problem. It is proved that the stationary solution of the problem, is a homogeneous stationary distribution of one of the components, is asymptotically stable. Any other stationary solution is not asymptotically stable and is even unstable if there are no close stationary solutions corresponding to a smaller energy level. Translated fromMatematicheskie Zametki, Vol. 62, No. 2, pp. 293–305, August, 1997. Translated by A. M. Chebotarev     

AN ENCOMPASSING APPROACH TO MODELING FISHERY DYNAMICS: MODELING DYNAMIC NONLINEAR SYSTEMS

Uncertainty about the nature and significance of nonlinearities and the manner in which dynamics affect future realizations makes model specification the most difficult aspect of modeling dynamic systems. By interpreting several popular fishery models as subcases of a nesting dynamic Taylor series approximation, we isolate the specification differences between these models in a way that accounts for commonalities. On the argument that the differences due to alternative nonlinear forms are likely to be small compared to more mundane considerations such as delay difference and general dynamic lag specification, we propose an alternative model that uses the terms from the first order approximation common to all models combined with a data-based determination of the appropriate lags using the methods of state space time series analysis. Finally, the success of the alternative models is judged in an application to Pacific halibut data.     

Single-server queues with spatially distributed arrivals

Consider a queueing system where customers arrive at a circle according to a homogeneous Poisson process. After choosing their positions on the circle, according to a uniform distribution, they wait for a single server who travels on the circle. The server's movement is modelled by a Brownian motion with drift. Whenever the server encounters a customer, he stops and serves this customer. The service times are independent, but arbitrarily distributed. The model generalizes the continuous cyclic polling system (the diffusion coefficient of the Brownian motion is zero in this case) and can be interpreted as a continuous version of a Markov polling system. Using Tweedie's lemma for positive recurrence of Markov chains with general state space, we show that the system is stable if and only if the traffic intensity is less than one. Moreover, we derive a stochastic decomposition result which leads to equilibrium equations for the stationary configuration of customers on the circle. Steady-state performance characteristics are determined, in particular the expected number of customers in the system as seen by a travelling server and at an arbitrary point in time.     

SOLVABILITY OF FORWARD-BACKWARD SDES AND THE NODAL SET OF HAMILTON-JACOBI-BELLMAN EQUATIONS

The solvability of a class of forward-backward stochastic differential equations (SDEs for short) over an arbitrarily prescribed time duration is studied. The authors design a stochastic relaxed control problem, with both drift and difftusion all being controlled, so that the solvability problem is converted to a problem of finding the nodal set of the viscosity solution to a certain Hamilton-Jacobi-Bellman equation. This method overcomes the fatal difficulty encountered in the traditional contraction mapping approach to the existence theorem of such SDEs.     

Low-energy ion-assisted control of interfacial structures in metallic multilayers

A molecular dynamics method has been used to simulate the argon ion-assisted deposition of Cu/Co/Cu multilayers and to explore ion beam assistance strategies that can be used during or after the growth of each layer to control interfacial structures. A low-argon ion energy of 5–10 eV was found to minimize a combination of interfacial roughness and interlayer mixing (alloying) during the ion-assisted deposition of multilayers. However, complete flattening with simultaneous ion assistance could not be achieved without some mixing between the layers when a constant ion energy approach was used. It was found that multilayers with lower interfacial roughness and intermixing could be grown either by modulating the ion energy during the growth of each metal layer or by utilizing ion assistance only after the completion of each layers deposition. In these latter approaches, relatively high-energy ions could be used since the interface is buried and less susceptible to intermixing. The interlayer mixing dependence upon the thickness of the over layer has been determined as a function of ion energy.