A marker-and-cell approach to free surface 2-D multiphase flows

This work describes a methodology to simulate free surface incompressible multiphase flows. This novel methodology allows the simulation of multiphase flows with an arbitrary number of phases, each of them having different densities and viscosities. Surface and interfacial tension effects are also included. The numerical technique is based on the GENSMAC front-tracking method. The velocity field is computed using a finite-difference discretization of a modification of the NavierStokes equations. These equations together with the continuity equation are solved for the two-dimensional multiphase flows, with different densities and viscosities in the different phases. The governing equations are solved on a regular Eulerian grid, and a Lagrangian mesh is employed to track free surfaces and interfaces. The method is validated by comparing numerical with analytic results for a number of simple problems; it was also employed to simulate complex problems for which no analytic solutions are available. The method presented in this paper has been shown to be robust and computationally efficient. Copyright (c) 2012 John Wiley & Sons, Ltd.

Virasoro Action on Imaginary Verma Modules and the Operator Form of the KZ-Equation

We define the Virasoro algebra action on imaginary Verma modules for affine and construct an analogue of the Knizhnik-Zamolodchikov equation in the operator form. Both these results are based on a realization of imaginary Verma modules in terms of sums of partial differential operators.

Equation of Motion Governing the Dynamics of Vertically Collapsing Buildings

The present paper aims at contributing to a discussion, opened by several authors, on the proper equation of motion that governs the vertical collapse of buildings. The most striking and tragic example is that of the World Trade Center Twin Towers, in New York City, about 10 years ago. This is a very complex problem and, besides dynamics, the analysis involves several areas of knowledge in mechanics, such as structural engineering, materials sciences, and thermodynamics, among others. Therefore, the goal of this work is far from claiming to deal with the problem in its completeness, leaving aside discussions about the modeling of the resistive load to collapse, for example. However, the following analysis, restricted to the study of motion, shows that the problem in question holds great similarity to the classic falling-chain problem, very much addressed in a number of different versions as the pioneering one, by von Buquoy or the one by Cayley. Following previous works, a simple single-degree-of-freedom model was readdressed and conceptually discussed. The form of Lagrange's equation, which leads to a proper equation of motion for the collapsing building, is a general and extended dissipative form, which is proper for systems with mass varying explicitly with position. The additional dissipative generalized force term, which was present in the extended form of the Lagrange equation, was shown to be derivable from a Rayleigh-like energy function. DOI: 10.1061/(ASCE)EM.1943-7889.0000453. (C) 2012 American Society of Civil Engineers.

Continuity of Dynamical Structures for Nonautonomous Evolution Equations Under Singular Perturbations

In this paper we study the continuity of invariant sets for nonautonomous infinite-dimensional dynamical systems under singular perturbations. We extend the existing results on lower-semicontinuity of attractors of autonomous and nonautonomous dynamical systems. This is accomplished through a detailed analysis of the structure of the invariant sets and its behavior under perturbation. We prove that a bounded hyperbolic global solutions persists under singular perturbations and that their nonlinear unstable manifold behave continuously. To accomplish this, we need to establish results on roughness of exponential dichotomies under these singular perturbations. Our results imply that, if the limiting pullback attractor of a nonautonomous dynamical system is the closure of a countable union of unstable manifolds of global bounded hyperbolic solutions, then it behaves continuously (upper and lower) under singular perturbations.

Three dimensional Lifshitz black hole and the Korteweg-de Vries equation

We consider a solution of three dimensional New Massive Gravity with a negative cosmological constant and use the AdS/CTF correspondence to inquire about the equivalent two dimensional model at the boundary. We conclude that there should be a close relation of the theory with the Korteweg-de Vries equation. (C) 2012 Elsevier B.V..All rights reserved.

Additional scaled solutions to Richards' equation for infiltration and drainage

Warrick and Hussen developed in the nineties of the last century a method to scale Richards' equation (RE) for similar soils. In this paper, new scaled solutions are added to the method of Warrick and Hussen considering a wider range of soils regardless of their dissimilarity. Gardner-Kozeny hydraulic functions are adopted instead of Brooks-Corey functions used originally by Warrick and Hussen. These functions allow to reduce the dependence of the scaled RE on the soil properties. To evaluate the proposed method (PM), the scaled RE was solved numerically using a finite difference method with a fully implicit scheme. Three cases were considered: constant-head infiltration, constant-flux infiltration, and drainage of an initially uniform wet soil. The results for five texturally different soils ranging from sand to clay (adopted from the literature) showed that the scaled solutions were invariant to a satisfactory degree. However, slight deviations were observed mainly for the sandy soil. Moreover, the scaled solutions deviated when the soil profile was initially wet in the infiltration case or when deeply wet in the drainage condition. Based on the PM, a Philip-type model was also developed to approximate RE solutions for the constant-head infiltration. The model showed a good agreement with the scaled RE for the same range of soils and conditions, however only for Gardner-Kozeny soils. Such a procedure reduces numerical calculations and provides additional opportunities for solving the highly nonlinear RE for unsaturated water flow in soils. (C) 2011 Elsevier B.V. All rights reserved.

