Integrating strong and weak discontinuities without integration subcells and example applications in an XFEM/GFEM framework

Partition of unity methods, such as the extended finite element method, allows discontinuities to be simulated independently of the mesh (Int. J. Numer. Meth. Engng. 1999; 45:601-620). This eliminates the need for the mesh to be aligned with the discontinuity or cumbersome re-meshing, as the discontinuity evolves. However, to compute the stiffness matrix of the elements intersected by the discontinuity, a subdivision of the elements into quadrature subcells aligned with the discontinuity is commonly adopted. In this paper, we use a simple integration technique, proposed for polygonal domains (Int. J. Nuttier Meth. Engng 2009; 80(1):103-134. DOI: 10.1002/nme.2589) to suppress the need for element subdivision. Numerical results presented for a few benchmark problems in the context of linear elastic fracture mechanics and a multi-material problem show that the proposed method yields accurate results. Owing to its simplicity, the proposed integration technique can be easily integrated in any existing code. Copyright (C) 2010 John Wiley & Sons, Ltd.

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.

Temperature response of trehalase from a mesophilic (Neurospora crassa) and a thermophilic (Thermomyces lanuginosus) fungus

The purified trehalases of the mesophilic fungus, Neurospora crassa, and the thermophilic fungus, Thermomyces lanuginosus, had similar temperature and pH optima for activity, but differed in molecular weight, electrophoretic mobility and Michaelis constant. At lower concentration, trehalases from both fungi were inactivated to similar extent at 60°C. While purified trehalase of T. lanuginosus was afforded protection against heat-inactivation by proteinaceous protective factor(s) present in mycelial extracts, by bovine serum albumin and by casein, these did not afford protection to N. crassa trehalase against heat inactivation. Both trehalases exhibited discontinuous Arrhenius plots with temperature of discontinuity at 40°C. The activation energy calculated from the slope of the Arrhenius plot was higher for the T. lanuginosus enzyme. The plots of apparent K m versus 1/T for trehalases of N. crassa and T. lanuginosus were linear from 30° to 60°C. The results show that purified trehalases of the mesophilic and the thermophilic fungus are distinct. Although, these exhibit similar thermostability of their catalytic function at low concentration, distinctive thermal stability characteristics of thermophilic enzyme become apparent at high protein concentration. This could be brought about in the cell by the enzyme itself, or by other proteins.

A Study of copper-compounds by x-Ray Absorption-spectroscopy

Chemical shifts, ΔE, of the K-absorption discontinuity in several compounds of copper possessing formal oxidation states between 0 and III have been measured. The shifts show a parabolic dependence on the formal oxidation state as well as on the effective atomic charge, q, on copper. Anomalous chemical shifts shown by some of the compounds are discussed in terms of the bonding in these compounds. The ΔE values have also been correlated with the core electron binding energies obtained from X-ray photoelectron spectroscopy.

High pressure study of Fe3O4 through the verwey transition

The thermopower (TEP) and electrical resistance of stoichiometric Fe3O4 crystals have been measured up to pressures of 6 GPa over the temperature range of 80-160 K. The resistance decreases markedly with increasing pressure below the Verwey transition temperature TV and TV decreases linearly with increasing pressure. The magnitude of the TEP as well as the discontinuity at TV decrease with increasing pressure. The thermopower of Fe3O4 shows an interesting upswing at low temperatures (lt;100 K) which is affected significantly by pressure.

Non-Fermi liquid behavior in a mean-field study of the t-J model: Role of zero-point spin fluctuations

We present analytic results to show that the Schwinger-boson hole-fermion mean-field state exhibits non-Fermi liquid behavior due to spin-charge separation. The physical electron Green's function consists of three additive components. (a) A Fermi-liquid component associated with the bose condensate. (b) A non-Fermi liquid component which has a logarithmic peak and a long tail that gives rise to a linear density of states that is symmetric about the Fermi level and a momentum distribution function with a logarithmic discontinuity at the Fermi surface. (c) A second non-Fermi liquid component associated with the thermal bosons which leads to a constant density of states. It is shown that zero-point fluctuations associated with the spin-degrees of freedom are responsible for the logarithmic instabilities and the restoration of particle-hole symmetry close to the Fermi surface.

Coherent phonons in pyrochlore titanates A(2)Ti(2)O(7) (A = Dy, Gd, Tb): A phase transition in Dy2Ti2O7 at 110 K

We study the generation of coherent optical phonons in spin-frustrated pyrochlore single crystals Dy2Ti2O7, Gd2Ti2O7, and Tb2Ti2O7 using femtosecond laser pulses (65 fs, 1.57 eV) in degenerate time-resolved transmission experiments as a function of temperature from 4 to 296 K. At 4 K, two coherent phonons are observed at similar to 5.3 THz (5.0 THz) and similar to 9.3 THz (9.4 THz) for Dy2Ti2O7 (Gd2Ti2O7), whereas three coherent phonons are generated at similar to 5.0, 8.6, and 9.7 THz for Tb2Ti2O7. In the case of spin-ice Dy2Ti2O7, a clear discontinuity is observed in the linewidths of both the coherent phonons as well as in the phase of lower-energy coherent phonon mode, indicating a subtle structural change at 110 K. Another important observation is a phase difference of pi between the modes in all the samples, thus suggesting that the driving forces behind the generation of these modes could be different in nature, unlike a purely impulsive or displacive mechanism.

