A Monte Carlo study of the anisotropic N=3 Ashkin-Teller model

The phase diagram of an asymmetric N = 3 Ashkin-Teller model is obtained by a numerical analysis which combines Monte Carlo renormalization group and reweighting techniques. Present results reveal several differences with those obtained by mean-field calculations and a Hamiltonian approach. In particular, we found Ising critical exponents along a line where Goldschmidt has located the Kosterlitz-Thouless multicritical point. On the other hand, we did find nonuniversal exponents along another transition line. Symmetry breaking in this case is very similar to the N = 2 case, since the symmetries associated to only two of the Ising variables are broken. However, for large values of the coupling constant ratio XW = W/K, when the only broken symmetry is of a hidden variable, we detected first-order phase transitions giving evidence supporting the existence of a multicritical point, as suggested by Goldschmidt, but in a different region of the phase diagram. © 2002 Elsevier Science B.V. All rights reserved.

Renormalization-group calculation of dynamical properties for impurity models

The numerical renormalization-group method was originally developed to calculate the thermodynamical properties of impurity Hamiltonians. A recently proposed generalization capable of computing dynamical properties is discussed. As illustrative applications, essentially exact results for the impurity specttral densities of the spin-degenerate Anderson model and of a model for electronic tunneling between two centers in a metal are presented. © 1991.

Joint economic design of X̄ and R control charts for processes subject to two independent assignable causes

A model for the joint economic design of X̄ and R control charts is developed. This model assumes that the process is subject to two assignable causes. One assignable cause shifts the process mean; the other shifts the process variance. The occurrence of the assignable cause of one kind does not block the occurrence of the assignable cause of another kind. Consequently, a second process parameter can go out-of-control after the first process parameter has gone out-of-control. A numerical study of the cost surface to the model considered has revealed that it is convex, at least in the interest region.

Análise numérica de uma geometria de coletor solar para sistemas de aquecimento de água

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

Análise numérica de um sistema integrado coletor solar/armazenador térmico

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

Estudo de sistemas não-lineares em mecânica quântica

The main goal of this work is to investigate the effects of a nonlinear cubic term inserted in the Schrödinger equation for one-dimensional potentials studied in Quantum Mechanics textbooks. Being the main tool the numerical analysis in a large number of works, the analysis of this effect by this term in the potential itself, in order to work with an analytical solution, can be considered something new. For the harmonic oscillator potential, the analysis was made from a numerical method, comparing the result with the known results in the literature. In the case of the infinite well potential and the step potential, hoping to work with an analytical solution, by construction we started with the known wavefunction for the linear case noting the effects in the other physical quantities. The coupling of the physical quantities involved in this work has yielded, besides many complications in the calculations, a series of conditions on the existence and validity of the solutions in regard to the system possible configurations

Contribuição ao estudo de lajes lisas unidirecionais de concreto armado com pilares retangulares ao puncionamento

São apresentadas análises das recomendações para estimar a resistência ao puncionamento de lajes lisas de acordo com a versão atual da norma brasileira, a NBR 6118:2003, que não considera o efeito da retangularidade de pilares internos, e sua predecessora, a NBR 6118:1980. Em seções transversais retangulares, o índice de polarização das tensões de cisalhamento é crescente. Isto altera a forma de ruptura da laje, já que a punção se inicia em torno das extremidades da seção do pilar, não se estendendo para os lados maiores, como prescreve a norma, comportamento este que reduz a resistência última das lajes frente às estimativas normativas. São apresentados os resultados de uma análise numérica por elementos finitos (M.E.F.) e uma proposta para melhorar a ductilidade das lajes lisas através das dimensões mínimas dos pilares. Observou-se que a NBR 6118:1980 mostrou-se conservadora e que as estimativas da NBR 6118: 2003 podem ser melhoradas.

Estudo de sistemas não-lineares em mecânica quântica

