Novel cationic polymers for use at biological interfaces

One of the main problems with the use of synthetic polymers as biomaterials is the invasion of micro-organisms causing infection. A study of the properties of polymeric antibacterial agents, in particular polyhexamethylene biguanide, has revealed that the essential components for the design of a novel polymeric antibacterial are a balance between hydrophilicity and hydrophobicity coupled with sites of cationicity. The effect of cation incorporation on the physical properties of hydrogels has been investigated. Hydrogel systems copolymerised with either N-vinyl imidazole or dimethylaminoethyl methacrylate have been characterised in terms of their water binding, mechanical and surface properties. It has been concluded that the incorporation of these monomers does not adversely affect the properties of such hydrogels and that these materials are potential candidates for further development for use in biomedical applications. It has been reported that hydro gels with ionic character may increase the deposition of biological material onto the hydrogel surface when it is in contact with body fluids. An investigation into the deposition characteristics of hydrogels containing the potentially cationic monomers has been carried out, using specific protein adsorption and in vitro spoilation techniques. The results suggest that at low levels of cationicity, the deposition of positively charged proteins is reduced without adversely affecting the uptake of the other proteins. The gross deposition characteristics were found to be comparable to some commercially available contact lens materials. A preliminary investigation into the development of novel antibacterial polymers has been completed and some novel methods of bacterial inhibition discussed. These methods include development of an hydrogel whose potential application is as a catheter coating.

Application of the finite element method to problems in fracture mechanics

Numerical techniques have been finding increasing use in all aspects of fracture mechanics, and often provide the only means for analyzing fracture problems. The work presented here, is concerned with the application of the finite element method to cracked structures. The present work was directed towards the establishment of a comprehensive two-dimensional finite element, linear elastic, fracture analysis package. Significant progress has been made to this end, and features which can now be studied include multi-crack tip mixed-mode problems, involving partial crack closure. The crack tip core element was refined and special local crack tip elements were employed to reduce the element density in the neighbourhood of the core region. The work builds upon experience gained by previous research workers and, as part of the general development, the program was modified to incorporate the eight-node isoparametric quadrilateral element. Also. a more flexible solving routine was developed, and provided a very compact method of solving large sets of simultaneous equations, stored in a segmented form. To complement the finite element analysis programs, an automatic mesh generation program has been developed, which enables complex problems. involving fine element detail, to be investigated with a minimum of input data. The scheme has proven to be versati Ie and reasonably easy to implement. Numerous examples are given to demonstrate the accuracy and flexibility of the finite element technique.

On Finite Element Methods for 2nd order (semi–) periodic Eigenvalue Problems

We deal with a class of elliptic eigenvalue problems (EVPs) on a rectangle Ω ⊂ R^2 , with periodic or semi–periodic boundary conditions (BCs) on ∂Ω. First, for both types of EVPs, we pass to a proper variational formulation which is shown to fit into the general framework of abstract EVPs for symmetric, bounded, strongly coercive bilinear forms in Hilbert spaces, see, e.g., [13, §6.2]. Next, we consider finite element methods (FEMs) without and with numerical quadrature. The aim of the paper is to show that well–known error estimates, established for the finite element approximation of elliptic EVPs with classical BCs, hold for the present types of EVPs too. Some attention is also paid to the computational aspects of the resulting algebraic EVP. Finally, the analysis is illustrated by two non-trivial numerical examples, the exact eigenpairs of which can be determined.

On a Variational Approach to some Quasilinear Problems

We prove some multiplicity results concerning quasilinear elliptic equations with natural growth conditions. Techniques of nonsmooth critical point theory are employed.

Ultrafast charge-transfer in organic photovoltaic interfaces: geometrical and functionalization effects

