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.

Field effect and light-assisted a-Si:H sensors for detection of ions in solution

In this paper we present an amorphous silicon device that can be used in two operation modes to measure the concentration of ions in solution. While crystalline devices present a higher sensitivity, their amorphous counterpart present a much lower fabrication cost, thus enabling the production of cheap disposable sensors for use, for example, in the food industry. The devices were fabricated on glass substrates by the PECVD technique in the top gate configuration, where the metallic gate is replaced by an electrolytic solution with an immersed Ag/AgCl reference electrode. Silicon nitride is used as gate dielectric enhancing the sensitivity and passivation layer used to avoid leakage and electrochemical reactions. In this article we report on the semiconductor unit, showing that the device can be operated in a light-assisted mode, where changes in the pH produce changes on the measured ac photocurrent. In alternative the device can be operated as a conventional ion selective field effect device where changes in the pH induce changes in the transistor's threshold voltage.

Enthalpies of Solution of 1-Butyl-3-methylimidazolium Tetrafluoroborate in 15 Solvents at 298.15 K

Enthalpies of solution of 1-butyl-3-methylimidazolium tetra fluoroborate, [BMIm]BF4, are reported at 298.15 K in a set of 15 hydrogen bond donor and hydrogen bond acceptor solvents, chosen by their diversity, namely, water, methanol, ethanol, 1,2-ethanediol, 2-choroethanol, 2-methoxyethanol, formamide, propylene carbonate, nitromethane, acetonitrile, dimethyl sulfoxide, acetone, N,N-dimethylformamide, N,N-dimethylacetamide, and aniline. These values are shown to be largely independent of [BMIm]BF4 concentration. The obtained enthalpies of solution vary from very endothermic to quite exothermic, thus showing a very high sensitivity of the enthalpies of solution of [BMIm]BF4 to solvent properties. Solvent effects on the solution process of this IL are analyzed by a quantitative structure-property relationship methodology, using the TAKA equation and a modified equation, which significantly improves the model's predictive ability. The observed differences in the enthalpies of solution are rationalized in terms of the solvent properties found to be relevant, that is, pi* and E-T(N).

Solvent effects on solution enthalpies of adamantyl derivatives

Solution enthalpies of 1-bromoadamantane, 1-adamantanol, and 2-adamantanone in a large set of protic and aprotic solvents are reported at 298.15 K. Solvent effects on the solution processes of these solutes are analyzed in terms of a modified TAKA equation, involving delta(cav) h (s) as the cavity term. The nature and magnitude of the major interactions which influence these processes are assessed and discussed in terms of the solutes' characteristics. New insights on the solution processes under scrutiny are presented.

Mathematical modeling of pressurized system behaviour with entrapped air

The presence of entrapped air in pressurized hydraulic systems is considered a critical condition for the infrastructure security, due to the transient pressure enhancement related with its dynamic behaviour, similar to non-linear spring action. A mathematical model for the assessment of hydraulic transients resulting from rapid pressurizations, under referred condition is presented. Water movement was modeled through the elastic column theory considering a moving liquid boundary and the entrapped air pocket as lumped gas mass, where the acoustic effects are negligible. The method of characteristics was used to obtain the numerical solution of the liquid flow. The resulting model is applied to an experimental set-up having entrapped air in the top of a vertical pipe section and the numerical results are analyzed.

Solution enthalpies of hydroxylic compounds

Solution enthalpies of adamantan-1-ol, 2-methyl- butan-2-ol, and 3-methylbutan-1-ol have been measured at 298.15 K, in a set of 16 protogenic and non-protogenic solvents. The identification and quantification of solvent effects on the solution processes under study were performed using quantitative-structure property relationships. The results are discussed in terms of solute-solvent-solvent interactions and also in terms of the influence of compound's size and position of its hydroxyl group.

A strategy to extend reactive distillation column performance under catalyst deactivation

This work addresses the effects of catalyst deactivation and investigates methods to reduce their impact on the reactive distillation columns performance. The use of variable feed quality and reboil ratio are investigated using a rigorous dynamic model developed in gPROMS and applied to an illustrative example, i.e., the olefin metathesis system, wherein 2-pentene reacts to form 2-butene and 3-hexene. Three designs and different strategies on column energy supply to tackle catalyst deactivation are investigated and the results compared.

