New inorganic nanotube polymer nanocomposites: Improved thermal, mechanical and tribological properties in isotactic polypropylene incorporating INT-MoS2.

Environmentally friendly molybdenum disulfide (INT-MoS2) inorganic nanotubes were introduced into an isotactic polypropylene (iPP) polymer matrix to generate novel nanocomposite materials through an advantageous melt-processing route. The effects of INT-MoS2 content on the thermal, mechanical and tribological properties were investigated. The incorporation of INT-MoS2 generates notable performance enhancements through reinforcement effects, highly efficient nucleation activity and excellent lubricating ability in comparison with other nanoparticle fillers such as nanoclays, carbon nanotubes, silicon nitrides and halloysite nanotubes. It was shown that these INT-MoS2 nanocomposites can provide an effective balance between performance, cost effectiveness and processability, and should be of some interest in the area of multifunctional polymer nanocomposite materials.

Monte Carlo Neutronics and Thermal-hydraulics coupling applied to Fast Reactors

Un escenario habitualmente considerado para el uso sostenible y prolongado de la energía nuclear contempla un parque de reactores rápidos refrigerados por metales líquidos (LMFR) dedicados al reciclado de Pu y la transmutación de actínidos minoritarios (MA). Otra opción es combinar dichos reactores con algunos sistemas subcríticos asistidos por acelerador (ADS), exclusivamente destinados a la eliminación de MA. El diseño y licenciamiento de estos reactores innovadores requiere herramientas computacionales prácticas y precisas, que incorporen el conocimiento obtenido en la investigación experimental de nuevas configuraciones de reactores, materiales y sistemas. A pesar de que se han construido y operado un cierto número de reactores rápidos a nivel mundial, la experiencia operacional es todavía reducida y no todos los transitorios se han podido entender completamente. Por tanto, los análisis de seguridad de nuevos LMFR están basados fundamentalmente en métodos deterministas, al contrario que las aproximaciones modernas para reactores de agua ligera (LWR), que se benefician también de los métodos probabilistas. La aproximación más usada en los estudios de seguridad de LMFR es utilizar una variedad de códigos, desarrollados a base de distintas teorías, en busca de soluciones integrales para los transitorios e incluyendo incertidumbres. En este marco, los nuevos códigos para cálculos de mejor estimación ("best estimate") que no incluyen aproximaciones conservadoras, son de una importancia primordial para analizar estacionarios y transitorios en reactores rápidos. Esta tesis se centra en el desarrollo de un código acoplado para realizar análisis realistas en reactores rápidos críticos aplicando el método de Monte Carlo. Hoy en día, dado el mayor potencial de recursos computacionales, los códigos de transporte neutrónico por Monte Carlo se pueden usar de manera práctica para realizar cálculos detallados de núcleos completos, incluso de elevada heterogeneidad material. Además, los códigos de Monte Carlo se toman normalmente como referencia para los códigos deterministas de difusión en multigrupos en aplicaciones con reactores rápidos, porque usan secciones eficaces punto a punto, un modelo geométrico exacto y tienen en cuenta intrínsecamente la dependencia angular de flujo. En esta tesis se presenta una metodología de acoplamiento entre el conocido código MCNP, que calcula la generación de potencia en el reactor, y el código de termohidráulica de subcanal COBRA-IV, que obtiene las distribuciones de temperatura y densidad en el sistema. COBRA-IV es un código apropiado para aplicaciones en reactores rápidos ya que ha sido validado con resultados experimentales en haces de barras con sodio, incluyendo las correlaciones más apropiadas para metales líquidos. En una primera fase de la tesis, ambos códigos se han acoplado en estado estacionario utilizando un método iterativo con intercambio de archivos externos. El principal problema en el acoplamiento neutrónico y termohidráulico en estacionario con códigos de Monte Carlo es la manipulación de las secciones eficaces para tener en cuenta el ensanchamiento Doppler cuando la temperatura del combustible aumenta. Entre todas las opciones disponibles, en esta tesis se ha escogido la aproximación de pseudo materiales, y se ha comprobado que proporciona resultados aceptables en su aplicación con reactores rápidos. Por otro lado, los cambios geométricos originados por grandes gradientes de temperatura en el núcleo de reactores rápidos resultan importantes para la neutrónica como consecuencia del elevado recorrido libre medio del neutrón en estos sistemas. Por tanto, se ha desarrollado un módulo adicional que simula la geometría del reactor en caliente y permite estimar la reactividad debido a la expansión del núcleo en un transitorio. éste módulo calcula automáticamente la longitud del combustible, el radio de la vaina, la separación de los elementos de combustible y el radio de la placa soporte en función de la temperatura. éste efecto es muy relevante en transitorios sin inserción de bancos de parada. También relacionado con los cambios geométricos, se ha implementado una herramienta que, automatiza el movimiento de las barras de control en busca d la criticidad del reactor, o bien calcula el valor de inserción axial las barras de control. Una segunda fase en la plataforma de cálculo que se ha desarrollado es la simulació dinámica. Puesto que MCNP sólo realiza cálculos estacionarios para sistemas críticos o supercríticos, la solución más directa que se propone sin modificar el código fuente de MCNP es usar la aproximación de factorización de flujo, que resuelve por separado la forma del flujo y la amplitud. En este caso se han estudiado en profundidad dos aproximaciones: adiabática y quasiestática. El método adiabático usa un esquema de acoplamiento que alterna en el tiempo los cálculos neutrónicos y termohidráulicos. MCNP calcula el modo fundamental de la distribución de neutrones y la reactividad al final de cada paso de tiempo, y COBRA-IV calcula las propiedades térmicas en el punto intermedio de los pasos de tiempo. La evolución de la amplitud de flujo se calcula resolviendo las ecuaciones de cinética puntual. Este método calcula la reactividad estática en cada paso de tiempo que, en general, difiere de la reactividad dinámica que se obtendría con la distribución de flujo exacta y dependiente de tiempo. No obstante, para entornos no excesivamente alejados de la criticidad ambas reactividades son similares y el método conduce a resultados prácticos aceptables. Siguiendo esta línea, se ha desarrollado después un método mejorado para intentar tener en cuenta el efecto de la fuente de neutrones retardados en la evolución de la forma del flujo durante el transitorio. El esquema consiste en realizar un cálculo cuasiestacionario por cada paso de tiempo con MCNP. La simulación cuasiestacionaria se basa EN la aproximación de fuente constante de neutrones retardados, y consiste en dar un determinado peso o importancia a cada ciclo computacial del cálculo de criticidad con MCNP para la estimación del flujo final. Ambos métodos se han verificado tomando como referencia los resultados del código de difusión COBAYA3 frente a un ejercicio común y suficientemente significativo. Finalmente, con objeto de demostrar la posibilidad de uso práctico del código, se ha simulado un transitorio en el concepto de reactor crítico en fase de diseño MYRRHA/FASTEF, de 100 MW de potencia térmica y refrigerado por plomo-bismuto. ABSTRACT Long term sustainable nuclear energy scenarios envisage a fleet of Liquid Metal Fast Reactors (LMFR) for the Pu recycling and minor actinides (MAs) transmutation or combined with some accelerator driven systems (ADS) just for MAs elimination. Design and licensing of these innovative reactor concepts require accurate computational tools, implementing the knowledge obtained in experimental research for new reactor configurations, materials and associated systems. Although a number of fast reactor systems have already been built, the operational experience is still reduced, especially for lead reactors, and not all the transients are fully understood. The safety analysis approach for LMFR is therefore based only on deterministic methods, different from modern approach for Light Water Reactors (LWR) which also benefit from probabilistic methods. Usually, the approach adopted in LMFR safety assessments is to employ a variety of codes, somewhat different for the each other, to analyze transients looking for a comprehensive solution and including uncertainties. In this frame, new best estimate simulation codes are of prime importance in order to analyze fast reactors steady state and transients. This thesis is focused on the development of a coupled code system for best estimate analysis in fast critical reactor. Currently due to the increase in the computational resources, Monte Carlo methods for neutrons transport can be used for detailed full core calculations. Furthermore, Monte Carlo codes are usually taken as reference for deterministic diffusion multigroups codes in fast reactors applications because they employ point-wise cross sections in an exact geometry model and intrinsically account for directional dependence of the ux. The coupling methodology presented here uses MCNP to calculate the power deposition within the reactor. The subchannel code COBRA-IV calculates the temperature and density distribution within the reactor. COBRA-IV is suitable for fast reactors applications because it has been validated against experimental results in sodium rod bundles. The proper correlations for liquid metal applications have been added to the thermal-hydraulics program. Both codes are coupled at steady state using an iterative method and external files exchange. The main issue in the Monte Carlo/thermal-hydraulics steady state coupling is the cross section handling to take into account Doppler broadening when temperature rises. Among every available options, the pseudo materials approach has been chosen in this thesis. This approach obtains reasonable results in fast reactor applications. Furthermore, geometrical changes caused by large temperature gradients in the core, are of major importance in fast reactor due to the large neutron mean free path. An additional module has therefore been included in order to simulate the reactor geometry in hot state or to estimate the reactivity due to core expansion in a transient. The module automatically calculates the fuel length, cladding radius, fuel assembly pitch and diagrid radius with the temperature. This effect will be crucial in some unprotected transients. Also related to geometrical changes, an automatic control rod movement feature has been implemented in order to achieve a just critical reactor or to calculate control rod worth. A step forward in the coupling platform is the dynamic simulation. Since MCNP performs only steady state calculations for critical systems, the more straight forward option without modifying MCNP source code, is to use the flux factorization approach solving separately the flux shape and amplitude. In this thesis two options have been studied to tackle time dependent neutronic simulations using a Monte Carlo code: adiabatic and quasistatic methods. The adiabatic methods uses a staggered time coupling scheme for the time advance of neutronics and the thermal-hydraulics calculations. MCNP computes the fundamental mode of the neutron flux distribution and the reactivity at the end of each time step and COBRA-IV the thermal properties at half of the the time steps. To calculate the flux amplitude evolution a solver of the point kinetics equations is used. This method calculates the static reactivity in each time step that in general is different from the dynamic reactivity calculated with the exact flux distribution. Nevertheless, for close to critical situations, both reactivities are similar and the method leads to acceptable practical results. In this line, an improved method as an attempt to take into account the effect of delayed neutron source in the transient flux shape evolutions is developed. The scheme performs a quasistationary calculation per time step with MCNP. This quasistationary simulations is based con the constant delayed source approach, taking into account the importance of each criticality cycle in the final flux estimation. Both adiabatic and quasistatic methods have been verified against the diffusion code COBAYA3, using a theoretical kinetic exercise. Finally, a transient in a critical 100 MWth lead-bismuth-eutectic reactor concept is analyzed using the adiabatic method as an application example in a real system.

