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.

Spiral-wave dynamics in ionically realistic mathematical models for human ventricular tissue: the effects of periodic deformation

We carry out an extensive numerical study of the dynamics of spiral waves of electrical activation, in the presence of periodic deformation (PD) in two-dimensional simulation domains, in the biophysically realistic mathematical models of human ventricular tissue due to (a) ten-Tusscher and Panfilov (the TP06 model) and (b) ten-Tusscher, Noble, Noble, and Panfilov (the TNNPO4 model). We first consider simulations in cable-type domains, in which we calculate the conduction velocity theta and the wavelength lambda of a plane wave; we show that PD leads to a periodic, spatial modulation of theta and a temporally periodic modulation of lambda; both these modulations depend on the amplitude and frequency of the PD. We then examine three types of initial conditions for both TP06 and TNNPO4 models and show that the imposition of PD leads to a rich variety of spatiotemporal patterns in the transmembrane potential including states with a single rotating spiral (RS) wave, a spiral-turbulence (ST) state with a single meandering spiral, an ST state with multiple broken spirals, and a state SA in which all spirals are absorbed at the boundaries of our simulation domain. We find, for both TP06 and TNNPO4 models, that spiral-wave dynamics depends sensitively on the amplitude and frequency of PD and the initial condition. We examine how these different types of spiral-wave states can be eliminated in the presence of PD by the application of low-amplitude pulses by square- and rectangular-mesh suppression techniques. We suggest specific experiments that can test the results of our simulations.

Study on current-random wave forces acting on a framework

The forces of random wave plus current acting on a simplified offshore platform (jacket) model have been studied numerically and experimentally. The numerical results are in good agreement with experiments. The mean force can be approximated as a function of equivalent velocity parameter and the root-mean-square force as a function of equivalent significant wave height parameter.

A Note On The Numerical Simulations Of Flow Past A Wavy Square-Section Cylinder

The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Karman vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.

Controlling wave propagation through nonlinear engineered granular systems

We study the fundamental dynamic behavior of a special class of ordered granular systems in order to design new, structured materials with unique physical properties. The dynamic properties of granular systems are dictated by the nonlinear, Hertzian, potential in compression and zero tensile strength resulting from the discrete material structure. Engineering the underlying particle arrangement of granular systems allows for unique dynamic properties, not observed in natural, disordered granular media. While extensive studies on 1D granular crystals have suggested their usefulness for a variety of engineering applications, considerably less attention has been given to higher-dimensional systems. The extension of these studies in higher dimensions could enable the discovery of richer physical phenomena not possible in 1D, such as spatial redirection and anisotropic energy trapping. We present experiments, numerical simulation (based on a discrete particle model), and in some cases theoretical predictions for several engineered granular systems, studying the effects of particle arrangement on the highly nonlinear transient wave propagation to develop means for controlling the wave propagation pathways. The first component of this thesis studies the stress wave propagation resulting from a localized impulsive loading for three different 2D particle lattice structures: square, centered square, and hexagonal granular crystals. By varying the lattice structure, we observe a wide range of properties for the propagating stress waves: quasi-1D solitary wave propagation, fully 2D wave propagation with tunable wave front shapes, and 2D pulsed wave propagation. Additionally the effects of weak disorder, inevitably present in real granular systems, are investigated. The second half of this thesis studies the solitary wave propagation through 2D and 3D ordered networks of granular chains, reducing the effective density compared to granular crystals by selectively placing wave guiding chains to control the acoustic wave transmission. The rapid wave front amplitude decay exhibited by these granular networks makes them highly attractive for impact mitigation applications. The agreement between experiments, numerical simulations, and applicable theoretical predictions validates the wave guiding capabilities of these engineered granular crystals and networks and opens a wide range of possibilities for the realization of increasingly complex granular material design.

Electromagnetic wave propagation and scattering in a randomly-inhomogeneous dielectric sphere

Electromagnetic wave propagation and scattering in a sphere composed of an inhomogeneous medium having random variations in its permittivity are studied by utilizing the Born approximation in solving the vector wave equation. The variations in the permittivity are taken to be isotropic and homogeneous, and are spatially characterized by a Gaussian correlation function. Temporal variations in the medium are not considered.

Two particular problems are considered: i) finding the far-zone electric field when an electric or magnetic dipole is situated at the center of the sphere, and ii) finding the electric field at the sphere's center when a linearly polarized plane wave is incident upon it. Expressions are obtained for the mean-square magnitudes of the scattered field components; it is found that the mean of the product of any two transverse components vanishes. The cases where the wavelength is much shorter than correlation distance of the medium and where it is much longer than it are both considered.

Scattering by two degenerate anisotropic modes in square-lattice dielectric photonic crystals

We systematically investigate the square-lattice dielectric photonic crystals that have been used to demonstrate flat slab imaging experimentally. A right-handed Bloch mode is found in the left-handed frequency region by using the plane wave expansion method to analyze the photonic band structure and equifrequency contours. Using the multiple scattering theory, numerical simulations demonstrate that the left-handed mode and the right-handed mode are excited simultaneously by a point source and result in two kinds of transmitted waves. Impacted by the evanescent waves, superposition of these transmitted waves brings on complicated near field distributions such as the so-called imaging and its disappearance.

Optical properties of a square-lattice photonic crystal within the partial bandgap

Optical properties of a two-dimensional square-lattice photonic crystal are systematically investigated within the partial bandgap through anisotropic characteristics analysis and numerical simulation of field pattern. Using the plane-wave expansion method and Hellmann-Feynman theorem, the relationships between the incident and refracted angles for both phase and group velocities are calculated to analyze light propagation from air to photonic crystals. Three kinds of flat slab focusing are summarized and demonstrated by numerical simulations using the multiple scattering method. (c) 2007 Optical Society of America