Fast optimum-predictive control and capacitor voltage balancing strategy for bipolar back-to-back NPC converters in high-voltage direct current transmission systems

Multilevel power converters have been introduced as the solution for high-power high-voltage switching applications where they have well-known advantages. Recently, full back-to-back connected multilevel neutral point diode clamped converters (NPC converter) have been used inhigh-voltage direct current (HVDC) transmission systems. Bipolar-connected back-to-back NPC converters have advantages in long-distance HVDCtransmission systems over the full back-to-back connection, but greater difficulty to balance the dc capacitor voltage divider on both sending and receiving end NPC converters. This study shows that power flow control and dc capacitor voltage balancing are feasible using fast optimum-predictive-based controllers in HVDC systems using bipolar back-to-back-connected five-level NPC multilevel converters. For both converter sides, the control strategytakes in account active and reactive power, which establishes ac grid currents in both ends, and guarantees the balancing of dc bus capacitor voltages inboth NPC converters. Additionally, the semiconductor switching frequency is minimised to reduce switching losses. The performance and robustness of the new fast predictive control strategy, and its capability to solve the DC capacitor voltage balancing problem of bipolar-connected back-to-back NPCconverters are evaluated.

Study of the one dimensional and transient bioheat transfer equation: Multi-layer solution development and applications

: In this work we derive an analytical solution given by Bessel series to the transient and one-dimensional (1D) bioheat transfer equation in a multi-layer region with spatially dependent heat sources. Each region represents an independent biological tissue characterized by temperature-invariant physiological parameters and a linearly temperature dependent metabolic heat generation. Moreover, 1D Cartesian, cylindrical or spherical coordinates are used to define the geometry and temperature boundary conditions of first, second and third kinds are assumed at the inner and outer surfaces. We present two examples of clinical applications for the developed solution. In the first one, we investigate two different heat source terms to simulate the heating in a tumor and its surrounding tissue, induced during a magnetic fluid hyperthermia technique used for cancer treatment. To obtain an accurate analytical solution, we determine the error associated with the truncated Bessel series that defines the transient solution. In the second application, we explore the potential of this model to study the effect of different environmental conditions in a multi-layered human head model (brain, bone and scalp). The convective heat transfer effect of a large blood vessel located inside the brain is also investigated. The results are further compared with a numerical solution obtained by the Finite Element Method and computed with COMSOL Multi-physics v4.1 (c). (c) 2013 Elsevier Ltd. All rights reserved.

Solution enthalpies of 1,4-dioxane: Study of solvent effects through quantitative structure-property relationships

Solution enthalpies of 1,4-dioxane have been obtained in 15 protic and aprotic solvents at 298.15 K. Breaking the overall process through the use of Solomonov's methodology the cavity term was calculated and interaction enthalpies (Delta H-int) were determined. Main factors involved in the interaction enthalpy have been identified and quantified using a QSPR approach based on the TAKA model equation. The relevant descriptors were found to be pi* and beta, which showed, respectively, exothermic and endothermic contributions. The magnitude of pi* coefficient points toward non-specific solute-solvent interactions playing a major role in the solution process. The positive value of the beta coefficient reflects the endothermic character of the solvents' hydrogen bond acceptor (HBA) basicity contribution, indicating that solvent molecules engaged in hydrogen bonding preferentially interact with each other rather than with 1,4-dioxane. (C) 2013 Elsevier B.V. All rights reserved.

Assessing the behaviour of RC beams subject to significant gravity loads under cyclic loads

Gravity loads can affect a reinforced concrete structure's response to seismic actions, however, traditional procedures for testing the beam behaviour do not take this effect into consideration. An experimental campaign was carried out in order to assess the influence of the gravity load on RC beam connection to the column subjected to cyclic loading. The experiments included the imposition of a conventional quasi-static test protocol based on the imposition of a reverse cyclic displacement history and of an alternative cyclic test procedure starting from the gravity load effects. The test results are presented, compared and analysed in this paper. The imposition of a cyclic test procedure that included the gravity loads effects on the RC beam ends reproduces the demands on the beams' critical zones more realistically than the traditional procedure. The consideration of the vertical load effects in the test procedure led to an accumulation of negative (hogging) deformation. This phenomenon is sustained with the behaviour of a portal frame system under cyclic loads subject to a significant level of the vertical load, leading to the formation of unidirectional plastic hinges. In addition, the hysteretic behaviour of the RC beam ends tested was simulated numerically using the nonlinear structural analysis software - OpenSees. The beam-column model simulates the global element behaviour very well, as there is a reasonable approximation to the hysteretic loops obtained experimentally. (C) 2013 Elsevier Ltd. All rights reserved.

