The most persistent soft-clique in a set of sampled graphs

When searching for characteristic subpatterns in potentially noisy graph data, it appears self-evident that having multiple observations would be better than having just one. However, it turns out that the inconsistencies introduced when different graph instances have different edge sets pose a serious challenge. In this work we address this challenge for the problem of finding maximum weighted cliques. We introduce the concept of most persistent soft-clique. This is subset of vertices, that 1) is almost fully or at least densely connected, 2) occurs in all or almost all graph instances, and 3) has the maximum weight. We present a measure of clique-ness, that essentially counts the number of edge missing to make a subset of vertices into a clique. With this measure, we show that the problem of finding the most persistent soft-clique problem can be cast either as: a) a max-min two person game optimization problem, or b) a min-min soft margin optimization problem. Both formulations lead to the same solution when using a partial Lagrangian method to solve the optimization problems. By experiments on synthetic data and on real social network data we show that the proposed method is able to reliably find soft cliques in graph data, even if that is distorted by random noise or unreliable observations. Copyright 2012 by the author(s)/owner(s).

On the existence of Zeno behavior in hybrid systems with non-isolated zeno equilibria

This paper presents proof-certificate based sufficient conditions for the existence of Zeno behavior in hybrid systems near non-isolated Zeno equilibria. To establish these conditions, we first prove sufficient conditions for Zeno behavior in a special class of hybrid systems termed first quadrant interval hybrid systems. The proof-certificate sufficient conditions are then obtained through a collection of functions that effectively "reduce" a general hybrid system to a first quadrant interval hybrid system. This paper concludes with an application of these ideas to Lagrangian hybrid systems, resulting in easily verifiable sufficient conditions for Zeno behavior. © 2008 IEEE.

A theoretical approach to the passive control of spiral vortex breakdown

Previous numerical simulations have shown that vortex breakdown starts with the formation of a steady axisymmetric bubble and that an unsteady spiralling mode then develops on top of this.We study how this spiral mode of vortex breakdown might be suppressed or promoted. We use a Lagrangian approach to identify regions of the flow which are sensitive to small open-loop steady and unsteady (harmonic) forces. We find these regions to be upstream of the vortex breakdown bubble. We investigate passive control using a small axisymmetric control ring. In this case, the steady and unsteady control forces are caused by the drag force on the control ring. We find a narrow region upstream of the bubble where the control ring will stabilise the flow and we verify this using numerical simulations. © 2012 IEEE.

Lyapunov theory for zeno stability

Zeno behavior is a dynamic phenomenon unique to hybrid systems in which an infinite number of discrete transitions occurs in a finite amount of time. This behavior commonly arises in mechanical systems undergoing impacts and optimal control problems, but its characterization for general hybrid systems is not completely understood. The goal of this paper is to develop a stability theory for Zeno hybrid systems that parallels classical Lyapunov theory; that is, we present Lyapunov-like sufficient conditions for Zeno behavior obtained by mapping solutions of complex hybrid systems to solutions of simpler Zeno hybrid systems defined on the first quadrant of the plane. These conditions are applied to Lagrangian hybrid systems, which model mechanical systems undergoing impacts, yielding simple sufficient conditions for Zeno behavior. Finally, the results are applied to robotic bipedal walking. © 2012 IEEE.

Time-resolved particle image velocimetry within the nozzle of a drop-on-demand printhead

Time-resolved particle image velocimetry (PIV) has been performed inside the nozzle of a commercially available inkjet print-head to obtain the time-dependent velocity waveform. A printhead with a single transparent nozzle 80 μm in orifice diameter was used to eject single droplets at a speed of 5 m/s. An optical microscope was used with an ultra-high-speed camera to capture the motion of particles suspended in a transparent liquid at the center of the nozzle and above the fluid meniscus at a rate of half a million frames per second. Time-resolved velocity fields were obtained from a fluid layer approximately 200 μm thick within the nozzle for a complete jetting cycle. A Lagrangian finite-element numerical model with experimental measurements as inputs was used to predict the meniscus movement. The model predictions showed good agreement with the experimental results. This work provides the first experimental verification of physical models and numerical simulations of flows within a drop-on-demand nozzle. © 2012 Society for Imaging Science and Technology.