Dynamic compressive response of composite square honeycombs

Carbon fibre-epoxy composite square honeycombs, and the parent composite material, were tested in quasi-static compression at a strain rate of 10 -3 s -1 and in dynamic compression at strain rates of 10 3-10 4 s -1 using an instrumented Kolsky bar arrangement. Taken together, these tests provide an assessment of the potential of this composite topology for use as a lightweight sandwich core. The honeycombs had two relative densities, 0.12 and 0.24, and two material orientations, ±45° and 0/90° with respect to the prismatic, loading direction of the honeycomb. Honeycomb manufacture was by slotting, assembling and bonding together carbon fibre/epoxy woven plies of composite sheets of 2 × 2 twill weave construction. The peak value of wall stress in the honeycombs was about one third that of the parent material, for all strain rates. An elastic finite element analysis was used to trace the source of this knock-down in strength: a stress concentration exists at the root of the slots and leads to premature failure by microbuckling. Shock-wave effects were evident at impact velocities exceeding 50 ms -1 for the honeycomb of relative density 0.12. This was traced to stubbing of the buckled cell walls against the face of the Kolsky bar. © 2011 Elsevier Ltd. All rights reserved.

Two-port InGaAsP/InP square resonator microlasers

Two-port InGaAsP/InP square resonator microlasers with a side length of 20 mm have been fabricated by the planar technology process, which have two 1 mu m-wide output ports connected to the vertices of the square resonator. Continuous-wave electrically injected microsquare lasers have been realised at room temperature with mode Q-factors of 1.75 x 10(4) at the threshold current.

The impact of imperfect symmetry on band edge modes of a two-dimensional photonic crystal with square lattice

A theoretical analysis has been performed by means of the plane-wave expansion method to examine the dispersion properties of photons at high symmetry points of an InP based two-dimensional photonic crystal with square lattice. The Q factors are compared qualitatively. The mechanism of surface-emitting is due to the photon manipulation by periodic dielectric materials in terms of Bragg diffraction. A surface-emitting photonic crystal resonator is designed based on the phenomenon of slow light. Photonic crystal slabs with different unit cells are utilized in the simulation. The results indicate that the change of the air holes can affect the polarization property of the modes. So we can find a way to improve the polarization by reducing the symmetry of the structure.

Equilateral-Triangle and Square Resonator Semiconductor Microlasers

The characteristics of equilateral-triangle resonator (ETR) and square resonator microlasers are reported, which are potential light sources in the photonic integrations. Based on the numerical simulations, we find that high-efficiency directional emission can be achieved for the triangle and square microlasers by directly connecting an output waveguide to the resonators. The electrically injected InP/InGaAsP ETR and square resonator microlasers with a 2-mu m-wide output waveguide were fabricated by standard photolithography and inductively coupled plasma etching techniques. Room-temperature continuous-wave (CW) operations were achieved for the ETR microlasers with the side length from 10 to 30 mu m and the square resonator microlasers with the side length of 20 mu m. The output power versus CW injection current and the laser spectra are presented for an ETR microlaser up to 310 K and a square resonator microlaser to 305 K. The lasing spectra with mode wavelength intervals as that of whispering-gallery-type modes and Fabry-Perot modes are observed for two square lasers, which can lase at low temperature and room temperature, respectively.

Modes in square resonators

Modes in square resonators are analyzed and classified according to the irreducible representations of the point group C-4v. If the mode numbers p and q that denote the number of wave nodes in the directions of two orthogonal square sides are unequal and have the same even-odd characteristics, the corresponding double modes are accidentally degenerate and can be combined into two new distributions with definite parities relative to the square diagonal mirror planes. The distributions with odd parities belong to the whispering-gallery-like modes in square resonators. The mode frequencies and quality factors are also calculated by the finite-difference time-domain technique and Pade approximation method. The numerically calculated mode frequencies agree with the theoretical ones very well and the whispering-gallery-like modes have quality factors much higher than other modes, including their accidentally degenerate counterparts in square resonators.

InGaAsP-InP Square Microlasers With a Vertex Output Waveguide

InGaAsP-InP square microlasers with a vertex output waveguide are fabricated by planar processes, and the etched sidewalls of the lasers are confined by insulating layer SiO2 and p-electrode TiAu metals. For a square microlaser with a side length of 30 mu m and a 2-mu m-wide output waveguide, a continuous-wave threshold current is 26 mA at room temperature and output power is 0.72 mW at 86 mA. The mode interval of 21 and 7.4 nm is observed for the microlasers with the side length of 10 and 30 mu m, respectively. Finite-difference time-domain (FDTD) simulations indicate that the lasing modes have incident angles of about 45 degrees at the boundaries of the resonator. In addition, square resonators surrounded by air, SiO2-Ti-Au, and SiO2-Au are compared by FDTD simulations.

Mechanism Responsible for the Complete Suppression of Karman Vortex in Flows Past a Wavy Square-Section Cylinder

The Karman vortex shedding is totally suppressed in flows past a wavy square-section cylinder at a Reynolds number of 100 and the wave steepness of 0.025. Such a phenomenon is illuminated by the numerical simulations. In the present study, the mechanism responsible for it is mainly attributed to the vertical vorticity. The geometric disturbance on the rear surface leads to the appearance of spanwise flow near the base. The specific vertical vorticity is generated on the rear surface and convecting into the near wake. The wake flow is recirculated with the appearance of the pair of recirculating cells. The interaction between the upper and lower shear layers is weakened by such cells, so that the vortex rolls could not be formed and the near wake flow becomes stable.