The extended unsymmetric frontal solution for multiple-point constraints

The purpose of this paper is to discuss the linear solution of equality constrained problems by using the Frontal solution method without explicit assembling. Design/methodology/approach - Re-written frontal solution method with a priori pivot and front sequence. OpenMP parallelization, nearly linear (in elimination and substitution) up to 40 threads. Constraints enforced at the local assembling stage. Findings - When compared with both standard sparse solvers and classical frontal implementations, memory requirements and code size are significantly reduced. Research limitations/implications - Large, non-linear problems with constraints typically make use of the Newton method with Lagrange multipliers. In the context of the solution of problems with large number of constraints, the matrix transformation methods (MTM) are often more cost-effective. The paper presents a complete solution, with topological ordering, for this problem. Practical implications - A complete software package in Fortran 2003 is described. Examples of clique-based problems are shown with large systems solved in core. Social implications - More realistic non-linear problems can be solved with this Frontal code at the core of the Newton method. Originality/value - Use of topological ordering of constraints. A-priori pivot and front sequences. No need for symbolic assembling. Constraints treated at the core of the Frontal solver. Use of OpenMP in the main Frontal loop, now quantified. Availability of Software.

Anti-Corrosion performance of a new silane coating for corrosion protection of AZ31 Magnesium alloy in Hank’s solution

Mg alloys can be used as bioresorsable metallic implants. However, the high corrosion rate of magnesium alloys has limited their biomedical applications. Although Mg ions are essential to the human body, an excess may cause undesirable health effects. Therefore, surface treatments are required to enhance the corrosion resistance of magnesium parts, decreasing its rate to biocompatible levels and allowing its safe application as bioresorbable metallic implants. The application of biocompatible silane coatings is envisaged as a suitable strategy for retarding the corrosion process of magnesium alloys. In the current work, a new glycidoxypropyltrimethoxysilane (GPTMS) based coating was tested on AZ31 magnesium substrates subjected to different surface conditioning procedures before coating deposition. The surface conditioning included a short etching with hydrofluoric acid (HF) or a dc polarisation in alkaline electrolyte. The silane coated samples were immersed in Hank's solution and the protective performance of the coating was studied through electrochemical impedance spectroscopy (EIS). The EIS data was treated by new equivalent circuit models and the results revealed that the surface conditioning process plays a key role in the effectiveness of the silane coating. The HF treated samples led to the highest impedance values and delayed the coating degradation, compared to the mechanically polished samples or to those submitted to dc polarisation.

Solution enthalpies of hydroxylic compounds

Solution enthalpies of adamantan-1-ol, 2-methyl- butan-2-ol, and 3-methylbutan-1-ol have been measured at 298.15 K, in a set of 16 protogenic and non-protogenic solvents. The identification and quantification of solvent effects on the solution processes under study were performed using quantitative-structure property relationships. The results are discussed in terms of solute-solvent-solvent interactions and also in terms of the influence of compound's size and position of its hydroxyl group.

Is the electric vehicle a solution for the wind power integration in the portuguese power system?

The integration of Plug-in electric vehicles in the transportation sector has a great potential to reduce oil dependency, the GHG emissions and to contribute for the integration of renewable sources into the electricity generation mix. Portugal has a high share of wind energy, and curtailment may occur, especially during the off-peak hours with high levels of hydro generation. In this context, the electric vehicles, seen as a distributed storage system, can help to reduce the potential wind curtailments and, therefore, increase the integration of wind power into the power system. In order to assess the energy and environmental benefits of this integration, a methodology based on a unit commitment and economic dispatch is adapted and implemented. From this methodology, the thermal generation costs, the CO2 emissions and the potential wind generation curtailment are computed. Simulation results show that a 10% penetration of electric vehicles in the Portuguese fleet would increase electrical load by 3% and reduce wind curtailment by only 26%. This results from the fact that the additional generation required to supply the electric vehicles is mostly thermal. The computed CO2 emissions of the EV are 92 g CO2/kWh which become closer to those of some new ICE engines.