Conservation Properties of the Trapezoidal Rule in Linear Time Domain Analysis of Acoustics and Structures

The trapezoidal rule, which is a special case of the Newmark family of algorithms, is one of the most widely used methods for transient hyperbolic problems. In this work, we show that this rule conserves linear and angular momenta and energy in the case of undamped linear elastodynamics problems, and an ``energy-like measure'' in the case of undamped acoustic problems. These conservation properties, thus, provide a rational basis for using this algorithm. In linear elastodynamics problems, variants of the trapezoidal rule that incorporate ``high-frequency'' dissipation are often used, since the higher frequencies, which are not approximated properly by the standard displacement-based approach, often result in unphysical behavior. Instead of modifying the trapezoidal algorithm, we propose using a hybrid finite element framework for constructing the stiffness matrix. Hybrid finite elements, which are based on a two-field variational formulation involving displacement and stresses, are known to approximate the eigenvalues much more accurately than the standard displacement-based approach, thereby either bypassing or reducing the need for high-frequency dissipation. We show this by means of several examples, where we compare the numerical solutions obtained using the displacement-based and hybrid approaches against analytical solutions.

Stability of the flow in a soft tube deformed due to an applied pressure gradient

A linear stability analysis is carried out for the flow through a tube with a soft wall in order to resolve the discrepancy of a factor of 10 for the transition Reynolds number between theoretical predictions in a cylindrical tube and the experiments of Verma and Kumaran J. Fluid Mech. 705, 322 (2012)]. Here the effect of tube deformation (due to the applied pressure difference) on the mean velocity profile and pressure gradient is incorporated in the stability analysis. The tube geometry and dimensions are reconstructed from experimental images, where it is found that there is an expansion and then a contraction of the tube in the streamwise direction. The mean velocity profiles at different downstream locations and the pressure gradient, determined using computational fluid dynamics, are found to be substantially modified by the tube deformation. The velocity profiles are then used in a linear stability analysis, where the growth rates of perturbations are calculated for the flow through a tube with the wall modeled as a neo-Hookean elastic solid. The linear stability analysis is carried out for the mean velocity profiles at different downstream locations using the parallel flow approximation. The analysis indicates that the flow first becomes unstable in the downstream converging section of the tube where the flow profile is more pluglike when compared to the parabolic flow in a cylindrical tube. The flow is stable in the upstream diverging section where the deformation is maximum. The prediction for the transition Reynolds number is in good agreement with experiments, indicating that the downstream tube convergence and the consequent modification in the mean velocity profile and pressure gradient could reduce the transition Reynolds number by an order of magnitude.

Optimal Linear Glauber Model

Contrary to the actual nonlinear Glauber model, the linear Glauber model (LGM) is exactly solvable, although the detailed balance condition is not generally satisfied. This motivates us to address the issue of writing the transition rate () in a best possible linear form such that the mean squared error in satisfying the detailed balance condition is least. The advantage of this work is that, by studying the LGM analytically, we will be able to anticipate how the kinetic properties of an arbitrary Ising system depend on the temperature and the coupling constants. The analytical expressions for the optimal values of the parameters involved in the linear are obtained using a simple Moore-Penrose pseudoinverse matrix. This approach is quite general, in principle applicable to any system and can reproduce the exact results for one dimensional Ising system. In the continuum limit, we get a linear time-dependent Ginzburg-Landau equation from the Glauber's microscopic model of non-conservative dynamics. We analyze the critical and dynamic properties of the model, and show that most of the important results obtained in different studies can be reproduced by our new mathematical approach. We will also show in this paper that the effect of magnetic field can easily be studied within our approach; in particular, we show that the inverse of relaxation time changes quadratically with (weak) magnetic field and that the fluctuation-dissipation theorem is valid for our model.

Tailored interface and enhanced elastic modulus in epoxy-based composites in presence of branched poly(ethyleneimine) grafted multiwall carbon nanotubes

In this study, branched poly(ethyleneimine), BPEI, was synthesized from carboxylic acid terminated multi-walled carbon nanotubes (c-MWNTs) and characterized using FTIR, TEM and TGA. The BPEI was then chemically grafted onto MWNTs to enhance the interfacial adhesion with the epoxy matrix. The epoxy composites with c-MWNTs and the BPEI-g-MWNTs were prepared using a sonication and mechanical stirring method, followed by curing at 100 degrees C and post-curing at 120 degrees C. The dynamic mechanical thermal analysis showed an impressive 49% increment in the storage elastic modulus in the composites. In addition, the nanoindentation on the composites exhibited significant improvement in the hardness and decrease in the plasticity index in the presence of the BPEI-g-MWNTs. Thus, epoxy composites with BPEI-g-MWNTs can be further explored as self-healing materials.

