Theory and Application of Characteristic Finite Element Domain Decomposition Procedures for Coupled System of Dynamics of Fluids in Porous Media

For a coupled system of multiplayer dynamics of fluids in porous media,the characteristic finiteelement domain decomposition procedures applicable to parallel arithmetic are put forward.Techniques suchas calculus of variations,domain decomposition,characteristic method,negative norm estimate,energy methodand the theory of prior estimates are adopted.Optimal order estimates in L~2 norm are derived for the error inthe approximate solution.     

质子中子星物质的热力学性质

在相对论σ-ω-ρ模型的平均场近似下, 研究了质子中子星物质在均熵状态下的组成、温度和物态方程. 如给定每一个重子的熵, 一些热力学量的值将随重子密度的增加而增加, 当考虑超子时, 这些值会减小. 给定重子密度, 中子在S=2时的组分比S=1时的小, 而质子、电子、μ子在S=2时的组分比S=1时的大, 特别是在低密度区域. S是每个重子的熵. 保持重子密度不变, 在低密度区域, 超子在S=2时的组分比S=1时的大, 在高密度区域则相反. 同样, 在同一重子密度处, S=2时的温度、能量密度及压强分别比S=1时的大. 另外, 有限熵对粒子组分和温度的影响比对质子中子星物质的物态方程的影响大. 还研究了反粒子的贡献, 他们确实很小.     

Electroencephalogram (EEG) time series classification: Applications in epilepsy

Epilepsy is among the most common brain disorders. Approximately 25–30% of epilepsy patients remain unresponsive to anti-epileptic drug treatment, which is the standard therapy for epilepsy. In this study, we apply optimization-based data mining techniques to classify the brain's normal and epilepsy activity using intracranial electroencephalogram (EEG), which is a tool for evaluating the physiological state of the brain. A statistical cross validation and support vector machines were implemented to classify the brain's normal and abnormal activities. The results of this study indicate that it may be possible to design and develop efficient seizure warning algorithms for diagnostic and therapeutic purposes. Research was partially supported by the Rutgers Research Council grant-202018, the NSF grants DBI-980821, CCF-0546574, IIS-0611998, and NIH grant R01-NS-39687-01A1.     

USING A SYSTEM DYNAMICS APPROACH TO MODEL SUSTAINABILITY INDICATORS FOR IRRIGATION SYSTEMS IN AUSTRALIA

ABSTRACT. An indicator defined as a function of the total water diversion through the Coleambally canal and the potential irrigation demand is selected to represent the sustainability of the irrigation water system in the Coleambally irrigation area, Australia. A simulation procedure using a system dynamics approach was developed to evaluate the indicator. The procedure includes water diversion assessment, potential crop water demand and total gross margin. Three cases of water supply options (surface water, ground water pumping and water trading), two cases of changes in the total agricultural area and three cropping pattern scenarios were simulated to better understand their impact on sustainability. The simulated results indicate that increasing the agricultural area reduces the sustainability of the irrigation system because the demand of water increased despite increase in the gross margin. The scenarios show that imposed water trading and ground water pumping would considerably increase the supply system having a positive impact on the sustainability. The paper concludes that a multi‐objective sustainability indicator taking account of economic and environmental issues could be more useful.     

Multiscale Finite Element Modeling of Arc Dynamics in a DC Plasma Torch

The dynamics of the electric arc inside a direct current non-transferred arc plasma torch are simulated using a three-dimensional, transient, equilibrium model. The fluid and electromagnetic equations are solved numerically in a fully coupled approach by a multiscale finite element method. Simulations of a torch operating with argon and argon–hydrogen under different operating conditions are presented. The model is able to predict the operation of the torch in steady and takeover modes without any further assumption on the reattachment process except for the use of an artificially high electrical conductivity near the electrodes, needed because of the equilibrium assumption. The results obtained indicate that the reattachment process in these operating modes may be driven by the movement of the arc rather than by a breakdown-like process. It is also found that, for a torch operating in these modes and using straight gas injection, the arc will tend to re-attach to the opposite side of its original attachment. This phenomenon seems to be produced by a net angular momentum on the arc due to the imbalance between magnetic and fluid drag forces.     