Stability of the viscous flow of a fluid through a flexible tube

The stability of Hagen-Poiseuille flow of a Newtonian fluid of viscosity eta in a tube of radius R surrounded by a viscoelastic medium of elasticity G and viscosity eta(s) occupying the annulus R < r < HR is determined using a linear stability analysis. The inertia of the fluid and the medium are neglected, and the mass and momentum conservation equations for the fluid and wall are linear. The only coupling between the mean flow and fluctuations enters via an additional term in the boundary condition for the tangential velocity at the interface, due to the discontinuity in the strain rate in the mean flow at the surface. This additional term is responsible for destabilizing the surface when the mean velocity increases beyond a transition value, and the physical mechanism driving the instability is the transfer of energy from the mean flow to the fluctuations due to the work done by the mean flow at the interface. The transition velocity Gamma(t) for the presence of surface instabilities depends on the wavenumber k and three dimensionless parameters: the ratio of the solid and fluid viscosities eta(r) = (eta(s)/eta), the capillary number Lambda = (T/GR) and the ratio of radii H, where T is the surface tension of the interface. For eta(r) = 0 and Lambda = 0, the transition velocity Gamma(t) diverges in the limits k much less than 1 and k much greater than 1, and has a minimum for finite k. The qualitative behaviour of the transition velocity is the same for Lambda > 0 and eta(r) = 0, though there is an increase in Gamma(t) in the limit k much greater than 1. When the viscosity of the surface is non-zero (eta(r) > 0), however, there is a qualitative change in the Gamma(t) vs. k curves. For eta(r) < 1, the transition velocity Gamma(t) is finite only when k is greater than a minimum value k(min), while perturbations with wavenumber k < k(min) are stable even for Gamma--> infinity. For eta(r) > 1, Gamma(t) is finite only for k(min) < k < k(max), while perturbations with wavenumber k < k(min) or k > k(max) are stable in the limit Gamma--> infinity. As H decreases or eta(r) increases, the difference k(max)- k(min) decreases. At minimum value H = H-min, which is a function of eta(r), the difference k(max)-k(min) = 0, and for H < H-min, perturbations of all wavenumbers are stable even in the limit Gamma--> infinity. The calculations indicate that H-min shows a strong divergence proportional to exp (0.0832 eta(r)(2)) for eta(r) much greater than 1.

An Efficient Robust Algorithm for the Surface-Potential Calculation of Independent DG MOSFET

Although the recently proposed single-implicit-equation-based input voltage equations (IVEs) for the independent double-gate (IDG) MOSFET promise faster computation time than the earlier proposed coupled-equations-based IVEs, it is not clear how those equations could be solved inside a circuit simulator as the conventional Newton-Raphson (NR)-based root finding method will not always converge due to the presence of discontinuity at the G-zero point (GZP) and nonremovable singularities in the trigonometric IVE. In this paper, we propose a unique algorithm to solve those IVEs, which combines the Ridders algorithm with the NR-based technique in order to provide assured convergence for any bias conditions. Studying the IDG MOSFET operation carefully, we apply an optimized initial guess to the NR component and a minimized solution space to the Ridders component in order to achieve rapid convergence, which is very important for circuit simulation. To reduce the computation budget further, we propose a new closed-form solution of the IVEs in the near vicinity of the GZP. The proposed algorithm is tested with different device parameters in the extended range of bias conditions and successfully implemented in a commercial circuit simulator through its Verilog-A interface.

The Wiener-Hopf solution of a class of mixed boundary value problems arising in surface water wave phenomena

Two mixed boundary value problems associated with two-dimensional Laplace equation, arising in the study of scattering of surface waves in deep water (or interface waves in two superposed fluids) in the linearised set up, by discontinuities in the surface (or interface) boundary conditions, are handled for solution by the aid of the Weiner-Hopf technique applied to a slightly more general differential equation to be solved under general boundary conditions and passing on to the limit in a manner so as to finally give rise to the solutions of the original problems. The first problem involves one discontinuity while the second problem involves two discontinuities. The reflection coefficient is obtained in closed form for the first problem and approximately for the second. The behaviour of the reflection coefficient for both the problems involving deep water against the incident wave number is depicted in a number of figures. It is observed that while the reflection coefficient for the first problem steadily increases with the wave number, that for the second problem exhibits oscillatory behaviour and vanishes at some discrete values of the wave number. Thus, there exist incident wave numbers for which total transmission takes place for the second problem. (C) 1999 Elsevier Science B.V. All rights reserved.

A fixed-grid finite element based enthalpy formulation for generalized phase change problems: role of superficial mushy region

