Geometric optimization methods for the analysis of gene expression data

DNA microarrays provide such a huge amount of data that unsupervised methods are required to reduce the dimension of the data set and to extract meaningful biological information. This work shows that Independent Component Analysis (ICA) is a promising approach for the analysis of genome-wide transcriptomic data. The paper first presents an overview of the most popular algorithms to perform ICA. These algorithms are then applied on a microarray breast-cancer data set. Some issues about the application of ICA and the evaluation of biological relevance of the results are discussed. This study indicates that ICA significantly outperforms Principal Component Analysis (PCA).

DISCUSSION OF MECHANISMS OF ACCELERATED AND RETARDED FATIGUE CRACK GROWTH.

The relevance of the effective stress intensity range to crack growth is considered for constant and for variable amplitude loading. The accelerated and retarded growth associated with simple programmed loadings is reported for two steels and an aluminium alloy. The load interaction effects are due to several competing mechanisms, and not due to the single, popular mechanism of crack closure.

Widening engineering horizons: Addressing the complexity of sustainable development

The complex, fragmented and diverse aspects of a sustainable development perspective are translated into an eight-point framework that defines a problem boundary larger than that traditionally adopted by civil engineers. This leads to practical questions intended to inform engineers who ask 'am I being sustainable?' during project implementation. The value of the questions is tested against a case history of a wastewater treatment project. This demonstrates the relevance of the questions to successive project delivery phases of defining the problem, choosing a solution and implementing that solution through design, construction and operation. The case history highlights that answers to several of the additional questions raised by considering this wider problem space are currently buried within government and clients' policies, regulations and standard practice; these answers may not be accessible to the professional engineer.

Task-dependent coordination of rapid bimanual motor responses.

Optimal feedback control postulates that feedback responses depend on the task relevance of any perturbations. We test this prediction in a bimanual task, conceptually similar to balancing a laden tray, in which each hand could be perturbed up or down. Single-limb mechanical perturbations produced long-latency reflex responses ("rapid motor responses") in the contralateral limb of appropriate direction and magnitude to maintain the tray horizontal. During bimanual perturbations, rapid motor responses modulated appropriately depending on the extent to which perturbations affected tray orientation. Specifically, despite receiving the same mechanical perturbation causing muscle stretch, the strongest responses were produced when the contralateral arm was perturbed in the opposite direction (large tray tilt) rather than in the same direction or not perturbed at all. Rapid responses from shortening extensors depended on a nonlinear summation of the sensory information from the arms, with the response to a bimanual same-direction perturbation (orientation maintained) being less than the sum of the component unimanual perturbations (task relevant). We conclude that task-dependent tuning of reflexes can be modulated online within a single trial based on a complex interaction across the arms.

Study of second generation, high-temperature superconducting coils: Determination of critical current

This paper presents the modeling of second generation (2 G) high-temperature superconducting (HTS) pancake coils using finite element method. The axial symmetric model can be used to calculate current and magnetic field distribution inside the coil. The anisotropic characteristics of 2 G tapes are included in the model by direct interpolation. The model is validated by comparing to experimental results. We use the model to study critical currents of 2 G coils and find that 100μV/m is too high a criterion to determine long-term operating current of the coils, because the innermost turns of a coil will, due to the effect of local magnetic field, reach their critical current much earlier than outer turns. Our modeling shows that an average voltage criterion of 20μV/m over the coil corresponds to the point at which the innermost turns' electric field exceeds 100μV/m. So 20μV/m is suggested to be the critical current criterion of the HTS coil. The influence of background field on the coil critical current is also studied in the paper. © 2012 American Institute of Physics.

Validation of a lattice Boltzmann model for gas-solid reactions with experiments

A lattice Boltzmann method is used to model gas-solid reactions where the composition of both the gas and solid phase changes with time, while the boundary between phases remains fixed. The flow of the bulk gas phase is treated using a multiple relaxation time MRT D3Q19 model; the dilute reactant is treated as a passive scalar using a single relaxation time BGK D3Q7 model with distinct inter- and intraparticle diffusivities. A first-order reaction is incorporated by modifying the method of Sullivan et al. [13] to include the conversion of a solid reactant. The detailed computational model is able to capture the multiscale physics encountered in reactor systems. Specifically, the model reproduced steady state analytical solutions for the reaction of a porous catalyst sphere (pore scale) and empirical solutions for mass transfer to the surface of a sphere at Re=10 (particle scale). Excellent quantitative agreement between the model and experiments for the transient reduction of a single, porous sphere of Fe 2O 3 to Fe 3O 4 in CO at 1023K and 10 5Pa is demonstrated. Model solutions for the reduction of a packed bed of Fe 2O 3 (reactor scale) at identical conditions approached those of experiments after 25 s, but required prohibitively long processor times. The presented lattice Boltzmann model resolved successfully mass transport at the pore, particle and reactor scales and highlights the relevance of LB methods for modelling convection, diffusion and reaction physics. © 2012 Elsevier Inc.

Structural and thermodynamic properties of different phases of supercooled liquid water.

