Development of a continuum mechanics model of passive skeletal muscle

Skeletal muscle force evaluation is difficult to implement in a clinical setting. Muscle force is typically assessed through either manual muscle testing, isokinetic/isometric dynamometry, or electromyography (EMG). Manual muscle testing is a subjective evaluation of a patient’s ability to move voluntarily against gravity and to resist force applied by an examiner. Muscle testing using dynamometers adds accuracy by quantifying functional mechanical output of a limb. However, like manual muscle testing, dynamometry only provides estimates of the joint moment. EMG quantifies neuromuscular activation signals of individual muscles, and is used to infer muscle function. Despite the abundance of work performed to determine the degree to which EMG signals and muscle forces are related, the basic problem remains that EMG cannot provide a quantitative measurement of muscle force. Intramuscular pressure (IMP), the pressure applied by muscle fibers on interstitial fluid, has been considered as a correlate for muscle force. Numerous studies have shown that an approximately linear relationship exists between IMP and muscle force. A microsensor has recently been developed that is accurate, biocompatible, and appropriately sized for clinical use. While muscle force and pressure have been shown to be correlates, IMP has been shown to be non-uniform within the muscle. As it would not be practicable to experimentally evaluate how IMP is distributed, computational modeling may provide the means to fully evaluate IMP generation in muscles of various shapes and operating conditions. The work presented in this dissertation focuses on the development and validation of computational models of passive skeletal muscle and the evaluation of their performance for prediction of IMP. A transversly isotropic, hyperelastic, and nearly incompressible model will be evaluated along with a poroelastic model.

Optimal game theoretic distributed control methodologies in small scale power systems

This dissertation presents the competitive control methodologies for small-scale power system (SSPS). A SSPS is a collection of sources and loads that shares a common network which can be isolated during terrestrial disturbances. Micro-grids, naval ship electric power systems (NSEPS), aircraft power systems and telecommunication system power systems are typical examples of SSPS. The analysis and development of control systems for small-scale power systems (SSPS) lacks a defined slack bus. In addition, a change of a load or source will influence the real time system parameters of the system. Therefore, the control system should provide the required flexibility, to ensure operation as a single aggregated system. In most of the cases of a SSPS the sources and loads must be equipped with power electronic interfaces which can be modeled as a dynamic controllable quantity. The mathematical formulation of the micro-grid is carried out with the help of game theory, optimal control and fundamental theory of electrical power systems. Then the micro-grid can be viewed as a dynamical multi-objective optimization problem with nonlinear objectives and variables. Basically detailed analysis was done with optimal solutions with regards to start up transient modeling, bus selection modeling and level of communication within the micro-grids. In each approach a detail mathematical model is formed to observe the system response. The differential game theoretic approach was also used for modeling and optimization of startup transients. The startup transient controller was implemented with open loop, PI and feedback control methodologies. Then the hardware implementation was carried out to validate the theoretical results. The proposed game theoretic controller shows higher performances over traditional the PI controller during startup. In addition, the optimal transient surface is necessary while implementing the feedback controller for startup transient. Further, the experimental results are in agreement with the theoretical simulation. The bus selection and team communication was modeled with discrete and continuous game theory models. Although players have multiple choices, this controller is capable of choosing the optimum bus. Next the team communication structures are able to optimize the players’ Nash equilibrium point. All mathematical models are based on the local information of the load or source. As a result, these models are the keys to developing accurate distributed controllers.

XML for Model Specification in Neuroscience

One of the main roles of the Neural Open Markup Language, NeuroML, is to facilitate cooperation in building, simulating, testing and publishing models of channels, neurons and networks of neurons. MorphML, which was developed as a common format for exchange of neural morphology data, is distributed as part of NeuroML but can be used as a stand-alone application. In this collection of tutorials and workshop summary, we provide an overview of these XML schemas and provide examples of their use in down-stream applications. We also summarize plans for the further development of XML specifications for modeling channels, channel distributions, and network connectivity.

Lag Camera: A Moving Multi-Camera Array for Scene-Acquisition

Many applications, such as telepresence, virtual reality, and interactive walkthroughs, require a three-dimensional(3D)model of real-world environments. Methods, such as lightfields, geometric reconstruction and computer vision use cameras to acquire visual samples of the environment and construct a model. Unfortunately, obtaining models of real-world locations is a challenging task. In particular, important environments are often actively in use, containing moving objects, such as people entering and leaving the scene. The methods previously listed have difficulty in capturing the color and structure of the environment while in the presence of moving and temporary occluders. We describe a class of cameras called lag cameras. The main concept is to generalize a camera to take samples over space and time. Such a camera, can easily and interactively detect moving objects while continuously moving through the environment. Moreover, since both the lag camera and occluder are moving, the scene behind the occluder is captured by the lag camera even from viewpoints where the occluder lies in between the lag camera and the hidden scene. We demonstrate an implementation of a lag camera, complete with analysis and captured environments.

