Geometrically non-linear dynamics of composite four-bar mechanisms

This work intends to demonstrate the importance of a geometrically nonlinear cross-sectional analysis of certain composite beam-based four-bar mechanisms in predicting system dynamic characteristics. All component bars of the mechanism are made of fiber reinforced laminates and have thin rectangular cross-sections. They could, in general, be pre-twisted and/or possess initial curvature, either by design or by defect. They are linked to each other by means of revolute joints. We restrict ourselves to linear materials with small strains within each elastic body (beam). Each component of the mechanism is modeled as a beam based on geometrically non-linear 3-D elasticity theory. The component problems are thus split into 2-D analyses of reference beam cross-sections and non-linear 1-D analyses along the three beam reference curves. For the thin rectangular cross-sections considered here, the 2-D cross-sectional non-linearity is also overwhelming. This can be perceived from the fact that such sections constitute a limiting case between thin-walled open and closed sections, thus inviting the non-linear phenomena observed in both. The strong elastic couplings of anisotropic composite laminates complicate the model further. However, a powerful mathematical tool called the Variational Asymptotic Method (VAM) not only enables such a dimensional reduction, but also provides asymptotically correct analytical solutions to the non-linear cross-sectional analysis. Such closed-form solutions are used here in conjunction with numerical techniques for the rest of the problem to predict multi-body dynamic responses more quickly and accurately than would otherwise be possible. The analysis methodology can be viewed as a three-step procedure: First, the cross-sectional properties of each bar of the mechanism is determined analytically based on an asymptotic procedure, starting from Classical Laminated Shell Theory (CLST) and taking advantage of its thin strip geometry. Second, the dynamic response of the non-linear, flexible four-bar mechanism is simulated by treating each bar as a 1-D beam, discretized using finite elements, and employing energy-preserving and -decaying time integration schemes for unconditional stability. Finally, local 3-D deformations and stresses in the entire system are recovered, based on the 1-D responses predicted in the previous step. With the model, tools and procedure in place, we identify and investigate a few four-bar mechanism problems where the cross-sectional non-linearities are significant in predicting better and critical system dynamic characteristics. This is carried out by varying stacking sequences (i.e. the arrangement of ply orientations within a laminate) and material properties, and speculating on the dominating diagonal and coupling terms in the closed-form non-linear beam stiffness matrix. A numerical example is presented which illustrates the importance of 2-D cross-sectional non-linearities and the behavior of the system is also observed by using commercial software (I-DEAS + NASTRAN + ADAMS). (C) 2012 Elsevier Ltd. All rights reserved.

NON-GLYCINE RESIDUES IN THE DISALLOWED REGIONS OF THE RAMACHANDRAN MAP AND THEIR NEIGHBORHOOD PREFERENCES

The allowed and the ``disallowed'' regions in the celebrated Ramachandran map (phi-psi] map) was elegantly deduced by Ramachandran, Ramakrishnan and Sasisekharan even before the protein crystal structures became available. This powerful map was derived based on rigid geometry of the peptide group and later several investigations on protein crystal structures reported the occurrence of a small fraction of the phi-psi] torsion angles in the disallowed region. The question is what factors make these residues adopt disallowed conformations? Is it driven by the necessity to maintain the overall topology or is it associated with function or is it just that the disallowed conformations are extreme limits of the allowed conformations? Today, with the availability of a large number of high resolution crystal structures, we have revisited this problem. Apart from validating some of the earlier findings such as residue propensities, preferred location in the secondary structure, we have explored their spatial neighborhood preferences using the protein structure network PSN] approach developed in our lab. Finally, the structural and functional implications of the disallowed conformations are examined.

Non-existence of tight neighborly triangulated manifolds with beta(1)=2