Multi-objective optimization of a fluid structure interaction benchmarking

The integration and application of a new multi-objective tabu search optimization algorithm for Fluid Structure Interaction (FSI) problems are presented. The aim is to enhance the computational design process for real world applications and to achieve higher performance of the whole system for the four considered objectives. The described system combines the optimizer with a well established FSI solver which is based on the fully implicit, monolithic formuFlation of the problem in the Arbitrary Lagrangian-Eulerian FEM approach. The proposed solver resolves the proposed uid-structure interaction benchmark which describes the self-induced elastic deformation of a beam attached to a cylinder in laminar channel ow. The optimized ow characteristics of the aforementioned geometrical arrangement illustrate the performance of the system in two dimensions. Special emphasis is given to the analysis of the simulation package, which is of high accuracy and is the core of application. The design process identifies the best combination of ow features for optimal system behavior and the most important objectives. In addition, the presented methodology has the potential to run in parallel, which will significantly speed-up the elapsed time. Finite Element Method (FEM), Fluid-Structure Interaction (FSI), Multi-Ojective Tabu search (MOTS2). Copyright © 2013 Tech Science Press.

Multi-scale modelling of granular avalanches

Avalanches, debris flows, and landslides are geophysical hazards, which involve rapid mass movement of granular solids, water and air as a single-phase system. The dynamics of a granular flow involve at least three distinct scales: the micro-scale, meso-scale, and the macro-scale. This study aims to understand the ability of continuum models to capture the micro-mechanics of dry granular collapse. Material Point Method (MPM), a hybrid Lagrangian and Eulerian approach, with Mohr-Coulomb failure criterion is used to describe the continuum behaviour of granular column collapse, while the micromechanics is captured using Discrete Element Method (DEM) with tangential contact force model. The run-out profile predicted by the continuum simulations matches with DEM simulations for columns with small aspect ratios ('h/r' < 2), however MPM predicts larger run-out distances for columns with higher aspect ratios ('h/r' > 2). Energy evolution studies in DEM simulations reveal higher collisional dissipation in the initial free-fall regime for tall columns. The lack of a collisional energy dissipation mechanism in MPM simulations results in larger run-out distances. Micro-structural effects, such as shear band formations, were observed both in DEM and MPM simulations. A sliding flow regime is observed above the distinct passive zone at the core of the column. Velocity profiles obtained from both the scales are compared to understand the reason for a slow flow run-out mobilization in MPM simulations. © 2013 AIP Publishing LLC.

The dynamic response of edge clamped plates loaded by spherically expanding sand shells

The dynamic deformation of both edge clamped stainless steel sandwich panels with a pyramidal truss core and equal mass monolithic plates loaded by spherically expanding shells of dry and water saturated sand has been investigated, both experimentally and via a particle based simulation methodology. The spherically expanding sand shell is generated by detonating a sphere of explosive surrounded by a shell of either dry or water saturated synthetic sand. The measurements show that the sandwich panel and plate deflections decrease with increasing stand-off between the center of the charge and the front of the test structures. Moreover, for the same charge and sand mass, the deflections of the plates are significantly higher in the water saturated sand case compared to that of dry sand. For a given stand-off, the mid-span deflection of the sandwich panel rear faces was substantially less than that of the corresponding monolithic plate for both the dry and water saturated sand cases. The experiments were simulated via a coupled discrete-particle/ finite element scheme wherein the high velocity impacting sand is modeled by interacting particles while the plate is modeled within a Lagrangian finite element setting. The simulations are in good agreement with the measurements for the dry sand impact of both the monolithic and sandwich structures. However, the simulations underestimate the effect of stand-off in the case of the water saturated sand explosion, i.e. the deflections decrease more sharply with increasing stand-off in the experiments compared to the simulations. The simulations reveal that the momentum transmitted into the sandwich and monolithic plate structures by the sand shell is approximately the same, consistent with a small fluid-structure interaction effect. The smaller deflection of the sandwich panels is therefore primarily due to the higher bending strength of sandwich structures. © 2013 The Authors. Published by Elsevier Ltd. All rights reserved.

Application of the characteristic CIP method to a shallow water model on the sphere

Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the efficiency of the characteristic constrained interpolation profile (CIP) method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.