Multichannel EEG fields during and without visual input: frequency domain model source locations and dimensional complexities

27-Channel EEG potential map series were recorded from 12 normals with closed and open eyes. Intracerebral dipole model source locations in the frequency domain were computed. Eye opening (visual input) caused centralization (convergence and elevation) of the source locations of the seven frequency bands, indicative of generalized activity; especially, there was clear anteriorization of α-2 (10.5–12 Hz) and β-2 (18.5–21 Hz) sources (α-2 also to the left). Complexity of the map series' trajectories in state space (assessed by Global Dimensional Complexity and Global OMEGA Complexity) increased significantly with eye opening, indicative of more independent, parallel, active processes. Contrary to PET and fMRI, these results suggest that brain activity is more distributed and independent during visual input than after eye closing (when it is more localized and more posterior).

SLA-driven predictive orchestration for distributed cloud-based mobile services

We describe a system for performing SLA-driven management and orchestration of distributed infrastructures composed of services supporting mobile computing use cases. In particular, we focus on a Follow-Me Cloud scenario in which we consider mobile users accessing cloud-enable services. We combine a SLA-driven approach to infrastructure optimization, with forecast-based performance degradation preventive actions and pattern detection for supporting mobile cloud infrastructure management. We present our system's information model and architecture including the algorithmic support and the proposed scenarios for system evaluation.

Involvement of heparin-like glycosaminoglycans and expression of a novel heparan sulfate binding protein (p24) during human placentation and in a model for human implantation

Previous studies in our laboratory have indicated that heparan sulfate proteoglycans (HSPGs) play an important role in murine embryo implantation. To investigate the potential function of HSPGs in human implantation, two human cell lines (RL95 and JAR) were selected to model uterine epithelium and embryonal trophectoderm, respectively. A heterologous cell-cell adhesion assay showed that initial binding between JAR and RL95 cells is mediated by cell surface glycosaminoglycans (GAG) with heparin-like properties, i.e., heparan sulfate and dermatan sulfate. Furthermore, a single class of highly specific, protease-sensitive heparin/heparan sulfate binding sites exist on the surface of RL95 cells. Three heparin binding, tryptic peptide fragments were isolated from RL95 cell surfaces and their amino termini partially sequenced. Reverse transcription-polymerase chain reaction (RT-PCR) generated 1 to 4 PCR products per tryptic peptide. Northern blot analysis of RNA from RL95 cells using one of these RT-PCR products identified a 1.2 Kb mRNA species (p24). The amino acid sequence predicted from the cDNA sequence contains a putative heparin-binding domain. A synthetic peptide representing this putative heparin binding domain was used to generate a rabbit polyclonal antibody (anti-p24). Indirect immunofluorescence studies on RL95 and JAR cells as well as binding studies of anti-p24 to intact RL95 cells demonstrate that p24 is distributed on the cell surface. Western blots of RL95 membrane preparations identify a 24 kDa protein (p24) highly enriched in the 100,000 g pellet plasma membrane-enriched fraction. p24 eluted from membranes with 0.8 M NaCl, but not 0.6 M NaCl, suggesting that it is a peripheral membrane component. Solubilized p24 binds heparin by heparin affinity chromatography and $\sp{125}$I-heparin binding assays. Furthermore, indirect immunofluorescence studies indicate that cytotrophoblast of floating and attached villi of the human fetal-maternal interface are recognized by anti-p24. The study also indicates that the HSPG, perlecan, accumulates where chorionic villi are attached to uterine stroma and where p24-expressing cytotrophoblast penetrate the stroma. Collectively, these data indicate that p24 is a cell surface membrane-associated heparin/heparan sulfate binding protein found in cytotrophoblast, but not many other cell types of the fetal-maternal interface. Furthermore, p24 colocalizes with HSPGs in regions of cytotrophoblast invasion. These observations are consistent with a role for HSPGs and HSPG binding proteins in human trophoblast-uterine cell interactions. ^

Determination of respiratory gas flow by electrical impedance tomography in an animal model of mechanical ventilation

