A Mass-Conserving Algorithm for Dynamical Lubrication Problems With Cavitation

A numerical algorithm for fully dynamical lubrication problems based on the Elrod-Adams formulation of the Reynolds equation with mass-conserving boundary conditions is described. A simple but effective relaxation scheme is used to update the solution maintaining the complementarity conditions on the variables that represent the pressure and fluid fraction. The equations of motion are discretized in time using Newmark`s scheme, and the dynamical variables are updated within the same relaxation process just mentioned. The good behavior of the proposed algorithm is illustrated in two examples: an oscillatory squeeze flow (for which the exact solution is available) and a dynamically loaded journal bearing. This article is accompanied by the ready-to-compile source code with the implementation of the proposed algorithm. [DOI: 10.1115/1.3142903]

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.

A Floating-point Extended Kalman Filter Implementation for Autonomous Mobile Robots

Localization and Mapping are two of the most important capabilities for autonomous mobile robots and have been receiving considerable attention from the scientific computing community over the last 10 years. One of the most efficient methods to address these problems is based on the use of the Extended Kalman Filter (EKF). The EKF simultaneously estimates a model of the environment (map) and the position of the robot based on odometric and exteroceptive sensor information. As this algorithm demands a considerable amount of computation, it is usually executed on high end PCs coupled to the robot. In this work we present an FPGA-based architecture for the EKF algorithm that is capable of processing two-dimensional maps containing up to 1.8 k features at real time (14 Hz), a three-fold improvement over a Pentium M 1.6 GHz, and a 13-fold improvement over an ARM920T 200 MHz. The proposed architecture also consumes only 1.3% of the Pentium and 12.3% of the ARM energy per feature.

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.

Boundary and internal conditions for adjoint fluid-flow problems

The ever-increasing robustness and reliability of flow-simulation methods have consolidated CFD as a major tool in virtually all branches of fluid mechanics. Traditionally, those methods have played a crucial role in the analysis of flow physics. In more recent years, though, the subject has broadened considerably, with the development of optimization and inverse design applications. Since then, the search for efficient ways to evaluate flow-sensitivity gradients has received the attention of numerous researchers. In this scenario, the adjoint method has emerged as, quite possibly, the most powerful tool for the job, which heightens the need for a clear understanding of its conceptual basis. Yet, some of its underlying aspects are still subject to debate in the literature, despite all the research that has been carried out on the method. Such is the case with the adjoint boundary and internal conditions, in particular. The present work aims to shed more light on that topic, with emphasis on the need for an internal shock condition. By following the path of previous authors, the quasi-1D Euler problem is used as a vehicle to explore those concepts. The results clearly indicate that the behavior of the adjoint solution through a shock wave ultimately depends upon the nature of the objective functional.

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

Reaction-diffusion systems coupled at the boundary and the Morse-Smale property

We study an one-dimensional nonlinear reaction-diffusion system coupled on the boundary. Such system comes from modeling problems of temperature distribution on two bars of same length, jointed together, with different diffusion coefficients. We prove the transversality property of unstable and stable manifolds assuming all equilibrium points are hyperbolic. To this end, we write the system as an equation with noncontinuous diffusion coefficient. We then study the nonincreasing property of the number of zeros of a linearized nonautonomous equation as well as the Sturm-Liouville properties of the solutions of a linear elliptic problem. (C) 2008 Elsevier Inc. All rights reserved.

On the Multiplicity of Orthogonal Geodesics in Riemannian Manifold With Concave Boundary. Applications to Brake Orbits and Homoclinics

Let (M, g) be a complete Riemannian Manifold, Omega subset of M an open subset whose closure is diffeomorphic to an annulus. If partial derivative Omega is smooth and it satisfies a strong concavity assumption, then it is possible to prove that there are at least two geometrically distinct geodesics in (Omega) over bar = Omega boolean OR partial derivative Omega starting orthogonally to one connected component of partial derivative Omega and arriving orthogonally onto the other one. The results given in [6] allow to obtain a proof of the existence of two distinct homoclinic orbits for an autonomous Lagrangian system emanating from a nondegenerate maximum point of the potential energy, and a proof of the existence of two distinct brake orbits for a. class of Hamiltonian systems. Under a further symmetry assumption, it is possible to show the existence of at least dim(M) pairs of geometrically distinct geodesics as above, brake orbits and homoclinics.

Essentialism i religionsundervisningen : Ett religionsdidaktiskt problem

Essentialist concepts of religion are common in the teaching of religion in schools and to a certain extent also in the academic discipline of religious studies. In this article, a number of problems with essentialist perceptions of religion are discussed. In the first part of the article a thesis is maintained, according to which essentialist conceptions of religion or specific religions are too limited to be of value in the teaching of religion. This is done through examples of essentialist expressions about religion. The examples are grouped according to a typology of different kinds of essentialism. Two main categories, each with two sub-categories are identified. Thus, the category of essentialism regarding the substance of religion is divided into transcendental or theological essentialism (which presupposes the existence of a sacred power of some kind, the experience of which is the basis for religion), and core essentialism (where it is presupposed that certain ideas or concepts constitute religion as a general category or specific religions). Likewise, the category of essentialism regarding the function of religion has two sub-categories: positive and negative essentialism. These kinds of essentialism presuppose that religion or specific religions are inherently good or harmful respectively to human beings. Examples from each of these categories are given and discussed. In the second part of the article, Benson Saler’s open concept of religion is presented as an alternative to essentialist or bounded perceptions. It is based on Ludwig Wittgenstein’s idea of family resemblances and on prototype theory. In connection with this, it is argued that a certain kind of conscious ethnocentrism is needed as a point of departure in the study and teaching of religion. The metaphor of education as a journey from the familiar out into the unfamiliar and back again is suggested as a possible pattern for such teaching. Finally,some examples of non-essentialist ways to introduce religions are offered.

Gravitational self-localization for spherical masses

In this work, I consider the center-of-mass wave function for a homogenous sphere under the influence of the self-interaction due to Newtonian gravity. I solve for the ground state numerically and calculate the average radius as a measure of its size. For small masses, M≲10−17 kg, the radial size is independent of density, and the ground state extends beyond the extent of the sphere. For masses larger than this, the ground state is contained within the sphere and to a good approximation given by the solution for an effective radial harmonic-oscillator potential. This work thus determines the limits of applicability of the point-mass Newton Schrödinger equations for spherical masses. In addition, I calculate the fringe visibility for matter-wave interferometry and find that in the low-mass case, interferometry can in principle be performed, whereas for the latter case, it becomes impossible. Based on this, I discuss this transition as a possible boundary for the quantum-classical crossover, independent of the usually evoked environmental decoherence. The two regimes meet at sphere sizes R≈10−7 m, and the density of the material causes only minor variations in this value.

The Nurse Professional Competence (NPC) scale : self-reported competence among nursing students on the point of graduation

BACKGROUND: International organisations, e.g. WHO, stress the importance of competent registered nurses (RN) for the safety and quality of healthcare systems. Low competence among RNs has been shown to increase the morbidity and mortality of inpatients. OBJECTIVES: To investigate self-reported competence among nursing students on the point of graduation (NSPGs), using the Nurse Professional Competence (NPC) Scale, and to relate the findings to background factors. METHODS AND PARTICIPANTS: The NPC Scale consists of 88 items within eight competence areas (CAs) and two overarching themes. Questions about socio-economic background and perceived overall quality of the degree programme were added. In total, 1086 NSPGs (mean age, 28.1 [20-56]years, 87.3% women) from 11 universities/university colleges participated. RESULTS: NSPGs reported significantly higher scores for Theme I "Patient-Related Nursing" than for Theme II "Organisation and Development of Nursing Care". Younger NSPGs (20-27years) reported significantly higher scores for the CAs "Medical and Technical Care" and "Documentation and Information Technology". Female NSPGs scored significantly higher for "Value-Based Nursing". Those who had taken the nursing care programme at upper secondary school before the Bachelor of Science in Nursing (BSN) programme scored significantly higher on "Nursing Care", "Medical and Technical Care", "Teaching/Learning and Support", "Legislation in Nursing and Safety Planning" and on Theme I. Working extra paid hours in healthcare alongside the BSN programme contributed to significantly higher self-reported scores for four CAs and both themes. Clinical courses within the BSN programme contributed to perceived competence to a significantly higher degree than theoretical courses (93.2% vs 87.5% of NSPGs). SUMMARY AND CONCLUSION: Mean scores reported by NSPGs were highest for the four CAs connected with patient-related nursing and lowest for CAs relating to organisation and development of nursing care. We conclude that the NPC Scale can be used to identify and measure aspects of self-reported competence among NSPGs.

Graphical models and point set matching

Point pattern matching in Euclidean Spaces is one of the fundamental problems in Pattern Recognition, having applications ranging from Computer Vision to Computational Chemistry. Whenever two complex patterns are encoded by two sets of points identifying their key features, their comparison can be seen as a point pattern matching problem. This work proposes a single approach to both exact and inexact point set matching in Euclidean Spaces of arbitrary dimension. In the case of exact matching, it is assured to find an optimal solution. For inexact matching (when noise is involved), experimental results confirm the validity of the approach. We start by regarding point pattern matching as a weighted graph matching problem. We then formulate the weighted graph matching problem as one of Bayesian inference in a probabilistic graphical model. By exploiting the existence of fundamental constraints in patterns embedded in Euclidean Spaces, we prove that for exact point set matching a simple graphical model is equivalent to the full model. It is possible to show that exact probabilistic inference in this simple model has polynomial time complexity with respect to the number of elements in the patterns to be matched. This gives rise to a technique that for exact matching provably finds a global optimum in polynomial time for any dimensionality of the underlying Euclidean Space. Computational experiments comparing this technique with well-known probabilistic relaxation labeling show significant performance improvement for inexact matching. The proposed approach is significantly more robust under augmentation of the sizes of the involved patterns. In the absence of noise, the results are always perfect.

Two additions to Lucas's "inflation and welfare"

This work adds to Lucas (2000) by providing analytical solutions to two problems that are solved only numerically by the author. The first part uses a theorem in control theory (Arrow' s sufficiency theorem) to provide sufficiency conditions to characterize the optimum in a shopping-time problem where the value function need not be concave. In the original paper the optimality of the first-order condition is characterized only by means of a numerical analysis. The second part of the paper provides a closed-form solution to the general-equilibrium expression of the welfare costs of inflation when the money demand is double logarithmic. This closed-form solution allows for the precise calculation of the difference between the general-equilibrium and Bailey's partial-equilibrium estimates of the welfare losses due to inflation. Again, in Lucas's original paper, the solution to the general-equilibrium-case underlying nonlinear differential equation is done only numerically, and the posterior assertion that the general-equilibrium welfare figures cannot be distinguished from those derived using Bailey's formula rely only on numerical simulations as well.

Análise de projetos de investimentos : uma abordagem teórico-prática

The objective of this work is to search a real case of capital budgeting, relating the practical technical aspects of the elaboration of project, with theoretical referential and following secondary objectives: (i) to analyze the relations established between the bibliographical material and the found practical technical problems of capital budgeting in the enterprise; (ii) to search and to describe the necessary pacing to the economic and financial elaboration of an project, from the prospecting of the demand, the projection of revenues and expenditures and the evaluation of the necessary investments to its development; (iii) to relate and to exemplify the influences of the restrictions presented for the methods of capital budgeting, correlating the practical theoretical referential with the enterprise; (iv) to analyze the yield of the investment project, (v) to verify the influence of the financing, on the yield of the project; and, finally, (vi) to demonstrate the choice process among some alternatives of supply, when used as tools of aid to the purchase decision, the methods of the Internal Tax of Return and the Net Present Value. To the end of the study one concluded that the methods of the Internal Tax of Return and the Net Present Value are powerful tools in the yield evaluation and viability of investments projects. However, to only understand the methods through what they teach in books is not enough for the daily practical of capital budgeting. Literature starts from two basic points: (i) the investments analyst dominates all the countable revenues, expenditures, and investments concepts.(ii) the numerical examples are simple and easy to understand, to infer its practical applications is a contouring question to be raised and passed by the analyst. This study intends to show the conjunction of the bibliography with the practical one, therefore, from the instant that demonstrates the countable concept of the prescription, it also explains as it was constituted from the calculation of the demand, until its inclusion in the project. Thus, searching concepts of revenues, expenditures, depreciation and capital assets, disclosing its constitution and, over all, the application inside of the project, it all takes the analyst to the final part of the process, that consists in the determination of the numerical calculations, allowing to dedicate more time to the difficult task to interpret the data. Finally, understood the analysis of the economic viability of the project, the study guides the purchase of the equipment under the economic-financial point of view.

A cultura do girassol como estratégia de competitividade para o agronegócio regional e nacional: importância para a agroenergia e a alimentação

Equacionar aquecimento global, escassez de alimentos e a crescente necessidade energética, tornou-se, atualmente, o grande desafio mundial. Existem diversas culturas agrícolas que podem ser exploradas de maneira estratégica e assim colaborar com a solução deste problema. Dentre elas, pode-se destacar a cultura do girassol (Helianthus Annuus). O girassol é uma das quatro maiores culturas oleaginosas no mundo, cultivado com sucesso nos cinco continentes, ocupando uma área de cultivo superior a 22 milhões de hectares. A participação do Brasil nesse montante é inferior a 1%. Acredita-se que essa pequena participação se deva a fatores sócio-econômicos e tecnológicos. Salienta-se, porém que o Brasil, por suas vantagens comparativas naturais e vantagens competitivas construídas possui condições favoráveis para seu desenvolvimento. Diante dos fatos, o objetivo deste trabalho é aprofundar o conhecimento da cadeia produtiva do girassol e através de sua utilização como estratégia de competitividade, avaliar de maneira sistêmica os impactos na matriz agrícola do País. Dentre as inúmeras vantagens dessa cultura, pode-se destacar: características agronômicas, físicas, químicas, organolépticas e versatilidade, que permitem a utilização e otimização dos fatores de produção já disponíveis; época de plantio (adaptabilidade a diferentes condições edafoclimáticas), podendo ser cultivado desde o Rio Grande do Sul até o Estado de Roraima; sistema radicular (raiz pode chegar a dois metros de profundidade), permitindo o melhor aproveitamento dos nutrientes e da água do solo e promovendo a reciclagem de nutrientes; alto teor de óleo nas sementes (30% a 55%) e; alto valor comercial dos co-produtos. Esse conjunto de características é analisado sob a ótica da teoria das vantagens competitivas e das economias de escala e escopo, demonstrando que com inteligência e pragmatismo, a cultura do girassol pode repetir, com vantagens, o que a soja representou para o agronegócio brasileiro.