A real symmetric Lanczos subspace implementation of quantum scattering using boundary inhomogeneities

We present an efficient and robust method for calculating state-to-state reaction probabilities utilising the Lanczos algorithm for a real symmetric Hamiltonian. The method recasts the time-independent Artificial Boundary Inhomogeneity technique recently introduced by Jang and Light (J. Chem. Phys. 102 (1995) 3262) into a tridiagonal (Lanczos) representation. The calculation proceeds at the cost of a single Lanczos propagation for each boundary inhomogeneity function and yields all state-to-state probabilities (elastic, inelastic and reactive) over an arbitrary energy range. The method is applied to the collinear H + H-2 reaction and the results demonstrate it is accurate and efficient in comparison with previous calculations. (C) 2002 Elsevier Science B.V. All rights reserved.

Outcome of stroke patients admitted to intensive care: Experience from an Australian teaching hospital

The objective of this study was to determine the mortality rate and the functional outcomes of stroke patients admitted to the intensive care unit (ICU) and to identify predictors of poor outcome in this population. The records of all patients admitted to the ICU with the diagnosis of stroke between January 1994 and December 1999 were reviewed. Patients with subarachnoid haemorrhage were excluded. Data were collected on clinical and biological variables, risk factors for stroke and the presence of comorbidities. Mortality (ICU, in-hospital and three-month) and functional outcome were used as end-points. In the six-year-period, 61 patients were admitted to the ICU with either haemorrhagic or ischaemic stroke. Medical records were available for only 58 patients. There were 23 ischaemic and 35 haemorrhagic strokes. The ICU, in-hospital and three-month mortality rates were 36%, 47% and 52% respectively. There were no significant differences in the prevalence of premorbid risk factors between survivors and non-survivors. The mean Barthel score was significantly different between the independent and dependent survivors (94 +/- 6 vs 45 +/- 26, P < 0.001). A substantial number of patients with good functional outcomes had lower Rankin scores (92% vs 11%, P < 0.001). Only 46% of those who were alive at three months were functionally independent. Intensive care admission was associated with a high mortality rate and a high likelihood of dependent lifestyle after hospital discharge. Haemorrhagic stroke, fixed dilated pupil(s) and GCS < 10 during assessment were associated with increased mortality and poor functional outcome.

Neural compensation for compliant loads during rhythmic movement

An experiment was performed to characterise the movement kinematics and the electromyogram (EMG) during rhythmic voluntary flexion and extension of the wrist against different compliant (elastic-viscous-inertial) loads. Three levels of each type of load, and an unloaded condition, were employed. The movements were paced at a frequency of I Hz by an auditory metronome, and visual feedback of wrist displacement in relation to a target amplitude of 100degrees was provided. Electro-myographic recordings were obtained from flexor carpi radialis (FCR) and extensor carpi radialis brevis (ECR). The movement profiles generated in the ten experimental conditions were indistinguishable, indicating that the CNS was able to compensate completely for the imposed changes in the task dynamics. When the level of viscous load was elevated, this compensation took the form of an increase in the rate of initial rise of the flexor and the extensor EMG burst. In response to increases in inertial load, the flexor and extensor EMG bursts commenced and terminated earlier in the movement cycle, and tended to be of greater duration. When the movements were performed in opposition to an elastic load, both the onset and offset of EMG activity occurred later than in the unloaded condition. There was also a net reduction in extensor burst duration with increases in elastic load, and an increase in the rate of initial rise of the extensor burst. Less pronounced alterations in the rate of initial rise of the flexor EMG burst were also observed. In all instances, increases in the magnitude of the external load led to elevations in the overall level of muscle activation. These data reveal that the elements of the central command that are modified in response to the imposition of a compliant load are contingent, not only upon the magnitude, but also upon the character of the load.

Investigation of ultrasonic properties of PAG and MAGIC polymer gel dosimeters

