A discontinuous-Galerkin-based immersed boundary method with non-homogeneous boundary conditions and its application to elasticity

We propose a discontinuous-Galerkin-based immersed boundary method for elasticity problems. The resulting numerical scheme does not require boundary fitting meshes and avoids boundary locking by switching the elements intersected by the boundary to a discontinuous Galerkin approximation. Special emphasis is placed on the construction of a method that retains an optimal convergence rate in the presence of non-homogeneous essential and natural boundary conditions. The role of each one of the approximations introduced is illustrated by analyzing an analog problem in one spatial dimension. Finally, extensive two- and three-dimensional numerical experiments on linear and nonlinear elasticity problems verify that the proposed method leads to optimal convergence rates under combinations of essential and natural boundary conditions. (C) 2009 Elsevier B.V. All rights reserved.

CONSTRUCTING QUANTUM OBSERVABLES AND SELF-ADJOINT EXTENSIONS OF SYMMETRIC OPERATORS. III. SELF-ADJOINT BOUNDARY CONDITIONS

This paper completes the review of the theory of self-adjoint extensions of symmetric operators for physicists as a basis for constructing quantum-mechanical observables. It contains a comparative presentation of the well-known methods and a newly proposed method for constructing ordinary self-adjoint differential operators associated with self-adjoint differential expressions in terms of self-adjoint boundary conditions. The new method has the advantage that it does not require explicitly evaluating deficient subspaces and deficiency indices (these latter are determined in passing) and that boundary conditions are of explicit character irrespective of the singularity of a differential expression. General assertions and constructions are illustrated by examples of well-known quantum-mechanical operators like momentum and Hamiltonian.

A Direct Numerov Sixth-order Numerical Scheme to Accurately Solve the Unidimensional Poisson Equation with Dirichlet Boundary Conditions

In this article, we present an analytical direct method, based on a Numerov three-point scheme, which is sixth order accurate and has a linear execution time on the grid dimension, to solve the discrete one-dimensional Poisson equation with Dirichlet boundary conditions. Our results should improve numerical codes used mainly in self-consistent calculations in solid state physics.

Algorithms for two-dimensional cutting stock and strip packing problems using dynamic programming and column generation

We investigate several two-dimensional guillotine cutting stock problems and their variants in which orthogonal rotations are allowed. We first present two dynamic programming based algorithms for the Rectangular Knapsack (RK) problem and its variants in which the patterns must be staged. The first algorithm solves the recurrence formula proposed by Beasley; the second algorithm - for staged patterns - also uses a recurrence formula. We show that if the items are not so small compared to the dimensions of the bin, then these algorithms require polynomial time. Using these algorithms we solved all instances of the RK problem found at the OR-LIBRARY, including one for which no optimal solution was known. We also consider the Two-dimensional Cutting Stock problem. We present a column generation based algorithm for this problem that uses the first algorithm above mentioned to generate the columns. We propose two strategies to tackle the residual instances. We also investigate a variant of this problem where the bins have different sizes. At last, we study the Two-dimensional Strip Packing problem. We also present a column generation based algorithm for this problem that uses the second algorithm above mentioned where staged patterns are imposed. In this case we solve instances for two-, three- and four-staged patterns. We report on some computational experiments with the various algorithms we propose in this paper. The results indicate that these algorithms seem to be suitable for solving real-world instances. We give a detailed description (a pseudo-code) of all the algorithms presented here, so that the reader may easily implement these algorithms. (c) 2007 Elsevier B.V. All rights reserved.

Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.

Taenia saginata: Differential diagnosis of human taeniasis by polymerase chain reaction-restriction fragment length polymorphism assay

Speciation of Taenia in human stool is important because of their different clinical and epidemiological features. DNA analysis has recently become possible which overcomes the problems of differentiating human taeniid cestodes morphologically. In the present study, we evaluated PCR coupled to restriction fragment length polymorphism to differentiate Taenia solium from Taenia saginata eggs present in fecal samples from naturally infected patients. A different Dral-RFLP pattern: a two-band pattern (421 and 100 bp) for T saginata and a three-band pattern (234, 188, and 99 bp) for T solium was observed allowing the two species to be separated.. The lower detection limit of the PCR-RFLP using a non-infected fecal sample prepared with a given number of T saginata eggs was 34 eggs in 2 g stool sediment. The 521 bp mtDNA fragment was detected in 8 out of 12 Taenia sp. carriers (66.6%). of these, three showed a T solium pattern and five a T saginata pattern. (c) 2005 Elsevier B.V. All rights reserved.