Background A recent method determines regional gas flow of the lung by electrical impedance tomography (EIT). The aim of this study is to show the applicability of this method in a porcine model of mechanical ventilation in healthy and diseased lungs. Our primary hypothesis is that global gas flow measured by EIT can be correlated with spirometry. Our secondary hypothesis is that regional analysis of respiratory gas flow delivers physiologically meaningful results. Methods In two sets of experiments n = 7 healthy pigs and n = 6 pigs before and after induction of lavage lung injury were investigated. EIT of the lung and spirometry were registered synchronously during ongoing mechanical ventilation. In-vivo aeration of the lung was analysed in four regions-of-interest (ROI) by EIT: 1) global, 2) ventral (non-dependent), 3) middle and 4) dorsal (dependent) ROI. Respiratory gas flow was calculated by the first derivative of the regional aeration curve. Four phases of the respiratory cycle were discriminated. They delivered peak and late inspiratory and expiratory gas flow (PIF, LIF, PEF, LEF) characterizing early or late inspiration or expiration. Results Linear regression analysis of EIT and spirometry in healthy pigs revealed a very good correlation measuring peak flow and a good correlation detecting late flow. PIFEIT = 0.702 · PIFspiro + 117.4, r2 = 0.809; PEFEIT = 0.690 · PEFspiro-124.2, r2 = 0.760; LIFEIT = 0.909 · LIFspiro + 27.32, r2 = 0.572 and LEFEIT = 0.858 · LEFspiro-10.94, r2 = 0.647. EIT derived absolute gas flow was generally smaller than data from spirometry. Regional gas flow was distributed heterogeneously during different phases of the respiratory cycle. But, the regional distribution of gas flow stayed stable during different ventilator settings. Moderate lung injury changed the regional pattern of gas flow. Conclusions We conclude that the presented method is able to determine global respiratory gas flow of the lung in different phases of the respiratory cycle. Additionally, it delivers meaningful insight into regional pulmonary characteristics, i.e. the regional ability of the lung to take up and to release air.

Service Level Agreements-Driven Management of Distributed Applications in Cloud Computing Environments

Cloud Computing enables provisioning and distribution of highly scalable services in a reliable, on-demand and sustainable manner. However, objectives of managing enterprise distributed applications in cloud environments under Service Level Agreement (SLA) constraints lead to challenges for maintaining optimal resource control. Furthermore, conflicting objectives in management of cloud infrastructure and distributed applications might lead to violations of SLAs and inefficient use of hardware and software resources. This dissertation focusses on how SLAs can be used as an input to the cloud management system, increasing the efficiency of allocating resources, as well as that of infrastructure scaling. First, we present an extended SLA semantic model for modelling complex service-dependencies in distributed applications, and for enabling automated cloud infrastructure management operations. Second, we describe a multi-objective VM allocation algorithm for optimised resource allocation in infrastructure clouds. Third, we describe a method of discovering relations between the performance indicators of services belonging to distributed applications and then using these relations for building scaling rules that a CMS can use for automated management of VMs. Fourth, we introduce two novel VM-scaling algorithms, which optimally scale systems composed of VMs, based on given SLA performance constraints. All presented research works were implemented and tested using enterprise distributed applications.

Evaluation of a Novel Screw Position in a Type III Distal Phalanx Fracture Model: An Ex Vivo Study

OBJECTIVE: Mechanical evaluation of a novel screw position used for repair in a type III distal phalanx fracture model and assessment of solar canal penetration (SCP). STUDY DESIGN: Experimental study. SAMPLE POPULATION: Disarticulated equine hooves (n = 24) and 24 isolated distal phalanges. METHODS: Hooves/distal phalanges cut in a sagittal plane were repaired with 1 of 2 different cortical screw placements in lag fashion. In group 1 (conventional screw placement), the screw was inserted halfway between the proximal border of the solar canal (SC) and the subchondral bone surface on a line parallel to the dorsal cortex, whereas in group 2, the screw was inserted more palmar/plantar, where a perpendicular line drawn from the group 1 position reached the palmar/plantar cortex. Construct strength was evaluated by 3-point bending to failure. SCP was assessed by CT imaging and macroscopically. RESULTS: Screws were significantly longer in group 2 and in forelimbs. Group 2 isolated distal phalanges had a significantly more rigid fixation compared with the conventional screw position (maximum point at failure 31%, bending stiffness 41% higher). Lumen reduction of the SC was observed in 13/52 specimens (all from group 2), of which 9 were forelimbs. CONCLUSIONS: More distal screw positioning compared with the conventionally recommended screw position for internal fixation of type III distal phalangeal fractures allows placement of a longer screw and renders a more rigid fracture fixation. The novel screw position, however, carries a higher risk of SCP

Assessing residents' readiness to screen for domestic violence: Utilizing the transtheoretical model - stages of change, decisional balance and self-efficacy

