An automatic algorithm for mixed mode crack growth rate based on drop potential method

A new automatic algorithm for the assessment of mixed mode crack growth rate characteristics is presented based on the concept of an equivalent crack. The residual ligament size approach is introduced to implementation this algorithm for identifying the crack tip position on a curved path with respect to the drop potential signal. The automatic algorithm accounting for the curvilinear crack trajectory and employing an electrical potential difference was calibrated with respect to the optical measurements for the growing crack under cyclic mixed mode loading conditions. The effectiveness of the proposed algorithm is confirmed by fatigue tests performed on ST3 steel compact tension-shear specimens in the full range of mode mixities from pure mode Ito pure mode II. (C) 2015 Elsevier Ltd. All rights reserved.

Learning-Based Fast Iterative Convergence of 3-D MoM via Eigen-AGMRES Method

Three-dimensional (3-D) full-wave electromagnetic simulation using method of moments (MoM) under the framework of fast solver algorithms like fast multipole method (FMM) is often bottlenecked by the speed of convergence of the Krylov-subspace-based iterative process. This is primarily because the electric field integral equation (EFIE) matrix, even with cutting-edge preconditioning techniques, often exhibits bad spectral properties arising from frequency or geometry-based ill-conditioning, which render iterative solvers slow to converge or stagnate occasionally. In this communication, a novel technique to expedite the convergence of MoMmatrix solution at a specific frequency is proposed, by extracting and applying Eigen-vectors from a previously solved neighboring frequency in an augmented generalized minimum residual (AGMRES) iterative framework. This technique can be applied in unison with any preconditioner. Numerical results demonstrate up to 40% speed-up in convergence using the proposed Eigen-AGMRES method.

Differential binding of ppGpp and pppGpp to E-coli RNA polymerase: photo-labeling and mass spectral studies

(p) ppGpp, a secondary messenger, is induced under stress and shows pleiotropic response. It binds to RNA polymerase and regulates transcription in Escherichia coli. More than 25 years have passed since the first discovery was made on the direct interaction of ppGpp with E. coli RNA polymerase. Several lines of evidence suggest different modes of ppGpp binding to the enzyme. Earlier cross-linking experiments suggested that the beta-subunit of RNA polymerase is the preferred site for ppGpp, whereas recent crystallographic studies pinpoint the interface of beta'/omega-subunits as the site of action. With an aim to validate the binding domain and to follow whether tetra-and pentaphosphate guanosines have different location on RNA polymerase, this work was initiated. RNA polymerase was photo-labeled with 8-azido-ppGpp/8-azido-pppGpp, and the product was digested with trypsin and subjected to mass spectrometry analysis. We observed three new peptides in the trypsin digest of the RNA polymerase labeled with 8-azido-ppGpp, of which two peptides correspond to the same pocket on beta'-subunit as predicted by X-ray structural analysis, whereas the third peptide was mapped on the beta-subunit. In the case of 8-azido-pppGpp-labeled RNA polymerase, we have found only one cross-linked peptide from the beta'-subunit. However, we were unable to identify any binding site of pppGpp on the beta-subunit. Interestingly, we observed that pppGpp at high concentration competes out ppGpp bound to RNA polymerase more efficiently, whereas ppGpp cannot titrate out pppGpp. The competition between tetraphosphate guanosine and pentaphosphate guanosine for E. coli RNA polymerase was followed by gel-based assay as well as by a new method known as DRaCALA assay.

An Effective Method of Determining the Residual Stress Gradients in a Micro-Cantilever

A method of determining the micro-cantilever residual stress gradients by studying its deflection and curvature is presented. The stress gradients contribute to both axial load and bending moment, which, in prebuckling regime, cause the structural stiffness change and curving up/down, respectively. As the axial load corresponds to the even polynomial terms of stress gradients and bending moment corresponds to the odd polynomial terms, the deflection itself is not enough to determine the axial load and bending moment. Curvature together with the deflection can uniquely determine these two parameters. Both linear analysis and nonlinear analysis of micro-cantilever deflection under axial load and bending moment are presented. Because of the stiffening effect due to the nonlinearity of (large) deformation, the difference between linear and nonlinear analyses enlarges as the micro-cantilever deflection increases. The model developed in this paper determines the resultant axial load and bending moment due to the stress gradients. Under proper assumptions, the stress gradients profile is obtained through the resultant axial load and bending moment.

Compact difference approximation with consistent boundary condition

For simulating multi-scale complex flow fields it should be noted that all the physical quantities we are interested in must be simulated well. With limitation of the computer resources it is preferred to use high order accurate difference schemes. Because of their high accuracy and small stencil of grid points computational fluid dynamics (CFD) workers pay more attention to compact schemes recently. For simulating the complex flow fields the treatment of boundary conditions at the far field boundary points and near far field boundary points is very important. According to authors' experience and published results some aspects of boundary condition treatment for far field boundary are presented, and the emphasis is on treatment of boundary conditions for the upwind compact schemes. The consistent treatment of boundary conditions at the near boundary points is also discussed. At the end of the paper are given some numerical examples. The computed results with presented method are satisfactory.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.

H4 exponential finite-difference scheme for convective diffusion-equations