Ultrasonic speed of propagation and attenuation were investigated as a function of absorbed radiation dose in PAG and MAGIC polymer gel dosimeters. Both PAG and MAGIC gel dosimeters displayed a dependence of ultrasonic parameters on absorbed dose with attenuation displaying significant changes in the dose range investigated. The ultrasonic attenuation dose sensitivity at 4 MHz in MAGIC gels was determined to be 4.7 +/- 0.3 dB m(-1) Gy(-1) and for PAG 3.9 +/- 0.3 dB m(-1) Gy(-1). Ultrasonic speed dose sensitivities were 0.178 +/- 0.006 m s(-1) Gy(-1) for MAGIC gel and -0.44 +/- 0.02 m s(-1) Gy(-1) for PAG. Density and compressional elastic modulus were investigated to explain the different sensitivities of ultrasonic speed to radiation for PAG and MAGIC gels. The different sensitivities were found to be due to differences in the compressional elastic modulus as a function of dose for the two formulations. To understand the physical phenomena underlying the increase in ultrasonic attenuation with dose, the viscoelastic properties of the gels were studied. Results suggest that at ultrasonic frequencies, attenuation in polymer gel dosimeters is primarily due to volume viscosity. It is concluded that ultrasonic attenuation significantly increases with absorbed dose. Also, the ultrasonic speed in polymer gel dosimeters is affected by changes in dosimeter elastic modulus that are likely to be a result of polymerization. It is suggested that ultrasound is a sufficiently sensitive technique for polymer gel dosimetry.

Simulation of the load-unload response ratio and critical sensitivity in the lattice solid model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

Leaf water relations of Eucalyptus cloeziana and Eucalyptus argophloia in response to water deficit

Leaf water relations responses to limited water supply were determined in 7-month-old plants of a dry inland provenance of Eucalyptus argophloia Blakely and in a humid coastal provenance (Gympie) and a dry inland provenance (Hungry Hills) of Eucalyptus cloeziana F. Muell. Each provenance of E. cloeziana exhibited a lower relative water content at the turgor loss point, a lower apoplastic water content, a smaller ratio of dry mass to turgid mass and a lower bulk modulus of elasticity than the single provenance of E. argophloia. Osmotic potential at full turgor and water potential at the turgor loss point were significantly lower in E. argophloia and the inland provenance of E. cloeziana than in the coastal provenance of E. cloeziana. There was limited osmotic adjustment in response to soil drying in E. cloeziana, but not in E. argophloia. Between-species differences in water relations parameters were larger than those between the E. cloeziana provenances. Both E. cloeziana provenances maintained turgor under moderate water stress through a combination of osmotic and elastic adjustments. Eucalyptus argophloia had more rigid cell walls and reached lower water potentials with less reduction in relative water content than either of the E. cloeziana provenances, thereby enabling it to extract water from dryer soils.

Difference equations in Banach spaces

Difference equations which discretely approximate boundary value problems for second-order ordinary differential equations are analysed. It is well known that the existence of solutions to the continuous problem does not necessarily imply existence of solutions to the discrete problem and, even if solutions to the discrete problem are guaranteed, they may be unrelated and inapplicable to the continuous problem. Analogues to theorems for the continuous problem regarding a priori bounds and existence of solutions are formulated for the discrete problem. Solutions to the discrete problem are shown to converge to solutions of the continuous problem in an aggregate sense. An example which arises in the study of the finite deflections of an elastic string under a transverse load is investigated. The earlier results are applied to show the existence of a solution; the sufficient estimates on the step size are presented. (C) 2003 Elsevier Science Ltd. All rights reserved.

Theoretical and experimental studies of unbraced tubular trusses allowing for torsional stiffness

This paper describes the buckling phenomenon of a tubular truss with unsupported length through a full-scale test and presents a practical computational method for the design of the trusses allowing for the contribution of torsional stiffness against buckling, of which the effect has never been considered previously by others. The current practice for the design of a planar truss has largely been based on the linear elastic approach which cannot allow for the contribution of torsional stiffness and tension members in a structural system against buckling. The over-simplified analytical technique is unable to provide a realistic and an economical design to a structure. In this paper the stability theory is applied to the second-order analysis and design of the structural form, with detailed allowance for the instability and second-order effects in compliance with design code requirements. Finally, the paper demonstrates the application of the proposed method to the stability design of a commonly adopted truss system used in support of glass panels in which lateral bracing members are highly undesirable for economical and aesthetic reasons.

Comparação de desempenho entre a formulação direta do Método dos Elementos de Contorno com Funções Radiais e o Método dos Elementos Finitos em problemas de Poisson e Helmholtz