All triangulated d-manifolds satisfy the inequality ((f0-d-1)(2)) >= ((d+2)(2))beta(1) for d >= 3. A triangulated d-manifold is called tight neighborly if it attains equality in this bound. For each d >= 3, a (2d + 3)-vertex tight neighborly triangulation of the Sd-1-bundle over S-1 with beta(1) = 1 was constructed by Kuhnel in 1986. In this paper, it is shown that there does not exist a tight neighborly triangulated manifold with beta(1) = 2. In other words, there is no tight neighborly triangulation of (Sd-1 x S-1)(#2) or (Sd-1 (sic) S-1)(#2) for d >= 3. A short proof of the uniqueness of K hnel's complexes for d >= 4 under the assumption beta(1) not equal 0 is also presented.

GRMHD formulation of highly super-Chandrasekhar magnetized white dwarfs: stable configurations of non-spherical white dwarfs

The topic of magnetized super-Chandrasekhar white dwarfs is in the limelight, particularly in the last few years, since our proposal of their existence. By full-scale general relativistic magnetohydrodynamic (GRMHD) numerical analysis, we confirm in this work the existence of stable, highly magnetized, significantly super-Chandrasekhar white dwarfs with mass more than 3 solar mass. While a poloidal field geometry renders the white dwarfs oblate, a toroidal field makes them prolate retaining an overall quasi-spherical shape, as speculated in our earlier work. These white dwarfs are expected to serve as the progenitors of over-luminous type Ia supernovae.

NON-DIMENSIONAL KINETOELASTOSTATIC MAPS FOR COMPLIANT MECHANISMS

How do we assess the capability of a compliant mechanism of given topology and shape? The kinetoelastostatic maps proposed in this paper help answer this question. These maps are drawn in 2D using two non-dimensional quantities, one capturing the nonlinear static response and the other the geometry, material, and applied forces. Geometrically nonlinear finite element analysis is used to create the maps for compliant mechanisms consisting of slender beams. In addition to the topology and shape, the overall proportions and the proportions of the cross-sections of the beam segments are kept fixed for a map. The finite region of the map is parameterized using a non-dimensional quantity defined as the slenderness ratio. The shape and size of the map and the parameterized curves inside it indicate the complete kinetoelastostatic capability of the corresponding compliant mechanism of given topology, shape, and fixed proportions. Static responses considered in this paper include input/output displacement, geometric amplification, mechanical advantage, maximum stress, etc. The maps can be used to compare mechanisms, to choose a suitable mechanism for an application, or re-design as may be needed. The usefulness of the non-dimensional maps is presented with multiple applications of different variety. Non-dimensional portrayal of snap-through mechanisms is one such example. The effect of the shape of the cross-section of the beam segments and the role of different segments in the mechanism as well as extension to 3D compliant mechanisms, the cases of multiple inputs and outputs, and moment loads are also explained. The effects of disproportionate changes on the maps are also analyzed.

Cold gas in cluster cores: global stability analysis and non-linear simulations of thermal instability

We perform global linear stability analysis and idealized numerical simulations in global thermal balance to understand the condensation of cold gas from hot/virial atmospheres (coronae), in particular the intracluster medium (ICM). We pay particular attention to geometry (e.g. spherical versus plane-parallel) and the nature of the gravitational potential. Global linear analysis gives a similar value for the fastest growing thermal instability modes in spherical and Cartesian geometries. Simulations and observations suggest that cooling in haloes critically depends on the ratio of the cooling time to the free-fall time (t(cool)/t(ff)). Extended cold gas condenses out of the ICM only if this ratio is smaller than a threshold value close to 10. Previous works highlighted the difference between the nature of cold gas condensation in spherical and plane-parallel atmospheres; namely, cold gas condensation appeared easier in spherical atmospheres. This apparent difference due to geometry arises because the previous plane-parallel simulations focused on in situ condensation of multiphase gas but spherical simulations studied condensation anywhere in the box. Unlike previous claims, our non-linear simulations show that there are only minor differences in cold gas condensation, either in situ or anywhere, for different geometries. The amount of cold gas depends on the shape of tcool/tff; gas has more time to condense if gravitational acceleration decreases towards the centre. In our idealized plane-parallel simulations with heating balancing cooling in each layer, there can be significant mass/energy/momentum transfer across layers that can trigger condensation and drive tcool/tff far beyond the critical value close to 10.

Boundary knot method for 2D and 3D Helmholtz and convection-diffusion problems under complicated geometry

The boundary knot method (BKM) of very recent origin is an inherently meshless, integration-free, boundary-type, radial basis function collocation technique for the numerical discretization of general partial differential equation systems. Unlike the method of fundamental solutions, the use of non-singular general solution in the BKM avoids the unnecessary requirement of constructing a controversial artificial boundary outside the physical domain. The purpose of this paper is to extend the BKM to solve 2D Helmholtz and convection-diffusion problems under rather complicated irregular geometry. The method is also first applied to 3D problems. Numerical experiments validate that the BKM can produce highly accurate solutions using a relatively small number of knots. For inhomogeneous cases, some inner knots are found necessary to guarantee accuracy and stability. The stability and convergence of the BKM are numerically illustrated and the completeness issue is also discussed.

Non-Linear Wave-Induced Transient Response of Soil Around a Trenched Pipeline

Submarine pipelines are always trenched within a seabed for reducing wave loads and thereby enhancing their stability. Based on Biot’s poroelastic theory, a two-dimensional finite element model is developed to investigate non-linear wave-induced responses of soil around a trenched pipeline, which is verified with the flume test results by Sudhan et al. [Sudhan, C.M., Sundar, V., Rao, S.N., 2002. Wave induced forces around buried pipeline. Ocean Engineering, 29, 533–544] and Turcotte et al. [Turcotte, B.R., Liu, P.L.F., Kulhawy, F.H., 1984. Laboratory evaluation of wave tank parameters for wave-sediment interaction. Joseph H. Defree Hydraulic Laboratory Report 84-1, School of Civil and Environmental Engineering, Cornell University]. Non-linear wave-induced transient pore pressure around pipeline at various phases of wave loading is examined firstly. Unlike most previous investigations, in which only a single sediment layer and linear wave loading were concerned, in this study, the influences of the non-linearity of wave loading, the physical properties of backfill materials and the geometry profile of trenches on the excess pore pressures within the soil around pipeline, respectively, were explored, taking into account the in situ conditions of buried pipeline in the shallow ocean zones. Based on the parametric study, it is concluded that the shear modulus and permeability of backfill soils significantly affect the wave-induced excess pore pressures around trenched pipeline, and that the effect of wave non-linearity becomes more pronounced and comparable with that of trench depth, especially at high wave steepness in shallow water.

Fast numerical methods for mixed, singular Helmholtz boundary value problems and Laplace eigenvalue problems - with applications to antenna design, sloshing, electromagnetic scattering and spectral geometry

This thesis presents a novel class of algorithms for the solution of scattering and eigenvalue problems on general two-dimensional domains under a variety of boundary conditions, including non-smooth domains and certain "Zaremba" boundary conditions - for which Dirichlet and Neumann conditions are specified on various portions of the domain boundary. The theoretical basis of the methods for the Zaremba problems on smooth domains concern detailed information, which is put forth for the first time in this thesis, about the singularity structure of solutions of the Laplace operator under boundary conditions of Zaremba type. The new methods, which are based on use of Green functions and integral equations, incorporate a number of algorithmic innovations, including a fast and robust eigenvalue-search algorithm, use of the Fourier Continuation method for regularization of all smooth-domain Zaremba singularities, and newly derived quadrature rules which give rise to high-order convergence even around singular points for the Zaremba problem. The resulting algorithms enjoy high-order convergence, and they can tackle a variety of elliptic problems under general boundary conditions, including, for example, eigenvalue problems, scattering problems, and, in particular, eigenfunction expansion for time-domain problems in non-separable physical domains with mixed boundary conditions.

Effects of atmosphere visibility on performances of non-line-of-sight ultraviolet communication systems

Dependence of performances of non-line-of-sight (NLOS) solar-blind ultraviolet (UV) communication systems on atmosphere visibility is investigated numerically by correlating the propagation of UV radiation with the visibility. A simplified solar-blind UV atmospheric propagation model is introduced, and the NLOS UV communication system model is constituted based on the single scattering assumption. Using the model, numerical simulation is conducted for two typical geometry configurations and different modulation formats. The results indicate that the performance of the NLOS UV communication system is insensitive to variation of visibility in quite a large range, and deteriorates significantly only in very low-visibility weather, and is also dependent on the geometry configuration of the system. The results also show that the pulse position modulation (PPM) is preferable due to its high-power efficiency to improve the system performance. (c) 2007 Elsevier GmbH. All rights reserved.

Investigation of the effect of combustor cooling geometry on acoustic energy absorption

To control combustion instabilities occurring in LPP gas turbine combustors, several active and passive systems have been developed in recent years. The combustion chamber cooling geometry has the potential to influence instability feedback loops by absorbing acoustical energy inside the combustor. The design of the cooling liner and the geometry of the cooling plenum and the cooling air flow rate have a significant influence on the absorption characteristics of the system. This paper presents the results of a cold flow study which was carried out in the course of a comprehensive study on the influence of the cooling geometry on combustor thermoacoustics. Absorption characteristics of three different cooling liner geometries and non-perforated plates were determined over a frequency range from 50 Hz to 600 Hz for different cooling flow rates and different cooling plenum volumes. The experimental results compared well with results from a low order thermoacoustic network model. The acoustic energy absorption spectrum of a cooling liner with 90°-hole configuration was found to be strongly dependent on cooling flow rate and cooling plenum volume, whereas the absorption spectrum of cooling liners with 25°-holes were found to be strongly dependent on the cooling plenum volume, but less dependent on the cooling air flow rate. All cooling liner setups with perforations were capable of increased acoustic absorption over a broad band of frequencies compared to the case of non-perforated combustor walls. © 2010 by Johannes Schmidt.

Effect of flow-geometry on turbulence-scalar interaction in premixed flames

Turbulent combustion of stoichiometric hydrogen-air mixture is simulated using direct numerical simulation methodology, employing complex chemical kinetics. Two flame configurations, freely propagating and V-flames stabilized behind a hot rod, are simulated. The results are analyzed to study the influence of flame configuration on the turbulence-scalar interaction, which is critical for the scalar gradient generation processes. The result suggests that this interaction process is not influenced by the flame configuration and the flame normal is found to align with the most extensive strain in the region of intense heat release. The combustion in the rod stabilized flame is found to be flamelet like in an average sense and the growth of flame-brush thickness with the downstream distance is represented well by Taylor theory of turbulent diffusion, when the flame-brushes are non-interacting. The thickness is observed to saturate when the flame-brushes interact, which is found to occur in the simulated rod stabilized flame with Taylor micro-scale Reynolds number of 97. © 2011 American Institute of Physics.

Nonsmooth geometry and collapse of flexible structures under smooth loads

While static equilibria of flexible strings subject to various load types (gravity, hydrostatic pressure, Newtonian wind) is well understood textbook material, the combinations of the very same loads can give rise to complex spatial behaviour at the core of which is the unilateral material constraint prohibiting compressive loads. While the effects of such constraints have been explored in optimisation problems involving straight cables, the geometric complexity of physical configurations has not yet been addressed. Here we show that flexible strings subject to combined smooth loads may not have smooth solutions in certain ranges of the load ratios. This non-smooth phenomenon is closely related to the collapse geometry of inflated tents. After proving the nonexistence of smooth solutions for a broad family of loadings we identify two alternative, critical geometries immediately preceding the collapse. We verify these analytical results by dynamical simulation of flexible chains as well as with simple table-top experiments with an inflated membrane.

Geometry and interaction of structures in homogeneous isotropic turbulence

A strategy to extract turbulence structures from direct numerical simulation (DNS) data is described along with a systematic analysis of geometry and spatial distribution of the educed structures. A DNS dataset of decaying homogeneous isotropic turbulence at Reynolds number Reλ = 141 is considered. A bandpass filtering procedure is shown to be effective in extracting enstrophy and dissipation structures with their smallest scales matching the filter width, L. The geometry of these educed structures is characterized and classified through the use of two non-dimensional quantities, planarity' and filamentarity', obtained using the Minkowski functionals. The planarity increases gradually by a small amount as L is decreased, and its narrow variation suggests a nearly circular cross-section for the educed structures. The filamentarity increases significantly as L decreases demonstrating that the educed structures become progressively more tubular. An analysis of the preferential alignment between the filtered strain and vorticity fields reveals that vortical structures of a given scale L are most likely to align with the largest extensional strain at a scale 3-5 times larger than L. This is consistent with the classical energy cascade picture, in which vortices of a given scale are stretched by and absorb energy from structures of a somewhat larger scale. The spatial distribution of the educed structures shows that the enstrophy structures at the 5η scale (where η is the Kolmogorov scale) are more concentrated near the ones that are 3-5 times larger, which gives further support to the classical picture. Finally, it is shown by analysing the volume fraction of the educed enstrophy structures that there is a tendency for them to cluster around a larger structure or clusters of larger structures. Copyright © 2012 Cambridge University Press.

Geometrical learning, descriptive geometry, and biomimetic pattern recognition

Studies on learning problems from geometry perspective have attracted an ever increasing attention in machine learning, leaded by achievements on information geometry. This paper proposes a different geometrical learning from the perspective of high-dimensional descriptive geometry. Geometrical properties of high-dimensional structures underlying a set of samples are learned via successive projections from the higher dimension to the lower dimension until two-dimensional Euclidean plane, under guidance of the established properties and theorems in high-dimensional descriptive geometry. Specifically, we introduce a hyper sausage like geometry shape for learning samples and provides a geometrical learning algorithm for specifying the hyper sausage shapes, which is then applied to biomimetic pattern recognition. Experimental results are presented to show that the proposed approach outperforms three types of support vector machines with either a three degree polynomial kernel or a radial basis function kernel, especially in the cases of high-dimensional samples of a finite size. (c) 2005 Elsevier B.V. All rights reserved.