An interior penalty method for a sixth-order elliptic equation

We derive and study a C(0) interior penalty method for a sixth-order elliptic equation on polygonal domains. The method uses the cubic Lagrange finite-element space, which is simple to implement and is readily available in commercial software. After introducing some notation and preliminary results, we provide a detailed derivation of the method. We then prove the well-posedness of the method as well as derive quasi-optimal error estimates in the energy norm. The proof is based on replacing Galerkin orthogonality with a posteriori analysis techniques. Using this approach, we are able to obtain a Cea-like lemma with minimal regularity assumptions on the solution. Numerical experiments are presented that support the theoretical findings.

A QUADRATIC C degrees INTERIOR PENALTY METHOD FOR LINEAR FOURTH ORDER BOUNDARY VALUE PROBLEMS WITH BOUNDARY CONDITIONS OF THE CAHN-HILLIARD TYPE

We develop a quadratic C degrees interior penalty method for linear fourth order boundary value problems with essential and natural boundary conditions of the Cahn-Hilliard type. Both a priori and a posteriori error estimates are derived. The performance of the method is illustrated by numerical experiments.

Analysis of an Interior Penalty Method for Fourth Order Problems on Polygonal Domains

Error analysis for a stable C (0) interior penalty method is derived for general fourth order problems on polygonal domains under minimal regularity assumptions on the exact solution. We prove that this method exhibits quasi-optimal order of convergence in the discrete H (2), H (1) and L (2) norms. L (a) norm error estimates are also discussed. Theoretical results are demonstrated by numerical experiments.

A C-0 Interior Penalty Method for a Fourth-order Variational Inequality of the Second Kind

In this article, we propose a C-0 interior penalty ((CIP)-I-0) method for the frictional plate contact problem and derive both a priori and a posteriori error estimates. We derive an abstract error estimate in the energy norm without additional regularity assumption on the exact solution. The a priori error estimate is of optimal order whenever the solution is regular. Further, we derive a reliable and efficient a posteriori error estimator. Numerical experiments are presented to illustrate the theoretical results. (c) 2015Wiley Periodicals, Inc.

Shock propagation in gas dynamics: explicit form of higher order compatibility conditions

With the use of tensor analysis and the method of singular surfaces, an infinite system of equations can be derived to study the propagation of curved shocks of arbitrary strength in gas dynamics. The first three of these have been explicitly given here. This system is further reduced to one involving scalars only. The choice of dependent variables in the infinite system is quite important, it leads to coefficients free from singularities for all values of the shock strength.

Unsteady laminar compressible stagnation-point boundary-layer flow over a three-dimensional body: Higher order effects

The unsteady laminar compressible three-dimensional stagnation-point boundary-layer flow with variable properties has been studied when the velocity of the incident stream, mass transfer and wall temperature vary arbitrarily with time. The second-order unsteady boundary-layer equations for all the effects have been derived by using the method of matched asymptotic expansions. Both nodal and saddle point flows as well as cold and hot wall cases have been considered. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. Computations have been carried out for an accelerating stream, a decelerating stream and a fluctuating stream. The results indicate that the unsteady free stream velocity distributions, the nature of the stagnation point, the mass transfer, the wall temperature and the variation of the density-viscosity product across the boundary significantly affect the skin friction and heat transfer. The variation of the wall temperature with time strongly affects the heat transfer whereas its effect is comparatively less on skin friction. Suction increases the skin friction and heat transfer but injection does the opposite. The skin friction in the x direction due to the combined effects of first- and second-order boundary layers is less than the skin-friction in the x direction due to the first-order boundary layers for all the parameters. The overall skin friction in the z direction and heat transfer are more or less than the first-order boundary layers depending upon the values of the various parameters.

Energy identities in water wave theory for free-surface boundary condition with higher-order derivatives

A modified form of Green's integral theorem is employed to derive the energy identity in any water wave diffraction problem in a single-layer fluid for free-surface boundary condition with higher-order derivatives. For a two-layer fluid with free-surface boundary condition involving higher-order derivatives, two forms of energy identities involving transmission and reflection coefficients for any wave diffraction problem are also derived here by the same method. Based on this modified Green's theorem, hydrodynamic relations such as the energy-conservation principle and modified Haskind–Hanaoka relation are derived for radiation and diffraction problems in a single as well as two-layer fluid.

On the Role of Higher-Order Evanescent Modes in End-Offset Inlet and End-Centered Outlet Elliptical Flow-Reversal Chamber Mufflers

This paper deals with the role of the higher-order evanescent modes generated at the area discontinuities in the acoustic attenuation characteristics of an elliptical end-chamber muffler with an end-offset inlet and end-centered outlet. It has been observed that with an increase in length, the muffler undergoes a transition from being acoustically short to acoustically long. Short end chambers and long end chambers are characterized by transverse plane waves and axial plane waves, respectively, in the low-frequency range. The nondimensional frequency limit k(0)(D-1/2) or k(0)R(0) as well as the chamber length to inlet/outlet pipe diameter ratio, i.e., L/d(0), up to which the muffler behaves like a short chamber and the corresponding limit beyond which the muffler is acoustically long are determined. The limits between which neither the transverse plane-wave model nor the conventional axial plane-wave model gives a satisfactory prediction have also been determined, the region being called the intermediate range. The end-correction expression for this muffler configuration in the acoustically long limit has been obtained using 3-D FEA carried on commercial software, covering most of the dimension range used in the design exercise. Development of a method of combining the transverse plane wave model with the axial plane wave model using the impedance Z] matrix is another noteworthy contribution of this work.