Understanding the microscopic mechanisms of electronic excitation in organic photovoltaic cells is a challenging problem in the design of efficient devices capable of performing sunlight harvesting. Here we develop and apply an ab initio approach based on time-dependent density functional theory and Ehrenfest dynamics to investigate photoinduced charge transfer in small organic molecules. Our calculations include mixed quantum–classical dynamics with ions moving classically and electrons quantum mechanically, where no experimental external parameter other than the material geometry is required. We show that the behavior of photocarriers in zinc phthalocyanine (ZnPc) and C60 systems, an effective prototype system for organic solar cells, is sensitive to the atomic orientation of the donor and the acceptor units as well as the functionalization of covalent molecules at the interface. In particular, configurations with the ZnPc molecules facing on C60 facilitate charge transfer between substrate and molecules that occurs within 200 fs. In contrast, configurations where ZnPc is tilted above C60 present extremely low carrier injection efficiency even at longer times as an effect of the larger interfacial potential level offset and higher energetic barrier between the donor and acceptor molecules. An enhancement of charge injection into C60 at shorter times is observed as binding groups connect ZnPc and C60 in a dyad system. Our results demonstrate a promising way of designing and controlling photoinduced charge transfer on the atomic level in organic devices that would lead to efficient carrier separation and maximize device performance.

DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.

Motion Planning and Controls with Safety and Temporal Constraints

Motion planning, or trajectory planning, commonly refers to a process of converting high-level task specifications into low-level control commands that can be executed on the system of interest. For different applications, the system will be different. It can be an autonomous vehicle, an Unmanned Aerial Vehicle(UAV), a humanoid robot, or an industrial robotic arm. As human machine interaction is essential in many of these systems, safety is fundamental and crucial. Many of the applications also involve performing a task in an optimal manner within a given time constraint. Therefore, in this thesis, we focus on two aspects of the motion planning problem. One is the verification and synthesis of the safe controls for autonomous ground and air vehicles in collision avoidance scenarios. The other part focuses on the high-level planning for the autonomous vehicles with the timed temporal constraints. In the first aspect of our work, we first propose a verification method to prove the safety and robustness of a path planner and the path following controls based on reachable sets. We demonstrate the method on quadrotor and automobile applications. Secondly, we propose a reachable set based collision avoidance algorithm for UAVs. Instead of the traditional approaches of collision avoidance between trajectories, we propose a collision avoidance scheme based on reachable sets and tubes. We then formulate the problem as a convex optimization problem seeking control set design for the aircraft to avoid collision. We apply our approach to collision avoidance scenarios of quadrotors and fixed-wing aircraft. In the second aspect of our work, we address the high level planning problems with timed temporal logic constraints. Firstly, we present an optimization based method for path planning of a mobile robot subject to timed temporal constraints, in a dynamic environment. Temporal logic (TL) can address very complex task specifications such as safety, coverage, motion sequencing etc. We use metric temporal logic (MTL) to encode the task specifications with timing constraints. We then translate the MTL formulae into mixed integer linear constraints and solve the associated optimization problem using a mixed integer linear program solver. We have applied our approach on several case studies in complex dynamical environments subjected to timed temporal specifications. Secondly, we also present a timed automaton based method for planning under the given timed temporal logic specifications. We use metric interval temporal logic (MITL), a member of the MTL family, to represent the task specification, and provide a constructive way to generate a timed automaton and methods to look for accepting runs on the automaton to find an optimal motion (or path) sequence for the robot to complete the task.

New Rigorous Decomposition Methods for Mixed-integer Linear and Nonlinear Programming

Process systems design, operation and synthesis problems under uncertainty can readily be formulated as two-stage stochastic mixed-integer linear and nonlinear (nonconvex) programming (MILP and MINLP) problems. These problems, with a scenario based formulation, lead to large-scale MILPs/MINLPs that are well structured. The first part of the thesis proposes a new finitely convergent cross decomposition method (CD), where Benders decomposition (BD) and Dantzig-Wolfe decomposition (DWD) are combined in a unified framework to improve the solution of scenario based two-stage stochastic MILPs. This method alternates between DWD iterations and BD iterations, where DWD restricted master problems and BD primal problems yield a sequence of upper bounds, and BD relaxed master problems yield a sequence of lower bounds. A variant of CD, which includes multiple columns per iteration of DW restricted master problem and multiple cuts per iteration of BD relaxed master problem, called multicolumn-multicut CD is then developed to improve solution time. Finally, an extended cross decomposition method (ECD) for solving two-stage stochastic programs with risk constraints is proposed. In this approach, a CD approach at the first level and DWD at a second level is used to solve the original problem to optimality. ECD has a computational advantage over a bilevel decomposition strategy or solving the monolith problem using an MILP solver. The second part of the thesis develops a joint decomposition approach combining Lagrangian decomposition (LD) and generalized Benders decomposition (GBD), to efficiently solve stochastic mixed-integer nonlinear nonconvex programming problems to global optimality, without the need for explicit branch and bound search. In this approach, LD subproblems and GBD subproblems are systematically solved in a single framework. The relaxed master problem obtained from the reformulation of the original problem, is solved only when necessary. A convexification of the relaxed master problem and a domain reduction procedure are integrated into the decomposition framework to improve solution efficiency. Using case studies taken from renewable resource and fossil-fuel based application in process systems engineering, it can be seen that these novel decomposition approaches have significant benefit over classical decomposition methods and state-of-the-art MILP/MINLP global optimization solvers.