Metastability and exponential approach to equilibrium for low-temperature stochastic Ising models

We present a proof of the exponential convergence to equilibrium of single-spin-flip stochastic dynamics for the two-dimensional Ising ferromagnet in the low-temperature case with not too small external magnetic fieldh uniformly in the volume and in the boundary conditions.     

Stochastic modeling of a billiard in a gravitational field: Power law behavior of Lyapunov exponents

We consider the motion of a point particle (billiard) in a uniform gravitational field constrained to move in a symmetric wedge-shaped region. The billiard is reflected at the wedge boundary. The phase space of the system naturally divides itself into two regions in which the tangent maps are respectively parabolic and hyperbolic. It is known that the system is integrable for two values of the wedge half-angle ₁ and ₂ and chaotic for ₁<< ₂. We study the system at three levels of approximation: first, where the deterministic dynamics is replaced by a random evolution; second, where, in addition, the tangent map in each region is, replaced by its average; and third, where the tangent map is replaced by a single global average. We show that at all three levels the Lyapunov exponent exhibits power law behavior near ₁ and ₂ with exponents 1/2 and 1, respectively. We indicate the origin of the exponent 1, which has not been observed in unaccelerated billiards.     

Computer simulation of some dynamical properties of the Lorentz gas

We carried out molecular dynamics simulations of a Lorentz gas, consisting of a lone hydrogen molecule moving in a sea of stationary argon atoms. A Lennard-Jones form was assumed for the H₂-Ar potential. The calculations were performed at a reduced temperatureK ^* =kT/_{H 2–Ar} = 4.64 and at reduced densities ^*= _{Ar Ar} ³ in the range 0.074–0.414. The placement of Ar atoms was assumed to be random rather than dictated by equilibrium considerations. We followed the trajectories of many H₂ molecules, each of which is assigned in turn a velocity given by the Maxwell-Boltzmann distribution at the temperature of the simulation. Solving the equations of motion classically, we obtained the translational part of the incoherent dynamic structure factor for the H₂ molecule,S _tr(q, ). This was convoluted with the rotational structure factorS _rot(q, ) calculated assuming unhindered rotation to obtain the total structure factorS(q, ). Our results agree well with experimental data on this function obtained by Egelstaffet al. At the highest density ( ^*=0.414) we studied the dependence ofS(q, ) on system size (number of Ar atoms), number of H₂ molecules for which trajectories are generated, and the length of time over which these trajectories are followed.     

Dynamics of the current-driven Josephson junction

The irreversible macroscopic dynamics of the Josephson junction coupled to external wires acting as a current source is derived rigorously from the underlying microscopic Hamiltonian quantum mechanics. The external systems are treated in the singular coupling limit. The use of this limit is explicitly justified via an interpretation of the singular coupling limit in terms of the relative magnitudes of system, reservoir, and coupling energies. The qualitative behavior of the macroscopic dynamical equations is shown to depend sensitively and crucially on the interaction between the wires and the superconductors and on the size of the wires: the dc Josephson effect only happens when one lets Cooper pairs be driven into the junction by collective (i.e., small) reservoirs.     

Ar—Kr溶液扩散系数的分子动力学模拟及其与温度的关系

用分子动力学模拟方法研究确定Ａｒ－Ｋｒ溶液的自扩散系数Ｄ１、Ｄ２和互扩散系数Ｄ１２以及它们随温度变化的规律。结果表明，分别用Ｇｒｅｅｎ－Ｋｕｂｏ法和Ｅｉｎｓｔｅｉｎ法得到的扩散系数在数值上一致；该溶液的３种扩散系数均满足Ｄ＝Ｄ０ｅ＾Ｅ／ＲＴ关系。