A spectral multiscale method for wave propagation analysis: Atomistic-continuum coupled simulation

In this paper, we present a new multiscale method which is capable of coupling atomistic and continuum domains for high frequency wave propagation analysis. The problem of non-physical wave reflection, which occurs due to the change in system description across the interface between two scales, can be satisfactorily overcome by the proposed method. We propose an efficient spectral domain decomposition of the total fine scale displacement along with a potent macroscale equation in the Laplace domain to eliminate the spurious interfacial reflection. We use Laplace transform based spectral finite element method to model the macroscale, which provides the optimum approximations for required dynamic responses of the outer atoms of the simulated microscale region very accurately. This new method shows excellent agreement between the proposed multiscale model and the full molecular dynamics (MD) results. Numerical experiments of wave propagation in a 1D harmonic lattice, a 1D lattice with Lennard-Jones potential, a 2D square Bravais lattice, and a 2D triangular lattice with microcrack demonstrate the accuracy and the robustness of the method. In addition, under certain conditions, this method can simulate complex dynamics of crystalline solids involving different spatial and/or temporal scales with sufficient accuracy and efficiency. (C) 2014 Elsevier B.V. All rights reserved.

Wave propagation analysis in adhesively bonded composite joints using the wavelet spectral finite element method

A wavelet spectral finite element (WSFE) model is developed for studying transient dynamics and wave propagation in adhesively bonded composite joints. The adherands are formulated as shear deformable beams using the first order shear deformation theory (FSDT) to obtain accurate results for high frequency wave propagation. Equations of motion governing wave motion in the bonded beams are derived using Hamilton's principle. The adhesive layer is modeled as a line of continuously distributed tension/compression and shear springs. Daubechies compactly supported wavelet scaling functions are used to transform the governing partial differential equations from time domain to frequency domain. The dynamic stiffness matrix is derived under the spectral finite element framework relating the nodal forces and displacements in the transformed frequency domain. Time domain results for wave propagation in a lap joint are validated with conventional finite element simulations using Abaqus. Frequency domain spectrum and dispersion relation results are presented and discussed. The developed WSFE model yields efficient and accurate analysis of wave propagation in adhesively-bonded composite joints. (C) 2014 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.

A New Version of Smoothed Point Method to Simulate Linear Elastic Wave Propagation

By using the kernel function of the smoothed particle hydrodynamics (SPH) and modification of statistical volumes of the boundary points and their kernel functions, a new version of smoothed point method is established for simulating elastic waves in solid. With the simplicity of SPH kept, the method is easy to handle stress boundary conditions, especially for the transmitting boundary condition. A result improving by de-convolution is also proposed to achieve high accuracy under a relatively large smooth length. A numerical example is given and compared favorably with the analytical solution.

Improved Analysis Method for Wave-Induced Pipeline Stability on Sandy Seabed

The existing Det Norske Veritas DNV Recommended Practice RP E305 for pipeline on-bottom stability is mainly based on the pipe–soil interaction model reported by Wagner et al. in 1987, and the wake model reported by Lambrakos et al. in 1987, to calculate the soil resistance and the hydrodynamic forces upon pipeline, respectively. Unlike the methods in the DNV Practice, in this paper, an improved analysis method is proposed for the on-bottom stability of a submarine pipeline, which is based on the relationships between Um/ gD 0.5 and Ws / D2 for various restraint conditions obtained by the hydrodynamic loading experiments, taking into account the coupling effects between wave, pipeline, and sandy seabed. The analysis procedure is illustrated with a detailed flow chart. A comparison is made between the submerged weights of pipeline predicted with the DNV Practice and those with the new method. The proposed analysis method may provide a helpful tool for the engineering practice of pipeline on-bottom stability design.

Optimization of tried wave-function in quantum monte-carlo method

A method for optimizing tried wave functions in quantum Monte Carlo method has been found and used to calculate the energies of molecules, such as H-2, Li-2, H-3+, H-3 and H-4. Good results were obtained.

A shock tube method in the kinetic study of nonequilibrium ionization-recombination by rarefaction wave cooling

A shock tube method is introduced to study the ionization–recombination kinetics of high temperature gas, in which a test gas is heated and ionized by a reflected shock wave and subsequently quenched by a strong rarefaction wave reflected on the end wall of the driver section as the main cooling wave associated with a rarefaction wave incident back into region 5 when the reflected shock wave interacts with the contact surface. As the quenching rate of the strong rarefaction wave reaches 106 K/s, a nonequilibrium ionization-recombination process occurs, during which the ion recombination with electrons dominates.

Wave-front analysis method of circular aperture sampling for collimation testing

A new type of wave-front analysis method for the collimation testing of laser beams is proposed. A concept of wave-front height is defined, and, on this basis, the wave-front analysis method of circular aperture sampling is introduced. The wave-front height of the tested noncollimated wave can be estimated from the distance between two identical fiducial diffraction planes of the sampled wave, and then the divergence is determined. The design is detailed, and the experiment is demonstrated. The principle and experiment results of the method are presented. Owing to the simplicity of the method and its low cost, it is a promising method for checking the collimation of a laser beam with a large divergence. © 2005 Optical Society of America.

Calculation of Valence Subband Structures for Strained Quantum-Wells by Plane Wave Expansion Method Within 6 * 6 Luttinger-Kohn Model

The valence subband energies and wave functions of a tensile strained quantum well are calculated by the plane wave expansion method within the 6 * 6 Luttinger-Kohn model. The effect of the number and period of plane-waves used for expansion on the stability of energy eigenvalues is examined. For practical calculation, it should choose the period large sufficiently to ensure the envelope functions vanish at the boundary and the number of plane waves large enough to ensure the energy eigenvalues keep unchanged within a prescribed range.

An augmented method for free boundary problems with moving contact lines

An augmented immersed interface method (IIM) is proposed for simulating one-phase moving contact line problems in which a liquid drop spreads or recoils on a solid substrate. While the present two-dimensional mathematical model is a free boundary problem, in our new numerical method, the fluid domain enclosed by the free boundary is embedded into a rectangular one so that the problem can be solved by a regular Cartesian grid method. We introduce an augmented variable along the free boundary so that the stress balancing boundary condition is satisfied. A hybrid time discretization is used in the projection method for better stability. The resultant Helmholtz/Poisson equations with interfaces then are solved by the IIM in an efficient way. Several numerical tests including an accuracy check, and the spreading and recoiling processes of a liquid drop are presented in detail. (C) 2010 Elsevier Ltd. All rights reserved.

A New High Order Accurate Shock Capture Method with Wave Booster

A new kind of shock capturing method is developed. Before applying the high order accurate traditional scheme which is called as base scheme in this paper the fluid parameters are preconditioned in order to control the group velocity. The newly constructed scheme is high order accurate, simple, has high resolution of the shock, and less computer time consumed.

A parametric method of retrieving ocean wave spectra from synthetic aperture radar images

A parametric method that extracts the ocean wave directional spectra from synthetic aperture radar (SAR) image is presented. The 180 degrees ambiguity of SAR image and the loss of information beyond the azimuthal cutoff can be overcome with this method. The ocean wave spectra can be obtained from SAR image directly by using iteration inversion mapping method with forward nonlinear mapping. Some numerical experiments have been made by using ERS-1 satellite SAR imagette data. The ocean wave direction retrieved from SAR imagette data is in agreement with the wind direction from the scatterometer data.