Electrolyte disorders in the human ventricle. A simulation study.

Ventricular cells are immersed in a bath of electrolytes and these ions are essential for a healthy heart and a regular rhythm. Maintaining physiological concentration of them is fundamental for reducing arrhythmias and risk of sudden cardiac death, especially in haemodialysis patients and in the heart diseases treatments. Models of electrically activity of the heart based on mathematical formulation are a part of the efforts to improve the understanding and prediction of heart behaviour. Modern models incorporate the extensive and ever increasing amounts of experimental data in incorporating biophysically detailed mechanisms to allow the detailed study of molecular and subcellular mechanisms of heart disease. The goal of this project was to simulate the effects of changes in potassium and calcium concentrations in the extracellular space between experimental data and and a description incorpored into two modern biophysically detailed models (Grandi et al. Model; O’Hara Rudy Model). Moreover the task was to analyze the changes in the ventricular electrical activity, in particular by studying the modifications on the simulated electrocardiographic signal. We used the cellular information obtained by the heart models in order to build a 1D tissue description. The fibre is composed by 165 cells, it is divided in four groups to differentiate the cell types that compound human ventricular tissue. The main results are the following: Grandi et al. (GBP) model is not even able to reproduce the correct action potential profile in hyperkalemia. Data from hospitalized patients indicates that the action potential duration (APD) should be shorter than physiological state but in this model we have the opposite. From the potassium point of view the results obtained by using O’Hara model (ORD) are in agreement with experimental data for the single cell action potential in hypokalemia and hyperkalemia, most of the currents follow the data from literature. In the 1D simulations we were able to reproduce ECGs signal in most the potassium concentrations we selected for this study and we collected data that can help physician in understanding what happens in ventricular cells during electrolyte disorder. However the model fails in the conduction of the stimulus under hyperkalemic conditions. The model emphasized the ECG modifications when the K+ is slightly more than physiological value. In the calcium setting using the ORD model we found an APD shortening in hypocalcaemia and an APD lengthening in hypercalcaemia, i.e. the opposite to experimental observation. This wrong behaviour is kept in one dimensional simulations bringing a longer QT interval in the ECG under higher [Ca2+]o conditions and vice versa. In conclusion it has highlighted that the actual ventricular models present in literature, even if they are useful in the original form, they need an improvement in the sensitivity of these two important electrolytes. We suggest an use of the GBP model with modifications introduced by Carro et al. who understood that the failure of this model is related to the Shannon et al. model (a rabbit model) from which the GBP model was built. The ORD model should be modified in the Ca2+ - dependent IcaL and in the influence of the Iks in the action potential for letting it him produce a correct action potential under different calcium concentrations. In the 1D tissue maybe a heterogeneity setting of intra and extracellular conductances for the different cell types should improve a reproduction of the ECG signal.

Crashworthiness and composite materials: development of an experimental test method for the energy absorption determination and implementation of the relative numerical model

