A Smooth Discretization Bridging Finite Element and Mesh-free Methods Using Polynomial Reproducing Simplex Splines

This work sets forth a `hybrid' discretization scheme utilizing bivariate simplex splines as kernels in a polynomial reproducing scheme constructed over a conventional Finite Element Method (FEM)-like domain discretization based on Delaunay triangulation. Careful construction of the simplex spline knotset ensures the success of the polynomial reproduction procedure at all points in the domain of interest, a significant advancement over its precursor, the DMS-FEM. The shape functions in the proposed method inherit the global continuity (Cp-1) and local supports of the simplex splines of degree p. In the proposed scheme, the triangles comprising the domain discretization also serve as background cells for numerical integration which here are near-aligned to the supports of the shape functions (and their intersections), thus considerably ameliorating an oft-cited source of inaccuracy in the numerical integration of mesh-free (MF) schemes. Numerical experiments show the proposed method requires lower order quadrature rules for accurate evaluation of integrals in the Galerkin weak form. Numerical demonstrations of optimal convergence rates for a few test cases are given and the method is also implemented to compute crack-tip fields in a gradient-enhanced elasticity model.

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.

Mode coupling theory analysis of electrolyte solutions: Time dependent diffusion, intermediate scattering function, and ion solvation dynamics

A self-consistent mode coupling theory (MCT) with microscopic inputs of equilibrium pair correlation functions is developed to analyze electrolyte dynamics. We apply the theory to calculate concentration dependence of (i) time dependent ion diffusion, (ii) intermediate scattering function of the constituent ions, and (iii) ion solvation dynamics in electrolyte solution. Brownian dynamics with implicit water molecules and molecular dynamics method with explicit water are used to check the theoretical predictions. The time dependence of ionic self-diffusion coefficient and the corresponding intermediate scattering function evaluated from our MCT approach show quantitative agreement with early experimental and present Brownian dynamic simulation results. With increasing concentration, the dispersion of electrolyte friction is found to occur at increasingly higher frequency, due to the faster relaxation of the ion atmosphere. The wave number dependence of intermediate scattering function, F(k, t), exhibits markedly different relaxation dynamics at different length scales. At small wave numbers, we find the emergence of a step-like relaxation, indicating the presence of both fast and slow time scales in the system. Such behavior allows an intriguing analogy with temperature dependent relaxation dynamics of supercooled liquids. We find that solvation dynamics of a tagged ion exhibits a power law decay at long times-the decay can also be fitted to a stretched exponential form. The emergence of the power law in solvation dynamics has been tested by carrying out long Brownian dynamics simulations with varying ionic concentrations. The solvation time correlation and ion-ion intermediate scattering function indeed exhibit highly interesting, non-trivial dynamical behavior at intermediate to longer times that require further experimental and theoretical studies. (c) 2015 AIP Publishing LLC.

A spectral element for wave propagation in honeycomb sandwich construction considering core flexibility

Spectral elements are found to be extremely resourceful to study the wave propagation characteristics of structures at high frequencies. Most of the aerospace structures use honeycomb sandwich constructions. The existing spectral elements use single layer theories for a sandwich construction wherein the two face sheets vibrate together and this model is sufficient for low frequency excitations. At high frequencies, the two face sheets vibrate independently. The Extended Higher order SAndwich Plate theory (EHSaPT) is suitable for representing the independent motion of the face sheets. A 1D spectral element based on EHSaPT is developed in this work. The wave number and the wave speed characteristics are obtained using the developed spectral element. It is shown that the developed spectral element is capable of representing independent wave motions of the face sheets. The propagation speeds of a high frequency modulated pulse in the face sheets and the core of a honeycomb sandwich are demonstrated. Responses of a typical honeycomb sandwich beam to high frequency shock loads are obtained using the developed spectral element and the response match very well with the finite element results. It is shown that the developed spectral element is able to represent the flexibility of the core resulting into independent wave motions in the face sheets, for which a finite element method needs huge degrees of freedom. (C) 2015 Elsevier Ltd. All rights reserved.

Multi-transform based spectral element to include first order shear deformation in plates

For obtaining dynamic response of structure to high frequency shock excitation spectral elements have several advantages over conventional methods. At higher frequencies transverse shear and rotary inertia have a predominant role. These are represented by the First order Shear Deformation Theory (FSDT). But not much work is reported on spectral elements with FSDT. This work presents a new spectral element based on the FSDT/Mindlin Plate Theory which is essential for wave propagation analysis of sandwich plates. Multi-transformation method is used to solve the coupled partial differential equations, i.e., Laplace transforms for temporal approximation and wavelet transforms for spatial approximation. The formulation takes into account the axial-flexure and shear coupling. The ability of the element to represent different modes of wave motion is demonstrated. Impact on the derived wave motion characteristics in the absence of the developed spectral element is discussed. The transient response using the formulated element is validated by the results obtained using Finite Element Method (FEM) which needs significant computational effort. Experimental results are provided which confirms the need to having the developed spectral element for the high frequency response of structures. (C) 2015 Elsevier Ltd. All rights reserved.

A Reliable Residual Based A Posteriori Error Estimator for a Quadratic Finite Element Method for the Elliptic Obstacle Problem

A residual based a posteriori error estimator is derived for a quadratic finite element method (FEM) for the elliptic obstacle problem. The error estimator involves various residuals consisting of the data of the problem, discrete solution and a Lagrange multiplier related to the obstacle constraint. The choice of the discrete Lagrange multiplier yields an error estimator that is comparable with the error estimator in the case of linear FEM. Further, an a priori error estimate is derived to show that the discrete Lagrange multiplier converges at the same rate as that of the discrete solution of the obstacle problem. The numerical experiments of adaptive FEM show optimal order convergence. This demonstrates that the quadratic FEM for obstacle problem exhibits optimal performance.

Effect of a mesh on boundary layer transitions induced by free-stream turbulence and an isolated roughness element

Streamwise streaks, their lift-up and streak instability are integral to the bypass transition process. An experimental study has been carried out to find the effect of a mesh placed normal to the flow and at different wall-normal locations in the late stages of two transitional flows induced by free-stream turbulence (FST) and an isolated roughness element. The mesh causes an approximately 30% reduction in the free-stream velocity, and mild acceleration, irrespective of its wall-normal location. Interestingly, when located near the wall, the mesh suppresses several transitional events leading to transition delay over a large downstream distance. The transition delay is found to be mainly caused by suppression of the lift-up of the high-shear layer and its distortion, along with modification of the spanwise streaky structure to an orderly one. However, with the mesh well away from the wall, the lifted-up shear layer remains largely unaffected, and the downstream boundary layer velocity profile develops an overshoot which is found to follow a plane mixing layer type profile up to the free stream. Reynolds stresses, and the size and strength of vortices increase in this mixing layer region. This high-intensity disturbance can possibly enhance transition of the accelerated flow far downstream, although a reduction in streamwise turbulence intensity occurs over a short distance downstream of the mesh. However, the shape of the large-scale streamwise structure in the wall-normal plane is found to be more or less the same as that without the mesh.

Multiobjective fuzzy optimization for sustainable groundwater management using particle swarm optimization and analytic element method

Groundwater management involves conflicting objectives as maximization of discharge contradicts the criteria of minimum pumping cost and minimum piping cost. In addition, available data contains uncertainties such as market fluctuations, variations in water levels of wells and variations of ground water policies. A fuzzy model is to be evolved to tackle the uncertainties, and a multiobjective optimization is to be conducted to simultaneously satisfy the contradicting objectives. Towards this end, a multiobjective fuzzy optimization model is evolved. To get at the upper and lower bounds of the individual objectives, particle Swarm optimization (PSO) is adopted. The analytic element method (AEM) is employed to obtain the operating potentio metric head. In this study, a multiobjective fuzzy optimization model considering three conflicting objectives is developed using PSO and AEM methods for obtaining a sustainable groundwater management policy. The developed model is applied to a case study, and it is demonstrated that the compromise solution satisfies all the objectives with adequate levels of satisfaction. Sensitivity analysis is carried out by varying the parameters, and it is shown that the effect of any such variation is quite significant. Copyright (c) 2015 John Wiley & Sons, Ltd.

Ab Initio Molecular Dynamics Simulations of Amino Acids in Aqueous Solutions: Estimating pK(a) Values from Metadynamics Sampling

Changes in the protonation and deprotonation of amino acid residues in proteins play a key role in many biological processes and pathways. Here, we report calculations of the free-energy profile for the protonation deprotonation reaction of the 20 canonical alpha amino acids in aqueous solutions using ab initio Car-Parrinello molecular dynamics simulations coupled with metad-ynamics sampling. We show here that the calculated change in free energy of the dissociation reaction provides estimates of the multiple pK(a) values of the amino acids that are in good agreement with experiment. We use the bond-length-dependent number of the protons coordinated to the hydroxyl oxygen of the carboxylic and the amine groups as the collective variables to explore the free-energy profiles of the Bronsted acid-base chemistry of amino acids in aqueous solutions. We ensure that the amino acid undergoing dissociation is solvated by at least three hydrations shells with all water molecules included in the simulations. The method works equally well for amino acids with neutral, acidic and basic side chains and provides estimates of the multiple pK(a) values with a mean relative error, with respect to experimental results, of 0.2 pK(a) units.

Bearing capacity of circular footings over rock mass by using axisymmetric quasi lower bound finite element limit analysis

The ultimate bearing capacity of a circular footing, placed over rock mass, is evaluated by using the lower bound theorem of the limit analysis in conjunction with finite elements and nonlinear optimization. The generalized Hoek-Brown (HB) failure criterion, but by keeping a constant value of the exponent, alpha = 0.5, was used. The failure criterion was smoothened both in the meridian and pi planes. The nonlinear optimization was carried out by employing an interior point method based on the logarithmic barrier function. The results for the obtained bearing capacity were presented in a non-dimensional form for different values of GSI, m(i), sigma(ci)/(gamma b) and q/sigma(ci). Failure patterns were also examined for a few cases. For validating the results, computations were also performed for a strip footing as well. The results obtained from the analysis compare well with the data reported in literature. Since the equilibrium conditions are precisely satisfied only at the centroids of the elements, not everywhere in the domain, the obtained lower bound solution will be approximate not true. (C) 2015 Elsevier Ltd. All rights reserved.

Stability analysis of thick piezoelectric metal based FGM plate using first order and higher order shear deformation theory

This paper presents the stability analysis of functionally graded plate integrated with piezoelectric actuator and sensor at the top and bottom face, subjected to electrical and mechanical loading. The finite element formulation is based on first order and higher order shear deformation theory, degenerated shell element, von-Karman hypothesis and piezoelectric effect. The equation for static analysis is derived by using the minimum energy principle and solutions for critical buckling load is obtained by solving eigenvalue problem. The material properties of the functionally graded plate are assumed to be graded along the thickness direction according to simple power law function. Two types of boundary conditions are used, such as SSSS (simply supported) and CSCS (simply supported along two opposite side perpendicular to the direction of compression and clamped along the other two sides). Sensor voltage is calculated using present analysis for various power law indices and FG (functionally graded) material gradations. The stability analysis of piezoelectric FG plate is carried out to present the effects of power law index, material variations, applied mechanical pressure and piezo effect on buckling and stability characteristics of FG plate.

A FRAMEWORK FOR THE ERROR ANALYSIS OF DISCONTINUOUS FINITE ELEMENT METHODS FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS AND APPLICATIONS TO C-0 IP METHODS

In this article, an abstract framework for the error analysis of discontinuous Galerkin methods for control constrained optimal control problems is developed. The analysis establishes the best approximation result from a priori analysis point of view and delivers a reliable and efficient a posteriori error estimator. The results are applicable to a variety of problems just under the minimal regularity possessed by the well-posedness of the problem. Subsequently, the applications of C-0 interior penalty methods for a boundary control problem as well as a distributed control problem governed by the biharmonic equation subject to simply supported boundary conditions are discussed through the abstract analysis. Numerical experiments illustrate the theoretical findings.

Ab Initio MD Simulations of the Bronsted Acidity of Glutathione in Aqueous Solutions: Predicting pK(a) Shifts of the Cysteine Residue

The tripeptide glutathione (GSH) is one of the most abundant peptides and the major repository for nonprotein sulfur in both animal and plant cells. It plays a critical role in intracellular oxidative stress management by the reversible formation of glutathione disulfide with the thiol-disulfide pair acting as a redox buffer. The state of charge of the ionizable groups of GSH can influence the redox couple, and hence the pK(a) value of the cysteine residue of GSH is critical to its functioning. Here we report ab initio Car-Parrinello molecular dynamics simulations of glutathione solvated by 200 water molecules, all of which are considered in the simulation. We show that the free-energy landscape for the protonation-deprotonation reaction of the cysteine residue of GSH computed using metadynamics sampling provides shift in the dissociation constant values as compared with the isolated accurate estimates of the pK(a) and correctly predicts the cysteine amino acid.

Spectral solutions to the Korteweg-de-Vries and nonlinear Schrodinger equations

In this paper, we present the solutions of 1-D and 2-D non-linear partial differential equations with initial conditions. We approach the solutions in time domain using two methods. We first solve the equations using Fourier spectral approximation in the spatial domain and secondly we compare the results with the approximation in the spatial domain using orthogonal functions such as Legendre or Chebyshev polynomials as their basis functions. The advantages and the applicability of the two different methods for different types of problems are brought out by considering 1-D and 2-D nonlinear partial differential equations namely the Korteweg-de-Vries and nonlinear Schrodinger equation with different potential function. (C) 2015 Elsevier Ltd. All rights reserved.