Unsaturated polyester-poly(epsilon-caprolactone) hybrid nanocomposites: Thermal-mechanical properties

This paper reports on the thermal behavior and mechanical properties of nanocomposites based on unsaturated polyester resin (UP) modified with poly(ɛ-caprolactone) (PCL) and reinforced with an organically modified clay (cloisite 30B). To optimize the dispersion of 30B and the mixing of PCL in the UP resin, two different methods were employed to prepare crosslinked UP–PCL-30B hybrid nanocomposites. Besides, two samples of poly(ɛ-caprolactone) of different molecular weight (PCL2: Mn = 2.103g.mol−1 and PCL50: Mn = 5.104g.mol−1) were used at several concentrations (4, 6, 10 wt%). The 30B concentration was 4 wt% in all the nanocomposites. The morphology of the samples was studied by scanning electron microscopy (SEM). The analysis of X-ray patterns reveals that intercalated structures have been found for all ternary nanocomposites, independently of the molecular weight, PCL concentration and the preparation method selected. A slight rise of the glass transition temperature, Tg, is observed in UP/PCL/4%30B ternary nanocomposites regarding to neat UP. The analysis of the tensile properties of the ternary (hybrid) systems indicates that UP/4%PCL2/4%30B nanocomposite improves the tensile strength and elongation at break respect to the neat UP while the Young modulus remains constant