Continuity of Lyapunov functions and of energy level for generalized gradient semigroup

The global attractor of a gradient-like semigroup has a Morse decomposition. Associated to this Morse decomposition there is a Lyapunov function (differentiable along solutions)-defined on the whole phase space- which proves relevant information on the structure of the attractor. In this paper we prove the continuity of these Lyapunov functions under perturbation. On the other hand, the attractor of a gradient-like semigroup also has an energy level decomposition which is again a Morse decomposition but with a total order between any two components. We claim that, from a dynamical point of view, this is the optimal decomposition of a global attractor; that is, if we start from the finest Morse decomposition, the energy level decomposition is the coarsest Morse decomposition that still produces a Lyapunov function which gives the same information about the structure of the attractor. We also establish sufficient conditions which ensure the stability of this kind of decomposition under perturbation. In particular, if connections between different isolated invariant sets inside the attractor remain under perturbation, we show the continuity of the energy level Morse decomposition. The class of Morse-Smale systems illustrates our results.

Average control of Markov decision processes with Feller transition probabilities and general action spaces

This paper studies the average control problem of discrete-time Markov Decision Processes (MDPs for short) with general state space, Feller transition probabilities, and possibly non-compact control constraint sets A(x). Two hypotheses are considered: either the cost function c is strictly unbounded or the multifunctions A(r)(x) = {a is an element of A(x) : c(x, a) <= r} are upper-semicontinuous and compact-valued for each real r. For these two cases we provide new results for the existence of a solution to the average-cost optimality equality and inequality using the vanishing discount approach. We also study the convergence of the policy iteration approach under these conditions. It should be pointed out that we do not make any assumptions regarding the convergence and the continuity of the limit function generated by the sequence of relative difference of the alpha-discounted value functions and the Poisson equations as often encountered in the literature. (C) 2012 Elsevier Inc. All rights reserved.

Invariant Solutions of Richards' Equation for Water Movement in Dissimilar Soils

Scaling methods allow a single solution to Richards' equation (RE) to suffice for numerous specific cases of water flow in unsaturated soils. During the past half-century, many such methods were developed for similar soils. In this paper, a new method is proposed for scaling RE for a wide range of dissimilar soils. Exponential-power (EP) functions are used to reduce the dependence of the scaled RE on the soil hydraulic properties. To evaluate the proposed method, the scaled RE was solved numerically considering two test cases: infiltration into relatively dry soils having initially uniform water content distributions, and gravity-dominant drainage occurring from initially wet soil profiles. Although the results for four texturally different soils ranging from sand to heavy clay (adopted from the UNSODA database) showed that the scaled solution were invariant for a wide range of flow conditions, slight deviations were observed when the soil profile was initially wet in the infiltration case or deeply wet in the drainage case. The invariance of the scaled RE makes it possible to generalize a single solution of RE to many dissimilar soils and conditions. Such a procedure reduces the numerical calculations and provides additional opportunities for solving the highly nonlinear RE for unsaturated water flow in soils.

Generalized Metropolis dynamics with a generalized master equation: An approach for time-independent and time-dependent Monte Carlo simulations of generalized spin systems

The extension of Boltzmann-Gibbs thermostatistics, proposed by Tsallis, introduces an additional parameter q to the inverse temperature beta. Here, we show that a previously introduced generalized Metropolis dynamics to evolve spin models is not local and does not obey the detailed energy balance. In this dynamics, locality is only retrieved for q = 1, which corresponds to the standard Metropolis algorithm. Nonlocality implies very time-consuming computer calculations, since the energy of the whole system must be reevaluated when a single spin is flipped. To circumvent this costly calculation, we propose a generalized master equation, which gives rise to a local generalized Metropolis dynamics that obeys the detailed energy balance. To compare the different critical values obtained with other generalized dynamics, we perform Monte Carlo simulations in equilibrium for the Ising model. By using short-time nonequilibrium numerical simulations, we also calculate for this model the critical temperature and the static and dynamical critical exponents as functions of q. Even for q not equal 1, we show that suitable time-evolving power laws can be found for each initial condition. Our numerical experiments corroborate the literature results when we use nonlocal dynamics, showing that short-time parameter determination works also in this case. However, the dynamics governed by the new master equation leads to different results for critical temperatures and also the critical exponents affecting universality classes. We further propose a simple algorithm to optimize modeling the time evolution with a power law, considering in a log-log plot two successive refinements.

