Flexible multibody approach in bone strain estimation during physical activity:

Bone strain plays a major role as the activation signal for the bone (re)modeling process, which is vital for keeping bones healthy. Maintaining high bone mineral density reduces the chances of fracture in the event of an accident. Numerous studies have shown that bones can be strengthened with physical exercise. Several hypotheses have asserted that a stronger osteogenic (bone producing) effect results from dynamic exercise than from static exercise. These previous studies are based on short-term empirical research, which provide the motivation for justifying the experimental results with a solid mathematical background. The computer simulation techniques utilized in this work allow for non-invasive bone strain estimation during physical activity at any bone site within the human skeleton. All models presented in the study are threedimensional and actuated by muscle models to replicate the real conditions accurately. The objective of this work is to determine and present loading-induced bone strain values resulting from physical activity. It includes a comparison of strain resulting from four different gym exercises (knee flexion, knee extension, leg press, and squat) and walking, with the results reported for walking and jogging obtained from in-vivo measurements described in the literature. The objective is realized primarily by carrying out flexible multibody dynamics computer simulations. The dissertation combines the knowledge of finite element analysis and multibody simulations with experimental data and information available from medical field literature. Measured subject-specific motion data was coupled with forward dynamics simulation to provide natural skeletal movement. Bone geometries were defined using a reverse engineering approach based on medical imaging techniques. Both computed tomography and magnetic resonance imaging were utilized to explore modeling differences. The predicted tibia bone strains during walking show good agreement with invivo studies found in the literature. Strain measurements were not available for gym exercises; therefore, the strain results could not be validated. However, the values seem reasonable when compared to available walking and running invivo strain measurements. The results can be used for exercise equipment design aimed at strengthening the bones as well as the muscles during workout. Clinical applications in post fracture recovery exercising programs could also be the target. In addition, the methodology introduced in this study, can be applied to investigate the effect of weightlessness on astronauts, who often suffer bone loss after long time spent in the outer space.

Multibody Simulation Model of the Roller Test Rig

Multibody simulation model of the roller test rig is presented in this work. The roller test rig consists of a paper machine’s tube roll supported with a hard bearing type balancing machine. The simulation model includes non-idealities that are measured from the physical structure. These non-idealities are the shell thickness variation of the roll and roundness errors of the shafts of the roll. These kinds of non-idealities are harmful since they can cause subharmonic resonances of the rotor system. In this case, the natural vibration mode of the rotor is excited when the rotation speed is a fraction of the natural frequency of the system. With the simulation model, the half critical resonance is studied in detail and a sensitivity analysis is performed by simulating several analyses with slightly different input parameters. The model is verified by comparing the simulation results with those obtained by measuring the real structure. Comparison shows that good accuracy is achieved, since equivalent responses are achieved within the error limit of the input parameters.

Fluid dynamic modelling and simulation of circulating fluidized bed reactors: importance of the interface turbulence transfer

The main objective of this work is to analyze the importance of the gas-solid interface transfer of the kinetic energy of the turbulent motion on the accuracy of prediction of the fluid dynamic of Circulating Fluidized Bed (CFB) reactors. CFB reactors are used in a variety of industrial applications related to combustion, incineration and catalytic cracking. In this work a two-dimensional fluid dynamic model for gas-particle flow has been used to compute the porosity, the pressure, and the velocity fields of both phases in 2-D axisymmetrical cylindrical co-ordinates. The fluid dynamic model is based on the two fluid model approach in which both phases are considered to be continuous and fully interpenetrating. CFB processes are essentially turbulent. The model of effective stress on each phase is that of a Newtonian fluid, where the effective gas viscosity was calculated from the standard k-epsilon turbulence model and the transport coefficients of the particulate phase were calculated from the kinetic theory of granular flow (KTGF). This work shows that the turbulence transfer between the phases is very important for a better representation of the fluid dynamics of CFB reactors, especially for systems with internal recirculation and high gradients of particle concentration. Two systems with different characteristics were analyzed. The results were compared with experimental data available in the literature. The results were obtained by using a computer code developed by the authors. The finite volume method with collocated grid, the hybrid interpolation scheme, the false time step strategy and SIMPLEC (Semi-Implicit Method for Pressure Linked Equations - Consistent) algorithm were used to obtain the numerical solution.

Gas-solid two-phase flow in the riser of circulating fluidized beds: mathematical modelling and numerical simulation