Novel poly(3-hydroxybutyrate) nanocomposites containing WS2 inorganic nanotubes with improved thermal, mechanical and tribological properties

Poly(3-hydroxybutyrate) (PHB) nanocomposites containing environmentally-friendly tungsten disulphide inorganic nanotubes (INTeWS2) have been successfully prepared by a simple solution blending method. The dynamic and isothermal crystallization studies by differential scanning calorimetry (DSC) demonstrated that the INTeWS2 exhibits much more prominent nucleation activity on the crystallization of PHB than specific nucleating agents or other nanoscale fillers. Both crystallization rate and crystallinity significantly increase in the nanocomposites compared to neat PHB. These changes occur without modifying the crystalline structure of PHB in the nanocomposites, as shown by wide-angle X-ray diffraction (WAXS) and infrared/Raman spectroscopy. Other parameters such as the Avrami exponent, the equilibrium melting temperature, global rate constant and the fold surface free energy of PHB chains in the nanocomposites were obtained from the calorimetric data in order to determine the influence of the INTeWS2 filler. The addition of INTeWS2 remarkably influences the energetics and kinetics of nucleation and growth of PHB, reducing the fold surface free energy by up to 20%. Furthermore, these nanocomposites also show an improvement in both tribological and mechanical (hardness and modulus) properties with respect to pure PHB evidenced by friction and nanoindentation tests, which is of important potential interest for industrial and medical applications.