Domestic violence is a major public health problem, yet most physicians do not effectively identify patients at risk. Medical students and residents are not routinely educated on this topic and little is known about the factors that influence their decisions to include screening for domestic violence in their subsequent practice. In order to assess the readiness of primary care residents to screen all patients for domestic violence, this study utilized a survey incorporating constructs from the Transtheoretical Model, including Stages of Change, Decisional Balance (Pros and Cons) and Self-Efficacy. The survey was distributed to residents at the University of Texas Health Science Center Medical School in Houston in: Internal Medicine, Medicine/Pediatrics, Pediatrics, Family Medicine, and Obstetrics and Gynecology. Data from the survey was analyzed to test the hypothesis that residents in the earlier Stages of Change report more costs and fewer benefits with regards to screening for domestic violence, and that those in the later stages exhibit higher Self-Efficacy scores. The findings from this study were consistent with the model in that benefits to screening (Pros) and Self-Efficacy were correlated with later Stages of Change, however reporting fewer costs (Cons) was not. Very few residents were ready to screen all of their patients.^

The number and timing of time -dependent covariate measurements for the Cox regression model

The problem of analyzing data with updated measurements in the time-dependent proportional hazards model arises frequently in practice. One available option is to reduce the number of intervals (or updated measurements) to be included in the Cox regression model. We empirically investigated the bias of the estimator of the time-dependent covariate while varying the effect of failure rate, sample size, true values of the parameters and the number of intervals. We also evaluated how often a time-dependent covariate needs to be collected and assessed the effect of sample size and failure rate on the power of testing a time-dependent effect.^ A time-dependent proportional hazards model with two binary covariates was considered. The time axis was partitioned into k intervals. The baseline hazard was assumed to be 1 so that the failure times were exponentially distributed in the ith interval. A type II censoring model was adopted to characterize the failure rate. The factors of interest were sample size (500, 1000), type II censoring with failure rates of 0.05, 0.10, and 0.20, and three values for each of the non-time-dependent and time-dependent covariates (1/4,1/2,3/4).^ The mean of the bias of the estimator of the coefficient of the time-dependent covariate decreased as sample size and number of intervals increased whereas the mean of the bias increased as failure rate and true values of the covariates increased. The mean of the bias of the estimator of the coefficient was smallest when all of the updated measurements were used in the model compared with two models that used selected measurements of the time-dependent covariate. For the model that included all the measurements, the coverage rates of the estimator of the coefficient of the time-dependent covariate was in most cases 90% or more except when the failure rate was high (0.20). The power associated with testing a time-dependent effect was highest when all of the measurements of the time-dependent covariate were used. An example from the Systolic Hypertension in the Elderly Program Cooperative Research Group is presented. ^

Model results of the impact of climate change on agriculture in China

This paper assesses the impact of climate change on China's agricultural production at a cross-provincial level using the Ricardian approach, incorporating a multilevel model with farm-level group data. The farm-level group data includes 13379 farm households, across 316 villages, distributed in 31 provinces. The empirical results show that, firstly, the marginal effects and elasticities of net crop revenue per hectare with respect to climate factors indicated that the annual impact of temperature on net crop revenue per hectare was positive, and the effect of increased precipitation was negative when looking at the national totals; secondly, the total impact of simulated climate change scenarios on net crop revenues per hectare at a Chinese national total level, was an increase of between 79 USD per hectare and 207 USD per hectare for the 2050s, and an increase from 140 USD per hectare to 355 USD per hectare for the 2080s. As a result, climate change may create a potential advantage for the development of Chinese agriculture, rather than a risk, especially for agriculture in the provinces of the Northeast, Northwest and North regions. However, the increased precipitation can lead to a loss of net crop revenue per hectare, especially for the provinces of the Southwest, Northwest, North and Northeast regions.

Using ground-penetrating radar, topography and classification of vegetation to model the sediment and active layer thickness in a periglacial lake catchment, Western Greenland, link to shapefiles

The geometries of a catchment constitute the basis for distributed physically based numerical modeling of different geoscientific disciplines. In this paper results from ground-penetrating radar (GPR) measurements, in terms of a 3D model of total sediment thickness and active layer thickness in a periglacial catchment in western Greenland, is presented. Using the topography, thickness and distribution of sediments is calculated. Vegetation classification and GPR measurements are used to scale active layer thickness from local measurements to catchment scale models. Annual maximum active layer thickness varies from 0.3 m in wetlands to 2.0 m in barren areas and areas of exposed bedrock. Maximum sediment thickness is estimated to be 12.3 m in the major valleys of the catchment. A method to correlate surface vegetation with active layer thickness is also presented. By using relatively simple methods, such as probing and vegetation classification, it is possible to upscale local point measurements to catchment scale models, in areas where the upper subsurface is relatively homogenous. The resulting spatial model of active layer thickness can be used in combination with the sediment model as a geometrical input to further studies of subsurface mass-transport and hydrological flow paths in the periglacial catchment through numerical modelling.