Exponential stability for a plate equation with p-Laplacian and memory terms

This paper is concerned with the energy decay for a class of plate equations with memory and lower order perturbation of p-Laplacian type, utt+?2u-?pu+?0tg(t-s)?u(s)ds-?ut+f(u)=0inOXR+, with simply supported boundary condition, where O is a bounded domain of RN, g?>?0 is a memory kernel that decays exponentially and f(u) is a nonlinear perturbation. This kind of problem without the memory term models elastoplastic flows.

Conditions for the validity of the quantum Langevin equation

From microscopic models, a Langevin equation can, in general, be derived only as an approximation. Two possible conditions to validate this approximation are studied. One is, for a linear Langevin equation, that the frequency of the Fourier transform should be close to the natural frequency of the system. The other is by the assumption of "slow" variables. We test this method by comparison with an exactly soluble model and point out its limitations. We base our discussion on two approaches. The first is a direct, elementary treatment of Senitzky. The second is via a generalized Langevin equation as an intermediate step.

Estimating total body fat using a skinfold prediction equation in Brazilian children

Background: The double burden of obesity and underweight is increasing in developing countries and simple methods for the assessment of fat mass in children are needed. Aim: To develop and validate a new anthropometric predication equation for assessment of fat mass in children. Subjects and methods: Body composition was assessed in 145 children aged 9.8 +/- 1.3 (SD) years from Sao Paulo, Brazil using dual energy X-ray absorptiometry (DEXA) and skinfold measurements. The study sample was divided into development and validation sub-sets to develop a new prediction equation for FM (PE). Results: Using multiple linear regression analyses, the best equation for predicting FM (R-2 - 0.77) included body weight, triceps skinfold, height, gender and age as independent variables. When cross-validated, the new PE was valid in this sample (R-2 = 0.80), while previously published equations were not. Conclusion: The PE was more valid for Brazilian children that existing equations, but further studies are needed to assess the validity of this PE in other populations.

New Equation for Prediction of Reverse Remodeling after Cardiac Resynchronization Therapy

Objectives: To integrate data from two-dimensional echocardiography (2D ECHO), three-dimensional echocardiography (3D ECHO), and tissue Doppler imaging (TDI) for prediction of left ventricular (LV) reverse remodeling (LVRR) after cardiac resynchronization therapy (CRT). It was also compared the evaluation of cardiac dyssynchrony by TDI and 3D ECHO. Methods: Twenty-four consecutive patients with heart failure, sinus rhythm, QRS = 120 msec, functional class III or IV and LV ejection fraction (LVEF) = 0.35 underwent CRT. 2D ECHO, 3D ECHO with systolic dyssynchrony index (SDI) analysis, and TDI were performed before, 3 and 6 months after CRT. Cardiac dyssynchrony analyses by TDI and SDI were compared with the Pearson's correlation test. Before CRT, a univariate analysis of baseline characteristics was performed for the construction of a logistic regression model to identify the best predictors of LVRR. Results: After 3 months of CRT, there was a moderate correlation between TDI and SDI (r = 0.52). At other time points, there was no strong correlation. Nine of twenty-four (38%) patients presented with LVRR 6 months after CRT. After logistic regression analysis, SDI (SDI > 11%) was the only independent factor in the prediction of LVRR 6 months of CRT (sensitivity = 0.89 and specificity = 0.73). After construction of receiver operator characteristic (ROC) curves, an equation was established to predict LVRR: LVRR =-0.4LVDD (mm) + 0.5LVEF (%) + 1.1SDI (%), with responders presenting values >0 (sensitivity = 0.67 and specificity = 0.87). Conclusions: In this study, there was no strong correlation between TDI and SDI. An equation is proposed for the prediction of LVRR after CRT. Although larger trials are needed to validate these findings, this equation may be useful to candidates for CRT. (Echocardiography 2012;29:678-687)

STRUCTURE AND BIFURCATION OF PULLBACK ATTRACTORS IN A NON-AUTONOMOUS CHAFEE-INFANTE EQUATION

The Chafee-Infante equation is one of the canonical infinite-dimensional dynamical systems for which a complete description of the global attractor is available. In this paper we study the structure of the pullback attractor for a non-autonomous version of this equation, u(t) = u(xx) + lambda(xx) - lambda u beta(t)u(3), and investigate the bifurcations that this attractor undergoes as A is varied. We are able to describe these in some detail, despite the fact that our model is truly non-autonomous; i.e., we do not restrict to 'small perturbations' of the autonomous case.