Mechanical behavior of alkali-cement as function of the temperature.

This investigation reports on a comparative study of the mechanical behavior at different temperatures of three different alkali-activated fly ash pastes chemically activated using sodium silicate. A control Portland cement (OPC) was used as a reference. In an attempt to simulate the conditions prevailing in the event of accidental fire, post-thermal mechanical tests were performed to determine the residual strength. It has therefore been established that FA based cements can be fabricated for construction purposes and these materials have great potential for fire resistance applications.

Influence of different organoclays on the curing, morphology, and dynamic mechanical properties of an epoxy adhesive

The thermal, mechanical, and adhesive properties of nanoclay-modified adhesives were investigated. Two organically modified montmorillonites: Cloisite 93A (C93A) and Nanomer I.30E (I.30E) were used as reinforcement of an epoxy adhesive. C93A and I.30E are modified with tertiary and primary alkyl ammonium cations, respectively. The aim was to study the influence of the organoclays on the curing, and on the mechanical and adhesive properties of the nanocomposites. A specific goal was to compare their behavior with that of Cloisite30B/epoxy and Cloisite15A/ epoxy nanocomposites that we have previously studied. Both C30B and C15A are modified with quaternary alkyl ammonium cations. Differential scanning calorimetry results showed that the clays accelerate the curing reaction, an effect that is related to the chemical structure of the ammonium cations. The three Cloisite/nanocomposites showed intercalated clay structures,the interlayer distance was independent of the clay content. The I.30E/epoxy nanocomposites presented exfoliated structure due to the catalytic effect of the organic modifier. Clay-epoxy nanocompo-sites showed lower glass transition temperature (Tg) and higher values of storage modulus than neat epoxy thermoset, with no significant differences between exfoliated or intercalated nanocom-posites. The shear strength of aluminum joints using clay/epoxy adhesives was lower than with the neat epoxy adhesive. The wáter aging was less damaging for joints with I.30E/epoxy adhesive.

Eigenstructure Assignment Based Controllers Applied to Flexible Spacecraft

The objective of this paper is to evaluate the behaviour of a controller designed using a parametric Eigenstructure Assignment method and to evaluate its suitability for use in flexible spacecraft. The challenge of this objective lies in obtaining a suitable controller that is specifically designated to alleviate the deflections and vibrations suffered by external appendages in flexible spacecraft while performing attitude manoeuvres. One of the main problems in these vehicles is the mechanical cross-coupling that exists between the rigid and flexible parts of the spacecraft. Spacecraft with fine attitude pointing requirements need precise control of the mechanical coupling to avoid undesired attitude misalignment. In designing an attitude controller, it is necessary to consider the possible vibration of the solar panels and how it may influence the performance of the rest of the vehicle. The nonlinear mathematical model of a flexible spacecraft is considered a close approximation to the real system. During the process of controller evaluation, the design process has also been taken into account as a factor in assessing the robustness of the system.

Tuning the properties of carbon fiber-reinforced poly (phenylene sulphide) laminates via incorporation of inorganic fullerene-like WS2 nanoparticles