In this PhD thesis the crashworthiness topic is studied with the perspective of the development of a small-scale experimental test able to characterize a material in terms of energy absorption. The material properties obtained are then used to validate a nu- merical model of the experimental test itself. Consequently, the numerical model, calibrated on the specific ma- terial, can be extended to more complex structures and used to simulate their energy absorption behavior. The experimental activity started at University of Washington in Seattle, WA (USA) and continued at Second Faculty of Engi- neering, University of Bologna, Forl`ı (Italy), where the numerical model for the simulation of the experimental test was implemented and optimized.

Use of VEGFR-2 targeted microbubbles (BR55, Bracco imaging) for the early ultrasound evaluation of response to antiangiogenic treatment in a xenograft model of hepatocarcinoma

Aim: To evaluate the early response to treatment to an antiangiogenetic drug (sorafenib) in a heterotopic murine model of hepatocellular carcinoma (HCC) using ultrasonographic molecular imaging. Material and Methods: the xenographt model was established injecting a suspension of HuH7 cells subcutaneously in 19 nude mice. When tumors reached a mean diameter of 5-10 mm, they were divided in two groups (treatment and vehicle). The treatment group received sorafenib (62 mg/kg) by daily oral gavage for 14 days. Molecular imaging was performed using contrast enhanced ultrasound (CEUS), by injecting into the mouse venous circulation a suspension of VEGFR-2 targeted microbubbles (BR55, kind gift of Bracco Swiss, Geneve, Switzerland). Video clips were acquired for 6 minutes, then microbubbles (MBs) were destroyed by a high mechanical index (MI) impulse, and another minute was recorded to evaluate residual circulating MBs. The US protocol was repeated at day 0,+2,+4,+7, and +14 from the beginning of treatment administration. Video clips were analyzed using a dedicated software (Sonotumor, Bracco Swiss) to quantify the signal of the contrast agent. Time/intensity curves were obtained and the difference of the mean MBs signal before and after high MI impulse (Differential Targeted Enhancement-dTE) was calculated. dTE represents a numeric value in arbitrary units proportional to the amount of bound MBs. At day +14 mice were euthanized and the tumors analyzed for VEGFR-2, pERK, and CD31 tissue levels using western blot analysis. Results: dTE values decreased from day 0 to day +14 both in treatment and vehicle groups, and they were statistically higher in vehicle group than in treatment group at day +2, at day +7, and at day +14. With respect to the degree of tumor volume increase, measured as growth percentage delta (GPD), treatment group was divided in two sub-groups, non-responders (GPD>350%), and responders (GPD<200%). In the same way vehicle group was divided in slow growth group (GPD<400%), and fast growth group (GPD>900%). dTE values at day 0 (immediately before treatment start) were higher in non-responders than in responders group, with statistical difference at day 2. While dTE values were higher in the fast growth group than in the slow growth group only at day 0. A significant positive correlation was found between VEGFR-2 tissue levels and dTE values, confirming that level of BR55 tissue enhancement reflects the amount of tissue VEGF receptor. Conclusions: the present findings show that, at least in murine experimental models, CEUS with BR55 is feasable and appears to be a useful tool in the prediction of tumor growth and response to sorafenib treatment in xenograft HCC.

Biodegradazione aerobica cometabolica di cloroformio: studio di fattibilità in reattori a biomassa adesa

Chlorinated Aliphatic Hydrocarbons (CAHs) are widespread wastewater and groundwater contaminants and represent a real danger for human health and environment. This research is related to the biodegradation technologies to treat chlorinated hydrocarbons. In particular the study of this thesis is focused on chloroform cometabolism by a butane-grown aerobic pure culture (Rhodococcus aetherovorans BCP1) in continuous-flow biofilm reactors, which are used for in-situ and on-site treatments. The work was divided in two parts: in the first one an experimental study has been conducted in two packed-bed reactors (PBRs) for a period of 370 days; in the second one a fluid dynamics and kinetic model has been developed in order to simulate the experimental data concerning a previous study made in a 2-m continuous-flow sand-filled reactor. The goals of the first study were to obtain preliminary information on the feasibility of chloroform biodegradation by BCP1 under attached-cell conditions and to evaluate the applicability of the pulsed injection of growth substrate and oxygen to biofilm reactors. The pulsed feeding represents a tool to control the clogging and to ensure a long bioreactive zone. The operational conditions implemented in the PBRs allowed the attainment of a 4-fold increase of the ratio of chloroform degraded to substrate consumed, in comparison with the phase of continuous substrate supply. The second study was aimed at identifying guidelines for optimizing the oxygen/substrate supply schedule, developing a reliable model of chloroform cometabolism in porous media. The tested model led to a suitable interpretation of the experimental data as long as the ratio of CF degraded to butane consumed was ≤ 0.27 mgchloroform /mgbutane. A long-term simulation of the best-performing schedule of pulsed oxygen/substrate supply indicated the attainment of a steady state condition characterized by unsatisfactory bioremediation performances, evidencing the need for a further optimization of the pulsed injection technique.

A Multilevel Model with Time Series Components for the Analysis of Tribal Art Prices

In the present work we perform an econometric analysis of the Tribal art market. To this aim, we use a unique and original database that includes information on Tribal art market auctions worldwide from 1998 to 2011. In Literature, art prices are modelled through the hedonic regression model, a classic fixed-effect model. The main drawback of the hedonic approach is the large number of parameters, since, in general, art data include many categorical variables. In this work, we propose a multilevel model for the analysis of Tribal art prices that takes into account the influence of time on artwork prices. In fact, it is natural to assume that time exerts an influence over the price dynamics in various ways. Nevertheless, since the set of objects change at every auction date, we do not have repeated measurements of the same items over time. Hence, the dataset does not constitute a proper panel; rather, it has a two-level structure in that items, level-1 units, are grouped in time points, level-2 units. The main theoretical contribution is the extension of classical multilevel models to cope with the case described above. In particular, we introduce a model with time dependent random effects at the second level. We propose a novel specification of the model, derive the maximum likelihood estimators and implement them through the E-M algorithm. We test the finite sample properties of the estimators and the validity of the own-written R-code by means of a simulation study. Finally, we show that the new model improves considerably the fit of the Tribal art data with respect to both the hedonic regression model and the classic multilevel model.

Parametric Sensitivity Analysis of the Most Recent Computational Models of Rabbit Cardiac Pacemaking

The cellular basis of cardiac pacemaking activity, and specifically the quantitative contributions of particular mechanisms, is still debated. Reliable computational models of sinoatrial nodal (SAN) cells may provide mechanistic insights, but competing models are built from different data sets and with different underlying assumptions. To understand quantitative differences between alternative models, we performed thorough parameter sensitivity analyses of the SAN models of Maltsev & Lakatta (2009) and Severi et al (2012). Model parameters were randomized to generate a population of cell models with different properties, simulations performed with each set of random parameters generated 14 quantitative outputs that characterized cellular activity, and regression methods were used to analyze the population behavior. Clear differences between the two models were observed at every step of the analysis. Specifically: (1) SR Ca2+ pump activity had a greater effect on SAN cell cycle length (CL) in the Maltsev model; (2) conversely, parameters describing the funny current (If) had a greater effect on CL in the Severi model; (3) changes in rapid delayed rectifier conductance (GKr) had opposite effects on action potential amplitude in the two models; (4) within the population, a greater percentage of model cells failed to exhibit action potentials in the Maltsev model (27%) compared with the Severi model (7%), implying greater robustness in the latter; (5) confirming this initial impression, bifurcation analyses indicated that smaller relative changes in GKr or Na+-K+ pump activity led to failed action potentials in the Maltsev model. Overall, the results suggest experimental tests that can distinguish between models and alternative hypotheses, and the analysis offers strategies for developing anti-arrhythmic pharmaceuticals by predicting their effect on the pacemaking activity.

Forecasts on the dark energy anisotropic stress for the esa euclid survey

The cosmological constant Λ seems to be a not satisfactory explanation of the late-time accelerated expansion of the Universe, for which a number of experimental evidences exist; therefore, it has become necessary in the last years to consider alternative models of dark energy, meant as cause of the accelerated expansion. In the study of dark energy models, it is important to understand which quantities can be determined starting from observational data, without assuming any hypothesis on the cosmological model; such quantities have been determined in Amendola, Kunz et al., 2012. In the same paper it has been further shown that it is possible to estabilish a relation between the model-independent parameters and the anisotropic stress η, which can be also expressed as a combination of the functions appearing in the most general Lagrangian for the scalar-tensor theories, the Horndeski Lagrangian. In the present thesis, the Fisher matrix formalism is used to perform a forecast on the constraints that will be possible to make on the anisotropic stress η in the future, starting from the estimated uncertainties for the galaxy clustering and weak lensing measurements which will be performed by the European Space Agency Euclid mission, to be launched in 2020. Further, constraints coming from supernovae-Ia observations are considered. The forecast is performed for two cases in which (a) η is considered as depending from redshift only and (b) η is constant and equal to one, as in the ΛCDM model.

Ruolo dell'elettrocardiogramma standard nella stratificazione prognostica della cardiomiopatia ipertrofica. Studio multicentrico

Lo scopo di questo studio è di valutare il significato prognostico dell'elettrocardiogramma standard in un'ampia casistica di pazienti affetti da cardiomiopatia ipertrofica. In questo studio multicentrico sono stati considerati 841 pazienti con cardiomiopatia ipertrofica (66% uomini, età media 48±17 anni) per un follow-up di 7.1±7.1 anni, per ognuno è stato analizzato il primo elettrocardiogramma disponibile. I risultati hanno dimostrato come fattori indipendentemente correlati a morte cardiaca improvvisa la sincope inspiegata (p 0.004), il sopraslivellamento del tratto ST e/o la presenza di onde T positive giganti (p 0.048), la durata del QRS >= 120 ms (p 0.017). Sono stati costruiti due modelli per predire il rischio di morte improvvisa: il primo basato sui fattori di rischio universalmente riconosciuti (spessore parietale >= 30 mm, tachicardie ventricolari non sostenute all'ECG Holter 24 ore, sincope e storia familiare di morte improvvisa) e il secondo con l'aggiunta delle variabili sopraslivellamento del tratto ST/onde T positive giganti e durata del QRS >= 120 ms. Entrambi i modelli stratificano i pazienti in base al numero dei fattori di rischio, ma il secondo modello risulta avere un valore predittivo maggiore (chi-square da 12 a 22, p 0.002). In conclusione nella cardiomiopatia ipertrofica l'elettrocardiogramma standard risulta avere un valore prognostico e migliora l'attuale modello di stratificazione per il rischio di morte improvvisa.

Carbonate conical mounds of the eastern Anti-Atlas(Morocco) and their relationship with potential analogues on Mars

Longstanding debates concerning the origin of the Kess Kess Emsian carbonate mounds exposed at Hamar Laghdad Ridge (eastern Anti-Atlas, Morocco) centre around the processes that induced precipitation of carbonate mud and the preservation of steep morphologies. Although in the last years an origin related to hydrothermalism seemed to be more likely, to date the Kess Kess are still considered controversial vent deposits. This study combines in updated research review information from previous work and new detailed field observations coupled with new analytical results to define a consistent framework and some new insight of current knowledge about the origin of these mounds. We obtain a complete minero-petrographic and palaeobiological data set and a detailed geochemical characterization of the different lithologies and facies of the Hamar Laghdad stratigraphic succession, including mounds, and we compared the results with the data from Maïder Basin mounds (Anti-Atlas, Morocco). Our data support the hydrothermal model proposed for the genesis and development of the Kess Kess mounds. The mechanisms linked to the mounds formation and growth are discussed in the light of the new finding of fluid-sediment interaction within a scenario driven by late magmatic fluids circulation. Conical mounds and other fluids related morphologies were also reported from Crommelin crater area (Arabia Terra, Mars). These mounds consist in meter-sized conical buildups hosted in the Equatorial Layered Deposits (ELDs) deposed during a regional groundwater fluid upwelling. Geometries and geological conditions that might have controlled the development of such morphologies were discussed. According to our data the morphological and stratigraphical characteristics of Crommelin area mounds are most consistent with a formation by fluids advection. Then we compare terrestrial and Martian data and examine the geological settings of hydrothermal mound occurrences on Earth in order to describe potential target areas for hydrothermal structures on Mars.

Mimicking the one-dimensional marginal distributions of stochastic diffusions

In the large maturity limit, we compute explicitly the Local Volatility surface for Heston, through Dupire’s formula, with Fourier pricing of the respective derivatives of the call price. Than we verify that the prices of European call options produced by the Heston model, concide with those given by the local volatility model where the Local Volatility is computed as said above.

"Influence of lower-limb joint models on subject-specific musculoskeletal model predictions during gait" ( "modelli muscoloscheletrici personalizzati dell'arto inferiore: Analisi dell'effetto della modellazione dei giunti sulla predizione dei carichi agenti sul sistema scheletrico durante il cammino").

The aim of the present thesis was to investigate the influence of lower-limb joint models on musculoskeletal model predictions during gait. We started our analysis by using a baseline model, i.e., the state-of-the-art lower-limb model (spherical joint at the hip and hinge joints at the knee and ankle) created from MRI of a healthy subject in the Medical Technology Laboratory of the Rizzoli Orthopaedic Institute. We varied the models of knee and ankle joints, including: knee- and ankle joints with mean instantaneous axis of rotation, universal joint at the ankle, scaled-generic-derived planar knee, subject-specific planar knee model, subject-specific planar ankle model, spherical knee, spherical ankle. The joint model combinations corresponding to 10 musculoskeletal models were implemented into a typical inverse dynamics problem, including inverse kinematics, inverse dynamics, static optimization and joint reaction analysis algorithms solved using the OpenSim software to calculate joint angles, joint moments, muscle forces and activations, joint reaction forces during 5 walking trials. The predicted muscle activations were qualitatively compared to experimental EMG, to evaluate the accuracy of model predictions. Planar joint at the knee, universal joint at the ankle and spherical joints at the knee and at the ankle produced appreciable variations in model predictions during gait trials. The planar knee joint model reduced the discrepancy between the predicted activation of the Rectus Femoris and the EMG (with respect to the baseline model), and the reduced peak knee reaction force was considered more accurate. The use of the universal joint, with the introduction of the subtalar joint, worsened the muscle activation agreement with the EMG, and increased ankle and knee reaction forces were predicted. The spherical joints, in particular at the knee, worsened the muscle activation agreement with the EMG. A substantial increase of joint reaction forces at all joints was predicted despite of the good agreement in joint kinematics with those of the baseline model. The introduction of the universal joint had a negative effect on the model predictions. The cause of this discrepancy is likely to be found in the definition of the subtalar joint and thus, in the particular subject’s anthropometry, used to create the model and define the joint pose. We concluded that the implementation of complex joint models do not have marked effects on the joint reaction forces during gait. Computed results were similar in magnitude and in pattern to those reported in literature. Nonetheless, the introduction of planar joint model at the knee had positive effect upon the predictions, while the use of spherical joint at the knee and/or at the ankle is absolutely unadvisable, because it predicted unrealistic joint reaction forces.

Palaeo wind system reconstruction of the last glacial period over Europe, using high resolution proxy data and model-data-comparison

The atmosphere is a global influence on the movement of heat and humidity between the continents, and thus significantly affects climate variability. Information about atmospheric circulation are of major importance for the understanding of different climatic conditions. Dust deposits from maar lakes and dry maars from the Eifel Volcanic Field (Germany) are therefore used as proxy data for the reconstruction of past aeolian dynamics.rnrnIn this thesis past two sediment cores from the Eifel region are examined: the core SM3 from Lake Schalkenmehren and the core DE3 from the Dehner dry maar. Both cores contain the tephra of the Laacher See eruption, which is dated to 12,900 before present. Taken together the cores cover the last 60,000 years: SM3 the Holocene and DE3 the marine isotope stages MIS-3 and MIS-2, respectively. The frequencies of glacial dust storm events and their paleo wind direction are detected by high resolution grain size and provenance analysis of the lake sediments. Therefore two different methods are applied: geochemical measurements of the sediment using µXRF-scanning and the particle analysis method RADIUS (rapid particle analysis of digital images by ultra-high-resolution scanning of thin sections).rnIt is shown that single dust layers in the lake sediment are characterized by an increased content of aeolian transported carbonate particles. The limestone-bearing Eifel-North-South zone is the most likely source for the carbonate rich aeolian dust in the lake sediments of the Dehner dry maar. The dry maar is located on the western side of the Eifel-North-South zone. Thus, carbonate rich aeolian sediment is most likely to be transported towards the Dehner dry maar within easterly winds. A methodology is developed which limits the detection to the aeolian transported carbonate particles in the sediment, the RADIUS-carbonate module.rnrnIn summary, during the marine isotope stage MIS-3 the storm frequency and the east wind frequency are both increased in comparison to MIS-2. These results leads to the suggestion that atmospheric circulation was affected by more turbulent conditions during MIS-3 in comparison to the more stable atmospheric circulation during the full glacial conditions of MIS-2.rnThe results of the investigations of the dust records are finally evaluated in relation a study of atmospheric general circulation models for a comprehensive interpretation. Here, AGCM experiments (ECHAM3 and ECHAM4) with different prescribed SST patterns are used to develop a synoptic interpretation of long-persisting east wind conditions and of east wind storm events, which are suggested to lead to an enhanced accumulation of sediment being transported by easterly winds to the proxy site of the Dehner dry maar.rnrnThe basic observations made on the proxy record are also illustrated in the 10 m-wind vectors in the different model experiments under glacial conditions with different prescribed sea surface temperature patterns. Furthermore, the analysis of long-persisting east wind conditions in the AGCM data shows a stronger seasonality under glacial conditions: all the different experiments are characterized by an increase of the relative importance of the LEWIC during spring and summer. The different glacial experiments consistently show a shift from a long-lasting high over the Baltic Sea towards the NW, directly above the Scandinavian Ice Sheet, together with contemporary enhanced westerly circulation over the North Atlantic.rnrnThis thesis is a comprehensive analysis of atmospheric circulation patterns during the last glacial period. It has been possible to reconstruct important elements of the glacial paleo climate in Central Europe. While the proxy data from sediment cores lead to a binary signal of the wind direction changes (east versus west wind), a synoptic interpretation using atmospheric circulation models is successful. This shows a possible distribution of high and low pressure areas and thus the direction and strength of wind fields which have the capacity to transport dust. In conclusion, the combination of numerical models, to enhance understanding of processes in the climate system, with proxy data from the environmental record is the key to a comprehensive approach to paleo climatic reconstruction.rn

Study of a wind-wave numerical model and its integration with ocean and oil-spill numerical models

The ability to represent the transport and fate of an oil slick at the sea surface is a formidable task. By using an accurate numerical representation of oil evolution and movement in seawater, the possibility to asses and reduce the oil-spill pollution risk can be greatly improved. The blowing of the wind on the sea surface generates ocean waves, which give rise to transport of pollutants by wave-induced velocities that are known as Stokes’ Drift velocities. The Stokes’ Drift transport associated to a random gravity wave field is a function of the wave Energy Spectra that statistically fully describe it and that can be provided by a wave numerical model. Therefore, in order to perform an accurate numerical simulation of the oil motion in seawater, a coupling of the oil-spill model with a wave forecasting model is needed. In this Thesis work, the coupling of the MEDSLIK-II oil-spill numerical model with the SWAN wind-wave numerical model has been performed and tested. In order to improve the knowledge of the wind-wave model and its numerical performances, a preliminary sensitivity study to different SWAN model configuration has been carried out. The SWAN model results have been compared with the ISPRA directional buoys located at Venezia, Ancona and Monopoli and the best model settings have been detected. Then, high resolution currents provided by a relocatable model (SURF) have been used to force both the wave and the oil-spill models and its coupling with the SWAN model has been tested. The trajectories of four drifters have been simulated by using JONSWAP parametric spectra or SWAN directional-frequency energy output spectra and results have been compared with the real paths traveled by the drifters.

The behaviour of osteocyte dendritic processes in bone under cyclic load (il comportamento dei processi dendritici degli osteociti nel tessuto osseo sottoposti a carico ciclico)

Bone is continually being removed and replaced through the actions of basic multicellular units (BMU). This constant upkeep is necessary to remove microdamage formed naturally due to fatigue and thus maintain the integrity of the bone. The repair process in bone is targeted, meaning that a BMU travels directly to the site of damage and repairs it. It is still unclear how targeted remodelling is stimulated and directed but it is highly likely that osteocytes play a role. A number of theories have been advanced to explain the microcrack osteocyte interaction but no complete mechanism has been demonstrated. Osteocytes are connected to each other by dendritic processes. The “scissors model" proposed that the rupture of these processes where they cross microcracks signals the degree of damage and the urgency of the necessary repair. In its original form it was proposed that under applied compressive loading, microcrack faces will be pressed together and undergo relative shear movement. If this movement is greater than the width of an osteocyte process, then the process will be cut in a “scissors like" motion, releasing RANKL, a cytokine known to be essential in the formation of osteoclasts from pre-osteoclasts. The main aim of this thesis was to investigate this theoretical model with a specific focus on microscopy and finite element modelling. Previous studies had proved that cyclic stress was necessary for osteocyte process rupture to occur. This was a divergence from the original “scissors model" which had proposed that the cutting of cell material occurred in one single action. The present thesis is the first study to show fatigue failure in cellular processes spanning naturally occurring cracks and it's the first study to estimate the cyclic strain range and relate it to the number of cycles to failure, for any type of cell. Rupture due to shear movement was ruled out as microcrack closing never occurred, as a result of plastic deformation of the bone. Fatigue failure was found to occur due to cyclic tensile stress in the locality of the damage. The strain range necessary for osteocyte process rupture was quantified. It was found that the lower the process strain range the greater the number of cycles to cell process failure. FEM modelling allowed to predict stress in the vicinity of an osteocyte process and to analyse its interaction with the bone surrounding it: simulations revealed evident creep effects in bone during cyclic loading. This thesis confirms and dismisses aspects of the “scissors model". The observations support the model as a viable mechanism of microcrack detection by the osteocyte network, albeit in a slightly modified form where cyclic loading is necessary and the method of rupture is fatigue failure due to cyclic tensile motion. An in depth study was performed focusing on microscopy analysis of naturally occurring cracks in bone and FEM simulation analysis of an osteocyte process spanning a microcrack in bone under cyclic load.

Fast and accurate numerical solutions in some problems of particle and radiation transport: synthetic acceleration for the method of short characteristics, Doppler-broadened scattering kernel, remote sensing of the cryosphere

The aim of this work is to present various aspects of numerical simulation of particle and radiation transport for industrial and environmental protection applications, to enable the analysis of complex physical processes in a fast, reliable, and efficient way. In the first part we deal with speed-up of numerical simulation of neutron transport for nuclear reactor core analysis. The convergence properties of the source iteration scheme of the Method of Characteristics applied to be heterogeneous structured geometries has been enhanced by means of Boundary Projection Acceleration, enabling the study of 2D and 3D geometries with transport theory without spatial homogenization. The computational performances have been verified with the C5G7 2D and 3D benchmarks, showing a sensible reduction of iterations and CPU time. The second part is devoted to the study of temperature-dependent elastic scattering of neutrons for heavy isotopes near to the thermal zone. A numerical computation of the Doppler convolution of the elastic scattering kernel based on the gas model is presented, for a general energy dependent cross section and scattering law in the center of mass system. The range of integration has been optimized employing a numerical cutoff, allowing a faster numerical evaluation of the convolution integral. Legendre moments of the transfer kernel are subsequently obtained by direct quadrature and a numerical analysis of the convergence is presented. In the third part we focus our attention to remote sensing applications of radiative transfer employed to investigate the Earth's cryosphere. The photon transport equation is applied to simulate reflectivity of glaciers varying the age of the layer of snow or ice, its thickness, the presence or not other underlying layers, the degree of dust included in the snow, creating a framework able to decipher spectral signals collected by orbiting detectors.