The application of interpolating MLS approximations to the analysis of MHD flows

The element-free Galerkin method (EFGM) is a very attractive technique for solutions of partial differential equations, since it makes use of nodal point configurations which do not require a mesh. Therefore, it differs from FEM-like approaches by avoiding the need of meshing, a very demanding task for complicated geometry problems. However, the imposition of boundary conditions is not straightforward, since the EFGM is based on moving-least-squares (MLS) approximations which are not necessarily interpolants. This feature requires, for instance, the introduction of modified functionals with additional unknown parameters such as Lagrange multipliers, a serious drawback which leads to poor conditionings of the matrix equations. In this paper, an interpolatory formulation for MLS approximants is presented: it allows the direct introduction of boundary conditions, reducing the processing time and improving the condition numbers. The formulation is applied to the study of two-dimensional magnetohydrodynamic flow problems, and the computed results confirm the accuracy and correctness of the proposed formulation. (C) 2002 Elsevier B.V. All rights reserved.

Raman scattering, differential scanning calorimetry and Nd3+ spectroscopy in alkali niobium tellurite glasses

Alkali niobium tellurite glasses have been prepared and some of their properties measured by differential scanning calorimetry and Raman scattering. The vitreous domain was established in the pseudo ternary phases diagram for the system TeO2-Nb2O5-(0.5K(2)O-0.5Li(2)O). Raman scattering shows that for samples in the TeO2 rich part of the phase diagram the vitreous structure is composed essentially of (TeO4) units connected by the vertices, as in the alpha-TeO2 crystal. The addition of alkali and niobium oxides causes depolymerization to occur with structures composed essentially of (TeO3) and (NbO6) units. Samples with the composition (mol%) 80TeO(2)-10Nb(2)O(5)-5K(2)O-5Li(2)O, stable against crystallization, were prepared containing up to 10% mol Nd3+. The addition of this oxide increases the rigidity of the vitreous network shifting characteristic temperatures to higher temperatures. For the 10% Nd3+ sample amorphous phase separation is assumed to exist from the observation of two glass transition temperatures. Spectroscopic properties such as Judd-Ofelt Omega(lambda) intensity parameters, radiative emission probabilities, and induced emission cross sections were calculated. From these results and also from the emission quenching observed as a function of Nd3+ concentration, we suggest that these glasses could be utilized in optical amplifying devices. (C) 1999 Elsevier B.V. B.V. All rights reserved.

Non-linear analysis of reinforced concrete plates by the boundary element method and continuum damage mechanics

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

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.

Necessary conditions of optimality for vector-valued impulsive control problems

A vector-valued impulsive control problem is considered whose dynamics, defined by a differential inclusion, are such that the vector fields associated with the singular term do not satisfy the so-called Frobenius condition. A concept of robust solution based on a new reparametrization procedure is adopted in order to derive necessary conditions of optimality. These conditions are obtained by taking a limit of those for an appropriate sequence of auxiliary standard optimal control problems approximating the original one. An example to illustrate the nature of the new optimality conditions is provided. © 2000 Elsevier Science B.V. All rights reserved.

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.

Limit cycles for singular perturbation problems via inverse integrating factor

In this paper singularly perturbed vector fields Xε defined in ℝ2 are discussed. The main results use the solutions of the linear partial differential equation XεV = div(Xε)V to give conditions for the existence of limit cycles converging to a singular orbit with respect to the Hausdorff distance. © SPM.

An exact series solution for the vibration analysis of cylindrical shells with arbitrary boundary conditions

In this paper, an exact series solution for the vibration analysis of circular cylindrical shells with arbitrary boundary conditions is obtained, using the elastic equations based on Flügge's theory. Each of the three displacements is represented by a Fourier series and auxiliary functions and sought in a strong form by letting the solution exactly satisfy both the governing differential equations and the boundary conditions on a point-wise basis. Since the series solution has to be truncated for numerical implementation, the term exactly satisfying should be understood as a satisfaction with arbitrary precision. One of the important advantages of this approach is that it can be universally applied to shells with a variety of different boundary conditions, without the need of making any corresponding modifications to the solution algorithms and implementation procedures as typically required in other techniques. Furthermore, the current method can be easily used to deal with more complicated boundary conditions such as point supports, partial supports, and non-uniform elastic restraints. Numerical examples are presented regarding the modal parameters of shells with various boundary conditions. The capacity and reliability of this solution method are demonstrated through these examples. © 2012 Elsevier Ltd. All rights reserved.

Funções de Green para problemas de valor de contorno com três pontos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)