COMPUTATIONAL MODELING OF A MULTI-LAYERED PIEZO-COMPOSITE BEAM MADE UP OF MFC

The Variational Asymptotic Method (VAM) is used for modeling a coupled non-linear electromechanical problem finding applications in aircrafts and Micro Aerial Vehicle (MAV) development. VAM coupled with geometrically exact kinematics forms a powerful tool for analyzing a complex nonlinear phenomena as shown previously by many in the literature 3 - 7] for various challenging problems like modeling of an initially twisted helicopter rotor blades, matrix crack propagation in a composite, modeling of hyper elastic plates and various multi-physics problems. The problem consists of design and analysis of a piezocomposite laminate applied with electrical voltage(s) which can induce direct and planar distributed shear stresses and strains in the structure. The deformations are large and conventional beam theories are inappropriate for the analysis. The behavior of an elastic body is completely understood by its energy. This energy must be integrated over the cross-sectional area to obtain the 1-D behavior as is typical in a beam analysis. VAM can be used efficiently to approximate 3-D strain energy as closely as possible. To perform this simplification, VAM makes use of thickness to width, width to length, width multiplied by initial twist and strain as small parameters embedded in the problem definition and provides a way to approach the exact solution asymptotically. In this work, above mentioned electromechanical problem is modeled using VAM which breaks down the 3-D elasticity problem into two parts, namely a 2-D non-linear cross-sectional analysis and a 1-D non-linear analysis, along the reference curve. The recovery relations obtained as a by-product in the cross-sectional analysis earlier are used to obtain 3-D stresses, displacements and velocity contours. The piezo-composite laminate which is chosen for an initial phase of computational modeling is made up of commercially available Macro Fiber Composites (MFCs) stacked together in an arbitrary lay-up and applied with electrical voltages for actuation. The expressions of sectional forces and moments as obtained from cross-sectional analysis in closed-form show the electro-mechanical coupling and relative contribution of electric field in individual layers of the piezo-composite laminate. The spatial and temporal constitutive law as obtained from the cross-sectional analysis are substituted into 1-D fully intrinsic, geometrically exact equilibrium equations of motion and 1-D intrinsic kinematical equations to solve for all 1-D generalized variables as function of time and an along the reference curve co-ordinate, x(1).

Temperature Dependence of Mechanical Properties in Molecular Crystals

Quantitative evaluation of the mechanical behavior of molecular materials by a nanoindentation technique has gained prominence recently. However, all the reported data have been on room-temperature properties despite many interesting phenomena observed in them with variations in temperature. In this paper, we report the results of nanoindentation experiments conducted as a function of temperature, T, between 283 and 343 K, on the major faces of three organic crystals: saccharin, sulfathiazole (form 2), and L-alanine, which are distinct in terms of the number and strength of intermolecular interactions in them. Results show that elastic modulus, E, and hardness, H, decrease markedly with increasing T. While E decreases linearly with T, the variations in H with T are not so, and were observed to drop by similar to 50% over the range of T investigated. The slope of the linear fits to E vs T for the organic crystals was found to be around 1, which is considerably higher than the values of 0.3-0.5 reported in the literature for metallic, ionic, and covalently bonded crystalline materials. Possible implications of the observed remarkable changes in H for pharmaceutical manufacturing are highlighted.

VAM Applied to Dimensional Reduction of Nonlinear Multifunctional Film-fabric Laminates