A mathematical model is developed for gas-solids flows in circulating fluidized beds. An Eulerian formulation is followed based on the two-fluids model approach where both the fluid and the particulate phases are treated as a continuum. The physical modelling is discussed, including the formulation of boundary conditions and the description of the numerical methodology. Results of numerical simulation are presented and discussed. The model is validated through comparison to experiment, and simulation is performed to investigate the effects on the flow hydrodynamics of the solids viscosity.

Numerical Simulation of the Cavitation in the Hydrodynamic Lubrication of Journal Bearings: a Parallel Algorithm

In this paper we present an algorithm for the numerical simulation of the cavitation in the hydrodynamic lubrication of journal bearings. Despite the fact that this physical process is usually modelled as a free boundary problem, we adopted the equivalent variational inequality formulation. We propose a two-level iterative algorithm, where the outer iteration is associated to the penalty method, used to transform the variational inequality into a variational equation, and the inner iteration is associated to the conjugate gradient method, used to solve the linear system generated by applying the finite element method to the variational equation. This inner part was implemented using the element by element strategy, which is easily parallelized. We analyse the behavior of two physical parameters and discuss some numerical results. Also, we analyse some results related to the performance of a parallel implementation of the algorithm.

Numerical Simulation of Stochastic Di erential Equations

Stochastic differential equation (SDE) is a differential equation in which some of the terms and its solution are stochastic processes. SDEs play a central role in modeling physical systems like finance, Biology, Engineering, to mention some. In modeling process, the computation of the trajectories (sample paths) of solutions to SDEs is very important. However, the exact solution to a SDE is generally difficult to obtain due to non-differentiability character of realizations of the Brownian motion. There exist approximation methods of solutions of SDE. The solutions will be continuous stochastic processes that represent diffusive dynamics, a common modeling assumption for financial, Biology, physical, environmental systems. This Masters' thesis is an introduction and survey of numerical solution methods for stochastic differential equations. Standard numerical methods, local linearization methods and filtering methods are well described. We compute the root mean square errors for each method from which we propose a better numerical scheme. Stochastic differential equations can be formulated from a given ordinary differential equations. In this thesis, we describe two kind of formulations: parametric and non-parametric techniques. The formulation is based on epidemiological SEIR model. This methods have a tendency of increasing parameters in the constructed SDEs, hence, it requires more data. We compare the two techniques numerically.

Transient current technique simulation of semiconductor detectors

Nowadays advanced simulation technologies of semiconductor devices occupies an important place in microelectronics production process. Simulation helps to understand devices internal processes physics, detect new effects and find directions for optimization. Computer calculation reduces manufacturing costs and time. Modern simulation suits such as Silcaco TCAD allow simulating not only individual semiconductor structures, but also these structures in the circuit. For that purpose TCAD include MixedMode tool. That tool can simulate circuits using compact circuit models including semiconductor structures with their physical models. In this work, MixedMode is used for simulating transient current technique setup, which include detector and supporting electrical circuit. This technique was developed by RD39 collaboration project for investigation radiation detectors radiation hard properties.

Hydration kinetics, physicochemical composition, and textural changes of transgenic corn kernels of flint, semi-flint, and dent varieties

The hydration kinetics of transgenic corn types flint DKB 245PRO, semi-flint DKB 390PRO, and dent DKB 240PRO was studied at temperatures of 30, 40, 50, and 67 °C. The concentrated parameters model was used, and it fits the experimental data well for all three cultivars. The chemical composition of the corn kernels was also evaluated. The corn cultivar influenced the initial rate of absorption and the water equilibrium concentration, and the dent corn absorbed more water than the other cultivars at the four temperatures analyzed. The effect of hydration on the kernel texture was also studied, and it was observed that there was no significant difference in the deformation force required for all three corn types analyzed with longer hydration period.

Experimentation on Wind Powered Hydraulic Heating System

It is common knowledge of the world’s dependency on fossil fuel for energy, its unsustainability on the long run and the changing trend towards renewable energy as an alternative energy source. This aims to cut down greenhouse gas emission and its impact on the rate of ecological and climatic change. Quite remarkably, wind energy has been one of many focus areas of renewable energy sources and has attracted lots of investment and technological advancement. The objective of this research is to explore wind energy and its application in household heating. This research aims at applying experimental approach in real time to study and verify a virtually simulated wind powered hydraulic house heating system. The hardware components comprise of an integrated hydraulic pump, flow control valve, hydraulic fluid and other hydraulic components. The system design and control applies hardware in-the-loop (HIL) simulation setup. Output signal from the semi-empirical turbine modelling controls the integrated motor to generate flow. Throttling the volume flow creates pressure drop across the valve and subsequently thermal power in the system to be outputted using a heat exchanger. Maximum thermal power is achieved by regulating valve orifice to achieve optimum system parameter. Savonius rotor is preferred for its low inertia, high starting torque and ease of design and maintenance characteristics, but lags in power efficiency. A prototype turbine design is used; with power output in range of practical Savonius turbine. The physical mechanism of the prototype turbine’s augmentation design is not known and will not be a focus in this study.