Novel carbon fiber (CF)-reinforced poly(phenylene sulphide) (PPS) laminates incorporating inorganic fullerene-like tungsten disulfide (IF-WS2) nanoparticles were prepared via melt-blending and hot-press processing. The influence of the IF-WS2 on the morphology, thermal, mechanical and tribological properties of PPS/CF composites was investigated. Efficient nanoparticle dispersion within the matrix was attained without using surfactants. A progressive rise in thermal stability was found with increasing IF-WS2 loading, as revealed by thermogravimetric analysis. The addition of low nanoparticle contents retarded the crystallization of the matrix, whereas concentrations equal or higher than 1.0 wt% increased both the crystallization temperature and degree of crystallinity compared to those of PPS/CF. Mechanical tests indicated that with only 1.0 wt% IF-WS2 the flexural modulus and strength of PPS/CF improved by 17 and 14%, respectively, without loss in toughness, ascribed to a synergistic effect between the two fillers. A significant enhancement in the storage modulus and glass transition temperature was also observed. Moreover, the wear rate and coefficient of friction strongly decreased, attributed to the lubricant role of the IF-WS2 combined with their reinforcing effect. These inorganic nanoparticles show great potential to improve the mechanical and tribological properties of conventional thermoplastic/CF composites for structural applications.

Opportunities and challenges in the use of inorganic fullerene-like nanoparticles to produce advanced polymer nanocomposites

Polymer/inorganic nanoparticle nanocomposites have garnered considerable academic and industrial interest over recent decades in the development of advanced materials for a wide range of applications. In this respect, the dispersion of so-called inorganic fullerene-like (IF) nanoparticles, e.g., tungsten disulfide (IF-WS2) or molybdenum disulfide (IF-MoS2), into polymeric matrices is emerging as a new strategy. The surprising properties of these layered metal dichalcogenides such as high impact resistance and superior tribological behavior, attributed to their nanoscale size and hollow quasi-spherical shape, open up a wide variety of opportunities for applications of these inorganic compounds. The present work presents a detailed overview on research in the area of IF-based polymer nanocomposites, with special emphasis on the use of IF-WS2 nanoparticles as environmentally friendly reinforcing fillers. The incorporation of IF particles has been shown to be efficient for improving thermal, mechanical and tribological properties of various thermoplastic polymers, such as polypropylene, nylon-6, poly(phenylene sulfide), poly(ether ether ketone), where nanocomposites were fabricated by simple melt-processing routes without the need for modifiers or surfactants. This new family of nanocomposites exhibits similar or enhanced performance when compared with nanocomposites that incorporate carbon nanotubes, carbon nanofibers or nanoclays, but are substantially more cost-effective, efficient and environmentally satisfactory. Most recently, innovative approaches have been described that exploit synergistic effects to produce new materials with enhanced properties, including the combined use of micro- and nanoparticles such as IF-WS2/nucleating agent or IF-WS2/carbon fiber, as well as dual nanoparticle systems such as SWCNT/IF-WS2 where each nanoparticle has different characteristics. The structure–property relationships of these nanocomposites are discussed and potential applications proposed ranging from medicine to the aerospace, automotive and electronics industries.

Large scale multiphysics modeling of neurites

