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).

An exact equation for the free surface of a fluid in a porous medium

We study the problem of the evolution of the free surface of a fluid in a saturated porous medium, bounded from below by a. at impermeable bottom, and described by the Laplace equation with moving-boundary conditions. By making use of a convenient conformal transformation, we show that the solution to this problem is equivalent to the solution of the Laplace equation on a fixed domain, with new variable coefficients, the boundary conditions. We use a kernel of the Laplace equation which allows us to write the Dirichlet-to-Neumann operator, and in this way we are able to find an exact differential-integral equation for the evolution of the free surface in one space dimension. Although not amenable to direct analytical solutions, this equation turns out to allow an easy numerical implementation. We give an explicit illustrative case at the end of the article.

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.

Dirac equation in magnetic-solenoid field

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

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.

Energy and the AdS/CFT correspondence

We consider a scalar field theory on AdS in both minimally and non-minimally coupled cases. We show that there exist constraints which arise in the quantization of the scalar field theory on AdS which cannot be reproduced through the usual AdS/CFT prescription. We argue that the usual energy, defined through the stress-energy tensor, is not the natural one to be considered in the context of the AdS/CFT correspondence. We analyze a new definition of the energy which makes use of the Noether current corresponding to time displacements in global coordinates. We compute the new energy for Dirichlet, Neumann and mixed boundary conditions on the scalar field and for both the minimally and non-minimally coupled cases. Then, we perform the quantization of the scalar field theory on AdS showing that, for 'regular' and 'irregular' modes, the new energy is conserved, positive and finite. We show that the quantization gives rise, in a natural way, to a generalized AdS/CFT prescription which maps to the boundary all the information contained in the bulk. In particular, we show that the divergent local terms of the on-shell action contain information about the Legendre transformed generating functional, and that the new constraints for which the irregular modes propagate in the bulk are the same constraints for which such divergent local terms cancel out. In this situation, the addition of counterterms is not required. We also show that there exist particular cases for which the unitarity bound is reached, and the conformai dimension becomes independent of the effective mass. This phenomenon has no bulk counterpart.

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.

Schrodinger and Dirac operators with the Aharonov-Bohm and magnetic-solenoid fields

We construct all self-adjoint Schrodinger and Dirac operators (Hamiltonians) with both the pure Aharonov-Bohm (AB) field and the so-called magnetic-solenoid field (a collinear superposition of the AB field and a constant magnetic field). We perform a spectral analysis for these operators, which includes finding spectra and spectral decompositions, or inversion formulae. In constructing the Hamiltonians and performing their spectral analysis, we follow, respectively, the von Neumann theory of self-adjoint extensions of symmetric operators and the Krein method of guiding functionals.

Edge states and zero modes in quadratic fermionic models on a 1-D lattice

In una formulazione rigorosa della teoria quantistica, la definizione della varietà Riemanniana spaziale su cui il sistema è vincolato gioca un ruolo fondamentale. La presenza di un bordo sottolinea l'aspetto quantistico del sistema: l'imposizione di condizioni al contorno determina la discretizzazione degli autovalori del Laplaciano, come accade con condizioni note quali quelle periodiche, di Neumann o di Dirichlet. Tuttavia, non sono le uniche possibili. Qualsiasi condizione al bordo che garantisca l'autoaggiunzione dell' operatore Hamiltoniano è ammissibile. Tutte le possibili boundary conditions possono essere catalogate a partire dalla richiesta di conservazione del flusso al bordo della varietà. Alcune possibili condizioni al contorno, permettono l'esistenza di stati legati al bordo, cioè autostati dell' Hamiltoniana con autovalori negativi, detti edge states. Lo scopo di questa tesi è quello di investigare gli effetti di bordo in sistemi unidimensionali implementati su un reticolo discreto, nella prospettiva di capire come simulare proprietà di edge in un reticolo ottico. Il primo caso considerato è un sistema di elettroni liberi. La presenza di edge states è completamente determinata dai parametri di bordo del Laplaciano discreto. Al massimo due edge states emergono, e possono essere legati all' estremità destra o sinistra della catena a seconda delle condizioni al contorno. Anche il modo in cui decadono dal bordo al bulk e completamente determinato dalla scelta delle condizioni. Ammettendo un' interazione quadratica tra siti primi vicini, un secondo tipo di stati emerge in relazione sia alle condizioni al contorno che ai parametri del bulk. Questi stati sono chiamati zero modes, in quanto esiste la possibilità che siano degeneri con lo stato fondamentale. Per implementare le più generali condizioni al contorno, specialmente nel caso interagente, è necessario utilizzare un metodo generale per la diagonalizzazione, che estende la tecnica di Lieb-Shultz-Mattis per Hamiltoniane quadratiche a matrici complesse.

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.

Compilation of measurements from two interdisciplinary STAR-shaped surveys in the North Sea: The ZISCH - Dataset