This work aims at asymptotically accurate dimensional reduction of non-linear multi-functional film-fabric laminates having specific application in design of envelopes for High Altitude Airships (HAA). The film-fabric laminate for airship envelope consists of a woven fabric core coated with thin films on each face. These films provide UV protection and Helium leakage prevention, while the core provides required structural strength. This problem is both geometrically and materially non-linear. To incorporate the geometric non-linearity, generalized warping functions are used and finite deformations are allowed. The material non-linearity is handled by using hyper-elastic material models for each layer. The development begins with three-dimensional (3-D) nonlinear elasticity and mathematically splits the analysis into a one-dimensional through-the-thickness analysis and a two-dimensional (2-D) plate analysis. The through-the-thickness analysis provides the 2-D constitutive law which is then given as an input to the 2-D reference surface analysis. The dimensional reduction is carried out using Variational Asymptotic Method (VAM) for moderate strains and very small thickness-to-wavelength ratio. It features the identification and utilization of additional small parameters such as ratio of thicknesses and stiffness coefficients of core and films. Closed form analytical expressions for warping functions and 2-D constitutive law of the film-fabric laminate are obtained.

A Polymatroid Approach to Separable Convex Optimization with Linear Ascending Constraints

We revisit a problem studied by Padakandla and Sundaresan SIAM J. Optim., August 2009] on the minimization of a separable convex function subject to linear ascending constraints. The problem arises as the core optimization in several resource allocation problems in wireless communication settings. It is also a special case of an optimization of a separable convex function over the bases of a specially structured polymatroid. We give an alternative proof of the correctness of the algorithm of Padakandla and Sundaresan. In the process we relax some of their restrictions placed on the objective function.

Alternative Formulation of the Discrete Linear Quadratic Regulator (DLQR)

In this paper, an alternative apriori and aposteriori formulation has been derived for the discrete linear quadratic regulator (DLQR) in a manner analogous to that used in the discrete Kalman filter. It has been shown that the formulation seamlessly fits into the available formulation of the DLQR and the equivalent terms in the existing formulation and the proposed formulation have been identified. Thereafter, the significance of this alternative formulation has been interpreted in terms of the sensitivity of the controller performances to any changes in the states or to changes in the control inputs. The implications of this alternative formulation to adaptive controller tuning have also been discussed.

Anisotropy in the mechanical properties of organic crystals: temperature dependence

The nanoindentation technique has recently been utilized for quantitative evaluation of the mechanical properties of molecular materials successfully, including their temperature (T) dependence. In this paper, we examine how the mechanical anisotropy varies with T in saccharin and L-alanine single crystals. Our results show that elastic modulus (E) decreases linearly in all the cases examined, with the T-dependence of E being anisotropic. Correspondence between directional dependence of the slopes of the E vs. T plots and the linear thermal expansion coefficients was found. The T-dependence of hardness (H), on the other hand, was found to be nonlinear and significant when (100) of saccharin and (001) of L-alanine are indented. While the anisotropies in E and H of saccharin and E of L-alanine enhance with T, the anisotropy in H of L-alanine was found to reduce with T. Possible mechanistic origins of these variations are discussed.

India-Madagascar paleo-fit based on flexural isostasy of their rifted margins

The present study contributes new constraints on, and definitions of, the reconstructed plate margins of India and Madagascar based on flexural isostasy along the Western Continental Margin of India (WCMI) and the Eastern Continental Margin of Madagascar (ECMM). We have estimated the nature of isostasy and crustal geometry along the two margins, and have examined their possible conjugate structure. Here we utilize elastic thickness (Te) and Moho depth data as the primary basis for the correlation of these passive margins. We employ the flexure inversion technique that operates in spatial domain in order to estimate the spatial variation of effective elastic thickness. Gravity inversion and flexure inversion techniques are used to estimate the configuration of the Moho/Crust-Mantle Interface that reveals regional correlations with the elastic thickness variations. These results correlate well with the continental and oceanic segments of the Indian and African plates. The present study has found a linear zone of anomalously low-Te (1-5 km) along the WCMI (similar to 1680 km), which correlates well with the low-Te patterns obtained all along the ECMM. We suggest that the low-Te zones along the WCMI and ECMM represent paleo-rift inception points of lithosphere thermally and mechanically weakened by the combined effects of the Marion hotspot and lithospheric extension due to rifting. We have produced an India-Madagascar paleo-fit representing the initial phase of separation based on the Te estimates of the rifted conjugate margins, which is confirmed by a close-fit correlation of Moho geometry and bathymetry of the shelf margins. The matching of tectonic lineaments, lithologies and geochronological belts between India and Madagascar provide an additional support for the present plate reconstruction. (C) 2014 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

On the Bounds of Certain Maximal Linear Codes in a Projective Space

