Monte Carlo thermodynamic and structural properties of the TIP4P water model: dependence on the computational conditions

Classical Monte Carlo simulations were carried out on the NPT ensemble at 25°C and 1 atm, aiming to investigate the ability of the TIP4P water model [Jorgensen, Chandrasekhar, Madura, Impey and Klein; J. Chem. Phys., 79 (1983) 926] to reproduce the newest structural picture of liquid water. The results were compared with recent neutron diffraction data [Soper; Bruni and Ricci; J. Chem. Phys., 106 (1997) 247]. The influence of the computational conditions on the thermodynamic and structural results obtained with this model was also analyzed. The findings were compared with the original ones from Jorgensen et al [above-cited reference plus Mol. Phys., 56 (1985) 1381]. It is notice that the thermodynamic results are dependent on the boundary conditions used, whereas the usual radial distribution functions g(O/O(r)) and g(O/H(r)) do not depend on them.

Dynamic stationary response of reinforced plates by the boundary element method

A direct version of the boundary element method (BEM) is developed to model the stationary dynamic response of reinforced plate structures, such as reinforced panels in buildings, automobiles, and airplanes. The dynamic stationary fundamental solutions of thin plates and plane stress state are used to transform the governing partial differential equations into boundary integral equations (BIEs). Two sets of uncoupled BIEs are formulated, respectively, for the in-plane state ( membrane) and for the out-of-plane state ( bending). These uncoupled systems are joined to formamacro-element, in which membrane and bending effects are present. The association of these macro-elements is able to simulate thin-walled structures, including reinforced plate structures. In the present formulation, the BIE is discretized by continuous and/or discontinuous linear elements. Four displacement integral equations are written for every boundary node. Modal data, that is, natural frequencies and the corresponding mode shapes of reinforced plates, are obtained from information contained in the frequency response functions (FRFs). A specific example is presented to illustrate the versatility of the proposed methodology. Different configurations of the reinforcements are used to simulate simply supported and clamped boundary conditions for the plate structures. The procedure is validated by comparison with results determined by the finite element method (FEM).

Numerical solution of the two-dimensional Gross-Pitaevskii equation for trapped interacting atoms

We present a numerical scheme for solving the time-independent nonlinear Gross-Pitaevskii equation in two dimensions describing the Bose-Einstein condensate of trapped interacting neutral atoms at zero temperature. The trap potential is taken to be of the harmonic-oscillator type and the interaction both attractive and repulsive. The Gross-Pitaevskii equation is numerically integrated consistent with the correct boundary conditions at the origin and in the asymptotic region. Rapid convergence is obtained in all cases studied. In the attractive case there is a limit Co the maximum number of atoms in the condensate. (C) 2000 Published by Elsevier B.V. B.V. All rights reserved.

Short-wave instabilities in the Benjamin-Bona-Mahoney-Peregrine equation: theory and numerics

In this paper we discuss the nonlinear propagation of waves of short wavelength in dispersive systems. We propose a family of equations that is likely to describe the asymptotic behaviour of a large class of systems. We then restrict our attention to the analysis of the simplest nonlinear short-wave dynamics given by U-0 xi tau, = U-0 - 3(U-0)(2). We integrate numerically this equation for periodic and non-periodic boundary conditions, and we find that short waves may exist only if the amplitude of the initial profile is not too large.

Geometrical properties of coupled oscillators at synchronization

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Dressing approach to the nonvanishing boundary value problem for the AKNS hierarchy

We propose an approach to the nonvanishing boundary value problem for integrable hierarchies based on the dressing method. Then we apply the method to the AKNS hierarchy. The solutions are found by introducing appropriate vertex operators that takes into account the boundary conditions.

Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.

Necessary conditions for optimal impulsive control problems

A Maximum Principle is derived for a class of optimal control problems arising in midcourse guidance, in which certain controls are represented by measures and, the state trajectories are functions of bounded variation. The optimality conditions improves on previous optimality conditions by allowing nonsmooth data, measurable time dependence, and a possibly time varying constraint set for the conventional controls.