The main goal of this work is to investigate the effects of a nonlinear cubic term inserted in the Schrödinger equation for one-dimensional potentials studied in Quantum Mechanics textbooks. Being the main tool the numerical analysis in a large number of works, the analysis of this effect by this term in the potential itself, in order to work with an analytical solution, can be considered something new. For the harmonic oscillator potential, the analysis was made from a numerical method, comparing the result with the known results in the literature. In the case of the infinite well potential and the step potential, hoping to work with an analytical solution, by construction we started with the known wavefunction for the linear case noting the effects in the other physical quantities. The coupling of the physical quantities involved in this work has yielded, besides many complications in the calculations, a series of conditions on the existence and validity of the solutions in regard to the system possible configurations

Localising and quantifying damage by means of a multi-chromosome genetic algorithm

This paper presents a structural damage detection methodology based on genetic algorithms and dynamic parameters. Three chromosomes are used to codify an individual in the population. The first and second chromosomes locate and quantify damage, respectively. The third permits the self-adaptation of the genetic parameters. The natural frequencies and mode shapes are used to formulate the objective function. A numerical analysis was performed for several truss structures under different damage scenarios. The results have shown that the methodology can reliably identify damage scenarios using noisy measurements and that it results in only a few misidentified elements. (C) 2012 Civil-Comp Ltd and Elsevier Ltd. All rights reserved.

Performance of polymeric reinforcements in vehicle structures submitted to frontal impact

Preformed structural reinforcements have shown good performance in crash tests, where the great advantage is their weight. These reinforcements are designed with the aim of increasing the rigidity of regions with large deformations, thus stabilising sections of the vehicle that work as load path during impact. The objective of this work is to show the application of structural reinforcements made of polymeric material PA66 in the field of vehicle safety, through finite element simulations. Simulations of frontal impact at 50 km/h and in ODB (offset deformable barrier) at 57 km/h configurations (standards such as ECE R-94 and ECE R-12) were performed in the software LS-DYNA R (R) and MADYMO (R). The simulations showed that the use of polymeric reinforcements leads to a 70% reduction in A-pillar intrusion, a 65% reduction in the displacement of the steering column and a 59% reduction in the deformation in the region of the occupant legs and feet. The level of occupant injuries was analysed by MADYMO (R) software, and a reduction of 23.5% in the chest compression and 80% in the tibia compression were verified. According to the standard, such conditions lead to an improvement in the occupant safety in a vehicle collision event.

BEM modeling of saturated porous media susceptible to damage

This paper deals with the numerical analysis of saturated porous media, taking into account the damage phenomena on the solid skeleton. The porous media is taken into poro-elastic framework, in full-saturated condition, based on Biot's Theory. A scalar damage model is assumed for this analysis. An implicit boundary element method (BEM) formulation, based on time-independent fundamental solutions, is developed and implemented to couple the fluid flow and two-dimensional elastostatic problems. The integration over boundary elements is evaluated using a numerical Gauss procedure. A semi-analytical scheme for the case of triangular domain cells is followed to carry out the relevant domain integrals. The non-linear problem is solved by a Newton-Raphson procedure. Numerical examples are presented, in order to validate the implemented formulation and to illustrate its efficacy. (C) 2011 Elsevier Ltd. All rights reserved.

Interpolation estimate for a finite-element space with embedded discontinuities

We consider a recently proposed finite-element space that consists of piecewise affine functions with discontinuities across a smooth given interface Γ (a curve in two dimensions, a surface in three dimensions). Contrary to existing extended finite element methodologies, the space is a variant of the standard conforming Formula space that can be implemented element by element. Further, it neither introduces new unknowns nor deteriorates the sparsity structure. It is proved that, for u arbitrary in Formula, the interpolant Formula defined by this new space satisfies Graphic where h is the mesh size, Formula is the domain, Formula, Formula, Formula and standard notation has been adopted for the function spaces. This result proves the good approximation properties of the finite-element space as compared to any space consisting of functions that are continuous across Γ, which would yield an error in the Formula-norm of order Graphic. These properties make this space especially attractive for approximating the pressure in problems with surface tension or other immersed interfaces that lead to discontinuities in the pressure field. Furthermore, the result still holds for interfaces that end within the domain, as happens for example in cracked domains.

CP asymmetries in the supersymmetric trilepton signal at the LHC