O presente trabalho objetiva avaliar o desempenho do MECID (Método dos Elementos de Contorno com Interpolação Direta) para resolver o termo integral referente à inércia na Equação de Helmholtz e, deste modo, permitir a modelagem do Problema de Autovalor assim como calcular as frequências naturais, comparando-o com os resultados obtidos pelo MEF (Método dos Elementos Finitos), gerado pela Formulação Clássica de Galerkin. Em primeira instância, serão abordados alguns problemas governados pela equação de Poisson, possibilitando iniciar a comparação de desempenho entre os métodos numéricos aqui abordados. Os problemas resolvidos se aplicam em diferentes e importantes áreas da engenharia, como na transmissão de calor, no eletromagnetismo e em problemas elásticos particulares. Em termos numéricos, sabe-se das dificuldades existentes na aproximação precisa de distribuições mais complexas de cargas, fontes ou sorvedouros no interior do domínio para qualquer técnica de contorno. No entanto, este trabalho mostra que, apesar de tais dificuldades, o desempenho do Método dos Elementos de Contorno é superior, tanto no cálculo da variável básica, quanto na sua derivada. Para tanto, são resolvidos problemas bidimensionais referentes a membranas elásticas, esforços em barras devido ao peso próprio e problemas de determinação de frequências naturais em problemas acústicos em domínios fechados, dentre outros apresentados, utilizando malhas com diferentes graus de refinamento, além de elementos lineares com funções de bases radiais para o MECID e funções base de interpolação polinomial de grau (um) para o MEF. São geradas curvas de desempenho através do cálculo do erro médio percentual para cada malha, demonstrando a convergência e a precisão de cada método. Os resultados também são comparados com as soluções analíticas, quando disponíveis, para cada exemplo resolvido neste trabalho.

Determinação do momento crítico de flambagem lateral com distorção em vigas mistas contínuas de aço e concreto com perfis de alma senoidal

As vigas mistas de aço e concreto estão sendo largamente utilizadas em construções de edifícios e pontes. Ao se combinar o aço com o concreto obtêm-se estruturas mais econômicas, uma vez que se tira proveito das melhores características de cada material. Nas regiões de momento negativo de uma viga mista contínua, a mesa inferior e parte da alma estão comprimidas, se a alma do perfil não tiver rigidez suficiente para evitar a flexão lateral, ela distorcerá gerando um deslocamento lateral e um giro na mesa comprimida, caracterizando um modo de flambagem denominado flambagem lateral com distorção (FLD). O procedimento de verificação à FLD da EN 1994-1-1:2004 originou o método de cálculo da ABNT NBR 8800:2008, entretanto a EN 1994-1-1:2004 não fornece expressão para o cálculo do momento crítico elástico, enquanto a ABNT NBR 8800:2008 prescreve uma formulação proposta por Roik, Hanswille e Kina (1990) desenvolvida para vigas mistas com perfis de alma plana. Embora as normas prescrevam um método de verificação à FLD para vigas mistas com perfis de alma plana, poucos estudos têm sido feitos sobre esse estado-limite. Além disso, tanto a ABNT NBR 8800:2008 quanto as normas internacionais não abordam perfis de alma senoidal. Neste trabalho, foram implementadas análises de flambagem elástica, com auxílio do software ANSYS 14.0 (2011), em modelos de elementos finitos que retratem o comportamento à FLD de vigas mistas de aço e concreto com perfis de alma plana e senoidal. Os modelos numéricos foram constituídos pelo perfil de aço, por uma mola rotacional que restringe parcialmente o giro da mesa superior e uma restrição ao deslocamento lateral, ao longo de todo o comprimento da viga. Os resultados numéricos são comparados com os obtidos pelas formulações de Roik, Hanswille e Kina (1990) e de Hanswille (2002), adaptadas para levar em consideração a corrugação da alma do perfil de aço. Para avaliação das formulações supracitadas e da consistência da modelagem numérica adotada, o momento crítico elástico foi determinado para vigas mistas com perfis de aço de alma plana. Como resultado, um método para o cálculo do momento crítico elástico de vigas mistas de alma senoidal é proposto.