Computer simulation results are reported for a realistic polarizable potential model of water in the supercooled region. Three states, corresponding to the low density amorphous ice, high density amorphous ice, and very high density amorphous ice phases are chosen for the analyses. These states are located close to the liquid-liquid coexistence lines already shown to exist for the considered model. Thermodynamic and structural quantities are calculated, in order to characterize the properties of the three phases. The results point out the increasing relevance of the interstitial neighbors, which clearly appear in going from the low to the very high density amorphous phases. The interstitial neighbors are found to be, at the same time, also distant neighbors along the hydrogen bonded network of the molecules. The role of these interstitial neighbors has been discussed in connection with the interpretation of recent neutron scattering measurements. The structural properties of the systems are characterized by looking at the angular distribution of neighboring molecules, volume and face area distribution of the Voronoi polyhedra, and order parameters. The cumulative analysis of all the corresponding results confirms the assumption that a close similarity between the structural arrangement of molecules in the three explored amorphous phases and that of the ice polymorphs I(h), III, and VI exists.

A PDE viewpoint on basic properties of coordination algorithms with symmetries

Several recent control applications consider the coordination of subsystems through local interaction. Often the interaction has a symmetry in state space, e.g. invariance with respect to a uniform translation of all subsystem values. The present paper shows that in presence of such symmetry, fundamental properties can be highlighted by viewing the distributed system as the discrete approximation of a partial differential equation. An important fact is that the symmetry on the state space differs from the popular spatial invariance property, which is not necessary for the present results. The relevance of the viewpoint is illustrated on two examples: (i) ill-conditioning of interaction matrices in coordination/consensus problems and (ii) the string instability issue. ©2009 IEEE.

Stabilization of periodic orbits in a wedge billiard

This paper introduces a stabilization problem for an elementary impact control system in the plane. The rich dynamical properties of the wedge billiard, combined to the relevance of the associated stabilization problem for feedback control issues in legged robotics make it a valuable benchmark for energy-based stabilization of impact control systems.

First-principles energetics of water clusters and ice: a many-body analysis.

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for cluster, solid, and liquid forms of water. Recent work has stressed the importance of DFT errors in describing dispersion, but we note that errors in other parts of the energy may also contribute. We obtain information about the nature of DFT errors by using a many-body separation of the total energy into its 1-body, 2-body, and beyond-2-body components to analyze the deficiencies of the popular PBE and BLYP approximations for the energetics of water clusters and ice structures. The errors of these approximations are computed by using accurate benchmark energies from the coupled-cluster technique of molecular quantum chemistry and from quantum Monte Carlo calculations. The systems studied are isomers of the water hexamer cluster, the crystal structures Ih, II, XV, and VIII of ice, and two clusters extracted from ice VIII. For the binding energies of these systems, we use the machine-learning technique of Gaussian Approximation Potentials to correct successively for 1-body and 2-body errors of the DFT approximations. We find that even after correction for these errors, substantial beyond-2-body errors remain. The characteristics of the 2-body and beyond-2-body errors of PBE are completely different from those of BLYP, but the errors of both approximations disfavor the close approach of non-hydrogen-bonded monomers. We note the possible relevance of our findings to the understanding of liquid water.

Modelling of settlement induced building damage

This thesis focuses on the modelling of settlement induced damage to masonry buildings. In densely populated areas, the need for new space is nowadays producing a rapid increment of underground excavations. Due to the construction of new metro lines, tunnelling activity in urban areas is growing. One of the consequences is a greater attention to the risk of damage on existing structures. Thus, the assessment of potential damage of surface buildings has become an essential stage in the excavation projects in urban areas (Chapter 1). The current damage risk assessment procedure is based on strong simplifications, which not always lead to conservative results. Object of this thesis is the development of an improved damage classification system, which takes into account the parameters influencing the structural response to settlement, like the non-linear behaviour of masonry and the soil-structure interaction. The methodology used in this research is based on experimental and numerical modelling. The design and execution of an experimental benchmark test representative of the problem allows to identify the principal factors and mechanisms involved. The numerical simulations enable to generalize the results to a broader range of physical scenarios. The methodological choice is based on a critical review of the currently available procedures for the assessment of settlement-induced building damage (Chapter 2). A new experimental test on a 1/10th masonry façade with a rubber base interface is specifically designed to investigate the effect of soil-structure interaction on the tunnelling-induced damage (Chapter 3). The experimental results are used to validate a 2D semi-coupled finite element model for the simulation of the structural response (Chapter 4). The numerical approach, which includes a continuum cracking model for the masonry and a non-linear interface to simulate the soil-structure interaction, is then used to perform a sensitivity study on the effect of openings, material properties, initial damage, initial conditions, normal and shear behaviour of the base interface and applied settlement profile (Chapter 5). The results assess quantitatively the major role played by the normal stiffness of the soil-structure interaction and by the material parameters defining the quasi-brittle masonry behaviour. The limitation of the 2D modelling approach in simulating the progressive 3D displacement field induced by the excavation and the consequent torsional response of the building are overcome by the development of a 3D coupled model of building, foundation, soil and tunnel (Chapter 6). Following the same method applied to the 2D semi-coupled approach, the 3D model is validated through comparison with the monitoring data of a literature case study. The model is then used to carry out a series of parametric analyses on geometrical factors: the aspect ratio of horizontal building dimensions with respect to the tunnel axis direction, the presence of adjacent structures and the position and alignment of the building with respect to the excavation (Chapter 7). The results show the governing effect of the 3D building response, proving the relevance of 3D modelling. Finally, the results from the 2D and 3D parametric analyses are used to set the framework of an overall damage model which correlates the analysed structural features with the risk for the building of being damaged by a certain settlement (Chapter 8). This research therefore provides an increased experimental and numerical understanding of the building response to excavation-induced settlements, and sets the basis for an operational tool for the risk assessment of structural damage (Chapter 9).