The set of all subspaces of F-q(n) is denoted by P-q(n). The subspace distance d(S)(X, Y) = dim(X) + dim(Y)-2dim(X boolean AND Y) defined on P-q(n) turns it into a natural coding space for error correction in random network coding. A subset of P-q(n) is called a code and the subspaces that belong to the code are called codewords. Motivated by classical coding theory, a linear coding structure can be imposed on a subset of P-q(n). Braun et al. conjectured that the largest cardinality of a linear code, that contains F-q(n), is 2(n). In this paper, we prove this conjecture and characterize the maximal linear codes that contain F-q(n).

In-plane elastic responses of circular cell honeycombs and bundled circular tubes in a diamond array structure

Assemblages of circular tubes and circular honeycombs in close packed arrangement are presently both competing and complementing regular honeycomb structures (HCS). The intrinsic isotropy of bundled tubes/rings in hexagonal arrays restricts their use to applications with isotopic need. With the aim of extending the utility of tubes/rings assemblages to anisotropic needs, this paper explores the prospects of bundled tubes and circular honeycombs in a general diamond array structure (DAS) to cater these needs. To this end, effective transverse Young's moduli and Poisson's ratio for thick/thin DAS are obtained theoretically. Analysis frameworks including thin ring theory (TRT), curved beam theory (CBT) and elasticity formulations are tested and corroborated by FEA employing contact elements. Results indicate that TRT and CBT are reasonable for thin tubes and honeycombs. Nevertheless, TRT yields compact formulae to study the anisotropy ratio, moduli spectrum and sensitivity of the assemblage as a function of thicknesses and array structure. These formulae supplement designers as a guide to tailor the structures. On the other hand, elasticity formulation can estimate over a larger range including very thick tubes/rings. In addition, this formulation offers to estimate refined transverse strengths of assemblages. (C) 2015 Elsevier Ltd. All rights reserved.

A fast algorithm for speech polarity detection using long-term linear prediction

Speech polarity detection is a crucial first step in many speech processing techniques. In this paper, an algorithm is proposed that improvises the existing technique using the skewness of the voice source (VS) signal. Here, the integrated linear prediction residual (ILPR) is used as the VS estimate, which is obtained using linear prediction on long-term frames of the low-pass filtered speech signal. This excludes the unvoiced regions from analysis and also reduces the computation. Further, a modified skewness measure is proposed for decision, which also considers the magnitude of the skewness of the ILPR along with its sign. With the detection error rate (DER) as the performance metric, the algorithm is tested on 8 large databases and its performance (DER=0.20%) is found to be comparable to that of the best technique (DER=0.06%) on both clean and noisy speech. Further, the proposed method is found to be ten times faster than the best technique.

Wave propagation in a 2D nonlinear structural acoustic waveguide using asymptotic expansions of wavenumbers

Nonlinear acoustic wave propagation in an infinite rectangular waveguide is investigated. The upper boundary of this waveguide is a nonlinear elastic plate, whereas the lower boundary is rigid. The fluid is assumed to be inviscid with zero mean flow. The focus is restricted to non-planar modes having finite amplitudes. The approximate solution to the acoustic velocity potential of an amplitude modulated pulse is found using the method of multiple scales (MMS) involving both space and time. The calculations are presented up to the third order of the small parameter. It is found that at some frequencies the amplitude modulation is governed by the Nonlinear Schrodinger equation (NLSE). The first objective here is to study the nonlinear term in the NLSE. The sign of the nonlinear term in the NLSE plays a role in determining the stability of the amplitude modulation. Secondly, at other frequencies, the primary pulse interacts with its higher harmonics, as do two or more primary pulses with their resultant higher harmonics. This happens when the phase speeds of the waves match and the objective is to identify the frequencies of such interactions. For both the objectives, asymptotic coupled wavenumber expansions for the linear dispersion relation are required for an intermediate fluid loading. The novelty of this work lies in obtaining the asymptotic expansions and using them for predicting the sign change of the nonlinear term at various frequencies. It is found that when the coupled wavenumbers approach the uncoupled pressure-release wavenumbers, the amplitude modulation is stable. On the other hand, near the rigid-duct wavenumbers, the amplitude modulation is unstable. Also, as a further contribution, these wavenumber expansions are used to identify the frequencies of the higher harmonic interactions. And lastly, the solution for the amplitude modulation derived through the MMS is validated using these asymptotic expansions. (C) 2015 Elsevier Ltd. All rights reserved.