Situado en el límite entre Ingeniería, Informática y Biología, la mecánica computacional de las neuronas aparece como un nuevo campo interdisciplinar que potencialmente puede ser capaz de abordar problemas clínicos desde una perspectiva diferente. Este campo es multiescala por naturaleza, yendo desde la nanoescala (como, por ejemplo, los dímeros de tubulina) a la macroescala (como, por ejemplo, el tejido cerebral), y tiene como objetivo abordar problemas que son complejos, y algunas veces imposibles, de estudiar con medios experimentales. La modelización computacional ha sido ampliamente empleada en aplicaciones Neurocientíficas tan diversas como el crecimiento neuronal o la propagación de los potenciales de acción compuestos. Sin embargo, en la mayoría de los enfoques de modelización hechos hasta ahora, la interacción entre la célula y el medio/estímulo que la rodea ha sido muy poco explorada. A pesar de la tremenda importancia de esa relación en algunos desafíos médicos—como, por ejemplo, lesiones traumáticas en el cerebro, cáncer, la enfermedad del Alzheimer—un puente que relacione las propiedades electrofisiológicas-químicas y mecánicas desde la escala molecular al nivel celular todavía no existe. Con ese objetivo, esta investigación propone un marco computacional multiescala particularizado para dos escenarios respresentativos: el crecimiento del axón y el acomplamiento electrofisiológicomecánico de las neuritas. En el primer caso, se explora la relación entre los constituyentes moleculares del axón durante su crecimiento y sus propiedades mecánicas resultantes, mientras que en el último, un estímulo mecánico provoca deficiencias funcionales a nivel celular como consecuencia de sus alteraciones electrofisiológicas-químicas. La modelización computacional empleada en este trabajo es el método de las diferencias finitas, y es implementada en un nuevo programa llamado Neurite. Aunque el método de los elementos finitos es también explorado en parte de esta investigación, el método de las diferencias finitas tiene la flexibilidad y versatilidad necesaria para implementar mode los biológicos, así como la simplicidad matemática para extenderlos a simulaciones a gran escala con un coste computacional bajo. Centrándose primero en el efecto de las propiedades electrofisiológicas-químicas sobre las propiedades mecánicas, una versión adaptada de Neurite es desarrollada para simular la polimerización de los microtúbulos en el crecimiento del axón y proporcionar las propiedades mecánicas como función de la ocupación de los microtúbulos. Después de calibrar el modelo de crecimiento del axón frente a resultados experimentales disponibles en la literatura, las características mecánicas pueden ser evaluadas durante la simulación. Las propiedades mecánicas del axón muestran variaciones dramáticas en la punta de éste, donde el cono de crecimiento soporta las señales químicas y mecánicas. Bansándose en el conocimiento ganado con el modelo de diferencias finitas, y con el objetivo de ir de 1D a 3D, este esquema preliminar pero de una naturaleza innovadora allana el camino a futuros estudios con el método de los elementos finitos. Centrándose finalmente en el efecto de las propiedades mecánicas sobre las propiedades electrofisiológicas- químicas, Neurite es empleado para relacionar las cargas mecánicas macroscópicas con las deformaciones y velocidades de deformación a escala microscópica, y simular la propagación de la señal eléctrica en las neuritas bajo carga mecánica. Las simulaciones fueron calibradas con resultados experimentales publicados en la literatura, proporcionando, por tanto, un modelo capaz de predecir las alteraciones de las funciones electrofisiológicas neuronales bajo cargas externas dañinas, y uniendo lesiones mecánicas con las correspondientes deficiencias funcionales. Para abordar simulaciones a gran escala, aunque otras arquitecturas avanzadas basadas en muchos núcleos integrados (MICs) fueron consideradas, los solvers explícito e implícito se implementaron en unidades de procesamiento central (CPU) y unidades de procesamiento gráfico (GPUs). Estudios de escalabilidad fueron llevados acabo para ambas implementaciones mostrando resultados prometedores para casos de simulaciones extremadamente grandes con GPUs. Esta tesis abre la vía para futuros modelos mecánicos con el objetivo de unir las propiedades electrofisiológicas-químicas con las propiedades mecánicas. El objetivo general es mejorar el conocimiento de las comunidades médicas y de bioingeniería sobre la mecánica de las neuronas y las deficiencias funcionales que aparecen de los daños producidos por traumatismos mecánicos, como lesiones traumáticas en el cerebro, o enfermedades neurodegenerativas como la enfermedad del Alzheimer. ABSTRACT Sitting at the interface between Engineering, Computer Science and Biology, Computational Neuron Mechanics appears as a new interdisciplinary field potentially able to tackle clinical problems from a new perspective. This field is multiscale by nature, ranging from the nanoscale (e.g., tubulin dimers) to the macroscale (e.g., brain tissue), and aims at tackling problems that are complex, and sometime impossible, to study through experimental means. Computational modeling has been widely used in different Neuroscience applications as diverse as neuronal growth or compound action potential propagation. However, in the majority of the modeling approaches done in this field to date, the interactions between the cell and its surrounding media/stimulus have been rarely explored. Despite of the tremendous importance of such relationship in several medical challenges—e.g., traumatic brain injury (TBI), cancer, Alzheimer’s disease (AD)—a bridge between electrophysiological-chemical and mechanical properties of neurons from the molecular scale to the cell level is still lacking. To this end, this research proposes a multiscale computational framework particularized for two representative scenarios: axon growth and electrophysiological-mechanical coupling of neurites. In the former case, the relation between the molecular constituents of the axon during its growth and its resulting mechanical properties is explored, whereas in the latter, a mechanical stimulus provokes functional deficits at cell level as a consequence of its electrophysiological-chemical alterations. The computational modeling approach chosen in this work is the finite difference method (FDM), and was implemented in a new program called Neurite. Although the finite element method (FEM) is also explored as part of this research, the FDM provides the necessary flexibility and versatility to implement biological models, as well as the mathematical simplicity to extend them to large scale simulations with a low computational cost. Focusing first on the effect of electrophysiological-chemical properties on the mechanical proper ties, an adaptation of Neurite was developed to simulate microtubule polymerization in axonal growth and provide the axon mechanical properties as a function of microtubule occupancy. After calibrating the axon growth model against experimental results available in the literature, the mechanical characteristics can be tracked during the simulation. The axon mechanical properties show dramatic variations at the tip of the axon, where the growth cone supports the chemical and mechanical signaling. Based on the knowledge gained from the FDM scheme, and in order to go from 1D to 3D, this preliminary yet novel scheme paves the road for future studies with FEM. Focusing then on the effect of mechanical properties on the electrophysiological-chemical properties, Neurite was used to relate macroscopic mechanical loading to microscopic strains and strain rates, and simulate the electrical signal propagation along neurites under mechanical loading. The simulations were calibrated against experimental results published in the literature, thus providing a model able to predict the alteration of neuronal electrophysiological function under external damaging load, and linking mechanical injuries to subsequent acute functional deficits. To undertake large scale simulations, although other state-of-the-art architectures based on many integrated cores (MICs) were considered, the explicit and implicit solvers were implemented for central processing units (CPUs) and graphics processing units (GPUs). Scalability studies were done for both implementations showing promising results for extremely large scale simulations with GPUs. This thesis opens the avenue for future mechanical modeling approaches aimed at linking electrophysiological- chemical properties to mechanical properties. Its overarching goal is to enhance the bioengineering and medical communities knowledge on neuronal mechanics and functional deficits arising from damages produced by direct mechanical insults, such as TBI, or neurodegenerative evolving illness, such as AD.