In the CP-violating Minimal Supersymmetric Standard Model, we study the production of a neutralino-chargino pair at the LHC. For their decays into three leptons, we analyze CP asymmetries which are sensitive to the CP phases of the neutralino and chargino sector. We present analytical formulas for the entire production and decay process, and identify the CP-violating contributions in the spin correlation terms. This allows us to define the optimal CP asymmetries. We present a detailed numerical analysis of the cross sections, branching ratios, and the CP observables. For light neutralinos, charginos, and squarks, the asymmetries can reach several 10%. We estimate the discovery potential for the LHC to observe CP violation in the trilepton channel.

A swirler stabilized combustion chamber for a micro-gas turbine fuelled with natural gas

Micro-gas turbines are a good alternative for on-site power generation, since their operation is very reliable. The possibility of operating with various fuels increases versatility and, as a result, the usage of these devices. Focusing on a performance improvement of a tri-fuel low-cost micro-gas turbine, this work presents investigations of the inner flow of its combustion chamber. The aim of this analysis was the characterization of the flame structure by the temperature field of the chamber inner flow. The chamber was fuelled with natural gas. In the current chamber, a swirler and a reversed flow configuration were utilized to provide flame stabilization. The inner flow investigations were done with numerical analysis, which were compared to experimental data. The analysis of the inner flow was done with numerical simulations, which used the RSM turbulence model. A β-PDF equilibrium model was adopted to account for the turbulent combustion process. Different models of heat transfer were compared. Thermal radiation and specially heat conduction in the liner walls played significant roles on results.

Development of musculoskeletal models for the design and the pre-clinical validation of hip resurfacing prosthesis

Background. The surgical treatment of dysfunctional hips is a severe condition for the patient and a costly therapy for the public health. Hip resurfacing techniques seem to hold the promise of various advantages over traditional THR, with particular attention to young and active patients. Although the lesson provided in the past by many branches of engineering is that success in designing competitive products can be achieved only by predicting the possible scenario of failure, to date the understanding of the implant quality is poorly pre-clinically addressed. Thus revision is the only delayed and reliable end point for assessment. The aim of the present work was to model the musculoskeletal system so as to develop a protocol for predicting failure of hip resurfacing prosthesis. Methods. Preliminary studies validated the technique for the generation of subject specific finite element (FE) models of long bones from Computed Thomography data. The proposed protocol consisted in the numerical analysis of the prosthesis biomechanics by deterministic and statistic studies so as to assess the risk of biomechanical failure on the different operative conditions the implant might face in a population of interest during various activities of daily living. Physiological conditions were defined including the variability of the anatomy, bone densitometry, surgery uncertainties and published boundary conditions at the hip. The protocol was tested by analysing a successful design on the market and a new prototype of a resurfacing prosthesis. Results. The intrinsic accuracy of models on bone stress predictions (RMSE < 10%) was aligned to the current state of the art in this field. The accuracy of prediction on the bone-prosthesis contact mechanics was also excellent (< 0.001 mm). The sensitivity of models prediction to uncertainties on modelling parameter was found below 8.4%. The analysis of the successful design resulted in a very good agreement with published retrospective studies. The geometry optimisation of the new prototype lead to a final design with a low risk of failure. The statistical analysis confirmed the minimal risk of the optimised design over the entire population of interest. The performances of the optimised design showed a significant improvement with respect to the first prototype (+35%). Limitations. On the authors opinion the major limitation of this study is on boundary conditions. The muscular forces and the hip joint reaction were derived from the few data available in the literature, which can be considered significant but hardly representative of the entire variability of boundary conditions the implant might face over the patients population. This moved the focus of the research on modelling the musculoskeletal system; the ongoing activity is to develop subject-specific musculoskeletal models of the lower limb from medical images. Conclusions. The developed protocol was able to accurately predict known clinical outcomes when applied to a well-established device and, to support the design optimisation phase providing important information on critical characteristics of the patients when applied to a new prosthesis. The presented approach does have a relevant generality that would allow the extension of the protocol to a large set of orthopaedic scenarios with minor changes. Hence, a failure mode analysis criterion can be considered a suitable tool in developing new orthopaedic devices.