The measurements were obtained during two North Sea wide STAR-shaped cruises during summer 1986 and winter 1987, which were performed to investigate the circulation induced transport and biologically induced pollutant transfer within the interdisciplinary research in the project "ZISCH - Zirkulation und Schadstoffumsatz in der Nordsee / Circulation and Contaminant Fluxes in the North Sea (1984-1989)". The inventory presents parameters measured on hydrodynamics, nutrient dynamics, ecosystem dynamics and pollutant dynamics in the pelagic and benthic realm. The research program had the objective of quantifying fluxes of major budgets, especially contaminants in the North Sea. In spring 1986, following the phytoplankton spring bloom, and in late winter 1987, at minimum primary production activity, the North Sea ecosystem was investigated on a station net covering the whole North Sea. The station net was shaped like a star. Sampling started in the centre, followed by the northwest section and moving counter clockwise around the North Sea following the residual currents. By this strategy, a time series was measured in the central North Sea and more synoptic data sets were obtained in the individual sections. Generally advection processes have to be considered when comparing the data from different stations. The entire sampling period lasted for more than six weeks in each cruise. Thus, a time-lag should be considered especially when comparing the data from the eastern and the western part of the central and northern North Sea, where samples were taken at the beginning and at the end of the campaign. The ZISCH investigations represented a qualitatively and quantitatively new approach to North Sea research in several respects. (1) The first simultaneous blanket coverage of all important biological, chemical and physical parameters in the entire North Sea ecosystem; (2) the first simultaneous measurements of major contaminants (metals and organohaline compounds) in the different ecosystem compartments; (3) simultaneous determinations of atmospheric inputs of momentum, energy and matter as important ecosystem boundary conditions; (4) performance of the complex measurement program during two seasons, namely the spring plankton bloom and the subsequent winter period of minimal biological activity; and (5) support of data analysis and interpretation by oceanographic and meteorological numerical models on the same scales.

Theoretical Analysis of the No-Slip Boundary Condition Enforcement in SPH Methods

The aim of the present work is to provide an in-depth analysis of the most representative mirroring techniques used in SPH to enforce boundary conditions (BC) along solid profiles. We specifically refer to dummy particles, ghost particles, and Takeda et al. [Prog. Theor. Phys. 92 (1994), 939] boundary integrals. The analysis has been carried out by studying the convergence of the first- and second-order differential operators as the smoothing length (that is, the characteristic length on which relies the SPH interpolation) decreases. These differential operators are of fundamental importance for the computation of the viscous drag and the viscous/diffusive terms in the momentum and energy equations. It has been proved that close to the boundaries some of the mirroring techniques leads to intrinsic inaccuracies in the convergence of the differential operators. A consistent formulation has been derived starting from Takeda et al. boundary integrals (see the above reference). This original formulation allows implementing no-slip boundary conditions consistently in many practical applications as viscous flows and diffusion problems.

On the non-slip boundary condition enforcement in SPH methods.

The implementation of boundary conditions is one of the points where the SPH methodology still has some work to do. The aim of the present work is to provide an in-depth analysis of the most representative mirroring techniques used in SPH to enforce boundary conditions (BC) along solid profiles. We specifically refer to dummy particles, ghost particles, and Takeda et al. [1] boundary integrals. A Pouseuille flow has been used as a example to gradually evaluate the accuracy of the different implementations. Our goal is to test the behavior of the second-order differential operator with the proposed boundary extensions when the smoothing length h and other dicretization parameters as dx/h tend simultaneously to zero. First, using a smoothed continuous approximation of the unidirectional Pouseuille problem, the evolution of the velocity profile has been studied focusing on the values of the velocity and the viscous shear at the boundaries, where the exact solution should be approximated as h decreases. Second, to evaluate the impact of the discretization of the problem, an Eulerian SPH discrete version of the former problem has been implemented and similar results have been monitored. Finally, for the sake of completeness, a 2D Lagrangian SPH implementation of the problem has been also studied to compare the consequences of the particle movement

Combining a Similar Coefficient Identification Algorithm with the Boundary Element Method

The boundary element method (BEM) has been applied successfully to many engineering problems during the last decades. Compared with domain type methods like the finite element method (FEM) or the finite difference method (FDM) the BEM can handle problems where the medium extends to infinity much easier than domain type methods as there is no need to develop special boundary conditions (quiet or absorbing boundaries) or infinite elements at the boundaries introduced to limit the domain studied. The determination of the dynamic stiffness of arbitrarily shaped footings is just one of these fields where the BEM has been the method of choice, especially in the 1980s. With the continuous development of computer technology and the available hardware equipment the size of the problems under study grew and, as the flop count for solving the resulting linear system of equations grows with the third power of the number of equations, there was a need for the development of iterative methods with better performance. In [1] the GMRES algorithm was presented which is now widely used for implementations of the collocation BEM. While the FEM results in sparsely populated coefficient matrices, the BEM leads, in general, to fully or densely populated ones, depending on the number of subregions, posing a serious memory problem even for todays computers. If the geometry of the problem permits the surface of the domain to be meshed with equally shaped elements a lot of the resulting coefficients will be calculated and stored repeatedly. The present paper shows how these unnecessary operations can be avoided reducing the calculation time as well as the storage requirement. To this end a similar coefficient identification algorithm (SCIA), has been developed and implemented in a program written in Fortran 90. The vertical dynamic stiffness of a single pile in layered soil has been chosen to test the performance of the implementation. The results obtained with the 3-d model may be compared with those obtained with an axisymmetric formulation which are considered to be the reference values as the mesh quality is much better. The entire 3D model comprises more than 35000 dofs being a soil region with 21168 dofs the biggest single region. Note that the memory necessary to store all coefficients of this single region is about 6.8 GB, an amount which is usually not available with personal computers. In the problem under study the interface zone between the two adjacent soil regions as well as the surface of the top layer may be meshed with equally sized elements. In this case the application of the SCIA leads to an important reduction in memory requirements. The maximum memory used during the calculation has been reduced to 1.2 GB. The application of the SCIA thus permits problems to be solved on personal computers which otherwise would require much more powerful hardware.

A code for direct numerical simulation of turbulent boundary layers at high Reynolds numbers in BG/P supercomputers

A new high-resolution code for the direct numerical simulation of a zero pressure gradient turbulent boundary layers over a flat plate has been developed. Its purpose is to simulate a wide range of Reynolds numbers from Reθ = 300 to 6800 while showing a linear weak scaling up to 32,768 cores in the BG/P architecture. Special attention has been paid to the generation of proper inflow boundary conditions. The results are in good agreement with existing numerical and experimental data sets.