Quasilinear dirichlet problems in RN with critical growth

In this study, a given quasilinear problem is solved using variational methods. In particular, the existence of nontrivial solutions for GP is examined using minimax methods. The main theorem on the existence of a nontrivial solution for GP is detailed.

Numerical study of the coupled time-dependent Gross-Pitaevskii equation: Application to Bose-Einstein condensation

The Bose-Einstein condensate of several types of trapped bosons at ultralow temperature was described using the coupled time dependent Gross-Pitaevskii equation. Both the stationary and time evolution problems were analyzed using this approach. The ground state stationary wave functions were found to be sharply peaked near the origin for attractive interatomic interaction for larger nonlinearity while for a repulsive interatomic interaction the wave function extends over a larger region of space.

Local dimension and finite time prediction in spatiotemporal chaotic systems

Predictability is related to the uncertainty in the outcome of future events during the evolution of the state of a system. The cluster weighted modeling (CWM) is interpreted as a tool to detect such an uncertainty and used it in spatially distributed systems. As such, the simple prediction algorithm in conjunction with the CWM forms a powerful set of methods to relate predictability and dimension.

On the role of boundary terms in the AdS/CFT correspondence

We consider a scalar field theory on AdS, and show that the usual AdS/CFT prescription is unable to map to the boundary a part of the information arising from the quantization in the bulk. We propose a solution to this problem by defining the energy of the theory in the bulk through the Noether current corresponding to time displacements, and, in addition, by introducing a proper generalized AdS/CFT prescription. We also show how this extended formulation could be used to consistently describe double-trace interactions in the boundary. The formalism is illustrated by focusing on the non-minimally coupled case using Dirichlet boundary conditions. © 2004 Published by Elsevier B.V.

One-dimensional superfluid Bose-Fermi mixture: Mixing, demixing, and bright solitons

We study an ultracold and dilute superfluid Bose-Fermi mixture confined in a strictly one-dimensional (1D) atomic waveguide by using a set of coupled nonlinear mean-field equations obtained from the Lieb-Liniger energy density for bosons and the Gaudin-Yang energy density for fermions. We consider a finite Bose-Fermi interatomic strength gbf and both periodic and open boundary conditions. We find that with periodic boundary conditions-i.e., in a quasi-1D ring-a uniform Bose-Fermi mixture is stable only with a large fermionic density. We predict that at small fermionic densities the ground state of the system displays demixing if gbf >0 and may become a localized Bose-Fermi bright soliton for gbf <0. Finally, we show, using variational and numerical solutions of the mean-field equations, that with open boundary conditions-i.e., in a quasi-1D cylinder-the Bose-Fermi bright soliton is the unique ground state of the system with a finite number of particles, which could exhibit a partial mixing-demixing transition. In this case the bright solitons are demonstrated to be dynamically stable. The experimental realization of these Bose-Fermi bright solitons seems possible with present setups. © 2007 The American Physical Society.

Interações modais não ressonantes em vigas cantilever flexíveis

Pós-graduação em Engenharia Mecânica - FEG

Boundary condition identification of a bridge pier using experimental data and FEM modal updating

Over the past twenty years, new technologies have required an increasing use of mathematical models in order to understand better the structural behavior: finite element method is the one mostly used. However, the reliability of this method applied to different situations has to be tried each time. Since it is not possible to completely model the reality, different hypothesis must be done: these are the main problems of FE modeling. The following work deals with this problem and tries to figure out a way to identify some of the unknown main parameters of a structure. This main research focuses on a particular path of study and development, but the same concepts can be applied to other objects of research. The main purpose of this work is the identification of unknown boundary conditions of a bridge pier using the data acquired experimentally with field tests and a FEM modal updating process. This work doesn’t want to be new, neither innovative. A lot of work has been done during the past years on this main problem and many solutions have been shown and published. This thesis just want to rework some of the main aspects of the structural optimization process, using a real structure as fitting model.