Real-time simulation of multibody systems with applications for working mobile vehicles

This dissertation describes an approach for developing a real-time simulation for working mobile vehicles based on multibody modeling. The use of multibody modeling allows comprehensive description of the constrained motion of the mechanical systems involved and permits real-time solving of the equations of motion. By carefully selecting the multibody formulation method to be used, it is possible to increase the accuracy of the multibody model while at the same time solving equations of motion in real-time. In this study, a multibody procedure based on semi-recursive and augmented Lagrangian methods for real-time dynamic simulation application is studied in detail. In the semirecursive approach, a velocity transformation matrix is introduced to describe the dependent coordinates into relative (joint) coordinates, which reduces the size of the generalized coordinates. The augmented Lagrangian method is based on usage of global coordinates and, in that method, constraints are accounted using an iterative process. A multibody system can be modelled as either rigid or flexible bodies. When using flexible bodies, the system can be described using a floating frame of reference formulation. In this method, the deformation mode needed can be obtained from the finite element model. As the finite element model typically involves large number of degrees of freedom, reduced number of deformation modes can be obtained by employing model order reduction method such as Guyan reduction, Craig-Bampton method and Krylov subspace as shown in this study The constrained motion of the working mobile vehicles is actuated by the force from the hydraulic actuator. In this study, the hydraulic system is modeled using lumped fluid theory, in which the hydraulic circuit is divided into volumes. In this approach, the pressure wave propagation in the hoses and pipes is neglected. The contact modeling is divided into two stages: contact detection and contact response. Contact detection determines when and where the contact occurs, and contact response provides the force acting at the collision point. The friction between tire and ground is modelled using the LuGre friction model, which describes the frictional force between two surfaces. Typically, the equations of motion are solved in the full matrices format, where the sparsity of the matrices is not considered. Increasing the number of bodies and constraint equations leads to the system matrices becoming large and sparse in structure. To increase the computational efficiency, a technique for solution of sparse matrices is proposed in this dissertation and its implementation demonstrated. To assess the computing efficiency, augmented Lagrangian and semi-recursive methods are implemented employing a sparse matrix technique. From the numerical example, the results show that the proposed approach is applicable and produced appropriate results within the real-time period.

Exploring the experiences of youth in a development program with integrated physical activity in Costa Rica

National governments, the United Nations, and other organizations have deemed sport and other means of physical activity such as recreation, games and play for development a useful means for addressing a wide range of problems in communities and more specifically, providing youth with an opportunity to experience the benefits of physical activity. There is a need for research that furthers our understanding of how participants experience these programs. Specifically, the purpose of this study, was to better understand the lived experiences of the participants in a YMCA camp program that integrated physical activity and play for the specific development of poor youth street workers. A phenomenological approach infonned by a critical perspective (Creswell, 2003; Rossman & Rallis, 2003) was used. The study took place through the Asociaci6n Cristiana de J6venes de Costa Rica (ACJ) in Central America. The focus was on a camp program and the lived experiences of six purposefully chosen, youth street workers between the ages of 13-17. Their experiences were explored through semi-structured interviews. Other data that fonn the study include: field notes, observations, a reflexive journal and document analysis. The findings that emerged from the data include main themes of relationships, poverty, personal change and empowennent. For many youth, the ACJ is a relatively safe place to play, to "detach," their minds, to "distract" and "disorient" themselves from their dysfunctional families, violent neighbourhood, the poverty they live in, and from the necessity of having to work in the street to supplement the family income. Although many studies have shown that programs that include physical activity, play and/or sport have a positive impact on youth with regard to healthy development and improvements in well-being, there has been little work done to address the voices and experiences of the youth that participate in these programs. Using an interpretive-critical approach, this study focused on the participants' personal backgrounds, their experiences within the program and their critical reflections on the program. This study draws from a phenomenological philosophy and method to report findings from participants in an ACJ program in Costa Rica. This research shows how these youth were given the opportunity to use the program and the ACJ property as a relatively safe place to play, to behave like the youth they are, to establish and maintain their friendship networks, and develop empathy and conflict resolution skills. The fmdings from this study reveal how by participating in the ACJ program they each described a personal change, wherein they felt empowered to learn they could positivel y control themselves and as a result positively affect their own futures. These fmdings contribute knowledge surrounding the lived experiences of youth in developmental programs that use physical activity.