Reaction time assessment in children with ADHD

Attention deficit, impulsivity and hyperactivity are the cardinal features of attention deficit hyperactivity disorder (ADHD) but executive function (EF) disorders, as problems with inhibitory control, working memory and reaction time, besides others EFs, may underlie many of the disturbs associated with the disorder. OBJECTIVE: To examine the reaction time in a computerized test in children with ADHD and normal controls. METHOD: Twenty-three boys (aged 9 to 12) with ADHD diagnosis according to Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition, 2000 (DSM-IV) criteria clinical, without comorbidities, Intelligence Quotient (IQ) >89, never treated with stimulant and fifteen normal controls, age matched were investigated during performance on a voluntary attention psychophysical test. RESULTS: Children with ADHD showed reaction time higher than normal controls. CONCLUSION: A slower reaction time occurred in our patients with ADHD. This findings may be related to problems with the attentional system, that could not maintain an adequate capacity of perceptual input processes and/or in motor output processes, to respond consistently during continuous or repetitive activity.

Finite element analysis of steady-state natural convection problems in fluid-saturated porous media heated from below

In this paper, a progressive asymptotic approach procedure is presented for solving the steady-state Horton-Rogers-Lapwood problem in a fluid-saturated porous medium. The Horton-Rogers-Lapwood problem possesses a bifurcation and, therefore, makes the direct use of conventional finite element methods difficult. Even if the Rayleigh number is high enough to drive the occurrence of natural convection in a fluid-saturated porous medium, the conventional methods will often produce a trivial non-convective solution. This difficulty can be overcome using the progressive asymptotic approach procedure associated with the finite element method. The method considers a series of modified Horton-Rogers-Lapwood problems in which gravity is assumed to tilt a small angle away from vertical. The main idea behind the progressive asymptotic approach procedure is that through solving a sequence of such modified problems with decreasing tilt, an accurate non-zero velocity solution to the Horton-Rogers-Lapwood problem can be obtained. This solution provides a very good initial prediction for the solution to the original Horton-Rogers-Lapwood problem so that the non-zero velocity solution can be successfully obtained when the tilted angle is set to zero. Comparison of numerical solutions with analytical ones to a benchmark problem of any rectangular geometry has demonstrated the usefulness of the present progressive asymptotic approach procedure. Finally, the procedure has been used to investigate the effect of basin shapes on natural convection of pore-fluid in a porous medium. (C) 1997 by John Wiley & Sons, Ltd.

Solution of adsorption problems involving steep moving profiles

The moving finite element collocation method proposed by Kill et al. (1995) Chem. Engng Sci. 51 (4), 2793-2799 for solution of problems with steep gradients is further developed to solve transient problems arising in the field of adsorption. The technique is applied to a model of adsorption in solids with bidisperse pore structures. Numerical solutions were found to match the analytical solution when it exists (i.e. when the adsorption isotherm is linear). The method is simple yet sufficiently accurate for use in adsorption problems, where global collocation methods fail. (C) 1998 Elsevier Science Ltd. All rights reserved.

The conserved quantity theory of causation and chance raising

In this paper I offer an 'integrating account' of singular causation, where the term 'integrating' refers to the following program for analysing causation. There are two intuitions about causation, both of which face serious counterexamples when used as the basis for an analysis of causation. The 'process' intuition, which says that causes and effects are linked by concrete processes, runs into trouble with cases of misconnections', where an event which serves to prevent another fails to do so on a particular occasion and yet the two events are linked by causal processes. The chance raising intuition, according to which causes raise the chance of their effects, easily accounts for misconnections but faces the problem of chance lowering causes, a problem easily accounted for by the process approach. The integrating program attempts to provide an analysis of singular causation by synthesising the two insights, so as to solve both problems. In this paper I show that extant versions of the integrating program due to Eells, Lewis, and Menzies fail to account for the chance-lowering counterexample. I offer a new diagnosis of the chance lowering case, and use that as a basis for an integrating account of causation which does solve both cases. In doing so, I accept various assumptions of the integrating program, in particular that there are no other problems with these two approaches. As an example of the process account, I focus on the recent CQ theory of Wesley Salmon (1997).