Hierarchical and hybrid polymer nanocomposites based on carbon nanotubes and inorganic fullerene type nanoparticles

The influence of singlewalled carbon nanotubes (SWCNT) and inorganic fullerenelike tungsten disulfide nanoparticles (IFWS2) on the morphology and thermal, mechanical and electrical performance of multifunctional fibrereinforced polymer composites has been investigated. Significant improvements were observed in stiffness, strength and toughness in poly (ether ether ketone) (PEEK) / (SWCNT) / glass fibre (GF) laminates when a compatibilizer was used for wrapping the CNTs. Hybrid poly(phenylene sulphide) (PPS)/IFWS2/ carbon fibre (CF) reinforced polymer composites showed improved mechanical and tribological properties attributed to a synergetic effect between the IF nanoparticles and CF.

Thermoplastic Polymer Nanocomposites Based on Inorganic Fullerene-like Nanoparticles and Inorganic Nanotubes

Using inorganic fullerene-like (IF) nanoparticles and inorganic nanotubes (INT) in organic-inorganic hybrid composite, materials provide the potential for improving thermal, mechanical, and tribological properties of conventional composites. The processing of such high-performance hybrid thermoplastic polymer nanocomposites is achieved via melt-blending without the aid of any modifier or compatibilizing agent. The incorporation of small quantities (0.1-4 wt.%) of IF/INTs (tungsten disulfide, IF-WS2 or molybdenum disulfide, MoS2) generates notable performance enhancements through reinforcement effects and excellent lubricating ability in comparison with promising carbon nanotubes or other inorganic nanoscale fillers. It was shown that these IF/INT nanocomposites can provide an effective balance between performance, cost effectiveness, and processability, which is of significant importance for extending the practical applications of diverse hierarchical thermoplastic-based composites.