耗散粒子动力学处理复杂固体壁面的一种有效方法

耗散粒子动力学(dissipative particle dynamics,DPD)作为一种介观尺度拉格朗日型粒子方法,已经成功地应用于微纳米流动和生化科技的研究中. 复杂固体壁面的处理和壁面边界条件的实施一直是DPD方法发展及应用的一个障碍. 提出了处理复杂固体壁面的一种新的方法. 复杂固体区域通过冻结随机分布并且达到平衡状态的DPD粒子代表;所冻结的DPD粒子位于临近流动区域的一个截距内;在靠近固体壁面的流动区域中设置流动反弹层,当流动DPD粒子进入此流动层后反弹回流动区域. 应用这种固体壁面处理方法.

Role of the N

New measurement by CELSIUS-WASA Collaboration on the pp →pnπ+ reaction reveals clear evidence for the presence of the Roper resonance N∗(1440) which has been ignored in previous theoretical calculations. In this article, based on an effective Lagrangian approach and available knowledge on the Roper resonance, we investigate the role of the Roper resonance for the pp→pnπ+ reaction. It is found that the contribution from the Roper resonance N∗(1440) becomes significant for kinetic energy above 1.1 GeV, consistent with the new experimental observation. The t -channel σ-meson exchange is dominant for the production of the Roper resonance.

Evidence of N-*(1535) resonance contribution in the pn -> d phi reaction

The N ∗(1535) resonance contributions to the pn → dφ reaction are evaluated in an effective Lagrangian model. The π-, η-, and ρ-meson exchange are considered. It is shown that the contributions from π- and ρ-meson exchange are dominant, while the contribution from η-meson exchange is negligibly small. Our theoretical results reproduce the experimental data of both total cross section and angular distribution well. This is more evidence that the N ∗(1535) resonance has a large s ¯s component leading to a large coupling to Nφ, which may be the real origin of the Okubo-Zweig-Iizuka rule violation in the πN and pN reactions.

N*(1535) contribution in pp -> pp eta ' and pn -> d phi

In an effective Lagrangian model we find that the N*(1535) resonance contribution might be important to the interpretation of the present data of the pp -> pp eta' and pn -> d phi reactions. The strong coupling strength of N*(1.535) to eta' and phi are indicated, and the possible implication to the intrinsic component of N*(1535) is explored. These results may provide hints to the real origin of the OZI rule violation. It is stressed that further measurements could be performed at the Cooling Storage Ring (CSR) at Lanzhou of China.

Symmetry energy and isovector giant dipole resonance in finite nuclei

We study the relationship between the properties of the isovector giant dipole resonance of finite nuclei and the symmetry energy in the framework of the relativistic mean field theory with six different parameter sets of nonlinear effective Lagrangian. A strong linear correlation of excited energies of the dipole resonance in finite nuclei and symmetry energy at and below the saturation density is found. This linear correlation leads to the symmetry energy at the saturation density at the interval 33.0MeV <= S(po) <= 37.0 MeV. The comparison to the present experimental data in the soft dipole mode of (132) Sn constrains approximately the symmetry energy at p = 0.1 fm(-3) at the interval 21.2MeV similar to 22.5 MeV. It is proposed that a precise measurement of the soft dipole mode in neutron rich nuclei could set up an important constraint on the equation of state for asymmetric nuclear matter.

Role of the N*(1535) in eta ' production

We study the near-threshold eta ' production mechanism in nucleon-nucleon and pi N collisions under the assumption that subthreshold resonance N*(1535) is predominant. In an effective Lagrangian approach that gives a reasonable description to the pN -> pN eta and pi(-) p -> n eta reactions, we find that the excitation of N*(1535) resonance from the t- channel p exchange makes the dominate contribution to the pN -> pN eta ' process, and a value of 6.5 for the ratio of s(pp -> pp eta ') to sigma (pp -> pp eta ') is predicted. A strongcoupling strength ofN*(1535) to eta ' N (g(eta ' NN*)(2)/4 pi = 1.1) is extracted from a combined analysis to pp -> pp eta ' and pi N -> N eta ', and the possible implication to the intrinsic component of N*(1535) is explored.