The enthalpy method is primarily developed for studying phase change in a multicomponent material, characterized by a continuous liquid volume fraction (phi(1)) vs temperature (T) relationship. Using the Galerkin finite element method we obtain solutions to the enthalpy formulation for phase change in 1D slabs of pure material, by assuming a superficial phase change region (linear (phi(1) vs T) around the discontinuity at the melting point. Errors between the computed and analytical solutions are evaluated for the fluxes at, and positions of, the freezing front, for different widths of the superficial phase change region and spatial discretizations with linear and quadratic basis functions. For Stefan number (St) varying between 0.1 and 10 the method is relatively insensitive to spatial discretization and widths of the superficial phase change region. Greater sensitivity is observed at St = 0.01, where the variation in the enthalpy is large. In general the width of the superficial phase change region should span at least 2-3 Gauss quadrature points for the enthalpy to be computed accurately. The method is applied to study conventional melting of slabs of frozen brine and ice. Regardless of the forms for the phi(1) vs T relationships, the thawing times were found to scale as the square of the slab thickness. The ability of the method to efficiently capture multiple thawing fronts which may originate at any spatial location within the sample, is illustrated with the microwave thawing of slabs and 2D cylinders. (C) 2002 Elsevier Science Ltd. All rights reserved.

Investigation on the Sources of Reflection at the Junction of Tall Grounded Object and Lightning Channel

During lightning strike to a tall grounded object (TGO), reflections of current waves are known to occur at either ends of the TGO. These reflection modify the channel current and hence, the lightning electromagnetic fields. This study aims to identify the possible contributing factors to reflection at a TGO-channel junction for the current waves ascending on the TGO. Possible sources of reflection identified are corona sheath and discontinuity of resistance and radius. For analyzing the contribution of corona sheath and discontinuity of resistance at the junction, a macroscopic physical model for the return stroke developed in our earlier work is employed. NEC-2D is used for assessing the contribution of abrupt change in radii at a TGO-channel junction. The wire-cage model adopted for the same is validated using laboratory experiments. Detailed investigation revealed the following. The main contributor for reflection at a TGO-channel junction is the difference between TGO and channel core radii. Also, the discontinuity of resistance at a TGO-channel junction can be of some relevance only for the first microsecond regime. Further, corona sheath does not play any significant role in the reflection.

Numerical modeling of the contaminant Transport through stratified soil

The transport processes of the dissolved chemicals in stratified or layered soils have been studied for several decades. In case of the solute transport through stratified layers, interface condition plays an important role in determining appropriate transport parameters. First‐ type and third‐ type interface conditions are generally used in the literature. A first‐type interface condition will result in a continuous concentration profile across the interface at the expense of solute mass balance. On the other hand, a discontinuity in concentration develops when a third‐ type interface condition is used. To overcome this problem, a combined first‐ and third‐ type condition at the interface has been widely employed which yields second‐ type condition. This results in a similar break‐through curve irrespective of the layering order, which is non‐physical. In this work, an interface condition is proposed which satisfies the mass balance implicitly and brings the distinction between the breakthrough curves for different layering sequence corroborating with the experimental observations. This is in disagreement with the earlier work by H. M. Selim and co‐workers but, well agreement with the hypothetical result by Bosma and van der Zee; and Van der Zee.

Wavelet transform analysis of truncated fringe patterns in 3-D surface profilometry

Wavelet transform analysis of projected fringe pattern for phase recovery in 3-D shape measurement of objects is investigated. The present communication specifically outlines and evaluates the errors that creep in to the reconstructed profiles when fringe images do not satisfy periodicity. Three specific cases that give raise to non-periodicity of fringe image are simulated and leakage effects caused by each one of them are analyzed with continuous complex Morlet wavelet transform. Same images are analyzed with FFT method to make a comparison of the reconstructed profiles with both methods. Simulation results revealed a significant advantage of wavelet transform profilometry (WTP), that the distortions that arise due to leakage are confined to the locations of discontinuity and do not spread out over the entire projection as in the case of Fourier transform profilometry (FTP).

Improvements in Efficiency of Surface Potential Computation for Independent DG MOSFET

A robust numerical solution of the input voltage equations (IVEs) for the independent-double-gate metal-oxide-semiconductor field-effect transistor requires root bracketing methods (RBMs) instead of the commonly used Newton-Raphson (NR) technique due to the presence of nonremovable discontinuity and singularity. In this brief, we do an exhaustive study of the different RBMs available in the literature and propose a single derivative-free RBM that could be applied to both trigonometric and hyperbolic IVEs and offers faster convergence than the earlier proposed hybrid NR-Ridders algorithm. We also propose some adjustments to the solution space for the trigonometric IVE that leads to a further reduction of the computation time. The improvement of computational efficiency is demonstrated to be about 60% for trigonometric IVE and about 15% for hyperbolic IVE, by implementing the proposed algorithm in a commercial circuit simulator through the Verilog-A interface and simulating a variety of circuit blocks such as ring oscillator, ripple adder, and twisted ring counter.