Processing and characterization of alpha-elastin electrospun membranes

Elastin isolated from fresh bovine ligaments was dissolved in a mixture of 1,1,1,3,3,3-Hexafluoro-2-propanol and water and electrospun into fiber membranes under different processing conditions. Fiber mats of randomly and aligned fibers were obtained with fixed and rotating ground collectors and fibrils were composed by thin ribbons whose width depends on electrospinning conditions; fibrils with 721 nm up to 2.12 m width were achieved. After cross-linking with glutaraldehyde, -elastin can uptake as much as 1700 % of PBS solution and a slight increase on fiber thickness was observed. The glass transition temperature of electrospun fiber mats was found to occur at ~ 80 ºC. Moreover, -Elastin showed to be a perfect elastomeric material, and no mechanical hysteresis was found in cycle mechanical measurements. The elastic modulus obtained for oriented and random fibers mats in a PBS solution was 330 ± 10 kPa and 732 ± 165 kPa, respectively. Finally, the electrospinning and cross-linking process does not inhibit MC-3T3-E1 cell adhesion. Cell culture results showed good cell adhesion and proliferation in the cross-linked elastin fiber mats.

Surgical correction of scoliosis: Numerical analysis and optimization of the procedure

A previously developed model is used to numerically simulate real clinical cases of the surgical correction of scoliosis. This model consists of one-dimensional finite elements with spatial deformation in which (i) the column is represented by its axis; (ii) the vertebrae are assumed to be rigid; and (iii) the deformability of the column is concentrated in springs that connect the successive rigid elements. The metallic rods used for the surgical correction are modeled by beam elements with linear elastic behavior. To obtain the forces at the connections between the metallic rods and the vertebrae geometrically, non-linear finite element analyses are performed. The tightening sequence determines the magnitude of the forces applied to the patient column, and it is desirable to keep those forces as small as possible. In this study, a Genetic Algorithm optimization is applied to this model in order to determine the sequence that minimizes the corrective forces applied during the surgery. This amounts to find the optimal permutation of integers 1, ... , n, n being the number of vertebrae involved. As such, we are faced with a combinatorial optimization problem isomorph to the Traveling Salesman Problem. The fitness evaluation requires one computing intensive Finite Element Analysis per candidate solution and, thus, a parallel implementation of the Genetic Algorithm is developed.

A two-Higgs doublet model with remarkable CP properties

We analyze generalized CP symmetries of two-Higgs doublet models, extending them from the scalar to the fermion sector of the theory. We show that, other than the usual CP transformation, there is only one of those symmetries which does not imply massless charged fermions. That single model which accommodates a fermionic mass spectrum compatible with experimental data possesses a remarkable feature. Through a soft breaking of the symmetry it displays a new type of spontaneous CP violation, which does not occur in the scalar sector responsible for the symmetry breaking mechanism but, rather, in the fermion sector.

Active and passive behaviors of soft tissues: Pelvic floor muscles

: A new active-contraction visco-elastic numerical model of the pelvic floor (skeletal) muscle is presented. Our model includes all elements that represent the muscle constitutive behavior, contraction and relaxation. In contrast with the previous models, the activation function can be null. The complete equations are shown and exactly linearized. Small verification and validation tests are performed and the pelvis is modeled using the data from the intra-abdominal pressure tests

Simulação numérica do comportamento dinâmico de uma laje

Nesta dissertação pretende-se simular o comportamento dinâmico de uma laje de betão armado aplicando o Método de Elementos Finitos através da sua implementação no programa FreeFEM++. Este programa permite-nos a análise do modelo matemático tridimensional da Teoria da Elasticidade Linear, englobando a Equação de Equilíbrio, Equação de Compatibilidade e Relações Constitutivas. Tratando-se de um problema dinâmico é necessário recorrer a métodos numéricos de Integração Directa de modo a obter a resposta em termos de deslocamento ao longo do tempo. Para este trabalho escolhemos o Método de Newmark e o Método de Euler para a discretização temporal, um pela sua popularidade e o outro pela sua simplicidade de implementação. Os resultados obtidos pelo FreeFEM++ são validados através da comparação com resultados adquiridos a partir do SAP2000 e de Soluções Teóricas, quando possível.