Causes of communication problems at work: A qualitative analysis of supervisor-subordinate differences

This study provided information about how individual workers perceive, describe and interpret episodes of problematic communication. Sixteen full-time workers (5 males, 11 females) were interviewed in depth about specific incidents of problematic communication within their workplace. Their descriptions of the attributed causes of the incidents were coded using a categorisation scheme developed from Coupland, Wieman, and Giles' (1991) model of sources of problematic communication. Communication problems were most commonly attributed to individual deficiency and group membership, although there were differences depending on the direction of communication. The most negative attributions (to personality flaws, to lack of skills, and to negative stereotypes of the outgroup) were most commonly applied by individuals to their supervisors, whilst attributions applied to co-workers and subordinates tended to be less negative, or even positive in some instances (where individuals attributed the fault to themselves). Overall, results highlighted distinctions between the perceptions of communication problems with supervisors and with subordinates, and are interpreted with reference to social identity theory.

Cephalometric evaluation of facial pattern and hyoid bone position in children with obstructive sleep apnea syndrome

Objectives: To assess the development of face and hyoid bone in children with obstructive sleep apnea syndrome (OSAS) through lateral cephalometries. Materials and methods: Children aged 7-10 years with mixed dentition and with no previous otorhinolaryngologic, orthodontic or speech therapy treatments were studied. Twenty nasal breathers were compared to 20 mouth breathing children diagnosed as OSAS patients. All children underwent otorhinolaryngologic evaluation and cephalometries; children with OSAS also underwent nocturnal polysomnography in a sleep laboratory. Results: Children with OSAS presented increase in total and lower anterior heights of the face when compared to nasal breathers. In addition, children with OSAS presented a significantly more anterior and inferior position of the hyoid bone than nasal breathers. No significant differences in upper, anterior or posterior heights of the face were observed between groups. Conclusion: The results suggest that there are evident and early changes in facial growth and development among children with OSAS, characterized by increased total and inferior anterior heights of the face, as well as more anterior and inferior position of the hyoid bone. (C) 2010 Elsevier Ireland Ltd. All rights reserved.

The impact of perinatal and socioeconomic factors on mental health problems of children from a poor Brazilian city: a longitudinal study

Low birth weight and preterm birth, and social disadvantage may negatively affect mental health of children, but findings have been inconsistent. To assess the influence of perinatal and social factors on mental health problems in children aged 7-9 years. A random sample of 805 births in So Luis, Brazil was studied in 1997/1998 and again in 2005/2006. Perinatal, socioeconomic and demographic variables were assessed within 24 h after delivery. The Strengths and Difficulties Questionnaire (SDQ) was used to assess mental health problems in the children. Simple and multiple Poisson regressions were used for statistical analysis. The overall prevalence of mental health problems in the total sample was 47.7%. The prevalences of emotional and conduct problems were 58.2 and 48.8%, respectively. Only paternal age (< 20 years) was associated with mental health problems as measured by the full SDQ scale (prevalence ratio PR = 1.27). Children born to single mothers (PR = 1.31) and those with birth weight from 1,500 to 2,499 g (PR = 1.18) and from 2,500 to 2,999 g (PR = 1.17) had a higher risk of emotional problems, but those from low income families had a lower risk (PR = 0.80). Children with a father of less than 20 years had a higher risk of having problems with their peers (PR = 1.75). A maternal education of 9 years or over was inversely associated with peer (PR = 0.70) and conduct problems (PR = 0.73). Girls had a lower risk of conduct (PR = 0.77) and hyperactivity problems (PR = 0.68). A maternal education of 4 years or less increased the risk of hyperactivity (PR = 1.48). Socioeconomic and demographic conditions were better predictors of mental health problems in children than birth weight or preterm birth. However, since most effect sizes were small most mental health problems were, unexplained by the variables in the study.