Energy-entropy-momentum time integration methods for coupled smooth dissipative problems

Esta tesis aborda la formulación, análisis e implementación de métodos numéricos de integración temporal para la solución de sistemas disipativos suaves de dimensión finita o infinita de manera que su estructura continua sea conservada. Se entiende por dichos sistemas aquellos que involucran acoplamiento termo-mecánico y/o efectos disipativos internos modelados por variables internas que siguen leyes continuas, de modo que su evolución es considerada suave. La dinámica de estos sistemas está gobernada por las leyes de la termodinámica y simetrías, las cuales constituyen la estructura que se pretende conservar de forma discreta. Para ello, los sistemas disipativos se describen geométricamente mediante estructuras metriplécticas que identifican claramente las partes reversible e irreversible de la evolución del sistema. Así, usando una de estas estructuras conocida por las siglas (en inglés) de GENERIC, la estructura disipativa de los sistemas es identificada del mismo modo que lo es la Hamiltoniana para sistemas conservativos. Con esto, métodos (EEM) con precisión de segundo orden que conservan la energía, producen entropía y conservan los impulsos lineal y angular son formulados mediante el uso del operador derivada discreta introducido para asegurar la conservación de la Hamiltoniana y las simetrías de sistemas conservativos. Siguiendo estas directrices, se formulan dos tipos de métodos EEM basados en el uso de la temperatura o de la entropía como variable de estado termodinámica, lo que presenta importantes implicaciones que se discuten a lo largo de esta tesis. Entre las cuales cabe destacar que las condiciones de contorno de Dirichlet son naturalmente impuestas con la formulación basada en la temperatura. Por último, se validan dichos métodos y se comprueban sus mejores prestaciones en términos de la estabilidad y robustez en comparación con métodos estándar. This dissertation is concerned with the formulation, analysis and implementation of structure-preserving time integration methods for the solution of the initial(-boundary) value problems describing the dynamics of smooth dissipative systems, either finite- or infinite-dimensional ones. Such systems are understood as those involving thermo-mechanical coupling and/or internal dissipative effects modeled by internal state variables considered to be smooth in the sense that their evolutions follow continuos laws. The dynamics of such systems are ruled by the laws of thermodynamics and symmetries which constitutes the structure meant to be preserved in the numerical setting. For that, dissipative systems are geometrically described by metriplectic structures which clearly identify the reversible and irreversible parts of their dynamical evolution. In particular, the framework known by the acronym GENERIC is used to reveal the systems' dissipative structure in the same way as the Hamiltonian is for conserving systems. Given that, energy-preserving, entropy-producing and momentum-preserving (EEM) second-order accurate methods are formulated using the discrete derivative operator that enabled the formulation of Energy-Momentum methods ensuring the preservation of the Hamiltonian and symmetries for conservative systems. Following these guidelines, two kind of EEM methods are formulated in terms of entropy and temperature as a thermodynamical state variable, involving important implications discussed throughout the dissertation. Remarkably, the formulation in temperature becomes central to accommodate Dirichlet boundary conditions. EEM methods are finally validated and proved to exhibit enhanced numerical stability and robustness properties compared to standard ones.

Thermal and mechanical effects of different excitation modes based on low frequency laser modulation in optical hyperthermia

The low frequency modulation of the laser source (menor que30KHz) allows the generation of a pulsed signal that intermittently excites the gold nanorods. The temperature curves obtained for different frequencies and duty cycles of modulation but with equal average power and identical laser parameters, show that the thermal behavior in continuous wave and modulation modes is the same. However, the cell death experiments suggest that the percentage of death is higher in the cases of modulation. This observation allows us to conclude that there are other effects in addition to temperature that contribute to the cellular death. The mechanical effects like sound or pressure waves are expected to be generated from thermal expansion of gold nanorods. In order to study the behavior and magnitude of these processes we have developed a measure device based on ultrasound piezoelectric receivers (25KHz) and a lock-in amplifier that is able to detect the sound waves generated in samples of gold nanorods during laser irradiation providing us a voltage result proportional to the pressure signal. The first results show that the pressure measurements are directly proportional to the concentration of gold nanorods and the laser power, therefore, our present work is focused on determine the real influence of these effects in the cell death process.