Perturbations are applied to the convective coefficients and source term of a convection-diffusion equation so that second-order corrections may be applied to a second-order exponential scheme. The basic Structure of the equations in the resulting fourth-order scheme is identical to that for the second order. Furthermore, the calculations are quite simple as the second-order corrections may be obtained in a single pass using a second-order scheme. For one to three dimensions, the fourth-order exponential scheme is unconditionally stable. As examples, the method is applied to Burgers' and other fluid mechanics problems. Compared with schemes normally used, the accuracies are found to be good and the method is applicable to regions with large gradients.

Space-resolved spectral diagnosis of line-shaped laser-produced magnesium plasmas

Space-resolved spectra of line-shaped laser-produced magnesium plasmas in the normal direction of the target have been obtained using a pinhole crystal spectrograph. These spectra are treated by a spectrum analyzing code for obtaining the true spectra and fine structures of overlapped lines. The spatial distributions of electron temperature and density along the normal direction of the target surface have been obtained with different spectral diagnostic techniques. Especially, the electron density plateaus beyond the critical surface in line-shaped magnesium plasmas have been obtained with a fitting technique applied to the Stark-broadened Ly-alpha wings of hydrogenic ions. The difference of plasma parameters between those obtained by different diagnostic techniques is discussed. Other phenomena, such as plasma satellites, population inversion, etc., which are observed in magnesium plasmas, are also presented.

Constrained high-order statistical blind deconvolution of spectral data

A constrained high-order statistical algorithm is proposed to blindly deconvolute the measured spectral data and estimate the response function of the instruments simultaneously. In this algorithm, no prior-knowledge is necessary except a proper length of the unit-impulse response. This length can be easily set to be the width of the narrowest spectral line by observing the measured data. The feasibility of this method has been demonstrated experimentally by the measured Raman and absorption spectral data.

Phase control of higher spectral components in the presence of a static electric field

We investigate the higher spectral component generations driven by a few-cycle laser pulse in a dense medium when a static electric field is present. Our results show that, when assisted by a static electric field, the dependence of the transmitted laser spectrum on the carrier-envelope phase (CEP) is significantly increased. Continuum and distinct peaks can be achieved by controlling the CEP of the few-cycle ultrashort laser pulse. Such a strong variation is due to the fact that the presence of the static electric field modifies the waveform of the combined electric field, which further affects the spectral distribution of the generated higher spectral components.

Propagation of a few-cycle laser pulse in a V-type three-level system

Propagation of a few-cycle laser pulse in a V-type three-level system (fine structure levels of rubidium) is investigated numerically. The full three-level Maxwell-Bloch equations without the rotating wave approximation and the standing slowly varying envelope approximation are solved by using a finite-difference time-domain method. It is shown that, when the usual unequal oscillator strengths are considered, self-induced transparency cannot be recovered and higher spectral components can be produced even for small-area pulses. (c) 2005 Pleiades Publishing, Inc.

Cross spectral purity and its influence on ghost imaging experiments

Pseudo-thermal light has been widely used in ghost imaging experiments. In order to understand the differences between the pseudo-thermal source and thermal source, we propose a method to investigate whether a light source has cross spectral purity (CSP), and experimentally measure the cross spectral properties of the pseudo-thermal light source in near-field and far-field zones. Moreover we present a theoretical analysis of the cross spectral influence on ghost imaging. (c) 2006 Elsevier B.V. All rights reserved.

Gauge Invariant Spectral Cauchy Characteristic Extraction of Gravitational Waves in Computational General Relativity

We present a complete system for Spectral Cauchy characteristic extraction (Spectral CCE). Implemented in C++ within the Spectral Einstein Code (SpEC), the method employs numerous innovative algorithms to efficiently calculate the Bondi strain, news, and flux.

Spectral CCE was envisioned to ensure physically accurate gravitational wave-forms computed for the Laser Interferometer Gravitational wave Observatory (LIGO) and similar experiments, while working toward a template bank with more than a thousand waveforms to span the binary black hole (BBH) problem’s seven-dimensional parameter space.

The Bondi strain, news, and flux are physical quantities central to efforts to understand and detect astrophysical gravitational wave sources within the Simulations of eXtreme Spacetime (SXS) collaboration, with the ultimate aim of providing the first strong field probe of the Einstein field equation.

In a series of included papers, we demonstrate stability, convergence, and gauge invariance. We also demonstrate agreement between Spectral CCE and the legacy Pitt null code, while achieving a factor of 200 improvement in computational efficiency.

Spectral CCE represents a significant computational advance. It is the foundation upon which further capability will be built, specifically enabling the complete calculation of junk-free, gauge-free, and physically valid waveform data on the fly within SpEC.

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.

Range of validity of the method of averaging

Sufficient conditions are derived for the validity of approximate periodic solutions of a class of second order ordinary nonlinear differential equations. An approximate solution is defined to be valid if an exact solution exists in a neighborhood of the approximation.

Two classes of validity criteria are developed. Existence is obtained using the contraction mapping principle in one case, and the Schauder-Leray fixed point theorem in the other. Both classes of validity criteria make use of symmetry properties of periodic functions, and both classes yield an upper bound on a norm of the difference between the approximate and exact solution. This bound is used in a procedure which establishes sufficient stability conditions for the approximated solution.

Application to a system with piecewise linear restoring force (bilinear system) reveals that the approximate solution obtained by the method of averaging is valid away from regions where the response exhibits vertical tangents. A narrow instability region is obtained near one-half the natural frequency of the equivalent linear system. Sufficient conditions for the validity of resonant solutions are also derived, and two term harmonic balance approximate solutions which exhibit ultraharmonic and subharmonic resonances are studied.