Integrated higher-order pulse-width modulation filter-transformer structure for single-phase static compensator

A power filter is necessary to connect the output of a power converter to the grid so as to reduce the harmonic distortion introduced in the line current and voltage by the power converter. Many a times, a transformer is also present before the point of common coupling. Magnetic components often constitute a significant part of the overall weight, size and cost of the grid interface scheme. So, a compact inexpensive design is desirable. A higher-order LCL-filter and a transformer are increasingly being considered for grid interconnection of the power converter. This study proposes a design method based on a three-winding transformer, that generates an integrated structure that behaves as an LCL-filter, with both the filter inductances and the transformer that are merged into a single electromagnetic component. The parameters of the transformer are derived analytically. It is shown that along with a filter capacitor, the transformer parameters provide the filtering action of an LCL-filter. A single-phase full-bridge power converter is operated as a static compensator for performance evaluation of the integrated filter transformer. A resonant integrator-based single-phase phase locked loop and stationary frame AC current controller are employed for grid frequency synchronisation and line current control, respectively.

Automatic finite element formulation and assembly of hyperelastic higher order structural models

The formulation of higher order structural models and their discretization using the finite element method is difficult owing to their complexity, especially in the presence of non-linearities. In this work a new algorithm for automating the formulation and assembly of hyperelastic higher-order structural finite elements is developed. A hierarchic series of kinematic models is proposed for modeling structures with special geometries and the algorithm is formulated to automate the study of this class of higher order structural models. The algorithm developed in this work sidesteps the need for an explicit derivation of the governing equations for the individual kinematic modes. Using a novel procedure involving a nodal degree-of-freedom based automatic assembly algorithm, automatic differentiation and higher dimensional quadrature, the relevant finite element matrices are directly computed from the variational statement of elasticity and the higher order kinematic model. Another significant feature of the proposed algorithm is that natural boundary conditions are implicitly handled for arbitrary higher order kinematic models. The validity algorithm is illustrated with examples involving linear elasticity and hyperelasticity. (C) 2013 Elsevier Inc. All rights reserved.

Higher-Order Nonlinearity In Accretion Disks: Quasi-Periodic Oscillations Of Black Hole And Neutron Star Sources And Their Spin

We propose a unified model to explain Quasi-Periodic Oscillation (QPO), particularly of high frequency, observed from black hole and neutron star systems globally. We consider accreting systems to be damped harmonic oscillators exhibiting epicyclic oscillations with higher-order nonlinear resonance to explain QPO. The resonance is expected to be driven by the disturbance from the compact object at its spin frequency. The model explains various properties parallelly for both types of the compact object. It describes QPOs successfully for ten different compact sources. Based on this, we predict the spin frequency of the neutron star Sco X-1 and specific angular momentum of black holes GRO J1655–40, XTE J1550–564, H1743–322, and GRS 1915+105.

Fracture Behavior of Multidirectional DCB Specimen: Higher-Order Beam Theories

Mathematical models, for the stress analysis of symmetric multidirectional double cantilever beam (DCB) specimen using classical beam theory, first and higher-order shear deformation beam theories, have been developed to determine the Mode I strain energy release rate (SERR) for symmetric multidirectional composites. The SERR has been calculated using the compliance approach. In the present study, both variationally and nonvariationally derived matching conditions have been applied at the crack tip of DCB specimen. For the unidirectional and cross-ply composite DCB specimens, beam models under both plane stress and plane strain conditions in the width direction are applicable with good performance where as for the multidirectional composite DCB specimen, only the beam model under plane strain condition in the width direction appears to be applicable with moderate performance. Among the shear deformation beam theories considered, the performance of higher-order shear deformation beam theory, having quadratic variation for transverse displacement over the thickness, is superior in determining the SERR for multidirectional DCB specimen.

On higher order shear deformation theory of laminated composite panels

In this paper we examine the suitability of higher order shear deformation theory based on cubic inplane displacements and parabolic normal displacements, for stress analysis of laminated composite plates including the interlaminar stresses. An exact solution of a symmetrical four layered infinite strip under static loading has been worked out and the results obtained by the present theory are compared with the exact solution. The present theory provides very good estimates of the deflections, and the inplane stresses and strains. Nevertheless, direct estimates of strains and stresses do not display the required interlaminar stress continuity and strain discontinuity across the interlaminar surface. On the other hand, ‘statically equivalent stresses and strains’ do display the required interlaminar stress continuity and strain discontinuity and agree very closely with the exact solution.

A simple method for potential flow simulation of cascades

A simple method using a combination of conformal mapping and vortex panel method to simulate potential flow in cascades is presented. The cascade is first transformed to a single body using a conformal mapping, and the potential flow over this body is solved using a simple higher order vortex panel method. The advantage of this method over existing methodologies is that it enables the use of higher order panel methods, as are used to solve flow past an isolated airfoil, to solve the cascade problem without the need for any numerical integrations or iterations. The fluid loading on the blades, such as the normal force and pitching moment, may be easily calculated from the resultant velocity field. The coefficient of pressure on cascade blades calculated with this methodology shows good agreement with previous numerical and experimental results.

Higher twist and higher order contributions to the pion electromagnetic form factor

The O(m(pi)4/(m(u) + (d))2Q2) and O(alpha(S)2) corrections to the leading term of the perturbative QCD calculation of the pion electromagnetic form factor are examined numerically. Both sets of terms provide significant corrections for values of Q2 between 1 and 15 GeV2/c2.