Overcoming victim behaviour through physical education

This study presents infonnation gathered during personal interviews in which participants were asked how they, as physical educators, might possibly assist the victim of bullying through their programs. The research is a qualitative study, using an inductive approach. Five participants were chosen, based on convenience sampling, with semi-structured interviews which were audio recorded. The theoretical research found that the most stable characteristics of victimized boys were lack of strength and lack of physical fitness. This suggested that Physical Education class might be the best place in which to empower victimized students to reduce their own victimization by addressing these areas of strength and fitness. From the interviews it became clear that, while these educators showed a willingness to help bullying victims through their programs, their adherence to a Physical Education model based primarily on elitism, as opposed to individual fitness, would make it difficult for them to do so effectively.

Understanding the Integration of Living Skills Through the Context of Health and Physical Education: A Case Study of Educators’ Experiences

The research presented is a qualitative case study of educators’ experiences in integrating living skills in the context of health and physical education (HPE). In using semi-structured interviews the study investigated HPE educators’ experiences and revealed their insights relative to three major themes; professional practice, challenges and support systems. Professional practice experiences detailed the use of progressive lesson planning, reflective and engaging activities, explicit student centered pedagogy as well as holistic teaching philosophies. Even further, the limited knowledge and awareness of living skills, conflicting teaching philosophies, competitive environments between subject areas and lack of time and accessibility were four major challenges that emerged throughout the data. Major supportive roles for HPE educators in the integration process included other educators, consultants, school administration, public health, parents, community programs and professional organizations. The study provides valuable discussion and suggestions for improvement of pedagogical practices in teaching living skills in the HPE setting.

Etudes ab initio des effets de la température sur le spectre optique des semi-conducteurs

Thèse réalisée en cotutelle avec l'Université Catholique de Louvain (Belgique)

Implémentation de la répartition de charge et du mode TOT pour la simulation d’un détecteur Timepix à pixels

Les détecteurs à pixels Medipix ont été développés par la collaboration Medipix et permettent de faire de l'imagerie en temps réel. Leur surface active de près de $2\cm^2$ est divisée en 65536~pixels de $55\times 55\um^2$ chacun. Seize de ces détecteurs, les Medipix2, sont installés dans l'expérience ATLAS au CERN afin de mesurer en temps réel les champs de radiation produits par les collisions de hadrons au LHC. Ils seront prochainement remplacés par des Timepix, la plus récente version de ces détecteurs, qui permettent de mesurer directement l'énergie déposée dans chaque pixel en mode \textit{time-over-threshold} (TOT) lors du passage d'une particule dans le semi-conducteur. En vue d'améliorer l'analyse des données recueillies avec ces détecteurs Timepix dans ATLAS, un projet de simulation Geant4 a été amorcé par John Id\'{a}rraga à l'Université de Montréal. Dans le cadre de l'expérience ATLAS, cette simulation pourra être utilisée conjointement avec Athena, le programme d'analyse d'ATLAS, et la simulation complète du détecteur ATLAS. Sous l'effet de leur propre répulsion, les porteurs de charge créés dans le semi-conducteur sont diffusés vers les pixels adjacents causant un dépôt d'énergie dans plusieurs pixels sous l'effet du partage de charges. Un modèle effectif de cette diffusion latérale a été développé pour reproduire ce phénomène sans résoudre d'équation différentielle de transport de charge. Ce modèle, ainsi que le mode TOT du Timepix, qui permet de mesurer l'énergie déposée dans le détecteur, ont été inclus dans la simulation afin de reproduire adéquatement les traces laissées par les particules dans le semi-conducteur. On a d'abord étalonné le détecteur pixel par pixel à l'aide d'une source de $\Am$ et de $\Ba$. Ensuite, on a validé la simulation à l'aide de mesures d'interactions de protons et de particules $\alpha$ produits au générateur Tandem van de Graaff du Laboratoire René-J.-A.-Lévesque de l'Université de Montréal.