Extended Jacobi Elliptic Function Rational Expansion Method and Its Application to (2+1)-Dimensional Stochastic Dispersive Long Wave System

International Nuclear Information System (INIS)

Song Lina; Zhang Hongqing

2007-01-01

In this work, by means of a generalized method and symbolic computation, we extend the Jacobi elliptic function rational expansion method to uniformly construct a series of stochastic wave solutions for stochastic evolution equations. To illustrate the effectiveness of our method, we take the (2+1)-dimensional stochastic dispersive long wave system as an example. We not only have obtained some known solutions, but also have constructed some new rational formal stochastic Jacobi elliptic function solutions.

Spherical-wave expansions of piston-radiator fields.

Science.gov (United States)

Wittmann, R C; Yaghjian, A D

1991-09-01

Simple spherical-wave expansions of the continuous-wave fields of a circular piston radiator in a rigid baffle are derived. These expansions are valid throughout the illuminated half-space and are useful for efficient numerical computation in the near-field region. Multipole coefficients are given by closed-form expressions which can be evaluated recursively.

Investigation of Solitary wave solutions for Vakhnenko-Parkes equation via exp-function and Exp(-ϕ(ξ))-expansion method.

Science.gov (United States)

Roshid, Harun-Or; Kabir, Md Rashed; Bhowmik, Rajandra Chadra; Datta, Bimal Kumar

2014-01-01

In this paper, we have described two dreadfully important methods to solve nonlinear partial differential equations which are known as exp-function and the exp(-ϕ(ξ)) -expansion method. Recently, there are several methods to use for finding analytical solutions of the nonlinear partial differential equations. The methods are diverse and useful for solving the nonlinear evolution equations. With the help of these methods, we are investigated the exact travelling wave solutions of the Vakhnenko- Parkes equation. The obtaining soliton solutions of this equation are described many physical phenomena for weakly nonlinear surface and internal waves in a rotating ocean. Further, three-dimensional plots of the solutions such as solitons, singular solitons, bell type solitary wave i.e. non-topological solitons solutions and periodic solutions are also given to visualize the dynamics of the equation.

Improved wave functions for large-N expansions

International Nuclear Information System (INIS)

Imbo, T.; Sukhatme, U.

1985-01-01

Existing large-N expansions of radial wave functions phi/sub n/,l(r) are only accurate near the minimum of the effective potential. Within the framework of the shifted 1/N expansion, we use known analytic results to motivate a simple modification so that the improved wave functions are accurate over a wide range of r and any choice of quantum numbers n and l. It is shown that these wave functions yield simple and accurate analytic expressions for certain quantities of interest in quarkonium physics

Expansion of continuum functions on resonance wave functions and amplitudes

International Nuclear Information System (INIS)

Bang, J.; Gareev, F.A.; Gizzatkulov, M.H.; Goncharov, S.A.

1978-01-01

To overcome difficulties encountered with wave functions of continuum spectrum (for example, in a shell model with continuum) the pole expansion (by the Mittag-Leffler theorem) of wave functions, scattering amplitudes and the Green functions with positive energies are considered. It is shown that resonance functions (the Gamov functions) form a complete set over which the continuum functions could be expanded. The general view of these expansions for final potentials and for the Coulomb repulsion potential are obtained and discussed. It is shown that the application of the method to nuclear structure calculations leads to simple algebraic equations

Travelling Wave Solutions of Coupled Burger’s Equations of Time-Space Fractional Order by Novel (Gʹ/G-Expansion Method

Directory of Open Access Journals (Sweden)

Rashida Hussain

2017-04-01

Full Text Available In this paper, Novel (Gʹ/G-expansion method is used to find new generalized exact travelling wave solutions of fractional order coupled Burger’s equations in terms of trigonometric functions, rational functions and hyperbolic functions with arbitrary parameters. For the conversion of the partial differential equation to the ordinary differential equation, complex transformation method is used. Novel (Gʹ/G-expansion method is very effective and provides a powerful mathematical tool to solve nonlinear equations. Moreover, for the representation of these exact solutions we have plotted graphs for different values of parameters which were in travelling waveform.

New solitary wave solutions of the time-fractional Cahn-Allen equation via the improved (G'/G)-expansion method

Science.gov (United States)

Batool, Fiza; Akram, Ghazala

2018-05-01

An improved (G'/G)-expansion method is proposed for extracting more general solitary wave solutions of the nonlinear fractional Cahn-Allen equation. The temporal fractional derivative is taken in the sense of Jumarie's fractional derivative. The results of this article are generalized and extended version of previously reported solutions.

On WKB expansions for Alfven waves in the solar wind

International Nuclear Information System (INIS)

Hollweg, J.V.

1990-01-01

The author reexamines the WKB expansion for toroidal Alfven waves in the solar wind, as described by equations (9) of Heinemann and Olbert (1980). His principal conclusions are as follows: (1) The WKB expansion used by Belcher (1971) and Hollweg (1973) is nonuniformly convergent. (2) Using the method of multiple scales (Nayfeh, 1981), he obtains an expansion which is uniform. (3) The uniform expansion takes into account the small modification to the Alfven wave phase speed due to spatial gradients of the background. (4) Both the uniform and nonuniform expansions reveal that each normal mode has both Elsaesser variables δz + ≠ 0 and δz - ≠ 0. Thus if δz - corresponds to the outgoing mode in a homogeneous background, an observation of δz + ≠ 0 does not necessarily imply the presence of the inward propagating mode, as is commonly assumed. (5) Even at the Alfven critical point (where V = υ A ) he finds that δz + ≠ 0. Thus incompressible MHD turbulence, which requires both δz + ≠ 0 and δz - ≠ 0, can proceed at the Alfven critical point (cf. Roberts, 1989). (6) With very few exceptions, the predictions of these calculations do not agree with recent observations (Marsch and Tu, 1990) of the power spectra of δz + and δz - in the solar wind. Thus the evolution of Alfven waves in the solar wind is governed by dynamics not included in the Heinemann and Olbert equations

An extended Jacobi elliptic function rational expansion method and its application to (2+1)-dimensional dispersive long wave equation

International Nuclear Information System (INIS)

Wang Qi; Chen Yong; Zhang Hongqing

2005-01-01

With the aid of computerized symbolic computation, a new elliptic function rational expansion method is presented by means of a new general ansatz, in which periodic solutions of nonlinear partial differential equations that can be expressed as a finite Laurent series of some of 12 Jacobi elliptic functions, is more powerful than exiting Jacobi elliptic function methods and is very powerful to uniformly construct more new exact periodic solutions in terms of rational formal Jacobi elliptic function solution of nonlinear partial differential equations. As an application of the method, we choose a (2+1)-dimensional dispersive long wave equation to illustrate the method. As a result, we can successfully obtain the solutions found by most existing Jacobi elliptic function methods and find other new and more general solutions at the same time. Of course, more shock wave solutions or solitary wave solutions can be gotten at their limit condition

Equation level matching: An extension of the method of matched asymptotic expansion for problems of wave propagation

Science.gov (United States)

Faria, Luiz; Rosales, Rodolfo

2017-11-01

We introduce an alternative to the method of matched asymptotic expansions. In the ``traditional'' implementation, approximate solutions, valid in different (but overlapping) regions are matched by using ``intermediate'' variables. Here we propose to match at the level of the equations involved, via a ``uniform expansion'' whose equations enfold those of the approximations to be matched. This has the advantage that one does not need to explicitly solve the asymptotic equations to do the matching, which can be quite impossible for some problems. In addition, it allows matching to proceed in certain wave situations where the traditional approach fails because the time behaviors differ (e.g., one of the expansions does not include dissipation). On the other hand, this approach does not provide the fairly explicit approximations resulting from standard matching. In fact, this is not even its aim, which to produce the ``simplest'' set of equations that capture the behavior. Ruben Rosales work was partially supported by NSF Grants DMS-1614043 and DMS-1719637.

New travelling wave solutions of the (1 + 1-dimensional cubic nonlinear Schrodinger equation using novel (G′/G-expansion method

Directory of Open Access Journals (Sweden)

M.G. Hafez

2016-06-01

Full Text Available In this paper, the novel (G′/G-expansion method is applied to construct exact travelling wave solutions of the cubic nonlinear Schrodinger equation. This technique is straightforward and simple to use, and gives more new general solutions than the other existing methods. Various types of solitary and periodic wave solutions of this equation are derived. The obtained results may be helpful to describe the wave propagation in soliton physics, such as soliton propagation in optical fibers, modulus instability in plasma physics, etc. and provided us the firm mathematical foundation in soliton physics or any varied instances. Furthermore, three-dimensional modules plot of the solutions are also given to visualize the dynamics of the equation.

Gabor Wave Packet Method to Solve Plasma Wave Equations

International Nuclear Information System (INIS)

Pletzer, A.; Phillips, C.K.; Smithe, D.N.

2003-01-01

A numerical method for solving plasma wave equations arising in the context of mode conversion between the fast magnetosonic and the slow (e.g ion Bernstein) wave is presented. The numerical algorithm relies on the expansion of the solution in Gaussian wave packets known as Gabor functions, which have good resolution properties in both real and Fourier space. The wave packets are ideally suited to capture both the large and small wavelength features that characterize mode conversion problems. The accuracy of the scheme is compared with a standard finite element approach

A double expansion method for the frequency response of finite-length beams with periodic parameters

Science.gov (United States)

Ying, Z. G.; Ni, Y. Q.

2017-03-01

A double expansion method for the frequency response of finite-length beams with periodic distribution parameters is proposed. The vibration response of the beam with spatial periodic parameters under harmonic excitations is studied. The frequency response of the periodic beam is the function of parametric period and then can be expressed by the series with the product of periodic and non-periodic functions. The procedure of the double expansion method includes the following two main steps: first, the frequency response function and periodic parameters are expanded by using identical periodic functions based on the extension of the Floquet-Bloch theorem, and the period-parametric differential equation for the frequency response is converted into a series of linear differential equations with constant coefficients; second, the solutions to the linear differential equations are expanded by using modal functions which satisfy the boundary conditions, and the linear differential equations are converted into algebraic equations according to the Galerkin method. The expansion coefficients are obtained by solving the algebraic equations and then the frequency response function is finally determined. The proposed double expansion method can uncouple the effects of the periodic expansion and modal expansion so that the expansion terms are determined respectively. The modal number considered in the second expansion can be reduced remarkably in comparison with the direct expansion method. The proposed double expansion method can be extended and applied to the other structures with periodic distribution parameters for dynamics analysis. Numerical results on the frequency response of the finite-length periodic beam with various parametric wave numbers and wave amplitude ratios are given to illustrate the effective application of the proposed method and the new frequency response characteristics, including the parameter-excited modal resonance, doubling-peak frequency response

Acoustic waves in shock tunnels and expansion tubes

Science.gov (United States)

Paull, A.; Stalker, R. J.

1992-01-01

It is shown that disturbances in shock and expansion tubes can be modelled as lateral acoustic waves. The ratio of sound speed across the driver-test gas interface is shown to govern the quantity of noise in the test gas. Frequency 'focusing' which is fundamental to centered unsteady expansions is discussed and displayed in centerline pitot pressure measurements.

Taylor-series method for four-nucleon wave functions

International Nuclear Information System (INIS)

Sandulescu, A.; Tarnoveanu, I.; Rizea, M.

1977-09-01

Taylor-series method for transforming the infinite or finite well two-nucleon wave functions from individual coordinates to relative and c.m. coordinates, by expanding the single particle shell model wave functions around c.m. of the system, is generalized to four-nucleon wave functions. Also the connections with the Talmi-Moshinsky method for two and four harmonic oscillator wave functions are deduced. For both methods Fortran IV programs for the expansion coefficients have been written and the equivalence of corresponding expressions numerically proved. (author)

A polynomial expansion method and its application in the coupled Zakharov-Kuznetsov equations

International Nuclear Information System (INIS)

Huang Wenhua

2006-01-01

A polynomial expansion method is presented to solve nonlinear evolution equations. Applying this method, the coupled Zakharov-Kuznetsov equations in fluid system are studied and many exact travelling wave solutions are obtained. These solutions include solitary wave solutions, periodic solutions and rational type solutions

Accurate expansion of cylindrical paraxial waves for its straightforward implementation in electromagnetic scattering

International Nuclear Information System (INIS)

Naserpour, Mahin; Zapata-Rodríguez, Carlos J.

2018-01-01

Highlights: • Paraxial beams are represented in a series expansion in terms of Bessel wave functions. • The coefficients of the series expansion can be analytically determined by using the pattern in the focal plane. • In particular, Gaussian beams and apertured wave fields have been critically examined. • This representation of the wave field is adequate for scattering problems with shaped beams. - Abstract: The evaluation of vector wave fields can be accurately performed by means of diffraction integrals, differential equations and also series expansions. In this paper, a Bessel series expansion which basis relies on the exact solution of the Helmholtz equation in cylindrical coordinates is theoretically developed for the straightforward yet accurate description of low-numerical-aperture focal waves. The validity of this approach is confirmed by explicit application to Gaussian beams and apertured focused fields in the paraxial regime. Finally we discuss how our procedure can be favorably implemented in scattering problems.

Critical point anomalies include expansion shock waves

Energy Technology Data Exchange (ETDEWEB)

Nannan, N. R., E-mail: ryan.nannan@uvs.edu [Mechanical Engineering Discipline, Anton de Kom University of Suriname, Leysweg 86, PO Box 9212, Paramaribo, Suriname and Process and Energy Department, Delft University of Technology, Leeghwaterstraat 44, 2628 CA Delft (Netherlands); Guardone, A., E-mail: alberto.guardone@polimi.it [Department of Aerospace Science and Technology, Politecnico di Milano, Via La Masa 34, 20156 Milano (Italy); Colonna, P., E-mail: p.colonna@tudelft.nl [Propulsion and Power, Delft University of Technology, Kluyverweg 1, 2629 HS Delft (Netherlands)

2014-02-15

From first-principle fluid dynamics, complemented by a rigorous state equation accounting for critical anomalies, we discovered that expansion shock waves may occur in the vicinity of the liquid-vapor critical point in the two-phase region. Due to universality of near-critical thermodynamics, the result is valid for any common pure fluid in which molecular interactions are only short-range, namely, for so-called 3-dimensional Ising-like systems, and under the assumption of thermodynamic equilibrium. In addition to rarefaction shock waves, diverse non-classical effects are admissible, including composite compressive shock-fan-shock waves, due to the change of sign of the fundamental derivative of gasdynamics.

Expansion shock waves in the implosion process from a time-reversible molecular-dynamics simulation of a dual explosion process

International Nuclear Information System (INIS)

Komatsu, Nobuyoshi; Abe, Takashi

2007-01-01

Why does not an expansion shock wave exist in a gaseous medium in nature? The reason has been widely believed to be the irreversibility in nature, while an obvious demonstration for this belief has not been accomplished yet. In order to resolve the question from a microscopic viewpoint, an implosion process dual to an explosion process was investigated by means of the molecular-dynamics method (MD). To this aim, we employed a ''bit-reversible algorithm (Bit MD)'' that was completely time-reversible in a microscopic viewpoint and was free from any round-off error. Here we show that, through a dual implosion simulation (i.e., a time-reversible simulation of the explosion), a kind of expansion shock wave is successfully formed in the Bit MD simulation. Furthermore, we show that when the controlled noise is intentionally added to the Bit MD, the expansion shock wave disappears dramatically and turns into an isentropic expansion wave, even if the noise is extremely small. Since the controlled noise gives rise to the irreversibility in the Bit MD simulation, it can be concluded that the irreversibility in the system prohibits the expansion shock wave from appearing in the system

Treatment of Ion-Atom Collisions Using a Partial-Wave Expansion of the Projectile Wavefunction

Science.gov (United States)

Wong, T. G.; Foster, M.; Colgan, J.; Madison, D. H.

2009-01-01

We present calculations of ion-atom collisions using a partial-wave expansion of the projectile wavefunction. Most calculations of ion-atom collisions have typically used classical or plane-wave approximations for the projectile wavefunction, since partial-wave expansions are expected to require prohibitively large numbers of terms to converge…

Experimental methods of shock wave research

CERN Document Server

Seiler, Friedrich

2016-01-01

This comprehensive and carefully edited volume presents a variety of experimental methods used in Shock Waves research. In 14 self contained chapters this 9th volume of the “Shock Wave Science and Technology Reference Library” presents the experimental methods used in Shock Tubes, Shock Tunnels and Expansion Tubes facilities. Also described is their set-up and operation. The uses of an arc heated wind tunnel and a gun tunnel are also contained in this volume. Whenever possible, in addition to the technical description some typical scientific results obtained using such facilities are described. Additionally, this authoritative book includes techniques for measuring physical properties of blast waves and laser generated shock waves. Information about active shock wave laboratories at different locations around the world that are not described in the chapters herein is given in the Appendix, making this book useful for every researcher involved in shock/blast wave phenomena.

The optimizied expansion method for wavefield extrapolation

KAUST Repository

Wu, Zedong

2013-01-01

Spectral methods are fast becoming an indispensable tool for wave-field extrapolation, especially in anisotropic media, because of its dispersion and artifact free, as well as highly accurate, solutions of the wave equation. However, for inhomogeneous media, we face difficulties in dealing with the mixed space-wavenumber domain operator.In this abstract, we propose an optimized expansion method that can approximate this operator with its low rank representation. The rank defines the number of inverse FFT required per time extrapolation step, and thus, a lower rank admits faster extrapolations. The method uses optimization instead of matrix decomposition to find the optimal wavenumbers and velocities needed to approximate the full operator with its low rank representation.Thus,we obtain more accurate wave-fields using lower rank representation, and thus cheaper extrapolations. The optimization operation to define the low rank representation depends only on the velocity model, and this is done only once, and valid for a full reverse time migration (many shots) or one iteration of full waveform inversion. Applications on the BP model yielded superior results than those obtained using the decomposition approach. For transversely isotopic media, the solutions were free of the shear wave artifacts, and does not require that eta>0.

An Enhanced Plane Wave Expansion Method to Solve Piezoelectric Phononic Crystal with Resonant Shunting Circuits

Directory of Open Access Journals (Sweden)

Ziyang Lian

2016-01-01

Full Text Available An enhanced plane wave expansion (PWE method is proposed to solve piezoelectric phononic crystal (PPC connected with resonant shunting circuits (PPC-C, which is named as PWE-PPC-C. The resonant shunting circuits can not only bring about the locally resonant (LR band gap for the PPC-C but also conveniently tune frequency and bandwidth of band gaps through adjusting circuit parameters. However, thus far, more than one-dimensional PPC-C has been studied just by Finite Element method. Compared with other methods, the PWE has great advantages in solving more than one-dimensional PC as well as various lattice types. Nevertheless, the conventional PWE cannot accurately solve coupling between the structure and resonant shunting circuits of the PPC-C since only taking one-way coupling from displacements to electrical parameters into consideration. A two-dimensional PPC-C model of orthorhombic lattice is established to demonstrate the whole solving process of PWE-PPC-C. The PWE-PPC-C method is validated by Transfer Matrix method as well as Finite Element method. The dependence of band gaps on circuit parameters has been investigated in detail by PWE-PPC-C. Its advantage in solving various lattice types is further illustrated by calculating the PPC-C of triangular and hexagonal lattices, respectively.

Rapid expansion method (REM) for time‐stepping in reverse time migration (RTM)

KAUST Repository

Pestana, Reynam C.

2009-01-01

We show that the wave equation solution using a conventional finite‐difference scheme, derived commonly by the Taylor series approach, can be derived directly from the rapid expansion method (REM). After some mathematical manipulation we consider an analytical approximation for the Bessel function where we assume that the time step is sufficiently small. From this derivation we find that if we consider only the first two Chebyshev polynomials terms in the rapid expansion method we can obtain the second order time finite‐difference scheme that is frequently used in more conventional finite‐difference implementations. We then show that if we use more terms from the REM we can obtain a more accurate time integration of the wave field. Consequently, we have demonstrated that the REM is more accurate than the usual finite‐difference schemes and it provides a wave equation solution which allows us to march in large time steps without numerical dispersion and is numerically stable. We illustrate the method with post and pre stack migration results.

Elliptic equation rational expansion method and new exact travelling solutions for Whitham-Broer-Kaup equations

International Nuclear Information System (INIS)

Chen Yong; Wang Qi; Li Biao

2005-01-01

Based on a new general ansatz and a general subepuation, a new general algebraic method named elliptic equation rational expansion method is devised for constructing multiple travelling wave solutions in terms of rational special function for nonlinear evolution equations (NEEs). We apply the proposed method to solve Whitham-Broer-Kaup equation and explicitly construct a series of exact solutions which include rational form solitary wave solution, rational form triangular periodic wave solutions and rational wave solutions as special cases. In addition, the links among our proposed method with the method by Fan [Chaos, Solitons and Fractals 2004;20:609], are also clarified generally

Interactions of solitary waves and compression/expansion waves in core-annular flows

Science.gov (United States)

Maiden, Michelle; Anderson, Dalton; El, Gennady; Franco, Nevil; Hoefer, Mark

2017-11-01

The nonlinear hydrodynamics of an initial step leads to the formation of rarefaction waves and dispersive shock waves in dispersive media. Another hallmark of these media is the soliton, a localized traveling wave whose speed is amplitude dependent. Although compression/expansion waves and solitons have been well-studied individually, there has been no mathematical description of their interaction. In this talk, the interaction of solitons and shock/rarefaction waves for interfacial waves in viscous, miscible core-annular flows are modeled mathematically and explored experimentally. If the interior fluid is continuously injected, a deformable conduit forms whose interfacial dynamics are well-described by a scalar, dispersive nonlinear partial differential equation. The main focus is on interactions of solitons with dispersive shock waves and rarefaction waves. Theory predicts that a soliton can either be transmitted through or trapped by the extended hydrodynamic state. The notion of reciprocity is introduced whereby a soliton interacts with a shock wave in a reciprocal or dual fashion as with the rarefaction. Soliton reciprocity, trapping, and transmission are observed experimentally and are found to agree with the modulation theory and numerical simulations. This work was partially supported by NSF CAREER DMS-1255422 (M.A.H.) and NSF GRFP (M.D.M.).

The linear potential propagator via wave function expansion

International Nuclear Information System (INIS)

Nassar, Antonio B.; Cattani, Mauro S.D.

2002-01-01

We evaluate the quantum propagator for the motion of a particle in a linear potential via a recently developed formalism [A.B. Nassar et al., Phys. Rev. E56, 1230, (1997)]. In this formalism, the propagator comes about as a type of expansion of the wave function over the space of the initial velocities. (author)

Efficient 3D frequency response modeling with spectral accuracy by the rapid expansion method

KAUST Repository

Chu, Chunlei

2012-07-01

Frequency responses of seismic wave propagation can be obtained either by directly solving the frequency domain wave equations or by transforming the time domain wavefields using the Fourier transform. The former approach requires solving systems of linear equations, which becomes progressively difficult to tackle for larger scale models and for higher frequency components. On the contrary, the latter approach can be efficiently implemented using explicit time integration methods in conjunction with running summations as the computation progresses. Commonly used explicit time integration methods correspond to the truncated Taylor series approximations that can cause significant errors for large time steps. The rapid expansion method (REM) uses the Chebyshev expansion and offers an optimal solution to the second-order-in-time wave equations. When applying the Fourier transform to the time domain wavefield solution computed by the REM, we can derive a frequency response modeling formula that has the same form as the original time domain REM equation but with different summation coefficients. In particular, the summation coefficients for the frequency response modeling formula corresponds to the Fourier transform of those for the time domain modeling equation. As a result, we can directly compute frequency responses from the Chebyshev expansion polynomials rather than the time domain wavefield snapshots as do other time domain frequency response modeling methods. When combined with the pseudospectral method in space, this new frequency response modeling method can produce spectrally accurate results with high efficiency. © 2012 Society of Exploration Geophysicists.

A new Riccati equation rational expansion method and its application to (2 + 1)-dimensional Burgers equation

International Nuclear Information System (INIS)

Wang Qi; Chen Yong; Zhang Hongqing

2005-01-01

In this paper, we present a new Riccati equation rational expansion method to uniformly construct a series of exact solutions for nonlinear evolution equations. Compared with most existing tanh methods and other sophisticated methods, the proposed method not only recover some known solutions, but also find some new and general solutions. The solutions obtained in this paper include rational triangular periodic wave solutions, rational solitary wave solutions and rational wave solutions. The efficiency of the method can be demonstrated on (2 + 1)-dimensional Burgers equation

Fully resolved simulations of expansion waves propagating into particle beds

Science.gov (United States)

Marjanovic, Goran; Hackl, Jason; Annamalai, Subramanian; Jackson, Thomas; Balachandar, S.

2017-11-01

There is a tremendous amount of research that has been done on compression waves and shock waves moving over particles but very little concerning expansion waves. Using 3-D direct numerical simulations, this study will explore expansion waves propagating into fully resolved particle beds of varying volume fractions and geometric arrangements. The objectives of these simulations are as follows: 1) To fully resolve all (1-way coupled) forces on the particles in a time varying flow and 2) to verify state-of-the-art drag models for such complex flows. We will explore a range of volume fractions, from very low ones that are similar to single particle flows, to higher ones where nozzling effects are observed between neighboring particles. Further, we will explore two geometric arrangements: body centered cubic and face centered cubic. We will quantify the effects that volume fraction and geometric arrangement plays on the drag forces and flow fields experienced by the particles. These results will then be compared to theoretical predictions from a model based on the generalized Faxen's theorem. This work was supported in part by the U.S. Department of Energy under the Predictive Science Academic Alliance Program, under Contract No. DE-NA0002378.

The modified alternative (G'/G)-expansion method to nonlinear evolution equation: application to the (1+1)-dimensional Drinfel'd-Sokolov-Wilson equation.

Science.gov (United States)

Akbar, M Ali; Mohd Ali, Norhashidah Hj; Mohyud-Din, Syed Tauseef

2013-01-01

Over the years, (G'/G)-expansion method is employed to generate traveling wave solutions to various wave equations in mathematical physics. In the present paper, the alternative (G'/G)-expansion method has been further modified by introducing the generalized Riccati equation to construct new exact solutions. In order to illustrate the novelty and advantages of this approach, the (1+1)-dimensional Drinfel'd-Sokolov-Wilson (DSW) equation is considered and abundant new exact traveling wave solutions are obtained in a uniform way. These solutions may be imperative and significant for the explanation of some practical physical phenomena. It is shown that the modified alternative (G'/G)-expansion method an efficient and advance mathematical tool for solving nonlinear partial differential equations in mathematical physics.

Teaching graphical simulations of Fourier series expansion of some periodic waves using spreadsheets

Science.gov (United States)

Singh, Iqbal; Kaur, Bikramjeet

2018-05-01

The present article demonstrates a way of programming using an Excel spreadsheet to teach Fourier series expansion in school/colleges without the knowledge of any typical programming language. By using this, a student learns to approximate partial sum of the n terms of Fourier series for some periodic signals such as square wave, saw tooth wave, half wave rectifier and full wave rectifier signals.

Gravity waves as a probe of the Hubble expansion rate during an electroweak scale phase transition

International Nuclear Information System (INIS)

Chung, Daniel J. H.; Zhou Peng

2010-01-01

Just as big bang nucleosynthesis allows us to probe the expansion rate when the temperature of the Universe was around 1 MeV, the measurement of gravity waves from electroweak scale first order phase transitions may allow us to probe the expansion rate when the temperature of the Universe was at the electroweak scale. We compute the simple transformation rule for the gravity wave spectrum under the scaling transformation of the Hubble expansion rate. We then apply this directly to the scenario of quintessence kination domination and show how gravity wave spectra would shift relative to Laser Interferometer Space Antenna and Big Bang Observer projected sensitivities.

Photonic band structures solved by a plane-wave-based transfer-matrix method.

Science.gov (United States)

Li, Zhi-Yuan; Lin, Lan-Lan

2003-04-01

Transfer-matrix methods adopting a plane-wave basis have been routinely used to calculate the scattering of electromagnetic waves by general multilayer gratings and photonic crystal slabs. In this paper we show that this technique, when combined with Bloch's theorem, can be extended to solve the photonic band structure for 2D and 3D photonic crystal structures. Three different eigensolution schemes to solve the traditional band diagrams along high-symmetry lines in the first Brillouin zone of the crystal are discussed. Optimal rules for the Fourier expansion over the dielectric function and electromagnetic fields with discontinuities occurring at the boundary of different material domains have been employed to accelerate the convergence of numerical computation. Application of this method to an important class of 3D layer-by-layer photonic crystals reveals the superior convergency of this different approach over the conventional plane-wave expansion method.

Photonic band structures solved by a plane-wave-based transfer-matrix method

International Nuclear Information System (INIS)

Li Zhiyuan; Lin Lanlan

2003-01-01

Transfer-matrix methods adopting a plane-wave basis have been routinely used to calculate the scattering of electromagnetic waves by general multilayer gratings and photonic crystal slabs. In this paper we show that this technique, when combined with Bloch's theorem, can be extended to solve the photonic band structure for 2D and 3D photonic crystal structures. Three different eigensolution schemes to solve the traditional band diagrams along high-symmetry lines in the first Brillouin zone of the crystal are discussed. Optimal rules for the Fourier expansion over the dielectric function and electromagnetic fields with discontinuities occurring at the boundary of different material domains have been employed to accelerate the convergence of numerical computation. Application of this method to an important class of 3D layer-by-layer photonic crystals reveals the superior convergency of this different approach over the conventional plane-wave expansion method

Accurate expansion of cylindrical paraxial waves for its straightforward implementation in electromagnetic scattering

Science.gov (United States)

Naserpour, Mahin; Zapata-Rodríguez, Carlos J.

2018-01-01

The evaluation of vector wave fields can be accurately performed by means of diffraction integrals, differential equations and also series expansions. In this paper, a Bessel series expansion which basis relies on the exact solution of the Helmholtz equation in cylindrical coordinates is theoretically developed for the straightforward yet accurate description of low-numerical-aperture focal waves. The validity of this approach is confirmed by explicit application to Gaussian beams and apertured focused fields in the paraxial regime. Finally we discuss how our procedure can be favorably implemented in scattering problems.

The (G′/G-Expansion Method and Its Application for Higher-Order Equations of KdV (III

Directory of Open Access Journals (Sweden)

Huizhang Yang

2014-01-01

Full Text Available New exact traveling wave solutions of a higher-order KdV equation type are studied by the (G′/G-expansion method, where G=G(ξ satisfies a second-order linear differential equation. The traveling wave solutions are expressed by the hyperbolic functions, the trigonometric functions, and the rational functions. The property of this method is that it is quite simple and understandable.

New exact solutions of the(2+1-dimensional Broer-Kaup equation by the consistent Riccati expansion method

Directory of Open Access Journals (Sweden)

Jiang Ying

2017-01-01

Full Text Available In this work, we study the (2+1-D Broer-Kaup equation. The composite periodic breather wave, the exact composite kink breather wave and the solitary wave solutions are obtained by using the coupled degradation technique and the consistent Riccati expansion method. These results may help us to investigate some complex dynamical behaviors and the interaction between composite non-linear waves in high dimensional models

Disjoint sum expansion method in FTA

International Nuclear Information System (INIS)

Ruan Keqiang

1987-01-01

An expansion formula for transforming boolean algebraic expressions into disjoint form was proved. Based on this expansion formula, a method for transforming system failure function into disjoint form was devised. The fact that the expansion can be done for several elements simulatneously makes the method flexible and fast. Some examples from fault tree analysis (FTA) and network analysis were examined by the new method to show its algorithm and its merit. Besides, by means of the proved expansion formula some boolean algebraic relations can proved very easily

Symmetrized partial-wave method for density-functional cluster calculations

International Nuclear Information System (INIS)

Averill, F.W.; Painter, G.S.

1994-01-01

The computational advantage and accuracy of the Harris method is linked to the simplicity and adequacy of the reference-density model. In an earlier paper, we investigated one way the Harris functional could be extended to systems outside the limits of weakly interacting atoms by making the charge density of the interacting atoms self-consistent within the constraints of overlapping spherical atomic densities. In the present study, a method is presented for augmenting the interacting atom charge densities with symmetrized partial-wave expansions on each atomic site. The added variational freedom of the partial waves leads to a scheme capable of giving exact results within a given exchange-correlation approximation while maintaining many of the desirable convergence and stability properties of the original Harris method. Incorporation of the symmetry of the cluster in the partial-wave construction further reduces the level of computational effort. This partial-wave cluster method is illustrated by its application to the dimer C 2 , the hypothetical atomic cluster Fe 6 Al 8 , and the benzene molecule

The rate of beneficial mutations surfing on the wave of a range expansion.

Directory of Open Access Journals (Sweden)

Rémi Lehe

Full Text Available Many theoretical and experimental studies suggest that range expansions can have severe consequences for the gene pool of the expanding population. Due to strongly enhanced genetic drift at the advancing frontier, neutral and weakly deleterious mutations can reach large frequencies in the newly colonized regions, as if they were surfing the front of the range expansion. These findings raise the question of how frequently beneficial mutations successfully surf at shifting range margins, thereby promoting adaptation towards a range-expansion phenotype. Here, we use individual-based simulations to study the surfing statistics of recurrent beneficial mutations on wave-like range expansions in linear habitats. We show that the rate of surfing depends on two strongly antagonistic factors, the probability of surfing given the spatial location of a novel mutation and the rate of occurrence of mutations at that location. The surfing probability strongly increases towards the tip of the wave. Novel mutations are unlikely to surf unless they enjoy a spatial head start compared to the bulk of the population. The needed head start is shown to be proportional to the inverse fitness of the mutant type, and only weakly dependent on the carrying capacity. The precise location dependence of surfing probabilities is derived from the non-extinction probability of a branching process within a moving field of growth rates. The second factor is the mutation occurrence which strongly decreases towards the tip of the wave. Thus, most successful mutations arise at an intermediate position in the front of the wave. We present an analytic theory for the tradeoff between these factors that allows to predict how frequently substitutions by beneficial mutations occur at invasion fronts. We find that small amounts of genetic drift increase the fixation rate of beneficial mutations at the advancing front, and thus could be important for adaptation during species invasions.

The optimized expansion based low-rank method for wavefield extrapolation

KAUST Repository

Wu, Zedong

2014-03-01

Spectral methods are fast becoming an indispensable tool for wavefield extrapolation, especially in anisotropic media because it tends to be dispersion and artifact free as well as highly accurate when solving the wave equation. However, for inhomogeneous media, we face difficulties in dealing with the mixed space-wavenumber domain extrapolation operator efficiently. To solve this problem, we evaluated an optimized expansion method that can approximate this operator with a low-rank variable separation representation. The rank defines the number of inverse Fourier transforms for each time extrapolation step, and thus, the lower the rank, the faster the extrapolation. The method uses optimization instead of matrix decomposition to find the optimal wavenumbers and velocities needed to approximate the full operator with its explicit low-rank representation. As a result, we obtain lower rank representations compared with the standard low-rank method within reasonable accuracy and thus cheaper extrapolations. Additional bounds set on the range of propagated wavenumbers to adhere to the physical wave limits yield unconditionally stable extrapolations regardless of the time step. An application on the BP model provided superior results compared to those obtained using the decomposition approach. For transversely isotopic media, because we used the pure P-wave dispersion relation, we obtained solutions that were free of the shear wave artifacts, and the algorithm does not require that n > 0. In addition, the required rank for the optimization approach to obtain high accuracy in anisotropic media was lower than that obtained by the decomposition approach, and thus, it was more efficient. A reverse time migration result for the BP tilted transverse isotropy model using this method as a wave propagator demonstrated the ability of the algorithm.

Teaching Graphical Simulations of Fourier Series Expansion of Some Periodic Waves Using Spreadsheets

Science.gov (United States)

Singh, Iqbal; Kaur, Bikramjeet

2018-01-01

The present article demonstrates a way of programming using an Excel spreadsheet to teach Fourier series expansion in school/colleges without the knowledge of any typical programming language. By using this, a student learns to approximate partial sum of the n terms of Fourier series for some periodic signals such as square wave, saw tooth wave,…

Treatment of divergent expansions in scattering theory

International Nuclear Information System (INIS)

Gersten, A.; Malin, S.

1978-01-01

One of the biggest obstacles in applying quantum field theory to realistic scattering problems are the divergencies of pertubation expansions for large coupling constants and the divergencies of partial wave expansions for massless particles exchanges. There exist, however, methods of summation of the divergent expansions which can lead to significant application in physics. In this paper we treat the problem of summing such expansions using three methods: (i) a generalization of the Pade approximation to the multivariable case. The suggested definition is unique and preserves unitarity. (ii) The summation of divergent partial waves for arbitrary spins. (iii) A successful application of a series inversion to the 3 P 1 nucleon-nucleon phase shift up to 200 MeV. (orig./WL) [de

Range expansions transition from pulled to pushed waves with increasing cooperativity in an experimental microbial population

Science.gov (United States)

Gandhi, Saurabh; Yurtsev, Eugene; Korolev, Kirill; Gore, Jeff

Range expansions are becoming more frequent due to environmental changes and rare long distance dispersal, often facilitated by anthropogenic activities. Simple models in theoretical ecology explain many emergent properties of range expansions, such as a constant expansion velocity, in terms of organism-level properties such as growth and dispersal rates. Testing these quantitative predictions in natural populations is difficult because of large environmental variability. Here, we used a controlled microbial model system to study range expansions of populations with and without intra-specific cooperativity. For non-cooperative growth, the expansion dynamics were dominated by population growth at the low-density front, which pulled the expansion forward. We found these expansions to be in close quantitative agreement with the classical theory of pulled waves by Fisher and Skellam, suitably adapted to our experimental system. However, as cooperativity increased, the expansions transitioned to being pushed, i.e. controlled by growth in the bulk as well as in the front. Although both pulled and pushed waves expand at a constant velocity and appear otherwise similar, their distinct dynamics leads to very different evolutionary consequences. Given the prevalence of cooperative growth in nature, understanding the effects of cooperativity is essential to managing invading species and understanding their evolution.

Comparative study of the expansion dynamics of laser-driven plasma and shock wave in in-air and underwater ablation regimes

Science.gov (United States)

Nguyen, Thao T. P.; Tanabe, Rie; Ito, Yoshiro

2018-03-01

We compared the expansion characteristics of the plasma plumes and shock waves generated in laser-induced shock process between the two ablation regimes: in air and under water. The observation was made from the initial moment when the laser pulse hit the target until 1.5 μs. The shock processes were driven by focusing a single laser pulse (1064 nm, FWHM = 13 ns) onto the surface of epoxy-resin blocks using a 40-mm focal length lens. The estimated laser intensity at the target plane is approximate to 9 ×109Wcm-2 . We used the fast-imaging technique to observe the expansion of the plasma plume and a custom-designed time-resolved photoelasticity imaging technique to observe the propagation of shock waves with the time resolution of nanoseconds. We found that at the same intensity of the laser beam, the plasma expansion during the laser pulse follows different mechanisms: the plasma plume that grows in air follows a radiation-wave model while a detonation-wave model can explain the expansion of the plasma plume induced in water. The ideal blast wave theory can be used to predict the decay of the shock wave in air but is not appropriate to describe the decay of the shock wave induced under water.

Time evolution of the wave equation using rapid expansion method

KAUST Repository

Pestana, Reynam C.; Stoffa, Paul L.

2010-01-01

Forward modeling of seismic data and reverse time migration are based on the time evolution of wavefields. For the case of spatially varying velocity, we have worked on two approaches to evaluate the time evolution of seismic wavefields. An exact solution for the constant-velocity acoustic wave equation can be used to simulate the pressure response at any time. For a spatially varying velocity, a one-step method can be developed where no intermediate time responses are required. Using this approach, we have solved for the pressure response at intermediate times and have developed a recursive solution. The solution has a very high degree of accuracy and can be reduced to various finite-difference time-derivative methods, depending on the approximations used. Although the two approaches are closely related, each has advantages, depending on the problem being solved. © 2010 Society of Exploration Geophysicists.

Time evolution of the wave equation using rapid expansion method

KAUST Repository

Pestana, Reynam C.

2010-07-01

Forward modeling of seismic data and reverse time migration are based on the time evolution of wavefields. For the case of spatially varying velocity, we have worked on two approaches to evaluate the time evolution of seismic wavefields. An exact solution for the constant-velocity acoustic wave equation can be used to simulate the pressure response at any time. For a spatially varying velocity, a one-step method can be developed where no intermediate time responses are required. Using this approach, we have solved for the pressure response at intermediate times and have developed a recursive solution. The solution has a very high degree of accuracy and can be reduced to various finite-difference time-derivative methods, depending on the approximations used. Although the two approaches are closely related, each has advantages, depending on the problem being solved. © 2010 Society of Exploration Geophysicists.

Molecular wave function and effective adiabatic potentials calculated by extended multi-configuration time-dependent Hartree-Fock method

Energy Technology Data Exchange (ETDEWEB)

Kato, Tsuyoshi; Ide, Yoshihiro; Yamanouchi, Kaoru [Department of Chemistry, School of Science, The University of Tokyo, 7-3-1, Hongo Bunkyo-ku, Tokyo, 113-0033 (Japan)

2015-12-31

We first calculate the ground-state molecular wave function of 1D model H{sub 2} molecule by solving the coupled equations of motion formulated in the extended multi-configuration time-dependent Hartree-Fock (MCTDHF) method by the imaginary time propagation. From the comparisons with the results obtained by the Born-Huang (BH) expansion method as well as with the exact wave function, we observe that the memory size required in the extended MCTDHF method is about two orders of magnitude smaller than in the BH expansion method to achieve the same accuracy for the total energy. Second, in order to provide a theoretical means to understand dynamical behavior of the wave function, we propose to define effective adiabatic potential functions and compare them with the conventional adiabatic electronic potentials, although the notion of the adiabatic potentials is not used in the extended MCTDHF approach. From the comparison, we conclude that by calculating the effective potentials we may be able to predict the energy differences among electronic states even for a time-dependent system, e.g., time-dependent excitation energies, which would be difficult to be estimated within the BH expansion approach.

Ultrasonic wave propagation through aberrating layers: experimental verification of the conjugate gradient Rayleigh method

NARCIS (Netherlands)

Ledoux, L.A.F.; Berkhoff, Arthur P.; Thijssen, J.M.

The Conjugate Gradient Rayleigh method for the calculation of acoustic reflection and transmission at a rough interface between two media was experimentally verified. The method is based on a continuous version of the conjugate gradient technique and plane-wave expansions. We measured the beam

Multipole expansion of acoustical Bessel beams with arbitrary order and location.

Science.gov (United States)

Gong, Zhixiong; Marston, Philip L; Li, Wei; Chai, Yingbin

2017-06-01

An exact solution of expansion coefficients for a T-matrix method interacting with acoustic scattering of arbitrary order Bessel beams from an obstacle of arbitrary location is derived analytically. Because of the failure of the addition theorem for spherical harmonics for expansion coefficients of helicoidal Bessel beams, an addition theorem for cylindrical Bessel functions is introduced. Meanwhile, an analytical expression for the integral of products including Bessel and associated Legendre functions is applied to eliminate the integration over the polar angle. Note that this multipole expansion may also benefit other scattering methods and expansions of incident waves, for instance, partial-wave series solutions.

A ''quadratized'' augmented plane wave method

International Nuclear Information System (INIS)

Smrcka, L.

1982-02-01

The exact radial solution inside the muffin-tin sphere is replaced by its Taylor expansion with respect to the energy, truncated after the quadratic term. Making use of it the energy independent augmented plane waves are formed which lead to the secular equations linear in energy. The method resembles the currently used linearized APW method but yields higher accuracy. The analysis of solution inside one muffin-tin sphere shows that the eigenvalue error is proportional to (E-E 0 ) 6 as compared with (E-E 0 ) 4 for LAPW. The error of eigenfunctions is (E-E 0 ) 3 ((E-E 0 ) 2 for LAPW). These conclusions are confirmed by direct numerical calculation of band structure of Cu and Al. (author)

Nonlocal symmetries, solitary waves and cnoidal periodic waves of the (2+1)-dimensional breaking soliton equation

Science.gov (United States)

Zou, Li; Tian, Shou-Fu; Feng, Lian-Li

2017-12-01

In this paper, we consider the (2+1)-dimensional breaking soliton equation, which describes the interaction of a Riemann wave propagating along the y-axis with a long wave along the x-axis. By virtue of the truncated Painlevé expansion method, we obtain the nonlocal symmetry, Bäcklund transformation and Schwarzian form of the equation. Furthermore, by using the consistent Riccati expansion (CRE), we prove that the breaking soliton equation is solvable. Based on the consistent tan-function expansion, we explicitly derive the interaction solutions between solitary waves and cnoidal periodic waves.

Parametric form of QCD travelling waves

OpenAIRE

Peschanski, R.

2005-01-01

We derive parametric travelling-wave solutions of non-linear QCD equations. They describe the evolution towards saturation in the geometric scaling region. The method, based on an expansion in the inverse of the wave velocity, leads to a solvable hierarchy of differential equations. A universal parametric form of travelling waves emerges from the first two orders of the expansion.

Covariant spectator theory of $np$ scattering:\\\\ Effective range expansions and relativistic deuteron wave functions

Energy Technology Data Exchange (ETDEWEB)

Franz Gross, Alfred Stadler

2010-09-01

We present the effective range expansions for the 1S0 and 3S1 scattering phase shifts, and the relativistic deuteron wave functions that accompany our recent high precision fits (with \\chi^2/N{data} \\simeq 1) to the 2007 world np data below 350 MeV. The wave functions are expanded in a series of analytical functions (with the correct asymptotic behavior at both large and small arguments) that can be Fourier-transformed from momentum to coordinate space and are convenient to use in any application. A fortran subroutine to compute these wave functions can be obtained from the authors.

A new generalized expansion method and its application in finding explicit exact solutions for a generalized variable coefficients KdV equation

International Nuclear Information System (INIS)

Sabry, R.; Zahran, M.A.; Fan Engui

2004-01-01

A generalized expansion method is proposed to uniformly construct a series of exact solutions for general variable coefficients non-linear evolution equations. The new approach admits the following types of solutions (a) polynomial solutions, (b) exponential solutions, (c) rational solutions, (d) triangular periodic wave solutions, (e) hyperbolic and solitary wave solutions and (f) Jacobi and Weierstrass doubly periodic wave solutions. The efficiency of the method has been demonstrated by applying it to a generalized variable coefficients KdV equation. Then, new and rich variety of exact explicit solutions have been found

Chemical Kinetics in the expansion flow field of a rotating detonation-wave engine

Science.gov (United States)

Kailasanath, Kazhikathra; Schwer, Douglas

2014-11-01

Rotating detonation-wave engines (RDE) are a form of continuous detonation-wave engines. They potentially provide further gains in performance than an intermittent or pulsed detonation-wave engine (PDE). The overall flow field in an idealized RDE, primarily consisting of two concentric cylinders, has been discussed in previous meetings. Because of the high pressures involved and the lack of adequate reaction mechanisms for this regime, previous simulations have typically used simplified chemistry models. However, understanding the exhaust species concentrations in propulsion devices is important for both performance considerations as well as estimating pollutant emissions. A key step towards addressing this need will be discussed in this talk. In this approach, an induction parameter model is used for simulating the detonation but a more detailed finite-chemistry model is used in the expansion flow region, where the pressures are lower and the uncertainties in the chemistry model are greatly reduced. Results show that overall radical concentrations in the exhaust flow are substantially lower than from earlier predictions with simplified models. The performance of a baseline hydrogen/air RDE increased from 4940 s to 5000 s with the expansion flow chemistry, due to recombination of radicals and more production of H2O, resulting in additional heat release.

Partial wave expansions for arbitrary spin and the role of non-central forces

International Nuclear Information System (INIS)

Johnson, R.C.

1976-09-01

The partial wave expansion of the amplitudes used by Hooton and Johnson for the scattering of particles of arbitrary spin is derived. A discussion is given of the extent to which effects arising from transition matrix elements that are diagonal and nondiagonal in orbital angular momentum can be distinguished in observables

New extended (G'/G)-expansion method to solve nonlinear evolution equation: the (3 + 1)-dimensional potential-YTSF equation.

Science.gov (United States)

Roshid, Harun-Or-; Akbar, M Ali; Alam, Md Nur; Hoque, Md Fazlul; Rahman, Nizhum

2014-01-01

In this article, a new extended (G'/G) -expansion method has been proposed for constructing more general exact traveling wave solutions of nonlinear evolution equations with the aid of symbolic computation. In order to illustrate the validity and effectiveness of the method, we pick the (3 + 1)-dimensional potential-YTSF equation. As a result, abundant new and more general exact solutions have been achieved of this equation. It has been shown that the proposed method provides a powerful mathematical tool for solving nonlinear wave equations in applied mathematics, engineering and mathematical physics.

Partial wave expansions for arbitrary spin and the role of non-central forces

International Nuclear Information System (INIS)

Johnson, R.C.

1977-01-01

The partial wave expansion of the amplitudes used by Hooton and Johnson for the scattering of particles of arbitrary spin is derived. A discussion is given of the extent to which effects arising from transition matrix elements that are diagonal and non-diagonal in orbital angular momentum can be distinguished in observables. (Auth.)

Low-Frequency Waves in the Near-Earth Magnetotail before Substorm Expansion Onsets

Science.gov (United States)

Miyashita, Y.; Saito, M. H.; Hiraki, Y.; Machida, S.

2013-12-01

Magnetic reconnection and dipolarization, which occur in the near-Earth magnetotail just before substorm expansion onsets, are important processes for the substorm triggering. To understand the triggering of these processes, we have investigated low-frequency waves that were observed in the near-Earth magnetotail before onsets, by performing statistical analysis based on Geotail observations and case studies based on multi-point THEMIS and Geotail observations. Here we focused our examination on ~10 min interval before onsets. We find that small-amplitude Alfven and slow-mode magnetosonic waves with a period of ~1 to 2 min continuously exist for more than 10 min before onsets. Such waves are seen not only in the initial dipolarization region but also midway between the magnetic reconnection and initial dipolarization regions. It seems that the amplitudes of the waves are larger in the off-equator plasma sheet and the plasma sheet boundary layer than at the magnetic equator and in the lobe. After onsets the waves considerably amplify in the plasma sheet. These results may imply that instabilities already begin to grow gradually in a wide region during the substorm growth phase, while their explosive growth begins in localized regions just before onsets.

Solitary wave and periodic wave solutions for Burgers, Fisher ...

Indian Academy of Sciences (India)

Home; Journals; Pramana – Journal of Physics; Volume 85; Issue 1. Solitary wave and periodic wave solutions for Burgers, Fisher, Huxley and combined forms of these equations by the (′/)-expansion method. Jalil Manafian Mehrdad Lakestani. Volume 85 Issue 1 July 2015 pp 31-52 ...

Assessment of the further improved (G'/G)-expansion method and the extended tanh-method in probing exact solutions of nonlinear PDEs.

Science.gov (United States)

Akbar, M Ali; Ali, Norhashidah Hj Mohd; Mohyud-Din, Syed Tauseef

2013-01-01

The (G'/G)-expansion method is one of the most direct and effective method for obtaining exact solutions of nonlinear partial differential equations (PDEs). In the present article, we construct the exact traveling wave solutions of nonlinear evolution equations in mathematical physics via the (2 + 1)-dimensional breaking soliton equation by using two methods: namely, a further improved (G'/G)-expansion method, where G(ξ) satisfies the auxiliary ordinary differential equation (ODE) [G'(ξ)](2) = p G (2)(ξ) + q G (4)(ξ) + r G (6)(ξ); p, q and r are constants and the well known extended tanh-function method. We demonstrate, nevertheless some of the exact solutions bring out by these two methods are analogous, but they are not one and the same. It is worth mentioning that the first method has not been exercised anybody previously which gives further exact solutions than the second one. PACS numbers 02.30.Jr, 05.45.Yv, 02.30.Ik.

Modeling guided wave excitation in plates with surface mounted piezoelectric elements: coupled physics and normal mode expansion

Science.gov (United States)

Ren, Baiyang; Lissenden, Cliff J.

2018-04-01

Guided waves have been extensively studied and widely used for structural health monitoring because of their large volumetric coverage and good sensitivity to defects. Effectively and preferentially exciting a desired wave mode having good sensitivity to a certain defect is of great practical importance. Piezoelectric discs and plates are the most common types of surface-mounted transducers for guided wave excitation and reception. Their geometry strongly influences the proportioning between excited modes as well as the total power of the excited modes. It is highly desirable to predominantly excite the selected mode while the total transduction power is maximized. In this work, a fully coupled multi-physics finite element analysis, which incorporates the driving circuit, the piezoelectric element and the wave guide, is combined with the normal mode expansion method to study both the mode tuning and total wave power. The excitation of circular crested waves in an aluminum plate with circular piezoelectric discs is numerically studied for different disc and adhesive thicknesses. Additionally, the excitation of plane waves in an aluminum plate, using a stripe piezoelectric element is studied both numerically and experimentally. It is difficult to achieve predominant single mode excitation as well as maximum power transmission simultaneously, especially for higher order modes. However, guidelines for designing the geometry of piezoelectric elements for optimal mode excitation are recommended.

Correlation expansion: a powerful alternative multiple scattering calculation method

International Nuclear Information System (INIS)

Zhao Haifeng; Wu Ziyu; Sebilleau, Didier

2008-01-01

We introduce a powerful alternative expansion method to perform multiple scattering calculations. In contrast to standard MS series expansion, where the scattering contributions are grouped in terms of scattering order and may diverge in the low energy region, this expansion, called correlation expansion, partitions the scattering process into contributions from different small atom groups and converges at all energies. It converges faster than MS series expansion when the latter is convergent. Furthermore, it takes less memory than the full MS method so it can be used in the near edge region without any divergence problem, even for large clusters. The correlation expansion framework we derive here is very general and can serve to calculate all the elements of the scattering path operator matrix. Photoelectron diffraction calculations in a cluster containing 23 atoms are presented to test the method and compare it to full MS and standard MS series expansion

Punctual Pade approximants as a regularization procedure for divergent and oscillatory partial wave expansions of the scattering amplitude

International Nuclear Information System (INIS)

Garibotti, C.R.; Grinstein, F.F.

1978-01-01

Previous theorems on the convergence of the [n,n+m] punctual Pade approximants to the scattering amplitude are extended. The new proofs include the cases of nonforward and backward scattering corresponding to potentials having 1/r and 1/r 2 long-range behaviors, for which the partial wave expansions are divergent and oscillatory, respectively. In this way, the ability of the approximation scheme as a summation method is established for all of the long-range potentials of interest in potential scattering

Punctual Pade Approximants as a regularization procedure for divergent and oscillatory partial wave expansions of the scattering amplitude

International Nuclear Information System (INIS)

Garibotti, C.R.; Grinstein, F.F.

1978-01-01

Previous theorems on the convergence of the [n, n+m] Punctual Pade Approximants to the scattering amplitude are extended. The new proofs include the cases of non-forward and backward scattering corresponding to potentials having 1/r and 1/r 2 long range behaviours, for which the partial wave expansions are divergent and oscillatory, respectively. In this way, the ability of the approximation scheme as a summation method is established for all of the long range potentials of interest in potential scattering [pt

New exact solutions of (2 + 1)-dimensional Gardner equation via the new sine-Gordon equation expansion method

International Nuclear Information System (INIS)

Chen Yong; Yan Zhenya

2005-01-01

In this paper (2 + 1)-dimensional Gardner equation is investigated using a sine-Gordon equation expansion method, which was presented via a generalized sine-Gordon reduction equation and a new transformation. As a consequence, it is shown that the method is more powerful to obtain many types of new doubly periodic solutions of (2 + 1)-dimensional Gardner equation. In particular, solitary wave solutions are also given as simple limits of doubly periodic solutions

An improved wave rotor refrigerator using an outside gas flow for recycling the expansion work

Science.gov (United States)

Zhao, J.; Hu, D.

2017-03-01

To overcome the bottleneck of traditional gas wave refrigeration, an improved wave rotor refrigerator (WRR) cycle has been proposed, in which the expansion work was recycled during the process of refrigeration. Thermodynamic analysis of the two cycles shows that the refrigeration efficiency of the improved WRR cycle has been greatly increased compared with the traditional WRR. The performance of an improved WRR was investigated by adjusting the major operational parameters, such as the rotational speed of the wave rotor, port size, and inflow overpressure. The experimental results show that pressure loss can be reduced by nearly 40 % in this improved refrigeration system. Meanwhile, a two-dimensional numerical simulation was performed to understand the wave interactions that take place inside the rotor channels.

Exact soliton solutions of the generalized Gross-Pitaevskii equation based on expansion method

Directory of Open Access Journals (Sweden)

Ying Wang

2014-06-01

Full Text Available We give a more generalized treatment of the 1D generalized Gross-Pitaevskii equation (GGPE with variable term coefficients. External harmonic trapping potential is fully considered and the nonlinear interaction term is of arbitrary polytropic index of superfluid wave function. We also eliminate the interdependence between variable coefficients of the equation terms avoiding the restrictions that occur in some other works. The exact soliton solutions of the GGPE are obtained through the delicate combined utilization of modified lens-type transformation and F-expansion method with dominant features like soliton type properties highlighted.

Travelling wave solutions for a surface wave equation in fluid mechanics

Directory of Open Access Journals (Sweden)

Tian Yi

2016-01-01

Full Text Available This paper considers a non-linear wave equation arising in fluid mechanics. The exact traveling wave solutions of this equation are given by using G'/G-expansion method. This process can be reduced to solve a system of determining equations, which is large and difficult. To reduce this process, we used Wu elimination method. Example shows that this method is effective.

Temporal quadratic expansion nodal Green's function method

International Nuclear Information System (INIS)

Liu Cong; Jing Xingqing; Xu Xiaolin

2000-01-01

A new approach is presented to efficiently solve the three-dimensional space-time reactor dynamics equation which overcomes the disadvantages of current methods. In the Temporal Quadratic Expansion Nodal Green's Function Method (TQE/NGFM), the Quadratic Expansion Method (QEM) is used for the temporal solution with the Nodal Green's Function Method (NGFM) employed for the spatial solution. Test calculational results using TQE/NGFM show that its time step size can be 5-20 times larger than that of the Fully Implicit Method (FIM) for similar precision. Additionally, the spatial mesh size with NGFM can be nearly 20 times larger than that using the finite difference method. So, TQE/NGFM is proved to be an efficient reactor dynamics analysis method

New generalized and improved (G′/G-expansion method for nonlinear evolution equations in mathematical physics

Directory of Open Access Journals (Sweden)

Hasibun Naher

2014-10-01

Full Text Available In this article, new extension of the generalized and improved (G′/G-expansion method is proposed for constructing more general and a rich class of new exact traveling wave solutions of nonlinear evolution equations. To demonstrate the novelty and motivation of the proposed method, we implement it to the Korteweg-de Vries (KdV equation. The new method is oriented toward the ease of utilize and capability of computer algebraic system and provides a more systematic, convenient handling of the solution process of nonlinear equations. Further, obtained solutions disclose a wider range of applicability for handling a large variety of nonlinear partial differential equations.

Solitary wave and periodic wave solutions for Burgers, Fisher ...

Indian Academy of Sciences (India)

The generalized (G′/G)-expansion method; Burgers equation; Fisher's equation; ... the travelling wave solutions plays an important role in nonlinear sciences. ... Burgers, Fisher, Huxley equations and combined forms of these equations will ...

Dispersive traveling wave solutions of the Equal-Width and Modified Equal-Width equations via mathematical methods and its applications

Science.gov (United States)

Lu, Dianchen; Seadawy, Aly R.; Ali, Asghar

2018-06-01

The Equal-Width and Modified Equal-Width equations are used as a model in partial differential equations for the simulation of one-dimensional wave transmission in nonlinear media with dispersion processes. In this article we have employed extend simple equation method and the exp(-varphi(ξ)) expansion method to construct the exact traveling wave solutions of equal width and modified equal width equations. The obtained results are novel and have numerous applications in current areas of research in mathematical physics. It is exposed that our method, with the help of symbolic computation, provides a effective and powerful mathematical tool for solving different kind nonlinear wave problems.

Hadronic wave functions at short distances and the operator product expansion

International Nuclear Information System (INIS)

Brodsky, S.J.; Lepage, G.P.

1980-01-01

The operator product expansion, of appropriate products of quark fields, is used to find the anamalous dimensions which control the short distance behavior of hadronic wave functions. This vehavior in turn controls the high-Q 2 limit of hadronic form factors. In particular, we relate each anamalous dimension of the nonsinglet structure functions to a corresponding logarithmic correction factor to the nominal αsub(s)(Q 2 )/Q 2 fall off of meson form factors. Unlike the case of deep inelastic lepton-hadron scattering, the operator product necessary here involves extra terms which do not contribute to forward matrix elements. (orig.)

On analyticity of linear waves scattered by a layered medium

Science.gov (United States)

Nicholls, David P.

2017-10-01

The scattering of linear waves by periodic structures is a crucial phenomena in many branches of applied physics and engineering. In this paper we establish rigorous analytic results necessary for the proper numerical analysis of a class of High-Order Perturbation of Surfaces methods for simulating such waves. More specifically, we prove a theorem on existence and uniqueness of solutions to a system of partial differential equations which model the interaction of linear waves with a multiply layered periodic structure in three dimensions. This result provides hypotheses under which a rigorous numerical analysis could be conducted for recent generalizations to the methods of Operator Expansions, Field Expansions, and Transformed Field Expansions.

Application of the Method of Auxiliary Sources for the Analysis of Plane-Wave Scattering by Impedance Spheres

DEFF Research Database (Denmark)

Karamehmedovic, Mirza; Breinbjerg, Olav

2002-01-01

The Method of Auxiliary Sources (MAS) is applied to 3D scattering problems involving spherical impedance scatterers. The MAS results are compared with the reference spherical wave expansion (SWE) solution. It is demonstrated that good agreement is achieved between the MAS and SWE results....

Semiclassical expansions on and near caustics

International Nuclear Information System (INIS)

Meetz, K.

1984-09-01

We show that the standard WKB expansion can be generalized so that it reproduces the behavior of the wave function on and near a caustic in two-dimensional space time. The expansion is related to the unfolding polynomials of the elementary catastrophes occurring in two dimensions: the fold and the cusp catastrophe. The method determines control parameters and transport coefficients in a self-consistent way from differential equations and does not refer to the asymptotic expansion of Feynman path integrals. The lowest order equations are solved explicitly in terms of the multivalued classical action. The result is a generalized semiclassical approximation on and beyond a caustic. (orig.)

Chemical graph-theoretic cluster expansions

International Nuclear Information System (INIS)

Klein, D.J.

1986-01-01

A general computationally amenable chemico-graph-theoretic cluster expansion method is suggested as a paradigm for incorporation of chemical structure concepts in a systematic manner. The cluster expansion approach is presented in a formalism general enough to cover a variety of empirical, semiempirical, and even ab initio applications. Formally such approaches for the utilization of chemical structure-related concepts may be viewed as discrete analogues of Taylor series expansions. The efficacy of the chemical structure concepts then is simply bound up in the rate of convergence of the cluster expansions. In many empirical applications, e.g., boiling points, chromatographic separation coefficients, and biological activities, this rate of convergence has been observed to be quite rapid. More note will be made here of quantum chemical applications. Relations to questions concerning size extensivity of energies and size consistency of wave functions are addressed

Planar ESPAR Array Design with Nonsymmetrical Pattern by Means of Finite-Element Method, Domain Decomposition, and Spherical Wave Expansion

Directory of Open Access Journals (Sweden)

Jesús García

2012-01-01

Full Text Available The application of a 3D domain decomposition finite-element and spherical mode expansion for the design of planar ESPAR (electronically steerable passive array radiator made with probe-fed circular microstrip patches is presented in this work. A global generalized scattering matrix (GSM in terms of spherical modes is obtained analytically from the GSM of the isolated patches by using rotation and translation properties of spherical waves. The whole behaviour of the array is characterized including all the mutual coupling effects between its elements. This procedure has been firstly validated by analyzing an array of monopoles on a ground plane, and then it has been applied to synthesize a prescribed radiation pattern optimizing the reactive loads connected to the feeding ports of the array of circular patches by means of a genetic algorithm.

Self-force calculations with matched expansions and quasinormal mode sums

International Nuclear Information System (INIS)

Casals, Marc; Dolan, Sam; Ottewill, Adrian C.; Wardell, Barry

2009-01-01

Accurate modeling of gravitational wave emission by extreme-mass ratio inspirals is essential for their detection by the LISA mission. A leading perturbative approach involves the calculation of the self-force acting upon the smaller orbital body. In this work, we present the first application of the Poisson-Wiseman-Anderson method of 'matched expansions' to compute the self-force acting on a point particle moving in a curved spacetime. The method employs two expansions for the Green function, which are, respectively, valid in the 'quasilocal' and 'distant past' regimes, and which may be matched together within the normal neighborhood. We perform our calculation in a static region of the spherically symmetric Nariai spacetime (dS 2 xS 2 ), in which scalar-field perturbations are governed by a radial equation with a Poeschl-Teller potential (frequently used as an approximation to the Schwarzschild radial potential) whose solutions are known in closed form. The key new ingredients in our study are (i) very high order quasilocal expansions and (ii) expansion of the distant past Green function in quasinormal modes. In combination, these tools enable a detailed study of the properties of the scalar-field Green function. We demonstrate that the Green function is singular whenever x and x ' are connected by a null geodesic, and apply asymptotic methods to determine the structure of the Green function near the null wave front. We show that the singular part of the Green function undergoes a transition each time the null wave front passes through a caustic point, following a repeating fourfold sequence δ(σ), 1/πσ, -δ(σ), -1/πσ, etc., where σ is Synge's world function. The matched-expansion method provides insight into the nonlocal properties of the self-force. We show that the self-force generated by the segment of the worldline lying outside the normal neighborhood is not negligible. We apply the matched-expansion method to compute the scalar self-force acting on

Modal analysis of wave propagation in dispersive media

Science.gov (United States)

Abdelrahman, M. Ismail; Gralak, B.

2018-01-01

Surveys on wave propagation in dispersive media have been limited since the pioneering work of Sommerfeld [Ann. Phys. 349, 177 (1914), 10.1002/andp.19143491002] by the presence of branches in the integral expression of the wave function. In this article a method is proposed to eliminate these critical branches and hence to establish a modal expansion of the time-dependent wave function. The different components of the transient waves are physically interpreted as the contributions of distinct sets of modes and characterized accordingly. Then, the modal expansion is used to derive a modified analytical expression of the Sommerfeld precursor improving significantly the description of the amplitude and the oscillating period up to the arrival of the Brillouin precursor. The proposed method and results apply to all waves governed by the Helmholtz equations.

Cylindrical and spherical space equivalents to the plane wave expansion technique of Maxwell's wave equations

Science.gov (United States)

Gauthier, Robert C.; Alzahrani, Mohammed A.; Jafari, Seyed Hamed

2015-02-01

The plane wave expansion (PWM) technique applied to Maxwell's wave equations provides researchers with a supply of information regarding the optical properties of dielectric structures. The technique is well suited for structures that display a linear periodicity. When the focus is directed towards optical resonators and structures that lack linear periodicity the eigen-process can easily exceed computational resources and time constraints. In the case of dielectric structures which display cylindrical or spherical symmetry, a coordinate system specific set of basis functions have been employed to cast Maxwell's wave equations into an eigen-matrix formulation from which the resonator states associated with the dielectric profile can be obtained. As for PWM, the inverse of the dielectric and field components are expanded in the basis functions (Fourier-Fourier-Bessel, FFB, in cylindrical and Fourier- Bessel-Legendre, BLF, in spherical) and orthogonality is employed to form the matrix expressions. The theoretical development details will be presented indicating how certain mathematical complications in the process have been overcome and how the eigen-matrix can be tuned to a specific mode type. The similarities and differences in PWM, FFB and BLF are presented. In the case of structures possessing axial cylindrical symmetry, the inclusion of the z axis component of propagation constant makes the technique applicable to photonic crystal fibers and other waveguide structures. Computational results will be presented for a number of different dielectric geometries including Bragg ring resonators, cylindrical space slot channel waveguides and bottle resonators. Steps to further enhance the computation process will be reported.

New travelling wave solutions for nonlinear stochastic evolution ...

Indian Academy of Sciences (India)

expansion method to look for travelling wave solutions of nonlinear partial differential equations. It is interesting to mention that, in this method the sign of the parameters can be used to judge the numbers and types of travelling wave solutions.

New exact travelling wave solutions of nonlinear physical models

International Nuclear Information System (INIS)

Bekir, Ahmet; Cevikel, Adem C.

2009-01-01

In this work, we established abundant travelling wave solutions for some nonlinear evolution equations. This method was used to construct travelling wave solutions of nonlinear evolution equations. The travelling wave solutions are expressed by the hyperbolic functions, the trigonometric functions and the rational functions. The ((G ' )/G )-expansion method presents a wider applicability for handling nonlinear wave equations.

Optimized t-expansion method for the Rabi Hamiltonian

International Nuclear Information System (INIS)

Travenec, Igor; Samaj, Ladislav

2011-01-01

A polemic arose recently about the applicability of the t-expansion method to the calculation of the ground state energy E 0 of the Rabi model. For specific choices of the trial function and very large number of involved connected moments, the t-expansion results are rather poor and exhibit considerable oscillations. In this Letter, we formulate the t-expansion method for trial functions containing two free parameters which capture two exactly solvable limits of the Rabi Hamiltonian. At each order of the t-series, E 0 is assumed to be stationary with respect to the free parameters. A high accuracy of E 0 estimates is achieved for small numbers (5 or 6) of involved connected moments, the relative error being smaller than 10 -4 (0.01%) within the whole parameter space of the Rabi Hamiltonian. A special symmetrization of the trial function enables us to calculate also the first excited energy E 1 , with the relative error smaller than 10 -2 (1%). -- Highlights: → We study the ground state energy of the Rabi Hamiltonian. → We use the t-expansion method with an optimized trial function. → High accuracy of estimates is achieved, the relative error being smaller than 0.01%. → The calculation of the first excited state energy is made. The method has a general applicability.

A matched expansion approach to practical self-force calculations

International Nuclear Information System (INIS)

Anderson, Warren G; Wiseman, Alan G

2005-01-01

We discuss a practical method of computing the self-force on a particle moving through a curved spacetime. This method involves two expansions to calculate the self-force, one arising from the particle's immediate past and the other from the more distant past. The expansion in the immediate past is a covariant Taylor series and can be carried out for all geometries. The more distant expansion is a mode sum, and may be carried out in those cases where the wave equation for the field mediating the self-force admits a mode expansion of the solution. In particular, this method can be used to calculate the gravitational self-force for a particle of mass μ orbiting a black hole of mass M to order μ 2 , provided μ/M << 1. We discuss how to use these two expansions to construct a full self-force, and in particular investigate criteria for matching the two expansions. As with all methods of computing self-forces for particles moving in black hole spacetimes, one encounters considerable technical difficulty in applying this method; nevertheless, it appears that the convergence of each series is good enough that a practical implementation may be plausible

The development of shock wave overpressure driven by channel expansion of high current impulse discharge arc

Science.gov (United States)

Xiong, Jia-ming; Li, Lee; Dai, Hong-yu; Wu, Hai-bo; Peng, Ming-yang; Lin, Fu-chang

2018-03-01

During the formation of a high current impulse discharge arc, objects near the discharge arc will be strongly impacted. In this paper, a high power, high current gas switch is used as the site of the impulse discharge arc. The explosion wave theory and the arc channel energy balance equation are introduced to analyze the development of the shock wave overpressure driven by the high current impulse discharge arc, and the demarcation point of the arc channel is given, from which the energy of the arc channel is no longer converted into shock waves. Through the analysis and calculation, it is found that the magnitude of the shock wave overpressure caused by impulse discharge arc expansion is closely related to the arc current rising rate. The arc shock wave overpressure will undergo a slow decay process and then decay rapidly. The study of this paper will perform the function of deepening the understanding of the physical nature of the impulse arc discharge, which can be used to explain the damage effect of the high current impulse discharge arc.

IRP methods for Environmental Impact Statements of utility expansion plans

International Nuclear Information System (INIS)

Cavallo, J.D.; Hemphill, R.C.; Veselka, T.D.

1992-01-01

Most large electric utilities and a growing number of gas utilities in the United States are using a planning method -- Integrated Resource Planning (IRP) - which incorporates demand-side management (DSM) programs whenever the marginal cost of the DSM programs are lower than the marginal cost of supply-side expansion options. Argonne National Laboratory has applied the IRP method in its socio-economic analysis of an Environmental Impact Statement (EIS) of power marketing for a system of electric utilities in the mountain and western regions of the United States. Applying the IRP methods provides valuable information to the participants in an EIS process involving capacity expansion of an electric or gas utility. The major challenges of applying the IRP method within an EIS are the time consuming and costly task of developing a least cost expansion path for each altemative, the detailed quantification of environmental damages associated with capacity expansion, and the explicit inclusion of societal-impacts to the region

2D XXZ model ground state properties using an analytic Lanczos expansion

International Nuclear Information System (INIS)

Witte, N.S.; Hollenberg, L.C.L.; Weihong Zheng

1997-01-01

A formalism was developed for calculating arbitrary expectation values for any extensive lattice Hamiltonian system using a new analytic Lanczos expansion, or plaquette expansion, and a recently proved exact theorem for ground state energies. The ground state energy, staggered magnetisation and the excited state gap of the 2D anisotropic antiferromagnetic Heisenberg Model are then calculated using this expansion for a range of anisotropy parameters and compared to other moment based techniques, such as the t-expansion, and spin-wave theory and series expansion methods. It was found that far from the isotropic point all moment methods give essentially very similar results, but near the isotopic point the plaquette expansion is generally better than the others. 20 refs., 6 tabs

Expansion Under Climate Change: The Genetic Consequences.

Science.gov (United States)

Garnier, Jimmy; Lewis, Mark A

2016-11-01

Range expansion and range shifts are crucial population responses to climate change. Genetic consequences are not well understood but are clearly coupled to ecological dynamics that, in turn, are driven by shifting climate conditions. We model a population with a deterministic reaction-diffusion model coupled to a heterogeneous environment that develops in time due to climate change. We decompose the resulting travelling wave solution into neutral genetic components to analyse the spatio-temporal dynamics of its genetic structure. Our analysis shows that range expansions and range shifts under slow climate change preserve genetic diversity. This is because slow climate change creates range boundaries that promote spatial mixing of genetic components. Mathematically, the mixing leads to so-called pushed travelling wave solutions. This mixing phenomenon is not seen in spatially homogeneous environments, where range expansion reduces genetic diversity through gene surfing arising from pulled travelling wave solutions. However, the preservation of diversity is diminished when climate change occurs too quickly. Using diversity indices, we show that fast expansions and range shifts erode genetic diversity more than slow range expansions and range shifts. Our study provides analytical insight into the dynamics of travelling wave solutions in heterogeneous environments.

Three-dimensional static and dynamic reactor calculations by the nodal expansion method

International Nuclear Information System (INIS)

Christensen, B.

1985-05-01

This report reviews various method for the calculation of the neutron-flux- and power distribution in an nuclear reactor. The nodal expansion method (NEM) is especially described in much detail. The nodal expansion method solves the diffusion equation. In this method the reactor core is divided into nodes, typically 10 to 20 cm in each direction, and the average flux in each node is calculated. To obtain the coupling between the nodes the local flux inside each node is expressed by use of a polynomial expansion. The expansion is one-dimensional, so inside each node such three expansions occur. To calculate the expansion coefficients it is necessary that the polynomial expansion is a solution to the one-dimensional diffusion equation. When the one-dimensional diffusion equation is established a term with the transversal leakage occur, and this term is expanded after the same polynomials. The resulting equation system with the expansion coefficients as the unknowns is solved with weigthed residual technique. The nodal expansion method is built into a computer program (also called NEM), which is divided into two parts, one part for steady-state calculations and one part for dynamic calculations. It is possible to take advantage of symmetry properties of the reactor core. The program is very flexible with regard to the number of energy groups, the node size, the flux expansion order and the transverse leakage expansion order. The boundary of the core is described by albedos. The program and input to it are described. The program is tested on a number of examples extending from small theoretical one up to realistic reactor cores. Many calculations are done on the wellknown IAEA benchmark case. The calculations have tested the accuracy and the computing time for various node sizes and polynomial expansions. In the dynamic examples various strategies for variation of the time step-length have been tested. (author)

Applications of exact traveling wave solutions of Modified Liouville and the Symmetric Regularized Long Wave equations via two new techniques

Science.gov (United States)

Lu, Dianchen; Seadawy, Aly R.; Ali, Asghar

2018-06-01

In this current work, we employ novel methods to find the exact travelling wave solutions of Modified Liouville equation and the Symmetric Regularized Long Wave equation, which are called extended simple equation and exp(-Ψ(ξ))-expansion methods. By assigning the different values to the parameters, different types of the solitary wave solutions are derived from the exact traveling wave solutions, which shows the efficiency and precision of our methods. Some solutions have been represented by graphical. The obtained results have several applications in physical science.

Dispersive optical soliton solutions for the hyperbolic and cubic-quintic nonlinear Schrödinger equations via the extended sinh-Gordon equation expansion method

Science.gov (United States)

Seadawy, Aly R.; Kumar, Dipankar; Chakrabarty, Anuz Kumar

2018-05-01

The (2+1)-dimensional hyperbolic and cubic-quintic nonlinear Schrödinger equations describe the propagation of ultra-short pulses in optical fibers of nonlinear media. By using an extended sinh-Gordon equation expansion method, some new complex hyperbolic and trigonometric functions prototype solutions for two nonlinear Schrödinger equations were derived. The acquired new complex hyperbolic and trigonometric solutions are expressed by dark, bright, combined dark-bright, singular and combined singular solitons. The obtained results are more compatible than those of other applied methods. The extended sinh-Gordon equation expansion method is a more powerful and robust mathematical tool for generating new optical solitary wave solutions for many other nonlinear evolution equations arising in the propagation of optical pulses.

Wave reflections from breakwaters

OpenAIRE

Dickson, William S.

1994-01-01

A new method is presented for estimating the reflection of a random, multi-directional sea from a coastal structure. The technique is applicable to an array of wave gauges of arbitrary geometry deployed seaward of the reflector. An expansion for small oblique wave incidence angles is used to derive an approximate relationship between measured array cross-spectra and a small number of parameters that describe the incident wave properties and the reflectivity of the structure. Model tests with ...

Adiabatic supernova expansion into the circumstellar medium

International Nuclear Information System (INIS)

Band, D.L.; Liang, E.P.

1987-01-01

We perform one dimensional numerical simulations with a Lagrangian hydrodynamics code of the adiabatic expansion of a supernova into the surrounding medium. The early expansion follows Chevalier's analytic self-similar solution until the reverse shock reaches the ejecta core. We follow the expansion as it evolves towards the adiabatic blast wave phase. Some memory of the earlier phases of expansion is retained in the interior even when the outer regions expand as a blast wave. We find the results are sensitive to the initial configuration of the ejecta and to the placement of gridpoints. 6 refs., 2 figs

The Green-function transform and wave propagation

Directory of Open Access Journals (Sweden)

Colin eSheppard

2014-11-01

Full Text Available Fourier methods well known in signal processing are applied to three-dimensional wave propagation problems. The Fourier transform of the Green function, when written explicitly in terms of a real-valued spatial frequency, consists of homogeneous and inhomogeneous components. Both parts are necessary to result in a pure out-going wave that satisfies causality. The homogeneous component consists only of propagating waves, but the inhomogeneous component contains both evanescent and propagating terms. Thus we make a distinction between inhomogeneous waves and evanescent waves. The evanescent component is completely contained in the region of the inhomogeneous component outside the k-space sphere. Further, propagating waves in the Weyl expansion contain both homogeneous and inhomogeneous components. The connection between the Whittaker and Weyl expansions is discussed. A list of relevant spherically symmetric Fourier transforms is given.

Unconditionally stable WLP-FDTD method for the modeling of electromagnetic wave propagation in gyrotropic materials.

Science.gov (United States)

Li, Zheng-Wei; Xi, Xiao-Li; Zhang, Jin-Sheng; Liu, Jiang-fan

2015-12-14

The unconditional stable finite-difference time-domain (FDTD) method based on field expansion with weighted Laguerre polynomials (WLPs) is applied to model electromagnetic wave propagation in gyrotropic materials. The conventional Yee cell is modified to have the tightly coupled current density components located at the same spatial position. The perfectly matched layer (PML) is formulated in a stretched-coordinate (SC) system with the complex-frequency-shifted (CFS) factor to achieve good absorption performance. Numerical examples are shown to validate the accuracy and efficiency of the proposed method.

Application of potential harmonic expansion method to BEC ...

Indian Academy of Sciences (India)

We adopt the potential harmonics expansion method for an ab initio solu- ... commonly adopted mean-field theories, our method is capable of handling ..... potentials in self-consistent mean-field calculation [7] gives wrong results as the.

High-frequency Rayleigh-wave method

Science.gov (United States)

Xia, J.; Miller, R.D.; Xu, Y.; Luo, Y.; Chen, C.; Liu, J.; Ivanov, J.; Zeng, C.

2009-01-01

High-frequency (???2 Hz) Rayleigh-wave data acquired with a multichannel recording system have been utilized to determine shear (S)-wave velocities in near-surface geophysics since the early 1980s. This overview article discusses the main research results of high-frequency surface-wave techniques achieved by research groups at the Kansas Geological Survey and China University of Geosciences in the last 15 years. The multichannel analysis of surface wave (MASW) method is a non-invasive acoustic approach to estimate near-surface S-wave velocity. The differences between MASW results and direct borehole measurements are approximately 15% or less and random. Studies show that simultaneous inversion with higher modes and the fundamental mode can increase model resolution and an investigation depth. The other important seismic property, quality factor (Q), can also be estimated with the MASW method by inverting attenuation coefficients of Rayleigh waves. An inverted model (S-wave velocity or Q) obtained using a damped least-squares method can be assessed by an optimal damping vector in a vicinity of the inverted model determined by an objective function, which is the trace of a weighted sum of model-resolution and model-covariance matrices. Current developments include modeling high-frequency Rayleigh-waves in near-surface media, which builds a foundation for shallow seismic or Rayleigh-wave inversion in the time-offset domain; imaging dispersive energy with high resolution in the frequency-velocity domain and possibly with data in an arbitrary acquisition geometry, which opens a door for 3D surface-wave techniques; and successfully separating surface-wave modes, which provides a valuable tool to perform S-wave velocity profiling with high-horizontal resolution. ?? China University of Geosciences (Wuhan) and Springer-Verlag GmbH 2009.

expansion method

Indian Academy of Sciences (India)

of a system under investigation is to model the system in terms of some ... The organization of the paper is as follows: In §2, a brief account of the (G /G)- expansion ...... It is interesting to note that from the general results, one can easily recover.

Discrete expansions of continuum functions. General concepts

International Nuclear Information System (INIS)

Bang, J.; Ershov, S.N.; Gareev, F.A.; Kazacha, G.S.

1979-01-01

Different discrete expansions of the continuum wave functions are considered: pole expansion (according to the Mittag-Lefler theorem), Weinberg states. The general property of these groups of states is their completeness in the finite region of space. They satisfy the Schroedinger type equations and are matched with free solutions of the Schroedinger equation at the boundary. Convergence of expansions for the S matrix, the Green functions and the continuous-spectrum wave functions is studied. A new group of states possessing the best convergence is introduced

Series solutions to partial differential equations. A study of the singularities, expansions, and solutions of Schroedinger's equation for the helium atom

International Nuclear Information System (INIS)

Mahlab, M.S.

1975-01-01

All the presently available techniques for solving Schroedinger's differential equation for helium-like atoms display poor convergence of the wave function in the neighborhood of the singularities of the Hamiltonian operator. In general most of the methods of solving this equation will converge in the appropriate limit to the exact wave function; however, convergence is slow, especially near the singularities of this differential equation. These difficulties become readily apparent from local energy studies. A technique is presented that avoids these difficulties. The wave function it produces is specifically most accurate at the singularities of the Hamiltonian. The novel aspect of this treatment is the subdivision of the space spanned by the wave function. Different expansions are picked such that they converge rapidly in each of the different subdivisions. These expansions may be constructed in such a way that they obey the boundary conditions in their respective subdivision. Most importantly, all the information available from the recursion relations associated with the differential equation may be incorporated into these expansions. A systematic procedure is presented such that these expansions may be brought together to form a wave function that satisfies all the continuity requirements. An S-state helium wave function was constructed to demonstrate that this method of treatment is feasible, and capable of indefinite systematic improvement. A discussion of several new asymptotic expansions that were constructed for the helium wave function, as well as an improved functional form for the small electron-nucleus wave function, is included in this presentation

A two-dimensional, semi-analytic expansion method for nodal calculations

International Nuclear Information System (INIS)

Palmtag, S.P.

1995-08-01

Most modern nodal methods used today are based upon the transverse integration procedure in which the multi-dimensional flux shape is integrated over the transverse directions in order to produce a set of coupled one-dimensional flux shapes. The one-dimensional flux shapes are then solved either analytically or by representing the flux shape by a finite polynomial expansion. While these methods have been verified for most light-water reactor applications, they have been found to have difficulty predicting the large thermal flux gradients near the interfaces of highly-enriched MOX fuel assemblies. A new method is presented here in which the neutron flux is represented by a non-seperable, two-dimensional, semi-analytic flux expansion. The main features of this method are (1) the leakage terms from the node are modeled explicitly and therefore, the transverse integration procedure is not used, (2) the corner point flux values for each node are directly edited from the solution method, and a corner-point interpolation is not needed in the flux reconstruction, (3) the thermal flux expansion contains hyperbolic terms representing analytic solutions to the thermal flux diffusion equation, and (4) the thermal flux expansion contains a thermal to fast flux ratio term which reduces the number of polynomial expansion functions needed to represent the thermal flux. This new nodal method has been incorporated into the computer code COLOR2G and has been used to solve a two-dimensional, two-group colorset problem containing uranium and highly-enriched MOX fuel assemblies. The results from this calculation are compared to the results found using a code based on the traditional transverse integration procedure

Experiences using DAKOTA stochastic expansion methods in computational simulations.

Energy Technology Data Exchange (ETDEWEB)

Templeton, Jeremy Alan; Ruthruff, Joseph R.

2012-01-01

Uncertainty quantification (UQ) methods bring rigorous statistical connections to the analysis of computational and experiment data, and provide a basis for probabilistically assessing margins associated with safety and reliability. The DAKOTA toolkit developed at Sandia National Laboratories implements a number of UQ methods, which are being increasingly adopted by modeling and simulation teams to facilitate these analyses. This report disseminates results as to the performance of DAKOTA's stochastic expansion methods for UQ on a representative application. Our results provide a number of insights that may be of interest to future users of these methods, including the behavior of the methods in estimating responses at varying probability levels, and the expansion levels for the methodologies that may be needed to achieve convergence.

Mass spectra and wave functions of meson systems and the covariant oscillator quark model as an expansion basis

International Nuclear Information System (INIS)

Oda, Ryuichi; Ishida, Shin; Wada, Hiroaki; Yamada, Kenji; Sekiguchi, Motoo

1999-01-01

We examine mass spectra and wave functions of the nn-bar, cc-bar and bb-bar meson systems within the framework of the covariant oscillator quark model with the boosted LS-coupling scheme. We solve nonperturbatively an eigenvalue problem for the squared-mass operator, which incorporates the four-dimensional color-Coulomb-type interaction, by taking a set of covariant oscillator wave functions as an expansion basis. We obtain mass spectra of these meson systems, which reproduce quite well their experimental behavior. The resultant manifestly covariant wave functions, which are applicable to analyses of various reaction phenomena, are given. Our results seem to suggest that the present model may be considered effectively as a covariant version of the nonrelativistic linear-plus-Coulomb potential quark model. (author)

The optimizied expansion method for wavefield extrapolation

KAUST Repository

Wu, Zedong; Alkhalifah, Tariq Ali

2013-01-01

, for inhomogeneous media, we face difficulties in dealing with the mixed space-wavenumber domain operator.In this abstract, we propose an optimized expansion method that can approximate this operator with its low rank representation. The rank defines the number

The exact rogue wave recurrence in the NLS periodic setting via matched asymptotic expansions, for 1 and 2 unstable modes

Science.gov (United States)

Grinevich, P. G.; Santini, P. M.

2018-04-01

The focusing Nonlinear Schrödinger (NLS) equation is the simplest universal model describing the modulation instability (MI) of quasi monochromatic waves in weakly nonlinear media, the main physical mechanism for the generation of rogue (anomalous) waves (RWs) in Nature. In this paper we investigate the x-periodic Cauchy problem for NLS for a generic periodic initial perturbation of the unstable constant background solution, in the case of N = 1 , 2 unstable modes. We use matched asymptotic expansion techniques to show that the solution of this problem describes an exact deterministic alternate recurrence of linear and nonlinear stages of MI, and that the nonlinear RW stages are described by the N-breather solution of Akhmediev type, whose parameters, different at each RW appearance, are always given in terms of the initial data through elementary functions. This paper is motivated by a preceding work of the authors in which a different approach, the finite gap method, was used to investigate periodic Cauchy problems giving rise to RW recurrence.

Analytical method for estimating the thermal expansion coefficient of metals at high temperature

International Nuclear Information System (INIS)

Takamoto, S; Izumi, S; Nakata, T; Sakai, S; Oinuma, S; Nakatani, Y

2015-01-01

In this paper, we propose an analytical method for estimating the thermal expansion coefficient (TEC) of metals at high-temperature ranges. Although the conventional method based on quasiharmonic approximation (QHA) shows good results at low temperatures, anharmonic effects caused by large-amplitude thermal vibrations reduces its accuracy at high temperatures. Molecular dynamics (MD) naturally includes the anharmonic effect. However, since the computational cost of MD is relatively high, in order to make an interatomic potential capable of reproducing TEC, an analytical method is essential. In our method, analytical formulation of the radial distribution function (RDF) at finite temperature realizes the estimation of the TEC. Each peak of the RDF is approximated by the Gaussian distribution. The average and variance of the Gaussian distribution are formulated by decomposing the fluctuation of interatomic distance into independent elastic waves. We incorporated two significant anharmonic effects into the method. One is the increase in the averaged interatomic distance caused by large amplitude vibration. The second is the variation in the frequency of elastic waves. As a result, the TECs of fcc and bcc crystals estimated by our method show good agreement with those of MD. Our method enables us to make an interatomic potential that reproduces the TEC at high temperature. We developed the GEAM potential for nickel. The TEC of the fitted potential showed good agreement with experimental data from room temperature to 1000 K. As compared with the original potential, it was found that the third derivative of the wide-range curve was modified, while the zeroth, first and second derivatives were unchanged. This result supports the conventional theory of solid state physics. We believe our analytical method and developed interatomic potential will contribute to future high-temperature material development. (paper)

The G′G-expansion method using modified Riemann–Liouville derivative for some space-time fractional differential equations

Directory of Open Access Journals (Sweden)

Ahmet Bekir

2014-09-01

Full Text Available In this paper, the fractional partial differential equations are defined by modified Riemann–Liouville fractional derivative. With the help of fractional derivative and traveling wave transformation, these equations can be converted into the nonlinear nonfractional ordinary differential equations. Then G′G-expansion method is applied to obtain exact solutions of the space-time fractional Burgers equation, the space-time fractional KdV-Burgers equation and the space-time fractional coupled Burgers’ equations. As a result, many exact solutions are obtained including hyperbolic function solutions, trigonometric function solutions and rational solutions. These results reveal that the proposed method is very effective and simple in performing a solution to the fractional partial differential equation.

Application of potential harmonic expansion method to BEC

Indian Academy of Sciences (India)

We adopt the potential harmonics expansion method for an ab initio solution of the many-body system in a Bose condensate containing interacting bosons. Unlike commonly adopted mean-field theories, our method is capable of handling two-body correlation properly. We disregard three- and higher-body correlations.

Design of materials with extreme thermal expansion using a three-phase topology optimization method

DEFF Research Database (Denmark)

Sigmund, Ole; Torquato, S.

1997-01-01

Composites with extremal or unusual thermal expansion coefficients are designed using a three-phase topology optimization method. The composites are made of two different material phases and a void phase. The topology optimization method consists in finding the distribution of material phases...... materials having maximum directional thermal expansion (thermal actuators), zero isotropic thermal expansion, and negative isotropic thermal expansion. It is shown that materials with effective negative thermal expansion coefficients can be obtained by mixing two phases with positive thermal expansion...

Optoheterodyne Doppler measurements of the ballistic expansion of the products of the shock wave-induced surface destruction: Experiment and theory

International Nuclear Information System (INIS)

Andriyash, A. V.; Astashkin, M. V.; Baranov, V. K.; Golubinskii, A. G.; Irinichev, D. A.; Kondrat’ev, A. N.; Kuratov, S. E.; Mazanov, V. A.; Rogozkin, D. B.; Stepushkin, S. N.; Khatunkin, V. Yu.

2016-01-01

The results of optoheterodyne Doppler measurements of the ballistic expansion of the products of surface destruction under shock-wave loading are presented. The possibility of determining the physical characteristics of a rapidly flying dust cloud, including the microparticle velocities, the microparticle sizes, and the areal density of the dust cloud, is shown. A compact stand for performing experiments on shock-wave loading of metallic samples is described. Shock-wave loading is performed by a 100-µm-thick tantalum flyer plate accelerated to a velocity of 2.8 km/s. As the samples, lead plates having various thicknesses and the same surface roughness are used. At a shock-wave pressure of 31.5 GPa, the destruction products are solid microparticles about 50 µm in size. At a pressure of 42 and 88 GPa, a liquid-drop dust cloud with a particle size of 10–15 µm is formed. To interpret the spectral data on the optoheterodyne Doppler measurements of the expansion of the surface destruction products (spalled fragments, dust microparticles), a transport equation for the function of mutual coherence of a multiply scattered field is used. The Doppler spectra of a backscattered signal are calculated with the model developed for the dust cloud that appears when a shock wave reaches the sample surface at the parameters that are typical of an experimental situation. Qualitative changes are found in the spectra, depending on the optical thickness of the dust cloud. The obtained theoretical results are in agreement with the experimental data.

Optoheterodyne Doppler measurements of the ballistic expansion of the products of the shock wave-induced surface destruction: Experiment and theory

Energy Technology Data Exchange (ETDEWEB)

Andriyash, A. V. [All-Russia Research Institute of Automatics (Russian Federation); Astashkin, M. V.; Baranov, V. K.; Golubinskii, A. G.; Irinichev, D. A. [Russian Federal Nuclear Center, All-Russia Research Institute of Experimental Physics (VNIIEF) (Russian Federation); Kondrat’ev, A. N., E-mail: an.kondratev@physics.msu.ru; Kuratov, S. E. [All-Russia Research Institute of Automatics (Russian Federation); Mazanov, V. A. [Russian Federal Nuclear Center, All-Russia Research Institute of Experimental Physics (VNIIEF) (Russian Federation); Rogozkin, D. B. [All-Russia Research Institute of Automatics (Russian Federation); Stepushkin, S. N.; Khatunkin, V. Yu. [Russian Federal Nuclear Center, All-Russia Research Institute of Experimental Physics (VNIIEF) (Russian Federation)

2016-06-15

The results of optoheterodyne Doppler measurements of the ballistic expansion of the products of surface destruction under shock-wave loading are presented. The possibility of determining the physical characteristics of a rapidly flying dust cloud, including the microparticle velocities, the microparticle sizes, and the areal density of the dust cloud, is shown. A compact stand for performing experiments on shock-wave loading of metallic samples is described. Shock-wave loading is performed by a 100-µm-thick tantalum flyer plate accelerated to a velocity of 2.8 km/s. As the samples, lead plates having various thicknesses and the same surface roughness are used. At a shock-wave pressure of 31.5 GPa, the destruction products are solid microparticles about 50 µm in size. At a pressure of 42 and 88 GPa, a liquid-drop dust cloud with a particle size of 10–15 µm is formed. To interpret the spectral data on the optoheterodyne Doppler measurements of the expansion of the surface destruction products (spalled fragments, dust microparticles), a transport equation for the function of mutual coherence of a multiply scattered field is used. The Doppler spectra of a backscattered signal are calculated with the model developed for the dust cloud that appears when a shock wave reaches the sample surface at the parameters that are typical of an experimental situation. Qualitative changes are found in the spectra, depending on the optical thickness of the dust cloud. The obtained theoretical results are in agreement with the experimental data.

Wave energy patterns of counterpulsation: a novel approach with wave intensity analysis.

Science.gov (United States)

Lu, Pong-Jeu; Yang, Chi-Fu Jeffrey; Wu, Meng-Yu; Hung, Chun-Hao; Chan, Ming-Yao; Hsu, Tzu-Cheng

2011-11-01

In counterpulsation, diastolic augmentation increases coronary blood flow and systolic unloading reduces left ventricular afterload. We present a new approach with wave intensity analysis to revisit and explain counterpulsation principles. In an acute porcine model, a standard intra-aortic balloon pump was placed in descending aorta in 4 pigs. We measured pressure and velocity with probes in left anterior descending artery and aorta during and without intra-aortic balloon pump assistance. Wave intensities of aortic and left coronary waves were derived from pressure and flow measurements with synchronization correction. We identified predominating waves in counterpulsation. In the aorta, during diastolic augmentation, intra-aortic balloon inflation generated a backward compression wave, with a "pushing" effect toward the aortic root that translated to a forward compression wave into coronary circulation. During systolic unloading, intra-aortic balloon pump deflation generated a backward expansion wave that "sucked" blood from left coronary bed into the aorta. While this backward expansion wave translated to reduced left ventricular afterload, the "sucking" effect resulted in left coronary blood steal, as demonstrated by a forward expansion wave in left anterior descending coronary flow. The waves were sensitive to inflation and deflation timing, with just 25 ms delay from standard deflation timing leading to weaker forward expansion wave and less coronary regurgitation. Intra-aortic balloon pumps generate backward-traveling waves that predominantly drive aortic and coronary blood flow during counterpulsation. Wave intensity analysis of arterial circulations may provide a mechanism to explain diastolic augmentation and systolic unloading of intra-aortic balloon pump counterpulsation. Copyright © 2011 The American Association for Thoracic Surgery. Published by Mosby, Inc. All rights reserved.

Heuristic method for determining outgoing waves in many-body wave functions

International Nuclear Information System (INIS)

Redish, E.F.; Tandy, P.C.; L'Huillier, M.

1975-12-01

A new and simple method is proposed for determining the kinds of outgoing waves present in a given many-body wave function. Whether any particular wave function contains ''hidden'' rearrangement components can be determined. 1 figure

Closed form solutions of two time fractional nonlinear wave equations

Directory of Open Access Journals (Sweden)

M. Ali Akbar

2018-06-01

Full Text Available In this article, we investigate the exact traveling wave solutions of two nonlinear time fractional wave equations. The fractional derivatives are described in the sense of conformable fractional derivatives. In addition, the traveling wave solutions are accomplished in the form of hyperbolic, trigonometric, and rational functions involving free parameters. To investigate such types of solutions, we implement the new generalized (G′/G-expansion method. The extracted solutions are reliable, useful and suitable to comprehend the optimal control problems, chaotic vibrations, global and local bifurcations and resonances, furthermore, fission and fusion phenomena occur in solitons, the relativistic energy-momentum relation, scalar electrodynamics, quantum relativistic one-particle theory, electromagnetic interactions etc. The results reveal that the method is very fruitful and convenient for exploring nonlinear differential equations of fractional order treated in theoretical physics. Keywords: Traveling wave solution, Soliton, Generalized (G′/G-expansion method, Time fractional Duffing equation, Time fractional Riccati equation

Amplitude reconstruction from complete photoproduction experiments and truncated partial-wave expansions

International Nuclear Information System (INIS)

Workman, R. L.; Tiator, L.; Wunderlich, Y.; Doring, M.; Haberzettl, H.

2017-01-01

Here, we compare the methods of amplitude reconstruction, for a complete experiment and a truncated partial-wave analysis, applied to the photoproduction of pseudoscalar mesons. The approach is pedagogical, showing in detail how the amplitude reconstruction (observables measured at a single energy and angle) is related to a truncated partial-wave analysis (observables measured at a single energy and a number of angles).

General method for eliminating wave reflection in 2D photonic crystal waveguides by introducing extra scatterers based on interference cancellation of waves

Science.gov (United States)

Huang, Hao; Ouyang, Zhengbiao

2018-01-01

We propose a general method for eliminating the reflection of waves in 2 dimensional photonic crystal waveguides (2D-PCWs), a kind of 2D material, by introducing extra scatterers inside the 2D-PCWs. The intrinsic reflection in 2D-PCWs is compensated by the backward-scattered waves from these scatterers, so that the overall reflection is greatly reduced and the insertion loss is improved accordingly. We first present the basic theory for the compensation method. Then, as a demonstration, we give four examples of extremely-low-reflection and high-transmission 90°bent 2D-PCWs created according to the method proposed. In the four examples, it is demonstrated by plane-wave expansion method and finite-difference time-domain method that the 90°bent 2D-PCWs can have high transmission ratio greater than 90% in a wide range of operating frequency, and the highest transmission ratio can be greater than 99.95% with a return loss higher than 43 dB, better than that in other typical 90°bent 2D-PCWs. With our method, the bent 2D-PCWs can be optimized to obtain high transmission ratio at different operating wavelengths. As a further application of this method, a waveguide-based optical bridge for light crossing is presented, showing an optimum return loss of 46.85 dB, transmission ratio of 99.95%, and isolation rates greater than 41.77 dB. The method proposed provides also a useful way for improving conventional waveguides made of cables, fibers, or metal walls in the optical, infrared, terahertz, and microwave bands.

The asymptotic expansion method via symbolic computation

OpenAIRE

Navarro, Juan F.

2012-01-01

This paper describes an algorithm for implementing a perturbation method based on an asymptotic expansion of the solution to a second-order differential equation. We also introduce a new symbolic computation system which works with the so-called modified quasipolynomials, as well as an implementation of the algorithm on it.

Higher order polynomial expansion nodal method for hexagonal core neutronics analysis

International Nuclear Information System (INIS)

Jin, Young Cho; Chang, Hyo Kim

1998-01-01

A higher-order polynomial expansion nodal(PEN) method is newly formulated as a means to improve the accuracy of the conventional PEN method solutions to multi-group diffusion equations in hexagonal core geometry. The new method is applied to solving various hexagonal core neutronics benchmark problems. The computational accuracy of the higher order PEN method is then compared with that of the conventional PEN method, the analytic function expansion nodal (AFEN) method, and the ANC-H method. It is demonstrated that the higher order PEN method improves the accuracy of the conventional PEN method and that it compares very well with the other nodal methods like the AFEN and ANC-H methods in accuracy

SCATTERING OF SPIN WAVES BY MAGNETIC DEFECTS

Energy Technology Data Exchange (ETDEWEB)

Callaway, Joseph

1962-12-15

The scattering of spin waves by magnetic point defects is considered using a Green's function method. A partial wave expansion for the scattering amplitude is derived. An expression for the cross section is determined that includes the effect of resonant states. Application is made to the calculation of the thermal conductivity of an insulating ferromagnet. (auth)

Modeling laser beam diffraction and propagation by the mode-expansion method.

Science.gov (United States)

Snyder, James J

2007-08-01

In the mode-expansion method for modeling propagation of a diffracted beam, the beam at the aperture can be expanded as a weighted set of orthogonal modes. The parameters of the expansion modes are chosen to maximize the weighting coefficient of the lowest-order mode. As the beam propagates, its field distribution can be reconstructed from the set of weighting coefficients and the Gouy phase of the lowest-order mode. We have developed a simple procedure to implement the mode-expansion method for propagation through an arbitrary ABCD matrix, and we have demonstrated that it is accurate in comparison with direct calculations of diffraction integrals and much faster.

Three-dimensional stability of solitary kinetic Alfven waves and ion-acoustic waves

International Nuclear Information System (INIS)

Ghosh, G.; Das, K.P.

1994-01-01

Starting from a set of equations that lead to a linear dispersion relation coupling kinetic Alfven waves and ion-acoustic waves, three-dimensional KdV equations are derived for these waves. These equations are then used to investigate the three-dimensional stability of solitary kinetic Alfven waves and ion-acoustic waves by the small-k perturbation expansion method of Rowlands and Infeld. For kinetic Alfven waves it is found that there is instability if the direction of the plane-wave perturbation lies inside a cone, and the growth rate of the instability attains a maximum when the direction of the perturbation lies in the plane containing the external magnetic field and the direction of propagation of the solitary wave. For ion-acoustic waves the growth rate of instability attains a maximum when the direction of the perturbation lies in a plane perpendicular to the direction of propagation of the solitary wave. (Author)

Analysis of Different Methods for Wave Generation and Absorption in a CFD-Based Numerical Wave Tank

Directory of Open Access Journals (Sweden)

Adria Moreno Miquel

2018-06-01

Full Text Available In this paper, the performance of different wave generation and absorption methods in computational fluid dynamics (CFD-based numerical wave tanks (NWTs is analyzed. The open-source CFD code REEF3D is used, which solves the Reynolds-averaged Navier–Stokes (RANS equations to simulate two-phase flow problems. The water surface is computed with the level set method (LSM, and turbulence is modeled with the k-ω model. The NWT includes different methods to generate and absorb waves: the relaxation method, the Dirichlet-type method and active wave absorption. A sensitivity analysis has been conducted in order to quantify and compare the differences in terms of absorption quality between these methods. A reflection analysis based on an arbitrary number of wave gauges has been adopted to conduct the study. Tests include reflection analysis of linear, second- and fifth-order Stokes waves, solitary waves, cnoidal waves and irregular waves generated in an NWT. Wave breaking over a sloping bed and wave forces on a vertical cylinder are calculated, and the influence of the reflections on the wave breaking location and the wave forces on the cylinder is investigated. In addition, a comparison with another open-source CFD code, OpenFOAM, has been carried out based on published results. Some differences in the calculated quantities depending on the wave generation and absorption method have been observed. The active wave absorption method is seen to be more efficient for long waves, whereas the relaxation method performs better for shorter waves. The relaxation method-based numerical beach generally results in lower reflected waves in the wave tank for most of the cases simulated in this study. The comparably better performance of the relaxation method comes at the cost of larger computational requirements due to the relaxation zones that have to be included in the domain. The reflections in the NWT in REEF3D are generally lower than the published results for

The Asymptotic Expansion Method via Symbolic Computation

Directory of Open Access Journals (Sweden)

Juan F. Navarro

2012-01-01

Full Text Available This paper describes an algorithm for implementing a perturbation method based on an asymptotic expansion of the solution to a second-order differential equation. We also introduce a new symbolic computation system which works with the so-called modified quasipolynomials, as well as an implementation of the algorithm on it.

The method of boson expansions in quantum theory

International Nuclear Information System (INIS)

Garbaczewski, P.

1977-06-01

A review is presented of boson expansion methods applied in quantum theory, e.g. expansions of spin, bifermion and fermion operators in cases of finite and infinite number of degrees of freedom. The basic purpose of the paper is to formulate the most general criterion allowing one to obtain the so-called finite spin approximation of any given Bose field theory and the class of fermion theories associated with it. On the other hand, we also need to be able to reconstruct the primary Bose field theory, when any finite spin or Fermi systems are given

Photoelectron wave function in photoionization: plane wave or Coulomb wave?

Science.gov (United States)

Gozem, Samer; Gunina, Anastasia O; Ichino, Takatoshi; Osborn, David L; Stanton, John F; Krylov, Anna I

2015-11-19

The calculation of absolute total cross sections requires accurate wave functions of the photoelectron and of the initial and final states of the system. The essential information contained in the latter two can be condensed into a Dyson orbital. We employ correlated Dyson orbitals and test approximate treatments of the photoelectron wave function, that is, plane and Coulomb waves, by comparing computed and experimental photoionization and photodetachment spectra. We find that in anions, a plane wave treatment of the photoelectron provides a good description of photodetachment spectra. For photoionization of neutral atoms or molecules with one heavy atom, the photoelectron wave function must be treated as a Coulomb wave to account for the interaction of the photoelectron with the +1 charge of the ionized core. For larger molecules, the best agreement with experiment is often achieved by using a Coulomb wave with a partial (effective) charge smaller than unity. This likely derives from the fact that the effective charge at the centroid of the Dyson orbital, which serves as the origin of the spherical wave expansion, is smaller than the total charge of a polyatomic cation. The results suggest that accurate molecular photoionization cross sections can be computed with a modified central potential model that accounts for the nonspherical charge distribution of the core by adjusting the charge in the center of the expansion.

Multimodal method for scattering of sound at a sudden area expansion in a duct with subsonic flow

Science.gov (United States)

Kooijman, G.; Testud, P.; Aurégan, Y.; Hirschberg, A.

2008-03-01

The scattering of sound at a sudden area expansion in a duct with subsonic mean flow has been modelled with a multimodal method. Technological applications are for instance internal combustion engine exhaust silencers and silencers in industrial duct systems. Both two-dimensional (2D) rectangular and 2D cylindrical geometry and uniform mean flow as well as non-uniform mean flow profiles are considered. Model results for the scattering of plane waves in case of uniform flow, in which case an infinitely thin shear layer is formed downstream of the area expansion, are compared to results obtained by other models in literature. Generally good agreement is found. Furthermore, model results for the scattering are compared to experimental data found in literature. Also here fairly good correspondence is observed. When employing a turbulent pipe flow profile in the model, instead of a uniform flow profile, the prediction for the downstream transmission- and upstream reflection coefficient is improved. However, worse agreement is observed for the upstream transmission and downstream reflection coefficient. On the contrary, employing a non-uniform jet flow profile, which represents a typical shear layer flow downstream of the expansion, gives worse agreement for the downstream transmission- and the upstream reflection coefficient, whereas prediction for the upstream transmission and downstream reflection coefficient improves.

Second-Order Perturbation Theory for Generalized Active Space Self-Consistent-Field Wave Functions.

Science.gov (United States)

Ma, Dongxia; Li Manni, Giovanni; Olsen, Jeppe; Gagliardi, Laura

2016-07-12

A multireference second-order perturbation theory approach based on the generalized active space self-consistent-field (GASSCF) wave function is presented. Compared with the complete active space (CAS) and restricted active space (RAS) wave functions, GAS wave functions are more flexible and can employ larger active spaces and/or different truncations of the configuration interaction expansion. With GASSCF, one can explore chemical systems that are not affordable with either CASSCF or RASSCF. Perturbation theory to second order on top of GAS wave functions (GASPT2) has been implemented to recover the remaining electron correlation. The method has been benchmarked by computing the chromium dimer ground-state potential energy curve. These calculations show that GASPT2 gives results similar to CASPT2 even with a configuration interaction expansion much smaller than the corresponding CAS expansion.

A new time–space domain high-order finite-difference method for the acoustic wave equation

KAUST Repository

Liu, Yang; Sen, Mrinal K.

2009-01-01

A new unified methodology was proposed in Finkelstein and Kastner (2007) [39] to derive spatial finite-difference (FD) coefficients in the joint time-space domain to reduce numerical dispersion. The key idea of this method is that the dispersion relation is completely satisfied at several designated frequencies. We develop this new time-space domain FD method further for 1D, 2D and 3D acoustic wave modeling using a plane wave theory and the Taylor series expansion. New spatial FD coefficients are frequency independent though they lead to a frequency dependent numerical solution. We prove that the modeling accuracy is 2nd-order when the conventional (2 M)th-order space domain FD and the 2nd-order time domain FD stencils are directly used to solve the acoustic wave equation. However, under the same discretization, the new 1D method can reach (2 M)th-order accuracy and is always stable. The 2D method can reach (2 M)th-order accuracy along eight directions and has better stability. Similarly, the 3D method can reach (2 M)th-order accuracy along 48 directions and also has better stability than the conventional FD method. The advantages of the new method are also demonstrated by the results of dispersion analysis and numerical modeling of acoustic wave equation for homogeneous and inhomogeneous acoustic models. In addition, we study the influence of the FD stencil length on numerical modeling for 1D inhomogeneous media, and derive an optimal FD stencil length required to balance the accuracy and efficiency of modeling. A new time-space domain high-order staggered-grid FD method for the 1D acoustic wave equation with variable densities is also developed, which has similar advantages demonstrated by dispersion analysis, stability analysis and modeling experiments. The methodology presented in this paper can be easily extended to solve similar partial difference equations arising in other fields of science and engineering. © 2009 Elsevier Inc.

A new time–space domain high-order finite-difference method for the acoustic wave equation

KAUST Repository

Liu, Yang

2009-12-01

A new unified methodology was proposed in Finkelstein and Kastner (2007) [39] to derive spatial finite-difference (FD) coefficients in the joint time-space domain to reduce numerical dispersion. The key idea of this method is that the dispersion relation is completely satisfied at several designated frequencies. We develop this new time-space domain FD method further for 1D, 2D and 3D acoustic wave modeling using a plane wave theory and the Taylor series expansion. New spatial FD coefficients are frequency independent though they lead to a frequency dependent numerical solution. We prove that the modeling accuracy is 2nd-order when the conventional (2 M)th-order space domain FD and the 2nd-order time domain FD stencils are directly used to solve the acoustic wave equation. However, under the same discretization, the new 1D method can reach (2 M)th-order accuracy and is always stable. The 2D method can reach (2 M)th-order accuracy along eight directions and has better stability. Similarly, the 3D method can reach (2 M)th-order accuracy along 48 directions and also has better stability than the conventional FD method. The advantages of the new method are also demonstrated by the results of dispersion analysis and numerical modeling of acoustic wave equation for homogeneous and inhomogeneous acoustic models. In addition, we study the influence of the FD stencil length on numerical modeling for 1D inhomogeneous media, and derive an optimal FD stencil length required to balance the accuracy and efficiency of modeling. A new time-space domain high-order staggered-grid FD method for the 1D acoustic wave equation with variable densities is also developed, which has similar advantages demonstrated by dispersion analysis, stability analysis and modeling experiments. The methodology presented in this paper can be easily extended to solve similar partial difference equations arising in other fields of science and engineering. © 2009 Elsevier Inc.

Tanh-travelling wave solutions, truncated Painleve expansion and reduction of Bullough-Dodd equation to a quadrature in magnetohydrodynamic equilibrium

International Nuclear Information System (INIS)

Ibrahim, R.S.

2003-01-01

The equations of magnetohydrodynamic (MHD) equilibria for a plasma in gravitational field are investigated. For equilibria with one ignorable spatial coordinate, the MHD equations are reduced to a single nonlinear elliptic equation for the magnetic potential u-tilde, known as the Grad-Shafranov equation. Specifying the arbitrary functions in this equation, the Bullough-Dodd equation can be obtained. The truncated Painleve expansion and reduction of the partial differential equation to a quadrature problem (RQ method) are described and applied to obtain the travelling wave solutions of the Bullough-Dodd equation for the case of isothermal magnetostatic atmosphere, in which the current density J is proportional to the exponentially of the magnetic flux and moreover falls off exponentially with distance vertical to the base, with an 'e-folding' distance equal to the gravitational scale height

Paraxial WKB solution of a scalar wave equation

International Nuclear Information System (INIS)

Pereverzev, G.V.

1993-04-01

An asymptotic method of solving a scalar wave equation in inhomogeneous media is developed. This method is an extension of the WKB method to the multidimensional case. It reduces a general wave equation to a set of ordinary differential equations similar to that of the eikonal approach and includes the latter as a particular case. However, the WKB method makes use of another kind of asymptotic expansion and, unlike the eikonal approach, describes the wave properties, i.e. diffraction and interference. At the same time, the three-dimensional WKB method is more simple for numerical treatment because the number of equations is less than in the eikonal approach. The method developed may be used for a calculation of wave fields in problems of RF heating, current drive and plasma diagnostics with microwave beams. (orig.)

Truncated Painleve expansion: Tanh-traveling wave solutions and reduction of sine-Poisson equation to a quadrature for stationary and nonstationary three-dimensional collisionless cold plasma

International Nuclear Information System (INIS)

Ibrahim, R. S.; El-Kalaawy, O. H.

2006-01-01

The relativistic nonlinear self-consistent equations for a collisionless cold plasma with stationary ions [R. S. Ibrahim, IMA J. Appl. Math. 68, 523 (2003)] are extended to 3 and 3+1 dimensions. The resulting system of equations is reduced to the sine-Poisson equation. The truncated Painleve expansion and reduction of the partial differential equation to a quadrature problem (RQ method) are described and applied to obtain the traveling wave solutions of the sine-Poisson equation for stationary and nonstationary equations in 3 and 3+1 dimensions describing the charge-density equilibrium configuration model

Further improved F-expansion method and new exact solutions of Konopelchenko-Dubrovsky equation

International Nuclear Information System (INIS)

Wang Dengshan; Zhang Hongqing

2005-01-01

In this paper, with the aid of the symbolic computation we improve the extended F-expansion method in [Chaos, Solitons and Fractals 2004; 22:111] and propose the further improved F-expansion method. Using this method, we have gotten many new exact solutions which we have never seen before within our knowledge of the (2 + 1)-dimensional Konopelchenko-Dubrovsky equation. In addition,the solutions we get are more general than the solutions that the extended F-expansion method gets.The solutions we get include Jacobi elliptic function solutions, soliton-like solutions, trigonometric function solutions and so on. Our method can also apply to other partial differential equations and can also get many new exact solutions

The focusing effect of electromagnetic waves in two-dimensional photonic crystals with gradually varying lattice constant

Directory of Open Access Journals (Sweden)

F Bakhshi Garmi

2016-02-01

Full Text Available In this paper we studied the focusing effect of electromagnetic wave in the two-dimensional graded photonic crystal consisting of Silicon rods in the air background with gradually varying lattice constant. The results showed that graded photonic crystal can focus wide beams on a narrow area at frequencies near the lower edge of the band gap, where equal frequency contours are not concave. For calculation of photonic band structure and equal frequency contours, we have used plane wave expansion method and revised plane wave expansion method, respectively. The calculation of the electric and magnetic fields was performed by finite difference time domain method.

Analysis shear wave velocity structure obtained from surface wave methods in Bornova, Izmir

Energy Technology Data Exchange (ETDEWEB)

Pamuk, Eren, E-mail: eren.pamuk@deu.edu.tr; Akgün, Mustafa, E-mail: mustafa.akgun@deu.edu.tr [Department of Geophysical Engineering, Dokuz Eylul University, Izmir (Turkey); Özdağ, Özkan Cevdet, E-mail: cevdet.ozdag@deu.edu.tr [Dokuz Eylul University Rectorate, Izmir (Turkey)

2016-04-18

Properties of the soil from the bedrock is necessary to describe accurately and reliably for the reduction of earthquake damage. Because seismic waves change their amplitude and frequency content owing to acoustic impedance difference between soil and bedrock. Firstly, shear wave velocity and depth information of layers on bedrock is needed to detect this changing. Shear wave velocity can be obtained using inversion of Rayleigh wave dispersion curves obtained from surface wave methods (MASW- the Multichannel Analysis of Surface Waves, ReMi-Refraction Microtremor, SPAC-Spatial Autocorrelation). While research depth is limeted in active source study, a passive source methods are utilized for deep depth which is not reached using active source methods. ReMi method is used to determine layer thickness and velocity up to 100 m using seismic refraction measurement systems.The research carried out up to desired depth depending on radius using SPAC which is utilized easily in conditions that district using of seismic studies in the city. Vs profiles which are required to calculate deformations in under static and dynamic loads can be obtained with high resolution using combining rayleigh wave dispersion curve obtained from active and passive source methods. In the this study, Surface waves data were collected using the measurements of MASW, ReMi and SPAC at the İzmir Bornova region. Dispersion curves obtained from surface wave methods were combined in wide frequency band and Vs-depth profiles were obtained using inversion. Reliability of the resulting soil profiles were provided by comparison with theoretical transfer function obtained from soil paremeters and observed soil transfer function from Nakamura technique and by examination of fitting between these functions. Vs values are changed between 200-830 m/s and engineering bedrock (Vs>760 m/s) depth is approximately 150 m.

Application of Rational Expansion Method for Differential-Difference Equation

International Nuclear Information System (INIS)

Wang Qi

2011-01-01

In this paper, we applied the rational formal expansion method to construct a series of soliton-like and period-form solutions for nonlinear differential-difference equations. Compared with most existing methods, the proposed method not only recovers some known solutions, but also finds some new and more general solutions. The efficiency of the method can be demonstrated on Toda Lattice and Ablowitz-Ladik Lattice. (general)

Adaptive Laguerre-Gaussian variant of the Gaussian beam expansion method.

Science.gov (United States)

Cagniot, Emmanuel; Fromager, Michael; Ait-Ameur, Kamel

2009-11-01

A variant of the Gaussian beam expansion method consists in expanding the Bessel function J0 appearing in the Fresnel-Kirchhoff integral into a finite sum of complex Gaussian functions to derive an analytical expression for a Laguerre-Gaussian beam diffracted through a hard-edge aperture. However, the validity range of the approximation depends on the number of expansion coefficients that are obtained by optimization-computation directly. We propose another solution consisting in expanding J0 onto a set of collimated Laguerre-Gaussian functions whose waist depends on their number and then, depending on its argument, predicting the suitable number of expansion functions to calculate the integral recursively.

New exact solutions to MKDV-Burgers equation and (2 + 1)-dimensional dispersive long wave equation via extended Riccati equation method

International Nuclear Information System (INIS)

Kong Cuicui; Wang Dan; Song Lina; Zhang Hongqing

2009-01-01

In this paper, with the aid of symbolic computation and a general ansaetz, we presented a new extended rational expansion method to construct new rational formal exact solutions to nonlinear partial differential equations. In order to illustrate the effectiveness of this method, we apply it to the MKDV-Burgers equation and the (2 + 1)-dimensional dispersive long wave equation, then several new kinds of exact solutions are successfully obtained by using the new ansaetz. The method can also be applied to other nonlinear partial differential equations.

Abundant closed form solutions of the conformable time fractional Sawada-Kotera-Ito equation using (G‧ / G) -expansion method

Science.gov (United States)

Al-Shawba, Altaf Abdulkarem; Gepreel, K. A.; Abdullah, F. A.; Azmi, A.

2018-06-01

In current study, we use the (G‧ / G) -expansion method to construct the closed form solutions of the seventh order time fractional Sawada-Kotera-Ito (TFSKI) equation based on conformable fractional derivative. As a result, trigonometric, hyperbolic and rational functions solutions with arbitrary constants are obtained. When the arbitrary constants are taken some special values, the periodic and soliton solutions are obtained from the travelling wave solutions. The obtained solutions are new and not found elsewhere. The effect of the fractional order on some of these solutions are represented graphically to illustrate the behavior of the exact solutions when the parameter take some special choose.

The essential theory of fast wave current drive with full wave method

International Nuclear Information System (INIS)

Liu Yan; Gong Xueyu; Yang Lei; Yin Chenyan; Yin Lan

2007-01-01

The full wave numerical method is developed for analyzing fast wave current drive in the range of ion cyclotron waves in tokamak plasmas, taking into account finite larmor radius effects and parallel dispersion. the physical model, the dispersion relation on the assumption of Finite Larmor Radius (FLR) effects and the form of full wave be used for computer simulation are developed. All of the work will contribute to further study of fast wave current drive. (authors)

Numerical simulation of stratified shear flow using a higher order Taylor series expansion method

Energy Technology Data Exchange (ETDEWEB)

Iwashige, Kengo; Ikeda, Takashi [Hitachi, Ltd. (Japan)

1995-09-01

A higher order Taylor series expansion method is applied to two-dimensional numerical simulation of stratified shear flow. In the present study, central difference scheme-like method is adopted for an even expansion order, and upwind difference scheme-like method is adopted for an odd order, and the expansion order is variable. To evaluate the effects of expansion order upon the numerical results, a stratified shear flow test in a rectangular channel (Reynolds number = 1.7x10{sup 4}) is carried out, and the numerical velocity and temperature fields are compared with experimental results measured by laser Doppler velocimetry thermocouples. The results confirm that the higher and odd order methods can simulate mean velocity distributions, root-mean-square velocity fluctuations, Reynolds stress, temperature distributions, and root-mean-square temperature fluctuations.

Surface wave velocity tracking by bisection method

International Nuclear Information System (INIS)

Maeda, T.

2005-01-01

Calculation of surface wave velocity is a classic problem dating back to the well-known Haskell's transfer matrix method, which contributes to solutions of elastic wave propagation, global subsurface structure evaluation by simulating observed earthquake group velocities, and on-site evaluation of subsurface structure by simulating phase velocity dispersion curves and/or H/V spectra obtained by micro-tremor observation. Recently inversion analysis on micro-tremor observation requires efficient method of generating many model candidates and also stable, accurate, and fast computation of dispersion curves and Raleigh wave trajectory. The original Haskell's transfer matrix method has been improved in terms of its divergence tendency mainly by the generalized transmission and reflection matrix method with formulation available for surface wave velocity; however, root finding algorithm has not been fully discussed except for the one by setting threshold to the absolute value of complex characteristic functions. Since surface wave number (reciprocal to the surface wave velocity multiplied by frequency) is a root of complex valued characteristic function, it is intractable to use general root finding algorithm. We will examine characteristic function in phase plane to construct two dimensional bisection algorithm with consideration on a layer to be evaluated and algorithm for tracking roots down along frequency axis. (author)

expansion method for the Burgers, Burgers–Huxley and modified

Indian Academy of Sciences (India)

expansion method; Burgers equation; Burgers–Huxley equation; modified. Burgers–KdV equation .... Substituting the solution set (12) and the corresponding solutions of (4) into (8), we have ..... During the past several years, many have done.

Study on evaluation methods for Rayleigh wave dispersion characteristic

Science.gov (United States)

Shi, L.; Tao, X.; Kayen, R.; Shi, H.; Yan, S.

2005-01-01

The evaluation of Rayleigh wave dispersion characteristic is the key step for detecting S-wave velocity structure. By comparing the dispersion curves directly with the spectra analysis of surface waves (SASW) method, rather than comparing the S-wave velocity structure, the validity and precision of microtremor-array method (MAM) can be evaluated more objectively. The results from the China - US joint surface wave investigation in 26 sites in Tangshan, China, show that the MAM has the same precision with SASW method in 83% of the 26 sites. The MAM is valid for Rayleigh wave dispersion characteristic testing and has great application potentiality for site S-wave velocity structure detection.

Conservative numerical methods for solitary wave interactions

Energy Technology Data Exchange (ETDEWEB)

Duran, A; Lopez-Marcos, M A [Departamento de Matematica Aplicada y Computacion, Facultad de Ciencias, Universidad de Valladolid, Paseo del Prado de la Magdalena s/n, 47005 Valladolid (Spain)

2003-07-18

The purpose of this paper is to show the advantages that represent the use of numerical methods that preserve invariant quantities in the study of solitary wave interactions for the regularized long wave equation. It is shown that the so-called conservative methods are more appropriate to study the phenomenon and provide a dynamic point of view that allows us to estimate the changes in the parameters of the solitary waves after the collision.

Rapid expansion and pseudo spectral implementation for reverse time migration in VTI media

KAUST Repository

Pestana, Reynam C

2012-04-24

In isotropic media, we use the scalar acoustic wave equation to perform reverse time migration (RTM) of the recorded pressure wavefield data. In anisotropic media, P- and SV-waves are coupled, and the elastic wave equation should be used for RTM. For computational efficiency, a pseudo-acoustic wave equation is often used. This may be solved using a coupled system of second-order partial differential equations. We solve these using a pseudo spectral method and the rapid expansion method (REM) for the explicit time marching. This method generates a degenerate SV-wave in addition to the P-wave arrivals of interest. To avoid this problem, the elastic wave equation for vertical transversely isotropic (VTI) media can be split into separate wave equations for P- and SV-waves. These separate wave equations are stable, and they can be effectively used to model and migrate seismic data in VTI media where |ε- δ| is small. The artifact for the SV-wave has also been removed. The independent pseudo-differential wave equations can be solved one for each mode using the pseudo spectral method for the spatial derivatives and the REM for the explicit time advance of the wavefield. We show numerically stable and high-resolution modeling and RTM results for the pure P-wave mode in VTI media. © 2012 Sinopec Geophysical Research Institute.

Rapid expansion and pseudo spectral implementation for reverse time migration in VTI media

KAUST Repository

Pestana, Reynam C; Ursin, Bjø rn; Stoffa, Paul L

2012-01-01

In isotropic media, we use the scalar acoustic wave equation to perform reverse time migration (RTM) of the recorded pressure wavefield data. In anisotropic media, P- and SV-waves are coupled, and the elastic wave equation should be used for RTM. For computational efficiency, a pseudo-acoustic wave equation is often used. This may be solved using a coupled system of second-order partial differential equations. We solve these using a pseudo spectral method and the rapid expansion method (REM) for the explicit time marching. This method generates a degenerate SV-wave in addition to the P-wave arrivals of interest. To avoid this problem, the elastic wave equation for vertical transversely isotropic (VTI) media can be split into separate wave equations for P- and SV-waves. These separate wave equations are stable, and they can be effectively used to model and migrate seismic data in VTI media where |ε- δ| is small. The artifact for the SV-wave has also been removed. The independent pseudo-differential wave equations can be solved one for each mode using the pseudo spectral method for the spatial derivatives and the REM for the explicit time advance of the wavefield. We show numerically stable and high-resolution modeling and RTM results for the pure P-wave mode in VTI media. © 2012 Sinopec Geophysical Research Institute.

An acoustic wave equation for pure P wave in 2D TTI media

KAUST Repository

Zhan, Ge; Pestana, Reynam C.; Stoffa, Paul L.

2011-01-01

In this paper, a pure P wave equation for an acoustic 2D TTI media is derived. Compared with conventional TTI coupled equations, the resulting equation is unconditionally stable due to the complete isolation of the SV wave mode. To avoid numerical dispersion and produce high quality images, the rapid expansion method REM is employed for numerical implementation. Synthetic results validate the proposed equation and show that it is a stable algorithm for modeling and reverse time migration RTM in a TTI media for any anisotropic parameter values. © 2011 Society of Exploration Geophysicists.

Bose Operator Expansions of Tensor Operators in the Theory of Magnetism

DEFF Research Database (Denmark)

Lindgård, Per-Anker; Kowalska, A.

1976-01-01

For pt.I see ibid., vol.7, p.1523 (1974). The matching of matrix element method is used to find a new self-consistent Bose operator expansion for tensor operators in spin systems with isotropic exchange interaction plus anisotropy. Tables are given for all tensor operators relevant for cubic...... and hexagonal symmetry. A discussion of renormalized spin-wave theory for a system with planar anisotropy shows that the Goldstone theorem is rigorously fulfilled to the considered order of perturbation. It is finally shown that the new expansion introduces wavevector-dependent terms from the single...

Iterative calculation of reflected and transmitted acoustic waves at a rough interface

NARCIS (Netherlands)

Berkhoff, Arthur P.; van den Berg, P.M.; Thijssen, J.M.

A rigorous iterative technique is described for calculating the acoustic wave reflection and transmission at an irregular interface between two different media. The method is based upon a plane-wave expansion technique in which the acoustic field equations and the radiation condition are satisfied

The periodic wave solutions for the (2 + 1)-dimensional Konopelchenko-Dubrovsky equations

International Nuclear Information System (INIS)

Sheng Zhang

2006-01-01

More periodic wave solutions expressed by Jacobi elliptic functions for the (2 + 1)-dimensional Konopelchenko-Dubrovsky equations are obtained by using the extended F-expansion method. In the limit cases, the solitary wave solutions and trigonometric function solutions for the equations are also obtained

Monte Carlo methods for flux expansion solutions of transport problems

International Nuclear Information System (INIS)

Spanier, J.

1999-01-01

Adaptive Monte Carlo methods, based on the use of either correlated sampling or importance sampling, to obtain global solutions to certain transport problems have recently been described. The resulting learning algorithms are capable of achieving geometric convergence when applied to the estimation of a finite number of coefficients in a flux expansion representation of the global solution. However, because of the nonphysical nature of the random walk simulations needed to perform importance sampling, conventional transport estimators and source sampling techniques require modification to be used successfully in conjunction with such flux expansion methods. It is shown how these problems can be overcome. First, the traditional path length estimators in wide use in particle transport simulations are generalized to include rather general detector functions (which, in this application, are the individual basis functions chosen for the flus expansion). Second, it is shown how to sample from the signed probabilities that arise as source density functions in these applications, without destroying the zero variance property needed to ensure geometric convergence to zero error

New Families of Rational Form Solitary Wave Solutions to (2+1)-Dimensional Broer-Kaup-Kupershmidt System

International Nuclear Information System (INIS)

Wang Qi; Li Biao; Zhang Hongqing; Chen Yong

2005-01-01

Taking the (2+1)-dimensional Broer-Kaup-Kupershmidt system as a simple example, some families of rational form solitary wave solutions, triangular periodic wave solutions, and rational wave solutions are constructed by using the Riccati equation rational expansion method presented by us. The method can also be applied to solve more nonlinear partial differential equation or equations.

An accurate optical design method for synchrotron radiation beamlines with wave-front aberration theory

Energy Technology Data Exchange (ETDEWEB)

Yu, Xiaojiang, E-mail: slsyxj@nus.edu.sg; Diao, Caozheng; Breese, Mark B. H. [Singapore Synchrotron Light Source, National University of Singapore, Singapore 117603 (Singapore)

2016-07-27

An aberration calculation method which was developed by Lu [1] can treat individual aberration term precisely. Spectral aberration is the linear sum of these aberration terms, and the aberrations of multi-element systems also can be calculated correctly when the stretching ratio, defined herein, is unity. Evaluation of focusing mirror-grating systems which are optimized according to Lu’s method, along with the Light Path Function (LPF) and the Spot Diagram method (SD) are discussed to confirm the advantage of Lu’s methodology. Lu’s aberration terms are derived from a precise wave-front treatment, whereas the terms of the power series expansion of the light path function do not yield an accurate sum of the aberrations. Moreover, Lu’s aberration terms can be individually optimized. This is not possible with the analytical spot diagram formulae.

Exact travelling wave solutions of the (3+1)-dimensional mKdV-ZK ...

Indian Academy of Sciences (India)

In this paper, the new generalized (′/)-expansion method is executed to find the travelling wave solutions of the (3+1)-dimensional mKdV-ZK equation and the (1+1)-dimensional compound KdVB equation. The efficiency of this method for finding exact and travelling wave solutions has been demonstrated. It is shown ...

Genetic traces of east-to-west human expansion waves in Eurasia.

Science.gov (United States)

Chaix, Raphaëlle; Austerlitz, Frédéric; Hegay, Tatyana; Quintana-Murci, Lluís; Heyer, Evelyne

2008-07-01

In this study, we describe the landscape of human demographic expansions in Eurasia using a large continental Y chromosome and mitochondrial DNA dataset. Variation at these two uniparentally-inherited genetic systems retraces expansions that occurred in the past 60 ky, and shows a clear decrease of expansion ages from east to west Eurasia. To investigate the demographic events at the origin of this westward decrease of expansion ages, the estimated divergence ages between Eurasian populations are compared with the estimated expansion ages within each population. Both markers suggest that the demographic expansion diffused from east to west in Eurasia in a demic way, i.e., through migrations of individuals (and not just through diffusion of new technologies), highlighting the prominent role of eastern regions within Eurasia during Palaeolithic times. (c) 2008 Wiley-Liss, Inc.

Abundant general solitary wave solutions to the family of KdV type equations

Directory of Open Access Journals (Sweden)

Md. Azmol Huda

2017-03-01

Full Text Available This work explores the construction of more general exact traveling wave solutions of some nonlinear evolution equations (NLEEs through the application of the (G′/G, 1/G-expansion method. This method is allied to the widely used (G′/G-method initiated by Wang et al. and can be considered as an extension of the (G′/G-expansion method. For effectiveness, the method is applied to the family of KdV type equations. Abundant general form solitary wave solutions as well as periodic solutions are successfully obtained through this method. Moreover, in the obtained wider set of solutions, if we set special values of the parameters, some previously known solutions are revived. The approach of this method is simple and elegantly standard. Having been computerized it is also powerful, reliable and effective.

A multi-stage stochastic transmission expansion planning method

International Nuclear Information System (INIS)

Akbari, Tohid; Rahimikian, Ashkan; Kazemi, Ahad

2011-01-01

Highlights: → We model a multi-stage stochastic transmission expansion planning problem. → We include available transfer capability (ATC) in our model. → Involving this criterion will increase the ATC between source and sink points. → Power system reliability will be increased and more money can be saved. - Abstract: This paper presents a multi-stage stochastic model for short-term transmission expansion planning considering the available transfer capability (ATC). The ATC can have a huge impact on the power market outcomes and the power system reliability. The transmission expansion planning (TEP) studies deal with many uncertainties, such as system load uncertainties that are considered in this paper. The Monte Carlo simulation method has been applied for generating different scenarios. A scenario reduction technique is used for reducing the number of scenarios. The objective is to minimize the sum of investment costs (IC) and the expected operation costs (OC). The solution technique is based on the benders decomposition algorithm. The N-1 contingency analysis is also done for the TEP problem. The proposed model is applied to the IEEE 24 bus reliability test system and the results are efficient and promising.

A direct method to transform between expansions in the configuration state function and Slater determinant bases

International Nuclear Information System (INIS)

Olsen, Jeppe

2014-01-01

A novel algorithm is introduced for the transformation of wave functions between the bases of Slater determinants (SD) and configuration state functions (CSF) in the genealogical coupling scheme. By modifying the expansion coefficients as each electron is spin-coupled, rather than performing a single many-electron transformation, the large transformation matrix that plagues previous approaches is avoided and the required number of operations is drastically reduced. As an example of the efficiency of the algorithm, the transformation for a configuration with 30 unpaired electrons and singlet spin is discussed. For this case, the 10 × 10 6 coefficients in the CSF basis is obtained from the 150 × 10 6 coefficients in the SD basis in 1 min, which should be compared with the seven years that the previously employed method is estimated to require

Optimization of nonlinear wave function parameters

International Nuclear Information System (INIS)

Shepard, R.; Minkoff, M.; Chemistry

2006-01-01

An energy-based optimization method is presented for our recently developed nonlinear wave function expansion form for electronic wave functions. This expansion form is based on spin eigenfunctions, using the graphical unitary group approach (GUGA). The wave function is expanded in a basis of product functions, allowing application to closed-shell and open-shell systems and to ground and excited electronic states. Each product basis function is itself a multiconfigurational function that depends on a relatively small number of nonlinear parameters called arc factors. The energy-based optimization is formulated in terms of analytic arc factor gradients and orbital-level Hamiltonian matrices that correspond to a specific kind of uncontraction of each of the product basis functions. These orbital-level Hamiltonian matrices give an intuitive representation of the energy in terms of disjoint subsets of the arc factors, they provide for an efficient computation of gradients of the energy with respect to the arc factors, and they allow optimal arc factors to be determined in closed form for subspaces of the full variation problem. Timings for energy and arc factor gradient computations involving expansion spaces of > 10 24 configuration state functions are reported. Preliminary convergence studies and molecular dissociation curves are presented for some small molecules

Problems of the orthogonalized plane wave method. 1

International Nuclear Information System (INIS)

Farberovich, O.V.; Kurganskii, S.I.; Domashevskaya, E.P.

1979-01-01

The main problems of the orthogonalized plane wave method are discussed including (a) consideration of core states; (b) effect of overlap of wave functions of external core states upon the band structure; (c) calculation of d-type states. The modified orthogonal plane wave method (MOPW method) of Deegan and Twose is applied in a general form to solve the problems of the usual OPW method. For the first time the influence on the spectrum of the main parameters of the MOPW method is studied systematically by calculating the electronic energy spectrum in the transition metals Nb and V. (author)

Lamb wave scattering by a surface-breaking crack in a plate

Science.gov (United States)

Datta, S. K.; Al-Nassar, Y.; Shah, A. H.

1991-01-01

An NDE method based on finite-element representation and modal expansion has been developed for solving the scattering of Lamb waves in an elastic plate waveguide. This method is very powerful for handling discontinuities of arbitrary shape, weldments of different orientations, canted cracks, etc. The advantage of the method is that it can be used to study the scattering of Lamb waves in anisotropic elastic plates and in multilayered plates as well.

Interfacial wave theory for dendritic structure of a growing needle crystal. I - Local instability mechanism. II - Wave-emission mechanism at the turning point

Science.gov (United States)

Xu, Jian-Jun

1989-01-01

The complicated dendritic structure of a growing needle crystal is studied on the basis of global interfacial wave theory. The local dispersion relation for normal modes is derived in a paraboloidal coordinate system using the multiple-variable-expansion method. It is shown that the global solution in a dendrite growth process incorporates the morphological instability factor and the traveling wave factor.

Breaking the Link between Environmental Degradation and Oil Palm Expansion: A Method for Enabling Sustainable Oil Palm Expansion

Science.gov (United States)

Smit, Hans Harmen; Meijaard, Erik; van der Laan, Carina; Mantel, Stephan; Budiman, Arif; Verweij, Pita

2013-01-01

Land degradation is a global concern. In tropical areas it primarily concerns the conversion of forest into non-forest lands and the associated losses of environmental services. Defining such degradation is not straightforward hampering effective reduction in degradation and use of already degraded lands for more productive purposes. To facilitate the processes of avoided degradation and land rehabilitation, we have developed a methodology in which we have used international environmental and social sustainability standards to determine the suitability of lands for sustainable agricultural expansion. The method was developed and tested in one of the frontiers of agricultural expansion, West Kalimantan province in Indonesia. The focus was on oil palm expansion, which is considered as a major driver for deforestation in tropical regions globally. The results suggest that substantial changes in current land-use planning are necessary for most new plantations to comply with international sustainability standards. Through visualizing options for sustainable expansion with our methodology, we demonstrate that the link between oil palm expansion and degradation can be broken. Application of the methodology with criteria and thresholds similar to ours could help the Indonesian government and the industry to achieve its pro-growth, pro-job, pro-poor and pro-environment development goals. For sustainable agricultural production, context specific guidance has to be developed in areas suitable for expansion. Our methodology can serve as a template for designing such commodity and country specific tools and deliver such guidance. PMID:24039700

New exact travelling wave solutions for the generalized nonlinear Schroedinger equation with a source

International Nuclear Information System (INIS)

Abdou, M.A.

2008-01-01

The generalized F-expansion method with a computerized symbolic computation is used for constructing a new exact travelling wave solutions for the generalized nonlinear Schrodinger equation with a source. As a result, many exact travelling wave solutions are obtained which include new periodic wave solution, trigonometric function solutions and rational solutions. The method is straightforward and concise, and it can also be applied to other nonlinear evolution equations in physics

Expansion of the cathode spot and generation of shock waves in the plasma of a volume discharge in atmospheric-pressure helium

International Nuclear Information System (INIS)

Omarov, O. A.; Kurbanismailov, V. S.; Arslanbekov, M. A.; Gadzhiev, M. Kh.; Ragimkhanov, G. B.; Al-Shatravi, Ali J. G.

2012-01-01

The expansion of the cathode spot and the generation of shock waves during the formation and development of a pulsed volume discharge in atmospheric-pressure helium were studied by analyzing the emission spectra of the cathode plasma and the spatiotemporal behavior of the plasma glow. The transition of a diffuse volume discharge in a centimeter-long gap into a high-current diffuse mode when the gas pressure increased from 1 to 5 atm and the applied voltage rose from the statistical breakdown voltage to a 100% overvoltage was investigated. Analytical expressions for the radius of the cathode spot and its expansion velocity obtained in the framework of a spherically symmetric model agree satisfactorily with the experimental data.

Density-functional expansion methods: Grand challenges.

Science.gov (United States)

Giese, Timothy J; York, Darrin M

2012-03-01

We discuss the source of errors in semiempirical density functional expansion (VE) methods. In particular, we show that VE methods are capable of well-reproducing their standard Kohn-Sham density functional method counterparts, but suffer from large errors upon using one or more of these approximations: the limited size of the atomic orbital basis, the Slater monopole auxiliary basis description of the response density, and the one- and two-body treatment of the core-Hamiltonian matrix elements. In the process of discussing these approximations and highlighting their symptoms, we introduce a new model that supplements the second-order density-functional tight-binding model with a self-consistent charge-dependent chemical potential equalization correction; we review our recently reported method for generalizing the auxiliary basis description of the atomic orbital response density; and we decompose the first-order potential into a summation of additive atomic components and many-body corrections, and from this examination, we provide new insights and preliminary results that motivate and inspire new approximate treatments of the core-Hamiltonian.

General method for designing wave shape transformers.

Science.gov (United States)

Ma, Hua; Qu, Shaobo; Xu, Zhuo; Wang, Jiafu

2008-12-22

An effective method for designing wave shape transformers (WSTs) is investigated by adopting the coordinate transformation theory. Following this method, the devices employed to transform electromagnetic (EM) wave fronts from one style with arbitrary shape and size to another style, can be designed. To verify this method, three examples in 2D spaces are also presented. Compared with the methods proposed in other literatures, this method offers the general procedure in designing WSTs, and thus is of great importance for the potential and practical applications possessed by such kinds of devices.

Quality of potential harmonics expansion method for dilute Bose ...

Indian Academy of Sciences (India)

Abstract. We present and examine an approximate but ab initio many-body approach, viz., potential harmonics expansion method (PHEM), which includes two-body correla- tions for dilute Bose–Einstein condensates. Comparing the total ground state energy for three trapped interacting bosons calculated in PHEM with the ...

Thin foil expansion into a vacuum

International Nuclear Information System (INIS)

Mora, P.

2005-01-01

Plasma expansion into a vacuum is an old problem which has been renewed recently in various contexts: expansion of ultra-cold plasmas, cluster expansion, of dust grains, expansion of thin foils. In this presentation I will first discuss the physics of the expansion of a thin foil irradiated by an ultra-short ultra-intense laser pulse. The expansion results in the formation of high energy ions. For an infinitely steep plasma-vacuum interface the fastest ions are located in the outer part of the expansion and their velocity is given by ν m ax∼ 2 C s (In ω p it) where c s (Zk B T e /m i )''1/2 is the ion-acoustic velocity ω p i=(n e 0Ze''2/m i e 0 )''1/2 is the ion plasma frequency, n e 0 is the electron density in the unperturbed plasma, Z is the ion charge number. In the above expression, t is either the pulse duration or the effective acceleration time (in particular t∼L/2c s , where L is the width of the foil, when the electron cooling is taken into account). A salient characteristic of the expansion is the occurrence of a double layer structure and a peak of the accelerating electric field at the ion front. I will explain the origin of the peak and predict its temporal behavior. This peak has been diagnosed in recent experiments. I will also discuss the effect of a 2-temperatures electron distribution function on the expansion, showing the dominant role of the hot electron component. Finally I will discuss the occurrence of ion spikes in the expansion when the initial density profile is smooth. The ion spike is due to a wave breaking which cannot be handled in a satisfactory way by a fluid code and requires a kinetic description. A. simple collisionless particle code has been used to treat the evolution of the spike after the wave breaking and the results will be shown. (Author)

Application of the N-quantum approximation method to bound state problems

International Nuclear Information System (INIS)

Raychaudhuri, A.

1977-01-01

The N-quantum approximation (NQA) method is examined in the light of its application to bound state problems. Bound state wave functions are obtained as expansion coefficients in a truncated Haag expansion. From the equations of motion for the Heisenberg field and the NQA expansion, an equation satisfied by the wave function is derived. Two different bound state systems are considered. In one case, the bound state problem of two identical scalars by scalar exchange is analyzed using the NQA. An integral equation satisfied by the wave function is derived. In the nonrelativistic limit, the equation is shown to reduce to the Schroedinger equation. The equation is solved numerically, and the results compared with those obtained for this system by other methods. The NQA method is also applied to the bound state of two spin 1/2 particles with electromagnetic interaction. The integral equation for the wave function is shown to agree with the corresponding Bethe Salpeter equation in the nonrelativistic limit. Using the Dirac (4 x 4) matrices the wave function is expanded in terms of structure functions and the equation for the wave function is reduced to two disjoint sets of coupled equation for the structure functions

An optimal implicit staggered-grid finite-difference scheme based on the modified Taylor-series expansion with minimax approximation method for elastic modeling

Science.gov (United States)

Yang, Lei; Yan, Hongyong; Liu, Hong

2017-03-01

Implicit staggered-grid finite-difference (ISFD) scheme is competitive for its great accuracy and stability, whereas its coefficients are conventionally determined by the Taylor-series expansion (TE) method, leading to a loss in numerical precision. In this paper, we modify the TE method using the minimax approximation (MA), and propose a new optimal ISFD scheme based on the modified TE (MTE) with MA method. The new ISFD scheme takes the advantage of the TE method that guarantees great accuracy at small wavenumbers, and keeps the property of the MA method that keeps the numerical errors within a limited bound at the same time. Thus, it leads to great accuracy for numerical solution of the wave equations. We derive the optimal ISFD coefficients by applying the new method to the construction of the objective function, and using a Remez algorithm to minimize its maximum. Numerical analysis is made in comparison with the conventional TE-based ISFD scheme, indicating that the MTE-based ISFD scheme with appropriate parameters can widen the wavenumber range with high accuracy, and achieve greater precision than the conventional ISFD scheme. The numerical modeling results also demonstrate that the MTE-based ISFD scheme performs well in elastic wave simulation, and is more efficient than the conventional ISFD scheme for elastic modeling.

Using nodal expansion method in calculation of reactor core with square fuel assemblies

International Nuclear Information System (INIS)

Abdollahzadeh, M. Y.; Boroushaki, M.

2009-01-01

A polynomial nodal method is developed to solve few-group neutron diffusion equations in cartesian geometry. In this article, the effective multiplication factor, group flux and power distribution based on the nodal polynomial expansion procedure is presented. In addition, by comparison of the results the superiority of nodal expansion method on finite-difference and finite-element are fully demonstrated. The comparison of the results obtained by these method with those of the well known benchmark problems have shown that they are in very good agreement.

A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

Energy Technology Data Exchange (ETDEWEB)

Sidler, Rolf, E-mail: rsidler@gmail.com [Center for Research of the Terrestrial Environment, University of Lausanne, CH-1015 Lausanne (Switzerland); Carcione, José M. [Istituto Nazionale di Oceanografia e di Geofisica Sperimentale (OGS), Borgo Grotta Gigante 42c, 34010 Sgonico, Trieste (Italy); Holliger, Klaus [Center for Research of the Terrestrial Environment, University of Lausanne, CH-1015 Lausanne (Switzerland)

2013-02-15

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge–Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid–solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently benchmarked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

Large band gaps of water waves through two-dimensional periodic topography

International Nuclear Information System (INIS)

Yang Shaohua; Wu Fugen; Zhong Huilin; Zhong Lanhua

2006-01-01

In this Letter, the band structures and band gaps of liquid surface waves propagating over two-dimensional periodic topography was investigated by plane-waves expansion method. The periodic topography drilled by square hollows with square lattice was considered. And the effects of the filling fraction and the orientation of bottom-hollows on the band gaps are investigated in detail

Solution of the agglomerate Brownian coagulation using Taylor-expansion moment method.

Science.gov (United States)

Yu, Mingzhou; Lin, Jianzhong

2009-08-01

The newly proposed Taylor-expansion moment method (TEMOM) is extended to solve agglomerate coagulation in the free-molecule regime and in the continuum regime, respectively. The moment equations with respect to fractal dimension are derived based on 3rd Taylor-series expansion technique. The validation of this method is done by comparing its result with the published data at each limited size regime. By comparing with analytical method, sectional method (SM) and quadrature method of moments (QMOMs), this new approach is shown to produce the most efficiency without losing much accuracy. At each limited size regime, the effect of fractal dimension on the decay of particle number and particle size growth is mainly investigated, and especially in the continuum regime the relation of mean diameters of size distributions with different fractal dimensions is first proposed. The agglomerate size distribution is found to be sensitive to the fractal dimension and the initial geometric mean deviation before the self-preserving size distribution is achieved in the continuum regime.

Method for calculating anisotropic neutron transport using scattering kernel without polynomial expansion

International Nuclear Information System (INIS)

Takahashi, Akito; Yamamoto, Junji; Ebisuya, Mituo; Sumita, Kenji

1979-01-01

A new method for calculating the anisotropic neutron transport is proposed for the angular spectral analysis of D-T fusion reactor neutronics. The method is based on the transport equation with new type of anisotropic scattering kernels formulated by a single function I sub(i) (μ', μ) instead of polynomial expansion, for instance, Legendre polynomials. In the calculation of angular flux spectra by using scattering kernels with the Legendre polynomial expansion, we often observe the oscillation with negative flux. But in principle this oscillation disappears by this new method. In this work, we discussed anisotropic scattering kernels of the elastic scattering and the inelastic scatterings which excite discrete energy levels. The other scatterings were included in isotropic scattering kernels. An approximation method, with use of the first collision source written by the I sub(i) (μ', μ) function, was introduced to attenuate the ''oscillations'' when we are obliged to use the scattering kernels with the Legendre polynomial expansion. Calculated results with this approximation showed remarkable improvement for the analysis of the angular flux spectra in a slab system of lithium metal with the D-T neutron source. (author)

Solitary Wave Solutions of the Boussinesq Equation and Its Improved Form

Directory of Open Access Journals (Sweden)

Reza Abazari

2013-01-01

Full Text Available This paper presents the general case study of previous works on generalized Boussinesq equations, (Abazari, 2011 and (Kılıcman and Abazari, 2012, that focuses on the application of G′/G-expansion method with the aid of Maple to construct more general exact solutions for the coupled Boussinesq equations. In this work, the mentioned method is applied to construct more general exact solutions of Boussinesq equation and improved Boussinesq equation, which the French scientist Joseph Valentin Boussinesq (1842–1929 described in the 1870s model equations for the propagation of long waves on the surface of water with small amplitude. Our work is motivated by the fact that the G′/G-expansion method provides not only more general forms of solutions but also periodic, solitary waves and rational solutions. The method appears to be easier and faster by means of a symbolic computation.

Applying the expansion method in hierarchical functions to the solution of Navier-Stokes equations for incompressible fluids

International Nuclear Information System (INIS)

Sabundjian, Gaiane

1999-01-01

This work presents a novel numeric method, based on the finite element method, applied for the solution of the Navier-Stokes equations for incompressible fluids in two dimensions in laminar flow. The method is based on the expansion of the variables in almost hierarchical functions. The used expansion functions are based on Legendre polynomials, adjusted in the rectangular elements in a such a way that corner, side and area functions are defined. The order of the expansion functions associated with the sides and with the area of the elements can be adjusted to the necessary or desired degree. This novel numeric method is denominated by Hierarchical Expansion Method. In order to validate the proposed numeric method three well-known problems of the literature in two dimensions are analyzed. The results show the method capacity in supplying precise results. From the results obtained in this thesis it is possible to conclude that the hierarchical expansion method can be applied successfully for the solution of fluid dynamic problems that involve incompressible fluids. (author)

Non-linear triangle-based polynomial expansion nodal method for hexagonal core analysis

International Nuclear Information System (INIS)

Cho, Jin Young; Cho, Byung Oh; Joo, Han Gyu; Zee, Sung Qunn; Park, Sang Yong

2000-09-01

This report is for the implementation of triangle-based polynomial expansion nodal (TPEN) method to MASTER code in conjunction with the coarse mesh finite difference(CMFD) framework for hexagonal core design and analysis. The TPEN method is a variation of the higher order polynomial expansion nodal (HOPEN) method that solves the multi-group neutron diffusion equation in the hexagonal-z geometry. In contrast with the HOPEN method, only two-dimensional intranodal expansion is considered in the TPEN method for a triangular domain. The axial dependence of the intranodal flux is incorporated separately here and it is determined by the nodal expansion method (NEM) for a hexagonal node. For the consistency of node geometry of the MASTER code which is based on hexagon, TPEN solver is coded to solve one hexagonal node which is composed of 6 triangular nodes directly with Gauss elimination scheme. To solve the CMFD linear system efficiently, stabilized bi-conjugate gradient(BiCG) algorithm and Wielandt eigenvalue shift method are adopted. And for the construction of the efficient preconditioner of BiCG algorithm, the incomplete LU(ILU) factorization scheme which has been widely used in two-dimensional problems is used. To apply the ILU factorization scheme to three-dimensional problem, a symmetric Gauss-Seidel Factorization scheme is used. In order to examine the accuracy of the TPEN solution, several eigenvalue benchmark problems and two transient problems, i.e., a realistic VVER1000 and VVER440 rod ejection benchmark problems, were solved and compared with respective references. The results of eigenvalue benchmark problems indicate that non-linear TPEN method is very accurate showing less than 15 pcm of eigenvalue errors and 1% of maximum power errors, and fast enough to solve the three-dimensional VVER-440 problem within 5 seconds on 733MHz PENTIUM-III. In the case of the transient problems, the non-linear TPEN method also shows good results within a few minute of

Polarization-dependent ponderomotive gradient force in a standing wave

NARCIS (Netherlands)

Smorenburg, P.W.; Kanters, J.H.M.; Lassise, A.; Brussaard, G.J.H.; Kamp, L.P.J.; Luiten, O.J.

2011-01-01

The ponderomotive force is derived for a relativistic charged particle entering an electromagnetic standing wave with a general three-dimensional field distribution and a nonrelativistic intensity, using a perturbation expansion method. It is shown that the well-known ponderomotive gradient force

Interacting wave fronts and rarefaction waves in a second order model of nonlinear thermoviscous fluids : Interacting fronts and rarefaction waves

DEFF Research Database (Denmark)

Rasmussen, Anders Rønne; Sørensen, Mads Peter; Gaididei, Yuri Borisovich

2011-01-01

A wave equation including nonlinear terms up to the second order for a thermoviscous Newtonian fluid is proposed. In the lossless case this equation results from an expansion to third order of the Lagrangian for the fundamental non-dissipative fluid dynamical equations. Thus it preserves the Hami...... is proposed. The dynamics of the rarefaction wave is approximated by a collective coordinate approach in the energy balance equation. © 2010 Springer Science+Business Media B.V.......A wave equation including nonlinear terms up to the second order for a thermoviscous Newtonian fluid is proposed. In the lossless case this equation results from an expansion to third order of the Lagrangian for the fundamental non-dissipative fluid dynamical equations. Thus it preserves...... the Hamiltonian structure, in contrast to the Kuznetsov equation, a model often used in nonlinear acoustics. An exact traveling wave front solution is derived from a generalized traveling wave assumption for the velocity potential. Numerical studies of the evolution of a number of arbitrary initial conditions...

Dispersive solitary wave solutions of Kadomtsev-Petviashvili and modified Kadomtsev-Petviashvili dynamical equations in unmagnetized dust plasma

Science.gov (United States)

Seadawy, A. R.; El-Rashidy, K.

2018-03-01

The Kadomtsev-Petviashvili (KP) and modified KP equations are two of the most universal models in nonlinear wave theory, which arises as a reduction of system with quadratic nonlinearity which admit weakly dispersive waves. The generalized extended tanh method and the F-expansion method are used to derive exact solitary waves solutions of KP and modified KP equations. The region of solutions are displayed graphically.

Numerical investigation of over expanded flow behavior in a single expansion ramp nozzle

Science.gov (United States)

Mousavi, Seyed Mahmood; Pourabidi, Reza; Goshtasbi-Rad, Ebrahim

2018-05-01

The single expansion ramp nozzle is severely over-expanded when the vehicle is at low speed, which hinders its ability to provide optimal configurations for combined cycle engines. The over-expansion leads to flow separation as a result of shock wave/boundary-layer interaction. Flow separation, and the presence of shocks themselves, result in a performance loss in the single expansion ramp nozzle, leading to reduced thrust and increased pressure losses. In the present work, the unsteady two dimensional compressible flow in an over expanded single expansion ramp nozzle has been investigated using finite volume code. To achieve this purpose, the Reynolds stress turbulence model and full multigrid initialization, in addition to the Smirnov's method for examining the errors accumulation, have been employed and the results are compared with available experimental data. The results show that the numerical code is capable of predicting the experimental data with high accuracy. Afterward, the effect of discontinuity jump in wall temperature as well as the length of straight ramp on flow behavior have been studied. It is concluded that variations in wall temperature and length of straight ramp change the shock wave boundary layer interaction, shock structure, shock strength as well as the distance between Lambda shocks.

A high-order perturbation of surfaces method for scattering of linear waves by periodic multiply layered gratings in two and three dimensions

Science.gov (United States)

Hong, Youngjoon; Nicholls, David P.

2017-09-01

The capability to rapidly and robustly simulate the scattering of linear waves by periodic, multiply layered media in two and three dimensions is crucial in many engineering applications. In this regard, we present a High-Order Perturbation of Surfaces method for linear wave scattering in a multiply layered periodic medium to find an accurate numerical solution of the governing Helmholtz equations. For this we truncate the bi-infinite computational domain to a finite one with artificial boundaries, above and below the structure, and enforce transparent boundary conditions there via Dirichlet-Neumann Operators. This is followed by a Transformed Field Expansion resulting in a Fourier collocation, Legendre-Galerkin, Taylor series method for solving the problem in a transformed set of coordinates. Assorted numerical simulations display the spectral convergence of the proposed algorithm.

A wave propagation matrix method in semiclassical theory

International Nuclear Information System (INIS)

Lee, S.Y.; Takigawa, N.

1977-05-01

A wave propagation matrix method is used to derive the semiclassical formulae of the multiturning point problem. A phase shift matrix and a barrier transformation matrix are introduced to describe the processes of a particle travelling through a potential well and crossing a potential barrier respectively. The wave propagation matrix is given by the products of phase shift matrices and barrier transformation matrices. The method to study scattering by surface transparent potentials and the Bloch wave in solids is then applied

Travelling wave solutions to the Kuramoto-Sivashinsky equation

International Nuclear Information System (INIS)

Nickel, J.

2007-01-01

Combining the approaches given by Baldwin [Baldwin D et al. Symbolic computation of exact solutions expressible in hyperbolic and elliptic functions for nonlinear PDEs. J Symbol Comput 2004;37:669-705], Peng [Peng YZ. A polynomial expansion method and new general solitary wave solutions to KS equation. Comm Theor Phys 2003;39:641-2] and by Schuermann [Schuermann HW, Serov VS. Weierstrass' solutions to certain nonlinear wave and evolution equations. Proc progress electromagnetics research symposium, 28-31 March 2004, Pisa. p. 651-4; Schuermann HW. Traveling-wave solutions to the cubic-quintic nonlinear Schroedinger equation. Phys Rev E 1996;54:4312-20] leads to a method for finding exact travelling wave solutions of nonlinear wave and evolution equations (NLWEE). The first idea is to generalize ansaetze given by Baldwin and Peng to find elliptic solutions of NLWEEs. Secondly, conditions used by Schuermann to find physical (real and bounded) solutions and to discriminate between periodic and solitary wave solutions are used. The method is shown in detail by evaluating new solutions of the Kuramoto-Sivashinsky equation

Spectral and partial-wave decomposition of time-dependent wave functions on a grid: Photoelectron spectra of H and H2+ in electromagnetic fields

International Nuclear Information System (INIS)

Nikolopoulos, L. A. A.; Kjeldsen, T. K.; Madsen, L. B.

2007-01-01

We present a method for spectral (bound and continuum) and partial-wave analysis of a three-dimensional time-dependent wave function, defined on a grid, without projecting onto the field-free eigenstates of the system. The method consists of propagating the time-dependent Schroedinger equation to obtain its autocorrelation function C(t)= after the end of the interaction, at time T, of the system with an external time-dependent field. The Fourier spectrum of this correlation function is directly related to the expansion coefficients of the wave function on the field-free bound and continuum energy eigenstates of the system. By expanding on a spherical harmonics basis we show how to calculate the contribution of the various partial waves to the total photoelectron energy spectrum

Precision die design by the die expansion method

CERN Document Server

Ibhadode, A O Akii

2009-01-01

This book presents a new method for the design of the precision dies used in cold-forging, extrusion and drawing processes. The method is based upon die expansion, and attempts to provide a clear-cut theoretical basis for the selection of critical die dimensions for this group of precision dies when the tolerance on product diameter (or thickness) is specified. It also presents a procedure for selecting the minimum-production-cost die from among a set of design alternatives. The mathematical content of the book is relatively simple and will present no difficulty to those who have taken basic c

Distorted wave method in reactions with composite particles

International Nuclear Information System (INIS)

Zelenskaya, N.S.; Teplov, I.B.

1980-01-01

The work deals with the distorbed wave method with a finite radius of interaction (DWBAFR) as applied to quantitative analysis of direct nuclear reactions with composite particles (including heavy ions) considering the reaction mechanisms other than the cluster stripping mechanism, in particular the exchange processes. The accurate equations of the distorbed-wave method in the three-body problem and the general formula dor calculating differential cross-sections of arbitrary binary reactions by DWBAFR are presented. Accurate and approximate methods allowing for finite interaction radius are discussed. Two main versions of exact account of recoil effects: separation of variables in wave functions of relative motion of particles and in interaction potentials and separation of variables in distorted waves are analysed. Given is a characteristic of the known calculated programs approximately and exactly taking account of recoil effects for direct and exchange processes [ru

Shock waves & explosions

CERN Document Server

Sachdev, PL

2004-01-01

Understanding the causes and effects of explosions is important to experts in a broad range of disciplines, including the military, industrial and environmental research, aeronautic engineering, and applied mathematics. Offering an introductory review of historic research, Shock Waves and Explosions brings analytic and computational methods to a wide audience in a clear and thorough way. Beginning with an overview of the research on combustion and gas dynamics in the 1970s and 1980s, the author brings you up to date by covering modeling techniques and asymptotic and perturbative methods and ending with a chapter on computational methods.Most of the book deals with the mathematical analysis of explosions, but computational results are also included wherever they are available. Historical perspectives are provided on the advent of nonlinear science, as well as on the mathematical study of the blast wave phenomenon, both when visualized as a point explosion and when simulated as the expansion of a high-pressure ...

New travelling wave solutions for nonlinear stochastic evolution

Indian Academy of Sciences (India)

The nonlinear stochastic evolution equations have a wide range of applications in physics, chemistry, biology, economics and finance from various points of view. In this paper, the (′/)-expansion method is implemented for obtaining new travelling wave solutions of the nonlinear (2 + 1)-dimensional stochastic ...

Bifurcations of solitary wave solutions for (two and three)-dimensional nonlinear partial differential equation in quantum and magnetized plasma by using two different methods

Science.gov (United States)

Khater, Mostafa M. A.; Seadawy, Aly R.; Lu, Dianchen

2018-06-01

In this research, we study new two techniques that called the extended simple equation method and the novel (G‧/G) -expansion method. The extended simple equation method depend on the auxiliary equation (dϕ/dξ = α + λϕ + μϕ2) which has three ways for solving depends on the specific condition on the parameters as follow: When (λ = 0) this auxiliary equation reduces to Riccati equation, when (α = 0) this auxiliary equation reduces to Bernoulli equation and when (α ≠ 0, λ ≠ 0, μ ≠ 0) we the general solutions of this auxiliary equation while the novel (G‧/G) -expansion method depends also on similar auxiliary equation (G‧/G)‧ = μ + λ(G‧/G) + (v - 1)(G‧/G) 2 which depend also on the value of (λ2 - 4 (v - 1) μ) and the specific condition on the parameters as follow: When (λ = 0) this auxiliary equation reduces to Riccati equation, when (μ = 0) this auxiliary equation reduces to Bernoulli equation and when (λ2 ≠ 4 (v - 1) μ) we the general solutions of this auxiliary equation. This show how both of these auxiliary equation are special cases of Riccati equation. We apply these methods on two dimensional nonlinear Kadomtsev-Petviashvili Burgers equation in quantum plasma and three-dimensional nonlinear modified Zakharov-Kuznetsov equation of ion-acoustic waves in a magnetized plasma. We obtain the exact traveling wave solutions of these important models and under special condition on the parameters, we get solitary traveling wave solutions. All calculations in this study have been established and verified back with the aid of the Maple package program. The executed method is powerful, effective and straightforward for solving nonlinear partial differential equations to obtain more and new solutions.

Determination of optimum "multi-channel surface wave method" field parameters.

Science.gov (United States)

2012-12-01

Multi-channel surface wave methods (especially the multi-channel analyses of surface wave method; MASW) are routinely used to : determine the shear-wave velocity of the subsurface to depths of 100 feet for site classification purposes. Users are awar...

Conformal operator product expansion in the Yukawa model

International Nuclear Information System (INIS)

Prati, M.C.

1983-01-01

Conformal techniques are applied to the Yukawa model, as an example of a theory with spinor fields. It is written the partial-wave analysis of the 4-point function of two scalars and two spinors in the channel phi psi → phi psi in terms of spinor tensor representations of the conformal group. Using this conformal expansion, it is diagonalized the Bethe-Salpeter equation, which is reduced to algebraic relations among the partial waves. It is shown that in the γ 5 -invariant model, but not in the general case, it is possible to derive dynamically from the expansions of the 4-point function the vacuum operator product phi psi>

Conformal invariance and pion wave functions of nonleading twist

International Nuclear Information System (INIS)

Braun, V.M.; Filyanov, I.E.

1989-01-01

The restrictions are studied for the general structure of pion wave functions of twist 3 and twist 4 imposed by the conformal symmetry and the equations of motion. A systematic expansion of wave functions in the conformal spin is built and the first order corrections to asymptotic formulae are calculated by the QCD sum rule method. In particular, we have found a multiplicatively renormalizable contribution into the two-particle wave function of twist 4 which cannot be expanded in a finite set of Gegenbauer polynomials. 19 refs.; 5 figs

A second-order shock-expansion method applicable to bodies of revolution near zero lift

Science.gov (United States)

1957-01-01

A second-order shock-expansion method applicable to bodies of revolution is developed by the use of the predictions of the generalized shock-expansion method in combination with characteristics theory. Equations defining the zero-lift pressure distributions and the normal-force and pitching-moment derivatives are derived. Comparisons with experimental results show that the method is applicable at values of the similarity parameter, the ratio of free-stream Mach number to nose fineness ratio, from about 0.4 to 2.

Design of Extended Depth-of-Focus Laser Beams Using Orthogonal Beam Expansions

Directory of Open Access Journals (Sweden)

Leonard Bergstein

2005-06-01

Full Text Available Laser beams with extended depth of focus have many practical applications, such as scanning printed bar codes. Previous work has concentrated on synthesizing such beams by approximating the nondiffracting Bessel beam solution to the wave equation. In this paper, we introduce an alternate novel synthesis method that is based on maintaining a minimum MTF value (contrast over the largest possible distance. To achieve this, the coefficients of an orthogonal beam expansion are sequentially optimized to this criterion. One of the main advantages of this method is that it can be easily generalized to noncircularly symmetrical beams by the appropriate choice of the beam expansion basis functions. This approach is found to be very useful for applications that involve scanning of the laser beam.

Two dimensional fully nonlinear numerical wave tank based on the BEM

Science.gov (United States)

Sun, Zhe; Pang, Yongjie; Li, Hongwei

2012-12-01

The development of a two dimensional numerical wave tank (NWT) with a rocker or piston type wavemaker based on the high order boundary element method (BEM) and mixed Eulerian-Lagrangian (MEL) is examined. The cauchy principle value (CPV) integral is calculated by a special Gauss type quadrature and a change of variable. In addition the explicit truncated Taylor expansion formula is employed in the time-stepping process. A modified double nodes method is assumed to tackle the corner problem, as well as the damping zone technique is used to absorb the propagation of the free surface wave at the end of the tank. A variety of waves are generated by the NWT, for example; a monochromatic wave, solitary wave and irregular wave. The results confirm the NWT model is efficient and stable.

Asymptotic expansion and statistical description of turbulent systems

International Nuclear Information System (INIS)

Hagan, W.K. III.

1986-01-01

A new approach to studying turbulent systems is presented in which an asymptotic expansion of the general dynamical equations is performed prior to the application of statistical methods for describing the evolution of the system. This approach has been applied to two specific systems: anomalous drift wave turbulence in plasmas and homogeneous, isotropic turbulence in fluids. For the plasma case, the time and length scales of the turbulent state result in the asymptotic expansion of the Vlasov/Poisson equations taking the form of nonlinear gyrokinetic theory. Questions regarding this theory and modern Hamiltonian perturbation methods are discussed and resolved. A new alternative Hamiltonian method is described. The Eulerian Direct Interaction Approximation (EDIA) is slightly reformulated and applied to the equations of nonlinear gyrokinetic theory. Using a similarity transformation technique, expressions for the thermal diffusivity are derived from the EDIA equations for various geometries, including a tokamak. In particular, the unique result for generalized geometry may be of use in evaluating fusion reactor designs and theories of anomalous thermal transport in tokamaks. Finally, a new and useful property of the EDIA is pointed out. For the fluid case, an asymptotic expansion is applied to the Navier-Stokes equation and the results lead to the speculation that such an approach may resolve the problem of predicting the Kolmogorov inertial range energy spectrum for homogeneous, isotropic turbulence. 45 refs., 3 figs

Expansion of a quantum wave packet in a one-dimensional disordered potential in the presence of a uniform bias force

Science.gov (United States)

Crosnier de Bellaistre, C.; Trefzger, C.; Aspect, A.; Georges, A.; Sanchez-Palencia, L.

2018-01-01

We study numerically the expansion dynamics of an initially confined quantum wave packet in the presence of a disordered potential and a uniform bias force. For white-noise disorder, we find that the wave packet develops asymmetric algebraic tails for any ratio of the force to the disorder strength. The exponent of the algebraic tails decays smoothly with that ratio and no evidence of a critical behavior on the wave density profile is found. Algebraic localization features a series of critical values of the force-to-disorder strength where the m th position moment of the wave packet diverges. Below the critical value for the m th moment, we find fair agreement between the asymptotic long-time value of the m th moment and the predictions of diagrammatic calculations. Above it, we find that the m th moment grows algebraically in time. For correlated disorder, we find evidence of systematic delocalization, irrespective to the model of disorder. More precisely, we find a two-step dynamics, where both the center-of-mass position and the width of the wave packet show transient localization, similar to the white-noise case, at short time and delocalization at sufficiently long time. This correlation-induced delocalization is interpreted as due to the decrease of the effective de Broglie wavelength, which lowers the effective strength of the disorder in the presence of finite-range correlations.

On-line reconstruction of in-core power distribution by harmonics expansion method

International Nuclear Information System (INIS)

Wang Changhui; Wu Hongchun; Cao Liangzhi; Yang Ping

2011-01-01

Highlights: → A harmonics expansion method for the on-line in-core power reconstruction is proposed. → A harmonics data library is pre-generated off-line and a code named COMS is developed. → Numerical results show that the maximum relative error of the reconstruction is less than 5.5%. → This method has a high computational speed compared to traditional methods. - Abstract: Fixed in-core detectors are most suitable in real-time response to in-core power distributions in pressurized water reactors (PWRs). In this paper, a harmonics expansion method is used to reconstruct the in-core power distribution of a PWR on-line. In this method, the in-core power distribution is expanded by the harmonics of one reference case. The expansion coefficients are calculated using signals provided by fixed in-core detectors. To conserve computing time and improve reconstruction precision, a harmonics data library containing the harmonics of different reference cases is constructed. Upon reconstruction of the in-core power distribution on-line, the two closest reference cases are searched from the harmonics data library to produce expanded harmonics by interpolation. The Unit 1 reactor of DayaBay Nuclear Power Plant (DayaBay NPP) in China is considered for verification. The maximum relative error between the measurement and reconstruction results is less than 5.5%, and the computing time is about 0.53 s for a single reconstruction, indicating that this method is suitable for the on-line monitoring of PWRs.

Simulation and Experimental Study of Arc Column Expansion After Ignition in Low-Voltage Circuit Breakers

Institute of Scientific and Technical Information of China (English)

MA Qiang; RONG Mingzhe; WU Yi; XU Tiejun; SUN Zhiqiang

2008-01-01

The dynamicprocess of arc pressure and corresponding arc column expansion, which is the main feature after arc ignition and has a significant effect on the breaking behaviour of low -voltage circuit breakers, is studied. By constructing a three dimensional mathematical model of air arc plasma and adopting the Control Volume Method, the parameters of arc plasma including temperature and pressure axe obtained. The variations of pressure field and temperature field with time are simulated. The result indicates that there are six stages for the process of arc column expansion according to the variation of pressure in arc chamber. In the first stage, the maximal pressure locates in the region close to cathode, and in the second stage the maximal pressure shifts to the region close to the anode. In the third stage, the pressure difference between the middle of arc column and the ambient gas is very large, so the arc column begins to expand apparently. In the fourth stage, the pressure wave propagates towards both ends and the maximal pressure appears at the two ends when the pressure wave reaches both sidewalls. In the fifth stage, the pressure wave is reflected and collides in the middle of the arc chamber. In the last stage, the propagation and reflection of pressure wave will repeat several times until a steady burning state is reached. In addition, the experimental results of arc column expansion, corresponding to the arc pressure variation, are presented to verify the simulation results.

Heat kernel expansion in the background field formalism

CERN Document Server

Barvinsky, Andrei

2015-01-01

Heat kernel expansion and background field formalism represent the combination of two calculational methods within the functional approach to quantum field theory. This approach implies construction of generating functionals for matrix elements and expectation values of physical observables. These are functionals of arbitrary external sources or the mean field of a generic configuration -- the background field. Exact calculation of quantum effects on a generic background is impossible. However, a special integral (proper time) representation for the Green's function of the wave operator -- the propagator of the theory -- and its expansion in the ultraviolet and infrared limits of respectively short and late proper time parameter allow one to construct approximations which are valid on generic background fields. Current progress of quantum field theory, its renormalization properties, model building in unification of fundamental physical interactions and QFT applications in high energy physics, gravitation and...

Selective generation of ultrasonic Lamb waves by electromagnetic acoustic transducers

International Nuclear Information System (INIS)

Li Ming-Liang; Deng Ming-Xi; Gao Guang-Jian

2016-01-01

In this paper, we describe a modal expansion approach for the analysis of the selective generation of ultrasonic Lamb waves by electromagnetic acoustic transducers (EMATs). With the modal expansion approach for waveguide excitation, an analytical expression of the Lamb wave’s mode expansion coefficient is deduced, which is related to the driving frequency and the geometrical parameters of the EMAT’s meander coil, and lays a theoretical foundation for exactly analyzing the selective generation of Lamb waves with EMATs. The influences of the driving frequency on the mode expansion coefficient of ultrasonic Lamb waves are analyzed when the EMAT’s geometrical parameters are given. The numerical simulations and experimental examinations show that the ultrasonic Lamb wave modes can be effectively regulated (strengthened or restrained) by choosing an appropriate driving frequency of EMAT, with the geometrical parameters given. This result provides a theoretical and experimental basis for selectively generating a single and pure Lamb wave mode with EMATs. (special topic)

Influence of the nozzle angle on refrigeration performance of a gas wave refrigerator

Science.gov (United States)

Liu, P.; Zhu, Y.; Wang, H.; Zhu, C.; Zou, J.; Wu, J.; Hu, D.

2017-05-01

A gas wave refrigerator (GWR) is a novel refrigerating device that refrigerates a medium by shock waves and expansion waves generated by gas pressure energy. In a typical GWR, the injection energy losses between the nozzle and the expansion tube are essential factors which influence the refrigeration efficiency. In this study, numerical simulations are used to analyze the underlying mechanism of the injection energy losses. The results of simulations show that the vortex loss, mixing energy loss, and oblique shock wave reflection loss are the main factors contributing to the injection energy losses in the expansion tube. Furthermore, the jet angle of the gas is found to dominate the injection energy losses. Therefore, the optimum jet angle is theoretically calculated based on the velocity triangle method. The value of the optimum jet angle is found to be 4^{circ }, 8^{circ }, and 12^{circ } when the refrigeration efficiency is the first-order, second-order, and third-order maximum value over all working ranges of jet frequency, respectively. Finally, a series of experiments are conducted with the jet angle ranging from -4^{circ } to 12^{circ } at a constant expansion ratio. The results indicate the optimal jet angle obtained by the experiments is in good agreement with the calculated value. The isentropic refrigeration efficiency increased by about 4 % after the jet angle was optimized.

Projector augmented wave method: ab initio molecular dynamics ...

Indian Academy of Sciences (India)

Unknown

kinetic energy is small and the wave function is smooth. However, the wave ... and various strategies have been developed. ... methods let us briefly review the history of augmented ..... alleviated by adding an intelligent zero: If an operator B.

About peculiarities of application of the method of fast expansions in the solution of the Navier-Stokes equations

Directory of Open Access Journals (Sweden)

A. D. Chernyshov

2017-01-01

Full Text Available The brief presentation of the method of fast expansions is given to solve nonlinear differential equations. Application  rules of the operator of fast expansions are specified for solving differential equations. According to the method of fast expansions, an unknown function can be represented as the sum of the boundary function and Fourier series sines and cosines for one variable. The special construction of the boundary functions leads to reasonably fast convergence of the Fourier series, so that for engineering calculations, it is sufficient to consider only the first three members. The method is applicable both to linear and nonlinear integro-differential systems. By means of applying the method of fast expansions to nonlinear Navier-Stokes equations the problem is reduced to a closed system of ordinary differential equations, which solution doesn't represent special difficulties. We can reapply the method of fast expansions to the resulting system of differential equations and reduce the original problem to a system of algebraic equations. If the problem is n-dimensional, then after n-fold application of the method of fast expansions the problem will be reduced to a closed algebraic system. Finally, we obtain an analytic-form solution of complicated boundary value problem in partial derivatives. The flow of an incompressible viscous fluid of Navier–Stokes is considered in a curvilinear pipe. The problem is reduced to solving a closed system of ordinary differential equations with boundary conditions by the method of fast expansions. The article considers peculiarities of finding the coefficients of boundary functions and Fourier coefficients for the zero-order and first-order operators of fast expansions. Obtaining the analytic-form solution is of great interest, because it allows to analyze and to investigate the influence of various factors on the properties of the viscous fluid in specific cases.

Nonlinear Wave Propagation and Solitary Wave Formation in Two-Dimensional Heterogeneous Media

KAUST Repository

Luna, Manuel

2011-05-01

Solitary wave formation is a well studied nonlinear phenomenon arising in propagation of dispersive nonlinear waves under suitable conditions. In non-homogeneous materials, dispersion may happen due to effective reflections between the material interfaces. This dispersion has been used along with nonlinearities to find solitary wave formation using the one-dimensional p-system. These solitary waves are called stegotons. The main goal in this work is to find two-dimensional stegoton formation. To do so we consider the nonlinear two-dimensional p-system with variable coefficients and solve it using finite volume methods. The second goal is to obtain effective equations that describe the macroscopic behavior of the variable coefficient system by a constant coefficient one. This is done through a homogenization process based on multiple-scale asymptotic expansions. We compare the solution of the effective equations with the finite volume results and find a good agreement. Finally, we study some stability properties of the homogenized equations and find they and one-dimensional versions of them are unstable in general.

A study on the fusion reactor - A study on wave physics of fast wave heating and the current drive in tokamaks

Energy Technology Data Exchange (ETDEWEB)

Cho, Su Won; Yeom, Hyun Ju [Kyonggi University, Suwon (Korea, Republic of); Hong, Sang Hee; Chung, Mo Se [Seoul National University, Seoul (Korea, Republic of)

1996-09-01

A full 3-dimensional code for fast wave heating and the current drive has been developed ant its results are compared with those of FASTWA for Phaedrus-T tokamak. The finite Larmour radius expansion and the order reduction method have been used to derive the wave equation in the toroidal coordinate from the Maxwell-Vlasov equations. By expanding the fields in poloidal Fourier series, the wave equations are reduced to the system of ordinary differential equations in the radial axis, which are then numerically integrated via the shooting method. In addition, the convergence of the solutions and energy conservation are discussed. Finally, and example calculation of the current drive is presented for the advanced superconducting tokamak which is in its conceptual design phase. 17 refs., 10 tabs., 31 figs. (author)

Inversion of residual stress profiles from ultrasonic Rayleigh wave dispersion data

Science.gov (United States)

Mora, P.; Spies, M.

2018-05-01

We investigate theoretically and with synthetic data the performance of several inversion methods to infer a residual stress state from ultrasonic surface wave dispersion data. We show that this particular problem may reveal in relevant materials undesired behaviors for some methods that could be reliably applied to infer other properties. We focus on two methods, one based on a Taylor-expansion, and another one based on a piecewise linear expansion regularized by a singular value decomposition. We explain the instabilities of the Taylor-based method by highlighting singularities in the series of coefficients. At the same time, we show that the other method can successfully provide performances which only weakly depend on the material.

Bose Operator Expansions of Tensor Operators in the Theory of Magnetism

DEFF Research Database (Denmark)

Kowalska, A.; Lindgård, Per-Anker

1977-01-01

A new Bose operator expansion is discussed for tensor operators in the spin systems with isotropic exchange interaction plus anisotropy. Spin wave theory for a system with planar anisotropy shows that the Goldstone theorem is fulfilled. The new expansion replaces the off diagonal single ion...

The verification of the Taylor-expansion moment method in solving aerosol breakage

Directory of Open Access Journals (Sweden)

Yu Ming-Zhou

2012-01-01

Full Text Available The combination of the method of moment, characterizing the particle population balance, and the computational fluid dynamics has been an emerging research issue in the studies on the aerosol science and on the multiphase flow science. The difficulty of solving the moment equation arises mainly from the closure of some fractal moment variables which appears in the transform from the non-linear integral-differential population balance equation to the moment equations. Within the Taylor-expansion moment method, the breakage-dominated Taylor-expansion moment equation is first derived here when the symmetric fragmentation mechanism is involved. Due to the high efficiency and the high precision, this proposed moment model is expected to become an important tool for solving population balance equations.

Optimal separable bases and series expansions

International Nuclear Information System (INIS)

Poirier, B.

1997-01-01

A method is proposed for the efficient calculation of the Green close-quote s functions and eigenstates for quantum systems of two or more dimensions. For a given Hamiltonian, the best possible separable approximation is obtained from the set of all Hilbert-space operators. It is shown that this determination itself, as well as the solution of the resultant approximation, is a problem of reduced dimensionality. Moreover, the approximate eigenstates constitute the optimal separable basis, in the sense of self-consistent field theory. The full solution is obtained from the approximation via iterative expansion. In the time-independent perturbation expansion for instance, all of the first-order energy corrections are zero. In the Green close-quote s function case, we have a distorted-wave Born series with optimized convergence properties. This series may converge even when the usual Born series diverges. Analytical results are presented for an application of the method to the two-dimensional shifted harmonic-oscillator system, in the course of which the quantum tanh 2 potential problem is solved exactly. The universal presence of bound states in the latter is shown to imply long-lived resonances in the former. In a comparison with other theoretical methods, we find that the reaction path Hamiltonian fails to predict such resonances. copyright 1997 The American Physical Society

Character expansion methods for matrix models of dually weighted graphs

International Nuclear Information System (INIS)

Kazakov, V.A.; Staudacher, M.; Wynter, T.

1996-01-01

We consider generalized one-matrix models in which external fields allow control over the coordination numbers on both the original and dual lattices. We rederive in a simple fashion a character expansion formula for these models originally due to Itzykson and Di Francesco, and then demonstrate how to take the large N limit of this expansion. The relationship to the usual matrix model resolvent is elucidated. Our methods give as a by-product an extremely simple derivation of the Migdal integral equation describing the large N limit of the Itzykson-Zuber formula. We illustrate and check our methods by analysing a number of models solvable by traditional means. We then proceed to solve a new model: a sum over planar graphs possessing even coordination numbers on both the original and the dual lattice. We conclude by formulating equations for the case of arbitrary sets of even, self-dual coupling constants. This opens the way for studying the deep problem of phase transitions from random to flat lattices. (orig.). With 4 figs

Feasibility of wavelet expansion methods to treat the energy variable

International Nuclear Information System (INIS)

Van Rooijen, W. F. G.

2012-01-01

This paper discusses the use of the Discrete Wavelet Transform (DWT) to implement a functional expansion of the energy variable in neutron transport. The motivation of the work is to investigate the possibility of adapting the expansion level of the neutron flux in a material region to the complexity of the cross section in that region. If such an adaptive treatment is possible, 'simple' material regions (e.g., moderator regions) require little effort, while a detailed treatment is used for 'complex' regions (e.g., fuel regions). Our investigations show that in fact adaptivity cannot be achieved. The most fundamental reason is that in a multi-region system, the energy dependence of the cross section in a material region does not imply that the neutron flux in that region has a similar energy dependence. If it is chosen to sacrifice adaptivity, then the DWT method can be very accurate, but the complexity of such a method is higher than that of an equivalent hyper-fine group calculation. The conclusion is thus that, unfortunately, the DWT approach is not very practical. (authors)

Comparison of matrix methods for elastic wave scattering problems

International Nuclear Information System (INIS)

Tsao, S.J.; Varadan, V.K.; Varadan, V.V.

1983-01-01

This article briefly describes the T-matrix method and the MOOT (method of optimal truncation) of elastic wave scattering as they apply to A-D, SH- wave problems as well as 3-D elastic wave problems. Two methods are compared for scattering by elliptical cylinders as well as oblate spheroids of various eccentricity as a function of frequency. Convergence, and symmetry of the scattering cross section are also compared for ellipses and spheroidal cavities of different aspect ratios. Both the T-matrix approach and the MOOT were programmed on an AMDHL 470 computer using double precision arithmetic. Although the T-matrix method and MOOT are not always in agreement, it is in no way implied that any of the published results using MOOT are in error

A robust and efficient stepwise regression method for building sparse polynomial chaos expansions

Energy Technology Data Exchange (ETDEWEB)

Abraham, Simon, E-mail: Simon.Abraham@ulb.ac.be [Vrije Universiteit Brussel (VUB), Department of Mechanical Engineering, Research Group Fluid Mechanics and Thermodynamics, Pleinlaan 2, 1050 Brussels (Belgium); Raisee, Mehrdad [School of Mechanical Engineering, College of Engineering, University of Tehran, P.O. Box: 11155-4563, Tehran (Iran, Islamic Republic of); Ghorbaniasl, Ghader; Contino, Francesco; Lacor, Chris [Vrije Universiteit Brussel (VUB), Department of Mechanical Engineering, Research Group Fluid Mechanics and Thermodynamics, Pleinlaan 2, 1050 Brussels (Belgium)

2017-03-01

Polynomial Chaos (PC) expansions are widely used in various engineering fields for quantifying uncertainties arising from uncertain parameters. The computational cost of classical PC solution schemes is unaffordable as the number of deterministic simulations to be calculated grows dramatically with the number of stochastic dimension. This considerably restricts the practical use of PC at the industrial level. A common approach to address such problems is to make use of sparse PC expansions. This paper presents a non-intrusive regression-based method for building sparse PC expansions. The most important PC contributions are detected sequentially through an automatic search procedure. The variable selection criterion is based on efficient tools relevant to probabilistic method. Two benchmark analytical functions are used to validate the proposed algorithm. The computational efficiency of the method is then illustrated by a more realistic CFD application, consisting of the non-deterministic flow around a transonic airfoil subject to geometrical uncertainties. To assess the performance of the developed methodology, a detailed comparison is made with the well established LAR-based selection technique. The results show that the developed sparse regression technique is able to identify the most significant PC contributions describing the problem. Moreover, the most important stochastic features are captured at a reduced computational cost compared to the LAR method. The results also demonstrate the superior robustness of the method by repeating the analyses using random experimental designs.

Simulation of breaking waves using the high-order spectral method with laboratory experiments: wave-breaking energy dissipation

Science.gov (United States)

Seiffert, Betsy R.; Ducrozet, Guillaume

2018-01-01

We examine the implementation of a wave-breaking mechanism into a nonlinear potential flow solver. The success of the mechanism will be studied by implementing it into the numerical model HOS-NWT, which is a computationally efficient, open source code that solves for the free surface in a numerical wave tank using the high-order spectral (HOS) method. Once the breaking mechanism is validated, it can be implemented into other nonlinear potential flow models. To solve for wave-breaking, first a wave-breaking onset parameter is identified, and then a method for computing wave-breaking associated energy loss is determined. Wave-breaking onset is calculated using a breaking criteria introduced by Barthelemy et al. (J Fluid Mech https://arxiv.org/pdf/1508.06002.pdf, submitted) and validated with the experiments of Saket et al. (J Fluid Mech 811:642-658, 2017). Wave-breaking energy dissipation is calculated by adding a viscous diffusion term computed using an eddy viscosity parameter introduced by Tian et al. (Phys Fluids 20(6): 066,604, 2008, Phys Fluids 24(3), 2012), which is estimated based on the pre-breaking wave geometry. A set of two-dimensional experiments is conducted to validate the implemented wave breaking mechanism at a large scale. Breaking waves are generated by using traditional methods of evolution of focused waves and modulational instability, as well as irregular breaking waves with a range of primary frequencies, providing a wide range of breaking conditions to validate the solver. Furthermore, adjustments are made to the method of application and coefficient of the viscous diffusion term with negligible difference, supporting the robustness of the eddy viscosity parameter. The model is able to accurately predict surface elevation and corresponding frequency/amplitude spectrum, as well as energy dissipation when compared with the experimental measurements. This suggests the model is capable of calculating wave-breaking onset and energy dissipation

Spin-Wave Wave Function for Quantum Spin Models : Condensed Matter and Statistical Physics

OpenAIRE

Franjo, FRANJIC; Sandro, SORELLA; Istituto Nazionale di Fisica della Materia International School for Advance Studies; Istituto Nazionale di Fisica della Materia International School for Advance Studies

1997-01-01

We present a new approach to determine an accurate variational wave function for general quantum spin models, completely defined by a consistency requirement with the simple and well-known linear spin-wave expansion. With this wave function, it is also possible to obtain the correct behavior of the long distance correlation functions for the 1D S=1/2 antiferromagnet. In 2D the proposed spin-wave wave function represents an excellent approximation to the exact ground state of the S=1.2 XY mode...

Multiple scattering and stop band characteristics of flexural waves on a thin plate with circular holes

Science.gov (United States)

Wang, Zuowei; Biwa, Shiro

2018-03-01

A numerical procedure is proposed for the multiple scattering analysis of flexural waves on a thin plate with circular holes based on the Kirchhoff plate theory. The numerical procedure utilizes the wave function expansion of the exciting as well as scattered fields, and the boundary conditions at the periphery of holes are incorporated as the relations between the expansion coefficients of exciting and scattered fields. A set of linear algebraic equations with respect to the wave expansion coefficients of the exciting field alone is established by the numerical collocation method. To demonstrate the applicability of the procedure, the stop band characteristics of flexural waves are analyzed for different arrangements and concentrations of circular holes on a steel plate. The energy transmission spectra of flexural waves are shown to capture the detailed features of the stop band formation of regular and random arrangements of holes. The increase of the concentration of holes is found to shift the dips of the energy transmission spectra toward higher frequencies as well as deepen them. The hexagonal hole arrangement can form a much broader stop band than the square hole arrangement for flexural wave transmission. It is also demonstrated that random arrangements of holes make the transmission spectrum more complicated.

A design of a mode converter for electron cyclotron heating by the method of normal mode expansion

International Nuclear Information System (INIS)

Hoshino, Katsumichi; Kawashima, Hisato; Hata, Kenichiro; Yamamoto, Takumi

1983-09-01

Mode conversion of electromagnetic wave propagating in the over-size circular waveguide is attained by giving a periodical perturbation in the guide wall. Mode coupling equation is expressed by ''generalized telegraphist's equations'' which are derived from the Maxwell's equations using a normal mode expansion. A computer code to solve the coupling equations is developed and mode amplitude, conversion efficiency, etc. of a particular type of mode converter for the 60 GHz electron cyclotron heating are obtained. (author)

Real-space grid implementation of the projector augmented wave method

DEFF Research Database (Denmark)

Mortensen, Jens Jørgen; Hansen, Lars Bruno; Jacobsen, Karsten Wedel

2005-01-01

A grid-based real-space implementation of the projector augmented wave sPAWd method of Blöchl fPhys. Rev. B 50, 17953 s1994dg for density functional theory sDFTd calculations is presented. The use of uniform three-dimensional s3Dd real-space grids for representing wave functions, densities...... valence wave functions that can be represented on relatively coarse grids. We demonstrate the accuracy of the method by calculating the atomization energies of 20 small molecules, and the bulk modulus and lattice constants of bulk aluminum. We show that the approach in terms of computational efficiency...... is comparable to standard plane-wave methods, but the memory requirements are higher....

Singular solitons and other solutions to a couple of nonlinear wave equations

International Nuclear Information System (INIS)

Inc Mustafa; Ulutaş Esma; Biswas Anjan

2013-01-01

This paper addresses the extended (G'/G)-expansion method and applies it to a couple of nonlinear wave equations. These equations are modified the Benjamin—Bona—Mahoney equation and the Boussinesq equation. This extended method reveals several solutions to these equations. Additionally, the singular soliton solutions are revealed, for these two equations, with the aid of the ansatz method

Nonlocal Symmetries, Consistent Riccati Expansion, and Analytical Solutions of the Variant Boussinesq System

Science.gov (United States)

Feng, Lian-Li; Tian, Shou-Fu; Zhang, Tian-Tian; Zhou, Jun

2017-07-01

Under investigation in this paper is the variant Boussinesq system, which describes the propagation of surface long wave towards two directions in a certain deep trough. With the help of the truncated Painlevé expansion, we construct its nonlocal symmetry, Bäcklund transformation, and Schwarzian form, respectively. The nonlocal symmetries can be localised to provide the corresponding nonlocal group, and finite symmetry transformations and similarity reductions are computed. Furthermore, we verify that the variant Boussinesq system is solvable via the consistent Riccati expansion (CRE). By considering the consistent tan-function expansion (CTE), which is a special form of CRE, the interaction solutions between soliton and cnoidal periodic wave are explicitly studied.

Reflection and diffraction of atomic de Broglie waves by evanescent laser waves. Bare-state method

International Nuclear Information System (INIS)

Feng, Xiaoping; Witte, N.S.; Hollenberg, C.L.; Opat, G.

1994-01-01

Two methods are presented for the investigation of the reflection and diffraction of atoms by gratings formed either by standing or travelling evanescent laser waves. Both methods use the bare-state rather than dressed-state picture. One method is based on the Born series, whereas the other is based on the Laplace transformation of the coupled differential equations. The two methods yield the same theoretical expressions for the reflected and diffracted atomic waves in the whole space including the interaction and the asymptotic regions. 1 ref., 1 fig

Symbolic computation and abundant travelling wave solutions to ...

Indian Academy of Sciences (India)

2016-12-09

Dec 9, 2016 ... Abstract. In this article, the novel (G /G)-expansion method is successfully applied to construct the abundant travelling wave solutions to the KdV–mKdV equation with the aid of symbolic computation. This equation is one of the most popular equation in soliton physics and appear in many practical scenarios ...

Analysis of efficient preconditioned defect correction methods for nonlinear water waves

DEFF Research Database (Denmark)

Engsig-Karup, Allan Peter

2014-01-01

Robust computational procedures for the solution of non-hydrostatic, free surface, irrotational and inviscid free-surface water waves in three space dimensions can be based on iterative preconditioned defect correction (PDC) methods. Such methods can be made efficient and scalable to enable...... prediction of free-surface wave transformation and accurate wave kinematics in both deep and shallow waters in large marine areas or for predicting the outcome of experiments in large numerical wave tanks. We revisit the classical governing equations are fully nonlinear and dispersive potential flow...... equations. We present new detailed fundamental analysis using finite-amplitude wave solutions for iterative solvers. We demonstrate that the PDC method in combination with a high-order discretization method enables efficient and scalable solution of the linear system of equations arising in potential flow...

Analysis of the validity of the asymptotic techniques in the lower hybrid wave equation solution for reactor applications

International Nuclear Information System (INIS)

Cardinali, A.; Morini, L.; Castaldo, C.; Cesario, R.; Zonca, F.

2007-01-01

Knowing that the lower hybrid (LH) wave propagation in tokamak plasmas can be correctly described with a full wave approach only, based on fully numerical techniques or on semianalytical approaches, in this paper, the LH wave equation is asymptotically solved via the Wentzel-Kramers-Brillouin (WKB) method for the first two orders of the expansion parameter, obtaining governing equations for the phase at the lowest and for the amplitude at the next order. The nonlinear partial differential equation (PDE) for the phase is solved in a pseudotoroidal geometry (circular and concentric magnetic surfaces) by the method of characteristics. The associated system of ordinary differential equations for the position and the wavenumber is obtained and analytically solved by choosing an appropriate expansion parameter. The quasilinear PDE for the WKB amplitude is also solved analytically, allowing us to reconstruct the wave electric field inside the plasma. The solution is also obtained numerically and compared with the analytical solution. A discussion of the validity limits of the WKB method is also given on the basis of the obtained results

Simulation of breaking waves using the high-order spectral method with laboratory experiments: Wave-breaking onset

Science.gov (United States)

Seiffert, Betsy R.; Ducrozet, Guillaume; Bonnefoy, Félicien

2017-11-01

This study investigates a wave-breaking onset criteria to be implemented in the non-linear potential flow solver HOS-NWT. The model is a computationally efficient, open source code, which solves for the free surface in a numerical wave tank using the High-Order Spectral (HOS) method. The goal of this study is to determine the best method to identify the onset of random single and multiple breaking waves over a large domain at the exact time they occur. To identify breaking waves, a breaking onset criteria based on the ratio of local energy flux velocity to the local crest velocity, introduced by Barthelemy et al. (2017) is selected. The breaking parameter is uniquely applied in the numerical model in that calculations of the breaking onset criteria ratio are not made only at the location of the wave crest, but at every point in the domain and at every time step. This allows the model to calculate the onset of a breaking wave the moment it happens, and without knowing anything about the wave a priori. The application of the breaking criteria at every point in the domain and at every time step requires the phase velocity to be calculated instantaneously everywhere in the domain and at every time step. This is achieved by calculating the instantaneous phase velocity using the Hilbert transform and dispersion relation. A comparison between more traditional crest-tracking techniques shows the calculation of phase velocity using Hilbert transform at the location of the breaking wave crest provides a good approximation of crest velocity. The ability of the selected wave breaking criteria to predict single and multiple breaking events in two dimensions is validated by a series of large-scale experiments. Breaking waves are generated by energy focusing and modulational instability methods, with a wide range of primary frequencies. Steep irregular waves which lead to breaking waves, and irregular waves with an energy focusing wave superimposed are also generated. This set of

Exchange splitting of the interaction energy and the multipole expansion of the wave function

Energy Technology Data Exchange (ETDEWEB)

Gniewek, Piotr, E-mail: pgniewek@tiger.chem.uw.edu.pl; Jeziorski, Bogumił, E-mail: jeziorsk@chem.uw.edu.pl [Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw (Poland)

2015-10-21

The exchange splitting J of the interaction energy of the hydrogen atom with a proton is calculated using the conventional surface-integral formula J{sub surf}[Φ], the volume-integral formula of the symmetry-adapted perturbation theory J{sub SAPT}[Φ], and a variational volume-integral formula J{sub var}[Φ]. The calculations are based on the multipole expansion of the wave function Φ, which is divergent for any internuclear distance R. Nevertheless, the resulting approximations to the leading coefficient j{sub 0} in the large-R asymptotic series J(R) = 2e{sup −R−1}R(j{sub 0} + j{sub 1}R{sup −1} + j{sub 2}R{sup −2} + ⋯) converge with the rate corresponding to the convergence radii equal to 4, 2, and 1 when the J{sub var}[Φ], J{sub surf}[Φ], and J{sub SAPT}[Φ] formulas are used, respectively. Additionally, we observe that also the higher j{sub k} coefficients are predicted correctly when the multipole expansion is used in the J{sub var}[Φ] and J{sub surf}[Φ] formulas. The symmetry adapted perturbation theory formula J{sub SAPT}[Φ] predicts correctly only the first two coefficients, j{sub 0} and j{sub 1}, gives a wrong value of j{sub 2}, and diverges for higher j{sub n}. Since the variational volume-integral formula can be easily generalized to many-electron systems and evaluated with standard basis-set techniques of quantum chemistry, it provides an alternative for the determination of the exchange splitting and the exchange contribution of the interaction potential in general.

Research of the elastic waves generated by a pulse laser. Excitation mechanism of elastic waves and application to nondestructive testing; Pulse laser de reikishita danseiha ni kansuru kenkyu. Danseiha reiki no mechanism to hihakai kensa eno oyo

Energy Technology Data Exchange (ETDEWEB)

Cho, H.; Takemoto, M. [Aoyama Gakuin University, Tokyo (Japan). College of Science and Engineering

1994-07-20

A bulk wave is generated when a pulse laser is irradiated to the material, and the characteristics of a Young`s modulus and Poisson`s ratio can be nondestructively estimated from the bulk wave. The generation mechanism of laser ultrasonic waves must be first clarified for such application. In this paper, fundamental research was conducted to study the generation mechanism of the elastic waves excited by a Q-switched Nd-YAG laser, and the generation method and characteristics of Rayleigh waves. The following result was obtained. A bulk wave is generated by the disk-like adiabatic expansion near the surface if the laser power is small when a spot-shape pulse laser was irradiated. A bulk wave is generated by the thin disk-like adiabatic expansion beneath the surface due to the thermal diffusion in the depth direction of a base material when the laser power becomes large. Moreover, a bulk wave is generated by the impact force due to abrasion and plasma when the power becomes still larger. The information on the bulk wave characteristics and Rayleigh wave was also obtained. 25 refs., 15 figs., 1 tab.

Damage detection in composite materials using Lamb wave methods

Science.gov (United States)

Kessler, Seth S.; Spearing, S. Mark; Soutis, Constantinos

2002-04-01

Cost-effective and reliable damage detection is critical for the utilization of composite materials. This paper presents part of an experimental and analytical survey of candidate methods for in situ damage detection of composite materials. Experimental results are presented for the application of Lamb wave techniques to quasi-isotropic graphite/epoxy test specimens containing representative damage modes, including delamination, transverse ply cracks and through-holes. Linear wave scans were performed on narrow laminated specimens and sandwich beams with various cores by monitoring the transmitted waves with piezoceramic sensors. Optimal actuator and sensor configurations were devised through experimentation, and various types of driving signal were explored. These experiments provided a procedure capable of easily and accurately determining the time of flight of a Lamb wave pulse between an actuator and sensor. Lamb wave techniques provide more information about damage presence and severity than previously tested methods (frequency response techniques), and provide the possibility of determining damage location due to their local response nature. These methods may prove suitable for structural health monitoring applications since they travel long distances and can be applied with conformable piezoelectric actuators and sensors that require little power.

Analytic function expansion nodal method for nuclear reactor core design

International Nuclear Information System (INIS)

Noh, Hae Man

1995-02-01

In most advanced nodal methods the transverse integration is commonly used to reduce the multi-dimensional diffusion equation into equivalent one- dimensional diffusion equations when derving the nodal coupling equations. But the use of the transverse integration results in some limitations. The first limitation is that the transverse leakage term which appears in the transverse integration procedure must be appropriately approximated. The second limitation is that the one-dimensional flux shapes in each spatial direction resulted from the nodal calculation are not accurate enough to be directly used in reconstructing the pinwise flux distributions. Finally the transverse leakage defined for a non-rectangular node such as a hexagonal node or a triangular node is too complicated to be easily handled and may contain non-physical singular terms of step-function and delta-function types. In this thesis, the Analytic Function Expansion Nodal (AFEN) method and its two variations : the Polynomial Expansion Nodal (PEN) method and the hybrid of the AFEN and PEN methods, have been developed to overcome the limitations of the transverse integration procedure. All of the methods solve the multidimensional diffusion equation without the transverse integration. The AFEN method which we believe is the major contribution of this study to the reactor core analysis expands the homogeneous flux distributions within a node in non-separable analytic basis functions satisfying the neutron diffusion equations at any point of the node and expresses the coefficients of the flux expansion in terms of the nodal unknowns which comprise a node-average flux, node-interface fluxes, and corner-point fluxes. Then, the nodal coupling equations composed of the neutron balance equations, the interface current continuity equations, and the corner-point leakage balance equations are solved iteratively to determine all the nodal unknowns. Since the AFEN method does not use the transverse integration in

Explicit and exact nontraveling wave solutions of the (3+1)-dimensional potential Yu-Toda-Sasa-Fukuyama equation

Science.gov (United States)

Yuan, Na

2018-04-01

With the aid of the symbolic computation, we present an improved ( G ‧ / G ) -expansion method, which can be applied to seek more types of exact solutions for certain nonlinear evolution equations. In illustration, we choose the (3 + 1)-dimensional potential Yu-Toda-Sasa-Fukuyama equation to demonstrate the validity and advantages of the method. As a result, abundant explicit and exact nontraveling wave solutions are obtained including two solitary waves solutions, nontraveling wave solutions and dromion soliton solutions. Some particular localized excitations and the interactions between two solitary waves are researched. The method can be also applied to other nonlinear partial differential equations.

An incident flux expansion transport theory method suitable for coupling to diffusion theory methods in hexagonal geometry

International Nuclear Information System (INIS)

Hayward, Robert M.; Rahnema, Farzad; Zhang, Dingkang

2013-01-01

Highlights: ► A new hybrid stochastic–deterministic transport theory method to couple with diffusion theory. ► The method is implemented in 2D hexagonal geometry. ► The new method produces excellent results when compared with Monte Carlo reference solutions. ► The method is fast, solving all test cases in less than 12 s. - Abstract: A new hybrid stochastic–deterministic transport theory method, which is designed to couple with diffusion theory, is presented. The new method is an extension of the incident flux response expansion method, and it combines the speed of diffusion theory with the accuracy of transport theory. With ease of use in mind, the new method is derived in such a way that it can be implemented with only minimal modifications to an existing diffusion theory method. A new angular expansion, which is necessary for the diffusion theory coupling, is developed in 2D and 3D. The method is implemented in 2D hexagonal geometry, and an HTTR benchmark problem is used to test its accuracy in a standalone configuration. It is found that the new method produces excellent results (with average relative error in partial current less than 0.033%) when compared with Monte Carlo reference solutions. Furthermore, the method is fast, solving all test cases in less than 12 s

Parabolic approximation method for fast magnetosonic wave propagation in tokamaks

International Nuclear Information System (INIS)

Phillips, C.K.; Perkins, F.W.; Hwang, D.Q.

1985-07-01

Fast magnetosonic wave propagation in a cylindrical tokamak model is studied using a parabolic approximation method in which poloidal variations of the wave field are considered weak in comparison to the radial variations. Diffraction effects, which are ignored by ray tracing mthods, are included self-consistently using the parabolic method since continuous representations for the wave electromagnetic fields are computed directly. Numerical results are presented which illustrate the cylindrical convergence of the launched waves into a diffraction-limited focal spot on the cyclotron absorption layer near the magnetic axis for a wide range of plasma confinement parameters

Numerical method for wave forces acting on partially perforated caisson

Science.gov (United States)

Jiang, Feng; Tang, Xiao-cheng; Jin, Zhao; Zhang, Li; Chen, Hong-zhou

2015-04-01

The perforated caisson is widely applied to practical engineering because of its great advantages in effectively wave energy consumption and cost reduction. The attentions of many scientists were paid to the fluid-structure interaction between wave and perforated caisson studies, but until now, most concerns have been put on theoretical analysis and experimental model set up. In this paper, interaction between the wave and the partial perforated caisson in a 2D numerical wave flume is investigated by means of the renewed SPH algorithm, and the mathematical equations are in the form of SPH numerical approximation based on Navier-Stokes equations. The validity of the SPH mathematical method is examined and the simulated results are compared with the results of theoretical models, meanwhile the complex hydrodynamic characteristics when the water particles flow in or out of a wave absorbing chamber are analyzed and the wave pressure distribution of the perforated caisson is also addressed here. The relationship between the ratio of total horizontal force acting on caisson under regular waves and its influence factors is examined. The data show that the numerical calculation of the ratio of total horizontal force meets the empirical regression equation very well. The simulations of SPH about the wave nonlinearity and breaking are briefly depicted in the paper, suggesting that the advantages and great potentiality of the SPH method is significant compared with traditional methods.

Expansions of general stationary stochastic optical fields: general formalism

International Nuclear Information System (INIS)

Martinez-Herrero, R.; Mejias, P.M.

1985-01-01

A new expansion of a general stationary stochastic optical field is derived. Each term of the series is seen to represent a recently defined new class of optical fields, the so-called spectrally quasi-factorizable fields. Alternative expansion in terms of nonstationary fields that obey the wave equation is also shown. A relationship between temporal and spatial features of stationary free optical fields is discussed

S-wave velocity measurements along levees in New Orleans using passive surface wave methods

Science.gov (United States)

Hayashi, K.; Lorenzo, J. M.; Craig, M. S.; Gostic, A.

2017-12-01

In order to develop non-invasive methods for levee inspection, geophysical investigations were carried out at four sites along levees in the New Orleans area: 17th Street Canal, London Avenue Canal, Marrero Levee, and Industrial Canal. Three of the four sites sustained damage from Hurricane Katrina in 2005 and have since been rebuilt. The geophysical methods used include active and passive surface wave methods, and capacitively coupled resistivity. This paper summarizes the acquisition and analysis of the 1D and 2D passive surface wave data. Twelve wireless seismic data acquisition units with 2 Hz vertical component geophones were used to record data. Each unit includes a GPS receiver so that all units can be synchronized over any distance without cables. The 1D passive method used L shaped arrays of three different sizes with geophone spacing ranging from 5 to 340 m. Ten minutes to one hour of ambient noise was recorded with each array, and total data acquisition took approximately two hours at each site. The 2D method used a linear array with a geophone spacing of 5m. Four geophones were moved forward every 10 minutes along 400 1000 m length lines. Data acquisition took several hours for each line. Recorded ambient noise was processed using the spatial autocorrelation method and clear dispersion curves were obtained at all sites (Figure 1a). Minimum frequencies ranged from 0.4 to 0.7 Hz and maximum frequencies ranged from 10 to 30 Hz depending on the site. Non-linear inversion was performed and 1D and 2D S-wave velocity models were obtained. The 1D method penetrated to depths ranging from 200 to 500 m depending on the site (Figure 1b). The 2D method penetrated to a depth of 40 60 m and provided 400 1000 m cross sections along the levees (Figure 2). The interpretation focused on identifying zones beneath the levees or canal walls having low S-wave velocities corresponding to saturated, unconsolidated sands, or low-rigidity clays. Resultant S-wave velocity profiles

Producing accurate wave propagation time histories using the global matrix method

International Nuclear Information System (INIS)

Obenchain, Matthew B; Cesnik, Carlos E S

2013-01-01

This paper presents a reliable method for producing accurate displacement time histories for wave propagation in laminated plates using the global matrix method. The existence of inward and outward propagating waves in the general solution is highlighted while examining the axisymmetric case of a circular actuator on an aluminum plate. Problems with previous attempts to isolate the outward wave for anisotropic laminates are shown. The updated method develops a correction signal that can be added to the original time history solution to cancel the inward wave and leave only the outward propagating wave. The paper demonstrates the effectiveness of the new method for circular and square actuators bonded to the surface of isotropic laminates, and these results are compared with exact solutions. Results for circular actuators on cross-ply laminates are also presented and compared with experimental results, showing the ability of the new method to successfully capture the displacement time histories for composite laminates. (paper)

International Franchising as a Method for Business Expansion

OpenAIRE

Karpushina, Darya Evgenjevna

2009-01-01

The present Master Thesis investigates the concept of international franchising from both business and legal standpoints. The actuality of the topic is obvious: Franchising becomes one of the most perspective and fast-developing method for business expansion, and this Diploma was written as a reflection of such tendency. In the meantime, Franchising is an extremely complex and arguable business issue and still causes a kind of confusion in people's mind. For this reason, my effort in this Wor...

Thermal Expansion and Magnetostriction Measurements at Cryogenic Temperature Using the Strain Gauge Method.

Science.gov (United States)

Wang, Wei; Liu, Huiming; Huang, Rongjin; Zhao, Yuqiang; Huang, Chuangjun; Guo, Shibin; Shan, Yi; Li, Laifeng

2018-01-01

Thermal expansion and magnetostriction, the strain responses of a material to temperature and a magnetic field, especially properties at low temperature, are extremely useful to study electronic and phononic properties, phase transitions, quantum criticality, and other interesting phenomena in cryogenic engineering and materials science. However, traditional dilatometers cannot provide magnetic field and ultra-low temperature (thermal expansion and magnetostriction at cryogenic temperature using the strain gauge method based on a Physical Properties Measurements System (PPMS). The interfacing software and automation were developed using LabVIEW. The sample temperature range can be tuned continuously between 1.8 and 400 K. With this PPMS-aided measuring system, we can observe temperature and magnetic field dependence of the linear thermal expansion of different solid materials easily and accurately.

Does query expansion limit our learning? A comparison of social-based expansion to content-based expansion for medical queries on the internet.

Science.gov (United States)

Pentoney, Christopher; Harwell, Jeff; Leroy, Gondy

2014-01-01

Searching for medical information online is a common activity. While it has been shown that forming good queries is difficult, Google's query suggestion tool, a type of query expansion, aims to facilitate query formation. However, it is unknown how this expansion, which is based on what others searched for, affects the information gathering of the online community. To measure the impact of social-based query expansion, this study compared it with content-based expansion, i.e., what is really in the text. We used 138,906 medical queries from the AOL User Session Collection and expanded them using Google's Autocomplete method (social-based) and the content of the Google Web Corpus (content-based). We evaluated the specificity and ambiguity of the expansion terms for trigram queries. We also looked at the impact on the actual results using domain diversity and expansion edit distance. Results showed that the social-based method provided more precise expansion terms as well as terms that were less ambiguous. Expanded queries do not differ significantly in diversity when expanded using the social-based method (6.72 different domains returned in the first ten results, on average) vs. content-based method (6.73 different domains, on average).

The Method of Lines Solution of the Regularized Long-Wave Equation Using Runge-Kutta Time Discretization Method

Directory of Open Access Journals (Sweden)

H. O. Bakodah

2013-01-01

Full Text Available A method of lines approach to the numerical solution of nonlinear wave equations typified by the regularized long wave (RLW is presented. The method developed uses a finite differences discretization to the space. Solution of the resulting system was obtained by applying fourth Runge-Kutta time discretization method. Using Von Neumann stability analysis, it is shown that the proposed method is marginally stable. To test the accuracy of the method some numerical experiments on test problems are presented. Test problems including solitary wave motion, two-solitary wave interaction, and the temporal evaluation of a Maxwellian initial pulse are studied. The accuracy of the present method is tested with and error norms and the conservation properties of mass, energy, and momentum under the RLW equation.

Application of structural symmetries in the plane-wave-based transfer-matrix method for three-dimensional photonic crystal waveguides

International Nuclear Information System (INIS)

Li Zhiyuan; Ho Kaiming

2003-01-01

The plane-wave-based transfer-matrix method (TMM) exhibits a peculiar advantage of being capable of solving eigenmodes involved in an infinite photonic crystal and electromagnetic (EM) wave propagation in finite photonic crystal slabs or even semi-infinite photonic crystal structures within the same theoretical framework. In addition, this theoretical approach can achieve much improved numerical convergency in solution of photonic band structures than the conventional plane-wave expansion method. In this paper we employ this TMM in combination with a supercell technique to handle two important kinds of three-dimensional (3D) photonic crystal waveguide structures. The first one is waveguides created in a 3D layer-by-layer photonic crystal that possesses a complete band gap, the other more popular one is waveguides built in a two-dimensional photonic crystal slab. These waveguides usually have mirror-reflection symmetries in one or two directions perpendicular to their axis. We have taken advantage of these structural symmetries to reduce the numerical burden of the TMM solution of the guided modes. The solution to the EM problems under these mirror-reflection symmetries in both the real space and the plane-wave space is discussed in a systematic way and in great detail. Both the periodic boundary condition and the absorbing boundary condition are employed to investigate structures with or without complete 3D optical confinement. The fact that the EM field components investigated in the TMM are collinear with the symmetric axes of the waveguide brings great convenience and clarity in exploring the eigenmode symmetry in both the real space and the plane-wave space. The classification of symmetry involved in the guided modes can help people to better understand the coupling of the photonic crystal waveguides with external channels such as dielectric slab or wire waveguides

Transient Analysis of Dispersive Power-Ground Plate Pairs With Arbitrarily Shaped Antipads by the DGTD Method With Wave Port Excitation

KAUST Repository

Li, Ping

2016-09-09

A discontinuous Galerkin time-domain (DGTD) method analyzing signal/power integrity on multilayered power-ground parallel plate pairs is proposed. The excitation is realized by introducing wave ports on the antipads where electric/magnetic current sources are represented in terms of the eigenmodes of the antipads. Since closed-forms solutions do not exist for the eigenmodes of the arbitrarily shaped antipads, they have to be calculated using numerical schemes. Spatial orthogonality of the eigenmodes permits determination of each mode\\'s temporal expansion coefficient by integrating the product of the electric field and the mode over the wave port. The temporal mode coefficients are then Fourier transformed to accurately calculate the S-parameters corresponding to different modes. Additionally, to generalize the DGTD to manipulate dispersive media, the auxiliary differential equation method is employed. This is done by introducing a time-dependent polarization volume current as an auxiliary unknown and the constitutive relation between this current and the electric field as an auxiliary equation. Consequently, computationally expensive temporal convolution is avoided. Various numerical examples, which demonstrate the applicability, robustness, and accuracy of the proposed method, are presented.

Asymptotic Expansions - Methods and Applications

International Nuclear Information System (INIS)

Harlander, R.

1999-01-01

Different viewpoints on the asymptotic expansion of Feynman diagrams are reviewed. The relations between the field theoretic and diagrammatic approaches are sketched. The focus is on problems with large masses or large external momenta. Several recent applications also for other limiting cases are touched upon. Finally, the pros and cons of the different approaches are briefly discussed. (author)

The Sturmian expansion: A well-depth-method for orbitals in a deformed potential

International Nuclear Information System (INIS)

Bang, J.M.; Vaagen, J.S.

1980-01-01

The Sturmian expansion method has over the years successfully been used to generate orbitals in a deformed potential. In this paper we review the method in detail including more recent extentions. The convergence properties are discussed in terms of examples of current interest for nucleon-transfer reactions. Comparisons with other methods are also made. (orig.)

Asymptotic expansions of Mathieu functions in wave mechanics

International Nuclear Information System (INIS)

Hunter, G.; Kuriyan, M.

1976-01-01

Solutions of the radial Schroedinger equation containing a polarization potential r -4 are expanded in a form appropriate for large values of r. These expansions of the Mathieu functions are used in association with the numerical solution of the Schroedinger equation to impose the asymptotic boundary condition in the case of bound states, and to extract phase shifts in the case of scattering states

New exact travelling wave solutions of generalised sinh- Gordon and (2 + 1-dimensional ZK-BBM equations

Directory of Open Access Journals (Sweden)

Sachin Kumar

2012-10-01

Full Text Available Exact travelling wave solutions have been established for generalised sinh-Gordon andgeneralised (2+1 dimensional ZK-BBM equations by using GG      expansion method whereG  G( satisfies a second-order linear ordinary differential equation. The travelling wave solutionsare expressed by hyperbolic, trigonometric and rational functions.

Solution of the Helmholtz-Poincare Wave Equation using the coupled boundary integral equations and optimal surface eigenfunctions

International Nuclear Information System (INIS)

Werby, M.F.; Broadhead, M.K.; Strayer, M.R.; Bottcher, C.

1992-01-01

The Helmholtz-Poincarf Wave Equation (H-PWE) arises in many areas of classical wave scattering theory. In particular it can be found for the cases of acoustical scattering from submerged bounded objects and electromagnetic scattering from objects. The extended boundary integral equations (EBIE) method is derived from considering both the exterior and interior solutions of the H-PWECs. This coupled set of expressions has the advantage of not only offering a prescription for obtaining a solution for the exterior scattering problem, but it also obviates the problem of irregular values corresponding to fictitious interior eigenvalues. Once the coupled equations are derived, they can be obtained in matrix form by expanding all relevant terms in partial wave expansions, including a bi-orthogonal expansion of the Green's function. However some freedom in the choice of the surface expansion is available since the unknown surface quantities may be expanded in a variety of ways so long as closure is obtained. Out of many possible choices, we develop an optimal method to obtain such expansions which is based on the optimum eigenfunctions related to the surface of the object. In effect, we convert part of the problem (that associated with the Fredholms integral equation of the first kind) an eigenvalue problem of a related Hermitian operator. The methodology will be explained in detail and examples will be presented

Thermodynamics of non-ideal QGP using Mayers cluster expansion method

International Nuclear Information System (INIS)

Prasanth, J.P; Simji, P.; Bannur, Vishnu M.

2013-01-01

The Quark gluon plasma (QGP) is the state in which the individual hadrons dissolve into a system of free (or almost free) quarks and gluons in strongly compressed system at high temperature. The present paper aims to calculate the critical temperature at which a non-ideal three quark plasma condenses into droplet of three quarks (i.e., into a liquid of baryons) using Mayers cluster expansion method

Wave propagation retrieval method for chiral metamaterials

DEFF Research Database (Denmark)

Andryieuski, Andrei; Malureanu, Radu; Lavrinenko, Andrei

2010-01-01

In this paper we present the wave propagation method for the retrieving of effective properties of media with circularly polarized eigenwaves, in particularly for chiral metamaterials. The method is applied for thick slabs and provides bulk effective parameters. Its strong sides are the absence...

Spectral analysis of surface waves method to assess shear wave velocity within centrifuge models

Science.gov (United States)

Murillo, Carol Andrea; Thorel, Luc; Caicedo, Bernardo

2009-06-01

The method of the spectral analysis of surface waves (SASW) is tested out on reduced scale centrifuge models, with a specific device, called the mini Falling Weight, developed for this purpose. Tests are performed on layered materials made of a mixture of sand and clay. The shear wave velocity VS determined within the models using the SASW is compared with the laboratory measurements carried out using the bender element test. The results show that the SASW technique applied to centrifuge testing is a relevant method to characterize VS near the surface.

Thermal expansion of coking coals

Energy Technology Data Exchange (ETDEWEB)

Orlik, M.; Klimek, J. (Vyzkumny a Zkusebni Ustav Nova Hut, Ostrava (Czechoslovakia))

1992-12-01

Analyzes expansion of coal mixtures in coke ovens during coking. Methods for measuring coal expansion on both a laboratory and pilot plant scale are comparatively evaluated. The method, developed, tested and patented in Poland by the Institute for Chemical Coal Processing in Zabrze (Polish standard PN-73/G-04522), is discussed. A laboratory device developed by the Institute for measuring coal expansion is characterized. Expansion of black coal from 10 underground mines in the Ostrava-Karvina coal district and from 9 coal mines in the Upper Silesia basin in Poland is comparatively evaluated. Investigations show that coal expansion reaches a maximum for coal types with a volatile matter ranging from 20 to 25%. With increasing volatile matter in coal, its expansion decreases. Coal expansion increases with increasing swelling index. Coal expansion corresponds with coal dilatation. With increasing coal density its expansion increases. Coal mixtures should be selected in such a way that their expansion does not cause a pressure exceeding 40 MPa. 11 refs.

Study of the asymptotic expansion of multiple integrals in mathematical physics

International Nuclear Information System (INIS)

Chako, N.

1968-01-01

We have applied the method of stationary phase to evaluate double and multiple integrals of the type: (A) U(k) = g(x)e ikφ(x) d(x), (x)=(x 1 ,..., x n ) for large values of the parameter k. In the first part we have established in a rigorous manner the stationary phase method to double and multiple integrals of type (A). Furthermore we have obtained an asymptotic expansion of (A), if the amplitude and phase functions can be developed in a canonical form near the vicinity of critical or stationary points of the integral. This development contains as particular cases all those which are important in physical applications, especially, to diffraction and scattering of electromagnetic and corpuscular waves by optical systems, diffracting bodies and potential scatterers. In the second part we have considered the problem of convergence of the expansion of the principal contribution to the integral in the asymptotic sense of Poincare. The proof is based on the increasing method used in mathematical analysis. The third part is devoted to the derivation of various asymptotic series due to different types of critical or stationary points associated with the amplitude and phase functions. In the fourth part we have generalized the method to multiple integrals and to the case where the parameter k enter implicitly in the phase function The latter type of integrals extend the scope of the former type to a number of important physical problems; for instance, to the propagation of waves in dispersive and absorbing media. In the last chapter we have made a study and compared the results obtained by the application of the stationary phase method to the integrals (double) of diffraction and the results derived by using the Young-Rubinowicz method. Result of our analysis shows the equivalence of the two methods of approach to the problems of diffraction based, on one hand, on the Fresnel-Kirchhoff theory and, on the other hand, the Young-Rubinowicz theory, provided one interprets in

Algebraic method for constructing singular steady solitary waves: a case study

Science.gov (United States)

Clamond, Didier; Dutykh, Denys; Galligo, André

2016-07-01

This article describes the use of algebraic methods in a phase plane analysis of ordinary differential equations. The method is illustrated by the study of capillary-gravity steady surface waves propagating in shallow water. We consider the (fully nonlinear, weakly dispersive) Serre-Green-Naghdi equation with surface tension, because it provides a tractable model that, at the same time, is not too simple, so interest in the method can be emphasized. In particular, we analyse a special class of solutions, the solitary waves, which play an important role in many fields of physics. In capillary-gravity regime, there are two kinds of localized infinitely smooth travelling wave solutions-solitary waves of elevation and of depression. However, if we allow the solitary waves to have an angular point, then the `zoology' of solutions becomes much richer, and the main goal of this study is to provide a complete classification of such singular localized solutions using the methods of the effective algebraic geometry.

Longitudinal propagation of nonlinear surface Alfven waves at a magnetic interface in a compressible atmosphere

Energy Technology Data Exchange (ETDEWEB)

Ruderman, M S

1988-08-01

Nonlinear Alfven surface wave propagation at a magnetic interface in a compressible fluid is considered. It is supposed that the magnetic field directions at both sides of the interface and the direction of wave propagation coincide. The equation governing time-evolution of nonlinear small-amplitude waves is derived by the method of multiscale expansions. This equation is similar to the equation for nonlinear Alfven surface waves in an incompressible fluid derived previously. The numerical solution of the equation shows that a sinusoidal disturbance overturns, i.e. infinite gradients arise.

Asymptotic expansion of unsteady gravity flow of a power-law fluid ...

African Journals Online (AJOL)

We present a paper on the asymptotic expansion of unsteady non-linear rheological effects of a power-law fluid under gravity. The fluid flows through a porous medium. The asymptotic expansion is employed to obtain solution of the nonlinear problem. The results show the existence of traveling waves. It is assumed that the ...

A robust absorbing layer method for anisotropic seismic wave modeling

Energy Technology Data Exchange (ETDEWEB)

Métivier, L., E-mail: ludovic.metivier@ujf-grenoble.fr [LJK, CNRS, Université de Grenoble, BP 53, 38041 Grenoble Cedex 09 (France); ISTerre, Université de Grenoble I, BP 53, 38041 Grenoble Cedex 09 (France); Brossier, R. [ISTerre, Université de Grenoble I, BP 53, 38041 Grenoble Cedex 09 (France); Labbé, S. [LJK, CNRS, Université de Grenoble, BP 53, 38041 Grenoble Cedex 09 (France); Operto, S. [Géoazur, Université de Nice Sophia-Antipolis, CNRS, IRD, OCA, Villefranche-sur-Mer (France); Virieux, J. [ISTerre, Université de Grenoble I, BP 53, 38041 Grenoble Cedex 09 (France)

2014-12-15

When applied to wave propagation modeling in anisotropic media, Perfectly Matched Layers (PML) exhibit instabilities. Incoming waves are amplified instead of being absorbed. Overcoming this difficulty is crucial as in many seismic imaging applications, accounting accurately for the subsurface anisotropy is mandatory. In this study, we present the SMART layer method as an alternative to PML approach. This method is based on the decomposition of the wavefield into components propagating inward and outward the domain of interest. Only outgoing components are damped. We show that for elastic and acoustic wave propagation in Transverse Isotropic media, the SMART layer is unconditionally dissipative: no amplification of the wavefield is possible. The SMART layers are not perfectly matched, therefore less accurate than conventional PML. However, a reasonable increase of the layer size yields an accuracy similar to PML. Finally, we illustrate that the selective damping strategy on which is based the SMART method can prevent the generation of spurious S-waves by embedding the source in a small zone where only S-waves are damped.

A robust absorbing layer method for anisotropic seismic wave modeling

International Nuclear Information System (INIS)

Métivier, L.; Brossier, R.; Labbé, S.; Operto, S.; Virieux, J.

2014-01-01

When applied to wave propagation modeling in anisotropic media, Perfectly Matched Layers (PML) exhibit instabilities. Incoming waves are amplified instead of being absorbed. Overcoming this difficulty is crucial as in many seismic imaging applications, accounting accurately for the subsurface anisotropy is mandatory. In this study, we present the SMART layer method as an alternative to PML approach. This method is based on the decomposition of the wavefield into components propagating inward and outward the domain of interest. Only outgoing components are damped. We show that for elastic and acoustic wave propagation in Transverse Isotropic media, the SMART layer is unconditionally dissipative: no amplification of the wavefield is possible. The SMART layers are not perfectly matched, therefore less accurate than conventional PML. However, a reasonable increase of the layer size yields an accuracy similar to PML. Finally, we illustrate that the selective damping strategy on which is based the SMART method can prevent the generation of spurious S-waves by embedding the source in a small zone where only S-waves are damped

Rapid expansion method (REM) for time‐stepping in reverse time migration (RTM)

KAUST Repository

Pestana, Reynam C.; Stoffa, Paul L.

2009-01-01

an analytical approximation for the Bessel function where we assume that the time step is sufficiently small. From this derivation we find that if we consider only the first two Chebyshev polynomials terms in the rapid expansion method we can obtain the second

Closed form solutions of two time fractional nonlinear wave equations

Science.gov (United States)

Akbar, M. Ali; Ali, Norhashidah Hj. Mohd.; Roy, Ripan

2018-06-01

In this article, we investigate the exact traveling wave solutions of two nonlinear time fractional wave equations. The fractional derivatives are described in the sense of conformable fractional derivatives. In addition, the traveling wave solutions are accomplished in the form of hyperbolic, trigonometric, and rational functions involving free parameters. To investigate such types of solutions, we implement the new generalized (G‧ / G) -expansion method. The extracted solutions are reliable, useful and suitable to comprehend the optimal control problems, chaotic vibrations, global and local bifurcations and resonances, furthermore, fission and fusion phenomena occur in solitons, the relativistic energy-momentum relation, scalar electrodynamics, quantum relativistic one-particle theory, electromagnetic interactions etc. The results reveal that the method is very fruitful and convenient for exploring nonlinear differential equations of fractional order treated in theoretical physics.

Nonlinear modulation of torsional waves in elastic rod. [Instability

Energy Technology Data Exchange (ETDEWEB)

Hirao, M; Sugimoto, N [Osaka Univ., Toyonaka (Japan). Faculty of Engineering Science

1977-06-01

Nonlinear Schroedinger equation, which describes the nonlinear modulation of dispersive torsional waves in an elastic rod of circular cross-section, is derived by the derivative expansion method. It is found, for the lowest dispersive mode, that the modulational instability occurs except in the range of the carrier wavenumber, 2.799waves can propagate simultaneously, the second-harmonic resonance takes place and then the nonlinear Schroedinger equation is no longer valid. In this case, another system of equations is derived, which governs both the wave amplitudes involved in this resonance between the fundamental torsional and its second-harmonic longitudinal modes.

A novel method for predicting the power outputs of wave energy converters

Science.gov (United States)

Wang, Yingguang

2018-03-01

This paper focuses on realistically predicting the power outputs of wave energy converters operating in shallow water nonlinear waves. A heaving two-body point absorber is utilized as a specific calculation example, and the generated power of the point absorber has been predicted by using a novel method (a nonlinear simulation method) that incorporates a second order random wave model into a nonlinear dynamic filter. It is demonstrated that the second order random wave model in this article can be utilized to generate irregular waves with realistic crest-trough asymmetries, and consequently, more accurate generated power can be predicted by subsequently solving the nonlinear dynamic filter equation with the nonlinearly simulated second order waves as inputs. The research findings demonstrate that the novel nonlinear simulation method in this article can be utilized as a robust tool for ocean engineers in their design, analysis and optimization of wave energy converters.

Classical scattering theory of waves from the view point of an eigenvalue problem and application to target identification

International Nuclear Information System (INIS)

Bottcher, C.; Strayer, M.R.; Werby, M.F.

1993-01-01

The Helmholtz-Poincare Wave Equation (H-PWE) arises in many areas of classical wave scattering theory. In particular it can be found for the cases of acoustical scattering from submerged bounded objects and electromagnetic scattering from objects. The extended boundary integral equations (EBIE) method is derived from considering both the exterior and interior solutions of the H-PWE's. This coupled set of expressions has the advantage of not only offering a prescription for obtaining a solution for the exterior scattering problem, but it also obviates the problem of irregular values corresponding to fictitious interior eigenvalues. Once the coupled equations are derived, they can by obtained in matrix form be expanding all relevant terms in partial wave expansions, including a biorthogonal expansion of the Green function. However some freedom of choice in the choice of the surface expansion is available since the unknown surface quantities may be expanded in a variety of ways to long as closure is obtained. Out of many possible choices, we develop an optimal method to obtain such expansions which is based on the optimum eigenfunctions related to the surface of the object. In effect, we convert part of the problem (that associated with the Fredholms integral equation of the first kind) an eigenvalue problem of a related Hermition operator. The methodology will be explained in detail and examples will be presented

Forecasting ocean wave energy: A Comparison of the ECMWF wave model with time series methods

DEFF Research Database (Denmark)

Reikard, Gordon; Pinson, Pierre; Bidlot, Jean

2011-01-01

Recently, the technology has been developed to make wave farms commercially viable. Since electricity is perishable, utilities will be interested in forecasting ocean wave energy. The horizons involved in short-term management of power grids range from as little as a few hours to as long as several...... days. In selecting a method, the forecaster has a choice between physics-based models and statistical techniques. A further idea is to combine both types of models. This paper analyzes the forecasting properties of a well-known physics-based model, the European Center for Medium-Range Weather Forecasts...... (ECMWF) Wave Model, and two statistical techniques, time-varying parameter regressions and neural networks. Thirteen data sets at locations in the Atlantic and Pacific Oceans and the Gulf of Mexico are tested. The quantities to be predicted are the significant wave height, the wave period, and the wave...

Detonative propagation and accelerative expansion of the Crab Nebula shock front.

Science.gov (United States)

Gao, Yang; Law, Chung K

2011-10-21

The accelerative expansion of the Crab Nebula's outer envelope is a mystery in dynamics, as a conventional expanding blast wave decelerates when bumping into the surrounding interstellar medium. Here we show that the strong relativistic pulsar wind bumping into its surrounding nebula induces energy-generating processes and initiates a detonation wave that propagates outward to form the current outer edge, namely, the shock front, of the nebula. The resulting detonation wave, with a reactive downstream, then provides the needed power to maintain propagation of the shock front. Furthermore, relaxation of the curvature-induced reduction of the propagation velocity from the initial state of formation to the asymptotic, planar state of Chapman-Jouguet propagation explains the observed accelerative expansion. Potential richness in incorporating reactive fronts in the description of various astronomical phenomena is expected. © 2011 American Physical Society

Mathematical Methods in Wave Propagation: Part 2--Non-Linear Wave Front Analysis

Science.gov (United States)

Jeffrey, Alan

1971-01-01

The paper presents applications and methods of analysis for non-linear hyperbolic partial differential equations. The paper is concluded by an account of wave front analysis as applied to the piston problem of gas dynamics. (JG)

Stochastic series expansion simulation of the t -V model

Science.gov (United States)

Wang, Lei; Liu, Ye-Hua; Troyer, Matthias

2016-04-01

We present an algorithm for the efficient simulation of the half-filled spinless t -V model on bipartite lattices, which combines the stochastic series expansion method with determinantal quantum Monte Carlo techniques widely used in fermionic simulations. The algorithm scales linearly in the inverse temperature, cubically with the system size, and is free from the time-discretization error. We use it to map out the finite-temperature phase diagram of the spinless t -V model on the honeycomb lattice and observe a suppression of the critical temperature of the charge-density-wave phase in the vicinity of a fermionic quantum critical point.

Macroscopic quantum waves in non local theories

International Nuclear Information System (INIS)

Ventura, I.

1979-01-01

By means of an expansion in the density, it is shown that Macroscopic Quantum Waves also appear in non local theories. This result reinforces the conjecture that these waves should exist in liquid 4 He [pt

Thermal Expansion and Magnetostriction Measurements at Cryogenic Temperature Using the Strain Gauge Method

Directory of Open Access Journals (Sweden)

Wei Wang

2018-03-01

Full Text Available Thermal expansion and magnetostriction, the strain responses of a material to temperature and a magnetic field, especially properties at low temperature, are extremely useful to study electronic and phononic properties, phase transitions, quantum criticality, and other interesting phenomena in cryogenic engineering and materials science. However, traditional dilatometers cannot provide magnetic field and ultra-low temperature (<77 K environment easily. This paper describes the design and test results of thermal expansion and magnetostriction at cryogenic temperature using the strain gauge method based on a Physical Properties Measurements System (PPMS. The interfacing software and automation were developed using LabVIEW. The sample temperature range can be tuned continuously between 1.8 and 400 K. With this PPMS-aided measuring system, we can observe temperature and magnetic field dependence of the linear thermal expansion of different solid materials easily and accurately.

Thermal expansion and magnetostriction measurements at cryogenic temperature using the strain gage method

Science.gov (United States)

Wang, Wei; Liu, Huiming; Huang, Rongjin; Zhao, Yuqiang; Huang, Chuangjun; Guo, Shibin; Shan, Yi; Li, Laifeng

2018-03-01

Thermal expansion and magnetostriction, the strain responses of a material to temperature and a magnetic field, especially properties at low temperature, are extremely useful to study electronic and phononic properties, phase transitions, quantum criticality, and other interesting phenomena in cryogenic engineering and materials science. However, traditional dilatometers cannot provide magnetic field and ultra low temperature (＜77 K) environment easily. This paper describes the design and test results of thermal expansion and magnetostriction at cryogenic temperature using the strain gage method based on a Physical Properties Measurements System (PPMS). The interfacing software and automation were developed using LabVIEW. The sample temperature range can be tuned continuously between 1.8 K and 400 K. With this PPMS-aided measuring system, we can observe temperature and magnetic field dependence of the linear thermal expansion of different solid materials easily and accurately.

Ground penetrating radar antenna measurements based on plane-wave expansions

DEFF Research Database (Denmark)

Lenler-Eriksen, Hans-Rudolph; Meincke, Peter

2005-01-01

The plane-wave transmitting spectrum of the system consisting of the ground penetrating radar (GPR) antenna and the air-soil interface is measured using a loop buried in the soil. The plane-wave spectrum is used to determine various parameters characterizing the radiation of the GPR antenna...

Ion-acoustic waves in ultracold neutral plasmas: Modulational instability and dissipative rogue waves

Energy Technology Data Exchange (ETDEWEB)

El-Tantawy, S.A., E-mail: samireltantawy@yahoo.com

2017-02-26

Progress is reported on the modulational instability (MI) of ion-acoustic waves (IAWs) and dissipative rogue waves (RWs) in ultracold neutral plasmas (UNPs). The UNPs consist of inertial ions fluid and Maxwellian inertialess hot electrons, and the presence of an ion kinematic viscosity is allowed. For this purpose, a modified nonlinear Schrödinger equation (NLSE) is derived and then solved analytically to show the occurrence of MI. It is found that the (in)stability regions of the wavepacks are dependent on time due to of the existence of the dissipative term. The existing regions of the MI of the IAWs are inventoried precisely. After that, we use a suitable transformation to convert the modified NLSE into the normal NLSE whose analytical solutions for rogue waves are known. The rogue wave propagation condition and its behavior are discussed. The impact of the relevant physical parameters on the profile of the RWs is examined. - Highlights: • UNPs are modeled by the phenomenological generalized hydrodynamic equations. • The derivative expansion method has been employed in order to derive a modified-NLSE. • A suitable transformation is used to transform the modified-NLSE into the standard NLSE. • The effect of the ion viscosity on the modulational instability and rogue waves is investigated.

Numerical simulation of electromagnetic wave propagation using time domain meshless method

International Nuclear Information System (INIS)

Ikuno, Soichiro; Fujita, Yoshihisa; Itoh, Taku; Nakata, Susumu; Nakamura, Hiroaki; Kamitani, Atsushi

2012-01-01

The electromagnetic wave propagation in various shaped wave guide is simulated by using meshless time domain method (MTDM). Generally, Finite Differential Time Domain (FDTD) method is applied for electromagnetic wave propagation simulation. However, the numerical domain should be divided into rectangle meshes if FDTD method is applied for the simulation. On the other hand, the node disposition of MTDM can easily describe the structure of arbitrary shaped wave guide. This is the large advantage of the meshless time domain method. The results of computations show that the damping rate is stably calculated in case with R < 0.03, where R denotes a support radius of the weight function for the shape function. And the results indicate that the support radius R of the weight functions should be selected small, and monomials must be used for calculating the shape functions. (author)

Separate P‐ and SV‐wave equations for VTI media

KAUST Repository

Pestana, Reynam C.; Ursin, Bjø rn; Stoffa, Paul L.

2011-01-01

In isotropic media we use the scalar acoustic wave equation to perform reverse time migration RTM of the recorded pressure wavefleld data. In anisotropic media P- and SV-waves are coupled and the elastic wave equation should be used for RTM. However, an acoustic anisotropic wave equation is often used instead. This results in significant shear wave energy in both modeling and RTM. To avoid this undesired SV-wave energy, we propose a different approach to separate P- and SV-wave components for vertical transversely isotropic VTI media. We derive independent pseudo-differential wave equations for each mode. The derived equations for P- and SV-waves are stable and reduce to the isotropic case. The equations presented here can be effectively used to model and migrate seismic data in VTI media where ε - δ is small. The SV-wave equation we develop is now well-posed and triplications in the SV wavefront are removed resulting in stable wave propagation. We show modeling and RTM results using the derived pure P-wave mode in complex VTI media and use the rapid expansion method REM to propagate the waveflelds in time. © 2011 Society of Exploration Geophysicists.

Macroscopic quantum waves in non local theories

International Nuclear Information System (INIS)

Ventura, I.

1979-01-01

By means of an expansion in the density, it is shown that Macroscopic Quantum Waves also apear in non local theories. This result reinforces the conjecture that these waves should exist in liquid 4 He. (Author) [pt

Algebraic internal wave solitons and the integrable Calogero--Moser--Sutherland N-body problem

International Nuclear Information System (INIS)

Chen, H.H.; Lee, Y.C.; Pereira, N.R.

1979-01-01

The Benjamin--Ono equation that describes nonlinear internal waves in a stratified fluid is solved by a pole expansion method. The dynamics of poles which characterize solitons is shown to be identical to the well-known integrable N-body problem of Calogero, Moser, and Sutherland

Design of materials with extreme thermal expansion using a three-phase topology optimization method

DEFF Research Database (Denmark)

Sigmund, Ole; Torquato, S.

1997-01-01

We show how composites with extremal or unusual thermal expansion coefficients can be designed using a numerical topology optimization method. The composites are composed of two different material phases and void. The optimization method is illustrated by designing materials having maximum therma...

Brillouin Corrosion Expansion Sensors for Steel Reinforced Concrete Structures Using a Fiber Optic Coil Winding Method

Directory of Open Access Journals (Sweden)

Xingjun Lv

2011-11-01

Full Text Available In this paper, a novel kind of method to monitor corrosion expansion of steel rebars in steel reinforced concrete structures named fiber optic coil winding method is proposed, discussed and tested. It is based on the fiber optical Brillouin sensing technique. Firstly, a strain calibration experiment is designed and conducted to obtain the strain coefficient of single mode fiber optics. Results have shown that there is a good linear relationship between Brillouin frequency and applied strain. Then, three kinds of novel fiber optical Brillouin corrosion expansion sensors with different fiber optic coil winding packaging schemes are designed. Sensors were embedded into concrete specimens to monitor expansion strain caused by steel rebar corrosion, and their performance was studied in a designed electrochemical corrosion acceleration experiment. Experimental results have shown that expansion strain along the fiber optic coil winding area can be detected and measured by the three kinds of sensors with different measurement range during development the corrosion. With the assumption of uniform corrosion, diameters of corrosion steel rebars were obtained using calculated average strains. A maximum expansion strain of 6,738 με was monitored. Furthermore, the uniform corrosion analysis model was established and the evaluation formula to evaluate mass loss rate of steel rebar under a given corrosion rust expansion rate was derived. The research has shown that three kinds of Brillouin sensors can be used to monitor the steel rebar corrosion expansion of reinforced concrete structures with good sensitivity, accuracy and monitoring range, and can be applied to monitor different levels of corrosion. By means of this kind of monitoring technique, quantitative corrosion expansion monitoring can be carried out, with the virtues of long durability, real-time monitoring and quasi-distribution monitoring.

Brillouin corrosion expansion sensors for steel reinforced concrete structures using a fiber optic coil winding method.

Science.gov (United States)

Zhao, Xuefeng; Gong, Peng; Qiao, Guofu; Lu, Jie; Lv, Xingjun; Ou, Jinping

2011-01-01

In this paper, a novel kind of method to monitor corrosion expansion of steel rebars in steel reinforced concrete structures named fiber optic coil winding method is proposed, discussed and tested. It is based on the fiber optical Brillouin sensing technique. Firstly, a strain calibration experiment is designed and conducted to obtain the strain coefficient of single mode fiber optics. Results have shown that there is a good linear relationship between Brillouin frequency and applied strain. Then, three kinds of novel fiber optical Brillouin corrosion expansion sensors with different fiber optic coil winding packaging schemes are designed. Sensors were embedded into concrete specimens to monitor expansion strain caused by steel rebar corrosion, and their performance was studied in a designed electrochemical corrosion acceleration experiment. Experimental results have shown that expansion strain along the fiber optic coil winding area can be detected and measured by the three kinds of sensors with different measurement range during development the corrosion. With the assumption of uniform corrosion, diameters of corrosion steel rebars were obtained using calculated average strains. A maximum expansion strain of 6,738 με was monitored. Furthermore, the uniform corrosion analysis model was established and the evaluation formula to evaluate mass loss rate of steel rebar under a given corrosion rust expansion rate was derived. The research has shown that three kinds of Brillouin sensors can be used to monitor the steel rebar corrosion expansion of reinforced concrete structures with good sensitivity, accuracy and monitoring range, and can be applied to monitor different levels of corrosion. By means of this kind of monitoring technique, quantitative corrosion expansion monitoring can be carried out, with the virtues of long durability, real-time monitoring and quasi-distribution monitoring.

Asymptotic expansions for solitary gravity-capillary waves in two and three dimensions

International Nuclear Information System (INIS)

Ablowitz, M J; Haut, T S

2010-01-01

High-order asymptotic series are obtained for gravity-capillary solitary waves, where the first term in the series is the well-known sech 2 solution of the KdV equation. The asymptotic series is used, with nine terms included, to investigate the effects of surface tension on the height and energy of large amplitude waves, and waves close to the solitary version of Stokes' extreme wave. In particular, for surface tension below a critical value, the solitary wave with the maximum energy is obtained. For large surface tension, the series is also used to study the energy related to the solitary waves of depression. Energy considerations suggest that, for large enough surface tension, there are solitary waves that can get close to the fluid bottom. Comparisons are also made with recent experiments.

Process Investigation of Tube Expansion by Gas Detonation

OpenAIRE

Bach, F.-W.; Beerwald, C.; Brosius, A.; Gershteyn, G.; Hermes, M.; Kleiner, M.; Olivier, H.; Weber, M.

2006-01-01

The present paper deals with the expansion of tubes by direct application of gas detonation waves, i.e. the gas is both pressure medium and energy source. After an introduction to gas detonation forming, measurements of the motion process and the internal pressures are presented. Results of free expansion and of forming into a die are thoroughly studied and compared to the results of quasi-static burst tests and hydroforming. Using pure aluminum Al99.5 and a medium strength alloy AlMgSi1, ...

Full-wave calculation of fast-wave current drive in tokamaks including kparallel upshifts

International Nuclear Information System (INIS)

Jaeger, E.F.; Batchelor, D.B.

1991-01-01

Numerical calculations of fast-wave current drive (FWCD) efficiency have generally been of two types: ray tracing or global wave calculations. Ray tracing shows that the projection of the wave number (k parallel) along the magnetic field can vary greatly over a ray trajectory, particularly when the launch point is above or below the equatorial plane. As the wave penetrates toward the center of the plasma, k parallel increases, causing a decrease in the parallel phase speed and a corresponding decrease in the current drive efficiency, γ. But the assumptions of geometrical optics, namely short wavelength and strong single-pass absorption, are not greatly applicable in FWCD scenarios. Eigenmode structure, which is ignored in ray tracing, can play an important role in determining electric field strength and Landau damping rates. In such cases, a full-wave or global solution for the wave fields is desirable. In full-wave calculations such as ORION k parallel appear as a differential operator (rvec B·∇) in the argument of the plasma dispersion function. Since this leads to a differential system of infinite order, such codes of necessity assume k parallel ∼ k var-phi = const, where k var-phi is the toroidal wave number. Thus, it is not possible to correctly include effects of the poloidal magnetic field on k parallel. The problem can be alleviated by expressing the electric field as a superposition of poloidal modes, in which case k parallel is purely algebraic. This paper describes a new full-wave calculation, Poloidal Ion Cyclotron Expansion Solution, which uses poloidal and toroidal mode expansions to solve the wave equation in general flux coordinates. The calculation includes a full solution for E parallel and uses a reduced-order form of the plasma conductivity tensor to eliminate numerical problems associated with resolution of the very short wavelength ion Bernstein wave

A nonlinear analytic function expansion nodal method for transient calculations

Energy Technology Data Exchange (ETDEWEB)

Joo, Han Gyn; Park, Sang Yoon; Cho, Byung Oh; Zee, Sung Quun [Korea Atomic Energy Research Institute, Taejon (Korea, Republic of)

1999-12-31

The nonlinear analytic function expansion nodal (AFEN) method is applied to the solution of the time-dependent neutron diffusion equation. Since the AFEN method requires both the particular solution and the homogeneous solution to the transient fixed source problem, the derivation of the solution method is focused on finding the particular solution efficiently. To avoid complicated particular solutions, the source distribution is approximated by quadratic polynomials and the transient source is constructed such that the error due to the quadratic approximation is minimized, In addition, this paper presents a new two-node solution scheme that is derived by imposing the constraint of current continuity at the interface corner points. The method is verified through a series of application to the NEACRP PWR rod ejection benchmark problems. 6 refs., 2 figs., 1 tab. (Author)

A nonlinear analytic function expansion nodal method for transient calculations

Energy Technology Data Exchange (ETDEWEB)

Joo, Han Gyn; Park, Sang Yoon; Cho, Byung Oh; Zee, Sung Quun [Korea Atomic Energy Research Institute, Taejon (Korea, Republic of)

1998-12-31

The nonlinear analytic function expansion nodal (AFEN) method is applied to the solution of the time-dependent neutron diffusion equation. Since the AFEN method requires both the particular solution and the homogeneous solution to the transient fixed source problem, the derivation of the solution method is focused on finding the particular solution efficiently. To avoid complicated particular solutions, the source distribution is approximated by quadratic polynomials and the transient source is constructed such that the error due to the quadratic approximation is minimized, In addition, this paper presents a new two-node solution scheme that is derived by imposing the constraint of current continuity at the interface corner points. The method is verified through a series of application to the NEACRP PWR rod ejection benchmark problems. 6 refs., 2 figs., 1 tab. (Author)

Wave field restoration using three-dimensional Fourier filtering method.

Science.gov (United States)

Kawasaki, T; Takai, Y; Ikuta, T; Shimizu, R

2001-11-01

A wave field restoration method in transmission electron microscopy (TEM) was mathematically derived based on a three-dimensional (3D) image formation theory. Wave field restoration using this method together with spherical aberration correction was experimentally confirmed in through-focus images of amorphous tungsten thin film, and the resolution of the reconstructed phase image was successfully improved from the Scherzer resolution limit to the information limit. In an application of this method to a crystalline sample, the surface structure of Au(110) was observed in a profile-imaging mode. The processed phase image showed quantitatively the atomic relaxation of the topmost layer.

Synthesis of Numerical Methods for Modeling Wave Energy Converter-Point Absorbers: Preprint

Energy Technology Data Exchange (ETDEWEB)

Li, Y.; Yu, Y. H.

2012-05-01

During the past few decades, wave energy has received significant attention among all ocean energy formats. Industry has proposed hundreds of prototypes such as an oscillating water column, a point absorber, an overtopping system, and a bottom-hinged system. In particular, many researchers have focused on modeling the floating-point absorber as the technology to extract wave energy. Several modeling methods have been used such as the analytical method, the boundary-integral equation method, the Navier-Stokes equations method, and the empirical method. However, no standardized method has been decided. To assist the development of wave energy conversion technologies, this report reviews the methods for modeling the floating-point absorber.

Engineered high expansion glass-ceramics having near linear thermal strain and methods thereof

Energy Technology Data Exchange (ETDEWEB)

Dai, Steve Xunhu; Rodriguez, Mark A.; Lyon, Nathanael L.

2018-01-30

The present invention relates to glass-ceramic compositions, as well as methods for forming such composition. In particular, the compositions include various polymorphs of silica that provide beneficial thermal expansion characteristics (e.g., a near linear thermal strain). Also described are methods of forming such compositions, as well as connectors including hermetic seals containing such compositions.

Dynamic Response of Underground Circular Lining Tunnels Subjected to Incident P Waves

Directory of Open Access Journals (Sweden)

Hua Xu

2014-01-01

Full Text Available Dynamic stress concentration in tunnels and underground structures during earthquakes often leads to serious structural damage. A series solution of wave equation for dynamic response of underground circular lining tunnels subjected to incident plane P waves is presented by Fourier-Bessel series expansion method in this paper. The deformation and stress fields of the whole medium of surrounding rock and tunnel were obtained by solving the equations of seismic wave propagation in an elastic half space. Based on the assumption of a large circular arc, a series of solutions for dynamic stress were deduced by using a wave function expansion approach for a circular lining tunnel in an elastic half space rock medium subjected to incident plane P waves. Then, the dynamic response of the circular lining tunnel was obtained by solving a series of algebraic equations after imposing its boundary conditions for displacement and stress of the circular lining tunnel. The effects of different factors on circular lining rock tunnels, including incident frequency, incident angle, buried depth, rock conditions, and lining stiffness, were derived and several application examples are presented. The results may provide a good reference for studies on the dynamic response and aseismic design of tunnels and underground structures.

Exact traveling wave solutions of modified KdV-Zakharov-Kuznetsov equation and viscous Burgers equation.

Science.gov (United States)

Islam, Md Hamidul; Khan, Kamruzzaman; Akbar, M Ali; Salam, Md Abdus

2014-01-01

Mathematical modeling of many physical systems leads to nonlinear evolution equations because most physical systems are inherently nonlinear in nature. The investigation of traveling wave solutions of nonlinear partial differential equations (NPDEs) plays a significant role in the study of nonlinear physical phenomena. In this article, we construct the traveling wave solutions of modified KDV-ZK equation and viscous Burgers equation by using an enhanced (G '/G) -expansion method. A number of traveling wave solutions in terms of unknown parameters are obtained. Derived traveling wave solutions exhibit solitary waves when special values are given to its unknown parameters. 35C07; 35C08; 35P99.

Analytical solution for the transient wave propagation of a buried cylindrical P-wave line source in a semi-infinite elastic medium with a fluid surface layer

Science.gov (United States)

Shan, Zhendong; Ling, Daosheng

2018-02-01

This article develops an analytical solution for the transient wave propagation of a cylindrical P-wave line source in a semi-infinite elastic solid with a fluid layer. The analytical solution is presented in a simple closed form in which each term represents a transient physical wave. The Scholte equation is derived, through which the Scholte wave velocity can be determined. The Scholte wave is the wave that propagates along the interface between the fluid and solid. To develop the analytical solution, the wave fields in the fluid and solid are defined, their analytical solutions in the Laplace domain are derived using the boundary and interface conditions, and the solutions are then decomposed into series form according to the power series expansion method. Each item of the series solution has a clear physical meaning and represents a transient wave path. Finally, by applying Cagniard's method and the convolution theorem, the analytical solutions are transformed into the time domain. Numerical examples are provided to illustrate some interesting features in the fluid layer, the interface and the semi-infinite solid. When the P-wave velocity in the fluid is higher than that in the solid, two head waves in the solid, one head wave in the fluid and a Scholte wave at the interface are observed for the cylindrical P-wave line source.

A Study on the Profile Change Measurement of Steam Generator Tubes with Tube Expansion Methods

International Nuclear Information System (INIS)

Kim, Young Kyu; Song Myung Ho; Choi, Myung Sik

2011-01-01

Steam generator tubes for nuclear power plants contain the local shape transitions on their inner or outer surface such as dent, bulge, over-expansion, eccentricity, deflection, and so on by the application of physical force during the tube manufacturing and steam generator assembling and by the sludge (that is, corrosion products) produced during the plant operation. The structural integrity of tubes will be degraded by generating the corrosive crack at that location. The profilometry using the traditional bobbin probes which are currently applied for measuring the profile change of tubes gives us basic information such as axial locations and average magnitudes of deformations. However, the three-dimensional quantitative evaluation on circumferential locations, distributional angle, and size of deformations will have to be conducted to understand the effects of residual stresses increased by local deformations on corrosive cracking of tubes. Steam generator tubes of Korean standard nuclear power plants expanded within their tube-sheets by the explosive expansion method and suffered from corrosive cracks in the early stage of power operation. Thus, local deformations of steam generator tubes at the top of tube-sheet were measured with an advanced rotating probe and a laser profiling system for the two cases where the tubes expanded by the explosive expansion method and hydraulic expansion. Also, the trends of eccentricity, deflection, and over-expansion of tubes were evaluated. The advanced eddy current profilometry was confirmed to provide accurate information of local deformations compared with laser profilometry

On the almost everywhere convergence of the eigenfunction expansions from Liouville classes L_1^\\alpha ({T^N})

Science.gov (United States)

Ahmedov, Anvarjon; Materneh, Ehab; Zainuddin, Hishamuddin

2017-09-01

The relevance of waves in quantum mechanics naturally implies that the decomposition of arbitrary wave packets in terms of monochromatic waves plays an important role in applications of the theory. When eigenfunction expansions does not converge, then the expansions of the functions with certain smoothness should be considered. Such functions gained prominence primarily through their application in quantum mechanics. In this work we study the almost everywhere convergence of the eigenfunction expansions from Liouville classes L_p^α ({T^N}), related to the self-adjoint extension of the Laplace operator in torus TN . The sufficient conditions for summability is obtained using the modified Poisson formula. Isomorphism properties of the elliptic differential operators is applied in order to obtain estimation for the Fourier series of the functions from the classes of Liouville L_p^α .

Microstrip natural wave spectrum mathematical model using partial inversion method

International Nuclear Information System (INIS)

Pogarsky, S.A.; Litvinenko, L.N.; Prosvirnin, S.L.

1995-01-01

It is generally agreed that both microstrip lines itself and different discontinuities based on microstrips are the most difficult problem for accurate electrodynamic analysis. Over the last years much has been published about principles and accurate (or full wave) methods of microstrip lines investigations. The growing interest for this problem may be explained by the microstrip application in the millimeter-wave range for purpose of realizing interconnects and a variety of passive components. At these higher operating rating frequencies accurate component modeling becomes more critical. A creation, examination and experimental verification of the accurate method for planar electrodynamical structures natural wave spectrum investigations are the objects of this manuscript. The moment method with partial inversion operator method using may be considered as a basical way for solving this problem. This method is outlook for accurate analysis of different planar discontinuities in microstrip: such as step discontinuities, microstrip turns, Y- and X-junctions and etc., substrate space steps dielectric constants and other anisotropy types

The relationship between the Johnson-Baranger time-dependent folded diagram expansion and the time-independent methods of perturbation theory

International Nuclear Information System (INIS)

Passos, E.M.J. de

1976-01-01

The relationship between the Johnson-Baranger time-dependent folded diagram (JBFD) expansion, and the time independent methods of perturbation theory, are investigated. In the nondegenerate case, the JBFD expansion and the Rayleigh-Schroedinger perturbation expansion, for the ground state energy, are identical. On the other hand, in the degenerate case, for the nonhermitian effective interaction considered, the JBFD expansion, of the effective interaction, is equal to the perturbative expansion of the effective interaction of the nonhermitian eigenvalue problem of Bloch and Brandow-Des Cloizeaux. For the two hermitian effective interactions, the JBFD expansion of the effective interaction differs from the perturbation expansion of the effective interaction of the hermitian eigenvalue problem of Des Cloizeaux [pt

Non-intrusive uncertainty quantification in structural-acoustic systems using polynomial chaos expansion method

Directory of Open Access Journals (Sweden)

Wang Mingjie

2017-01-01

Full Text Available A framework of non-intrusive polynomial chaos expansion method (PC was proposed to investigate the statistic characteristics of the response of structural-acoustic system containing random uncertainty. The PC method does not need to reformulate model equations, and the statistics of the response can be evaluated directly. The results show that compared to the direct Monte Carlo method (MCM based on the original numerical model, the PC method is effective and more efficient.

Wave resistance calculation method combining Green functions based on Rankine and Kelvin source

Directory of Open Access Journals (Sweden)

LI Jingyu

2017-12-01

Full Text Available [Ojectives] At present, the Boundary Element Method(BEM of wave-making resistance mostly uses a model in which the velocity distribution near the hull is solved first, and the pressure integral is then calculated using the Bernoulli equation. However,the process of this model of wave-making resistance is complex and has low accuracy.[Methods] To address this problem, the present paper deduces a compound method for the quick calculation of ship wave resistance using the Rankine source Green function to solve the hull surface's source density, and combining the Lagally theorem concerning source point force calculation based on the Kelvin source Green function so as to solve the wave resistance. A case for the Wigley model is given.[Results] The results show that in contrast to the thin ship method of the linear wave resistance theorem, this method has higher precision, and in contrast to the method which completely uses the Kelvin source Green function, this method has better computational efficiency.[Conclusions] In general, the algorithm in this paper provides a compromise between precision and efficiency in wave-making resistance calculation.

Application of the generalized multi structural (GMS) wave function to photoelectron spectra and electron scattering processes

International Nuclear Information System (INIS)

Nascimento, M.A.C. do

1992-01-01

A Generalized Multi Structural (GMS) wave function is presented which combines the advantages of the SCF-MO and VB models, preserving the classical chemical structures but optimizing the orbitals in a self-consistent way. This wave function is particularly suitable to treat situations where the description of the molecular state requires localized wave functions. It also provides a very convenient way of treating the electron correlation problem, avoiding large CI expansions. The final wave functions are much more compact and easier to interpret than the ones obtained by the conventional methods, using orthogonal orbitals. Applications of the GMS wave function to the study of the photoelectron spectra of the trans-glyoxal molecule and to electron impact excitation processes in the nitrogen molecule are presented as an illustration of the method. (author)

Electromagnetic waves in irregular multilayered spheroidal structures of finite conductivity: full wave solutions

International Nuclear Information System (INIS)

Bahar, E.

1976-01-01

The propagation of electromagnetic waves excited by electric dipoles oriented along the axis of multilayered spheroidal structures of finite conductivity is investigated. The electromagnetic parameters and the thickness of the layers of the structure are assumed to be functions of the latitude. In the analysis, electric and magnetic field transforms that constitute a discrete and a continuous spectrum of spherical waves are used to provide a suitable basis for the expansion of the electromagnetic fields at any point in the irregular spheroidal structure. For spheroidal structures with good conducting cores, the terms in the solutions associated with the continuous part of the wave spectrum vanish. In general, however, when the skin depth for the core is large compared to its dimensions or when the sources are located in the core of the structure and propagation in the core is of special interest, the contribution from the continuous part of the wave spectrum cannot be neglected. At each interface between the layers of the irregular spheroidal structure, exact boundary conditions are imposed. Since the terms of the field expansions in the irregular structure do not individually satisfy the boundary conditions, Maxwell's equations are reduced to sets of coupled ordinary first-order differential equations for the wave amplitudes. The solutions are shown to satisfy the reciprocity relationships in electromagnetic theory. The analysis may be applied to problems of radio wave propagation in a nonuniform model of the earth-ionosphere waveguide, particularly when focusing effects at the antipodes are important

Stress wave focusing transducers

Energy Technology Data Exchange (ETDEWEB)

Visuri, S.R., LLNL

1998-05-15

Conversion of laser radiation to mechanical energy is the fundamental process behind many medical laser procedures, particularly those involving tissue destruction and removal. Stress waves can be generated with laser radiation in several ways: creation of a plasma and subsequent launch of a shock wave, thermoelastic expansion of the target tissue, vapor bubble collapse, and ablation recoil. Thermoelastic generation of stress waves generally requires short laser pulse durations and high energy density. Thermoelastic stress waves can be formed when the laser pulse duration is shorter than the acoustic transit time of the material: {tau}{sub c} = d/c{sub s} where d = absorption depth or spot diameter, whichever is smaller, and c{sub s} = sound speed in the material. The stress wave due to thermoelastic expansion travels at the sound speed (approximately 1500 m/s in tissue) and leaves the site of irradiation well before subsequent thermal events can be initiated. These stress waves, often evolving into shock waves, can be used to disrupt tissue. Shock waves are used in ophthalmology to perform intraocular microsurgery and photodisruptive procedures as well as in lithotripsy to fragment stones. We have explored a variety of transducers that can efficiently convert optical to mechanical energy. One such class of transducers allows a shock wave to be focused within a material such that the stress magnitude can be greatly increased compared to conventional geometries. Some transducer tips could be made to operate regardless of the absorption properties of the ambient media. The size and nature of the devices enable easy delivery, potentially minimally-invasive procedures, and precise tissue- targeting while limiting thermal loading. The transducer tips may have applications in lithotripsy, ophthalmology, drug delivery, and cardiology.

Microtubules Nonlinear Models Dynamics Investigations through the exp(−Φ(ξ-Expansion Method Implementation

Directory of Open Access Journals (Sweden)

Nur Alam

2016-02-01

Full Text Available In this research article, we present exact solutions with parameters for two nonlinear model partial differential equations(PDEs describing microtubules, by implementing the exp(−Φ(ξ-Expansion Method. The considered models, describing highly nonlinear dynamics of microtubules, can be reduced to nonlinear ordinary differential equations. While the first PDE describes the longitudinal model of nonlinear dynamics of microtubules, the second one describes the nonlinear model of dynamics of radial dislocations in microtubules. The acquired solutions are then graphically presented, and their distinct properties are enumerated in respect to the corresponding dynamic behavior of the microtubules they model. Various patterns, including but not limited to regular, singular kink-like, as well as periodicity exhibiting ones, are detected. Being the method of choice herein, the exp(−Φ(ξ-Expansion Method not disappointing in the least, is found and declared highly efficient.

Vibration isolation in a free-piston driven expansion tube facility

Science.gov (United States)

Gildfind, D. E.; Jacobs, P. A.; Morgan, R. G.

2013-09-01

The stress waves produced by rapid piston deceleration are a fundamental feature of free-piston driven expansion tubes, and wave propagation has to be considered in the design process. For lower enthalpy test conditions, these waves can traverse the tube ahead of critical flow processes, severely interfering with static pressure measurements of the passing flow. This paper details a new device which decouples the driven tube from the free-piston driver, and thus prevents transmission of stress waves. Following successful incorporation of the concept in the smaller X2 facility, it has now been applied to the larger X3 facility, and results for both facilities are presented.

Exact solutions of nonlinear fractional differential equations by (G′/G)-expansion method

International Nuclear Information System (INIS)

Bekir Ahmet; Güner Özkan

2013-01-01

In this paper, we use the fractional complex transform and the (G′/G)-expansion method to study the nonlinear fractional differential equations and find the exact solutions. The fractional complex transform is proposed to convert a partial fractional differential equation with Jumarie's modified Riemann—Liouville derivative into its ordinary differential equation. It is shown that the considered transform and method are very efficient and powerful in solving wide classes of nonlinear fractional order equations

Nonlinear self-modulation of ion-acoustic waves

International Nuclear Information System (INIS)

Ikezi, H.; Schwarzenegger, K.; Simons, A.L.; Ohsawa, Y.; Kamimura, T.

1978-01-01

The nonlinear evolution of an ion-acoustic wave packet is studied. Experimentally, it is found that (i) nonlinear phase modulation develops in the wave packet; (ii) the phase modulation, together with the dispersion effect, causes expansion and breaking of the wave packet; (iii) the ions trapped in the troughs of the wave potential introduce self-phase modulation; and (iv) the ion-acoustic wave is stable with respect to the modulational instability. Computer simulations have reproduced the experimental results. The physical picture and the model equation describing the wave evolution are discussed

Wave function of the Universe in the early stage of its evolution

International Nuclear Information System (INIS)

Maydanyuk, Sergei P.

2008-01-01

In quantum cosmological models, constructed in the framework of Friedmann-Robertson-Walker metrics, a nucleation of the Universe with its further expansion is described as a tunneling transition through an effective barrier between regions with small and large values of the scale factor a at non-zero (or zero) energy. The approach for describing this tunneling consists of constructing a wave function satisfying an appropriate boundary condition. There are various ways for defining the boundary condition that lead to different estimates of the barrier penetrability and the tunneling time. In order to describe the escape from the tunneling region as accurately as possible and to construct the total wave function on the basis of its two partial solutions unambiguously, we use the tunneling boundary condition that the total wave function must represent only the outgoing wave at the point of escape from the barrier, where the following definition for the wave is introduced: the wave is represented by the wave function whose modulus changes minimally under a variation of the scale factor a. We construct a new method for a direct non-semiclassical calculation of the total stationary wave function of the Universe, analyze the behavior of this wave function in the tunneling region, near the escape point and in the asymptotic region, and estimate the barrier penetrability. We observe oscillations of the modulus of the wave function in the external region starting from the turning point which decrease with increasing of a and which are not shown in semiclassical calculations. The period of such an oscillation decreases uniformly with increasing a and can be used as a fully quantum dynamical characteristic of the expansion of the Universe. (orig.)

Computational and theoretical study of the wave-particle interaction of protons and waves

Directory of Open Access Journals (Sweden)

P. S. Moya

2012-09-01

Full Text Available We study the wave-particle interaction and the evolution of electromagnetic waves propagating through a plasma composed of electrons and protons, using two approaches. First, a quasilinear kinetic theory has been developed to study the energy transfer between waves and particles, with the subsequent acceleration and heating of protons. Second, a one-dimensional hybrid numerical simulation has been performed, with and without including an expanding-box model that emulates the spherical expansion of the solar wind, to investigate the fully nonlinear evolution of this wave-particle interaction. Numerical results of both approaches show that there is an anisotropic evolution of proton temperature.

The finite product method in the theory of linear wave propagation

DEFF Research Database (Denmark)

Sorokin, Sergey; Chapman, John

2012-01-01

of the method are presented for several non-trivial examples, that of symmetric/anti-symmetric elastic waves in a layer and in a thin plate. In each case, the method gives a sequence of polynomial approximations to the dispersion relation of remarkable accuracy over a broad range of frequencies and wave numbers...

Finite element and discontinuous Galerkin methods for transient wave equations

CERN Document Server

Cohen, Gary

2017-01-01

This monograph presents numerical methods for solving transient wave equations (i.e. in time domain). More precisely, it provides an overview of continuous and discontinuous finite element methods for these equations, including their implementation in physical models, an extensive description of 2D and 3D elements with different shapes, such as prisms or pyramids, an analysis of the accuracy of the methods and the study of the Maxwell’s system and the important problem of its spurious free approximations. After recalling the classical models, i.e. acoustics, linear elastodynamics and electromagnetism and their variational formulations, the authors present a wide variety of finite elements of different shapes useful for the numerical resolution of wave equations. Then, they focus on the construction of efficient continuous and discontinuous Galerkin methods and study their accuracy by plane wave techniques and a priori error estimates. A chapter is devoted to the Maxwell’s system and the important problem ...

Spectral analysis of surface waves method to assess shear wave velocity within centrifuge models

OpenAIRE

MURILLO, Carol Andrea; THOREL, Luc; CAICEDO, Bernardo

2009-01-01

The method of the spectral analysis of surface waves (SASW) is tested out on reduced scale centrifuge models, with a specific device, called the mini Falling Weight, developed for this purpose. Tests are performed on layered materials made of a mixture of sand and clay. The shear wave velocity VS determined within the models using the SASW is compared with the laboratory measurements carried out using the bender element test. The results show that the SASW technique applied to centrifuge test...

Nonlinear wave forces on large ocean structures

Science.gov (United States)

Huang, Erick T.

1993-04-01

This study explores the significance of second-order wave excitations on a large pontoon and tests the feasibility of reducing a nonlinear free surface problem by perturbation expansions. A simulation model has been developed based on the perturbation expansion technique to estimate the wave forces. The model uses a versatile finite element procedure for the solution of the reduced linear boundary value problems. This procedure achieves a fair compromise between computation costs and physical details by using a combination of 2D and 3D elements. A simple hydraulic model test was conducted to observe the wave forces imposed on a rectangle box by Cnoidal waves in shallow water. The test measurements are consistent with the numerical predictions by the simulation model. This result shows favorable support to the perturbation approach for estimating the nonlinear wave forces on shallow draft vessels. However, more sophisticated model tests are required for a full justification. Both theoretical and experimental results show profound second-order forces that could substantially impact the design of ocean facilities.

The Taylor-expansion method of moments for the particle system with bimodal distribution

Directory of Open Access Journals (Sweden)

Liu Yan-Hua

2013-01-01

Full Text Available This paper derives the multipoint Taylor expansion method of moments for the bimodal particle system. The collision effects are modeled by the internal and external coagulation terms. Simple theory and numerical tests are performed to prove the effect of the current model.

A tuning method for nonuniform traveling-wave accelerating structures

International Nuclear Information System (INIS)

Gong Cunkui; Zheng Shuxin; Shao Jiahang; Jia Xiaoyu; Chen Huaibi

2013-01-01

The tuning method of uniform traveling-wave structures based on non-resonant perturbation field distribution measurement has been widely used in tuning both constant-impedance and constant-gradient structures. In this paper, the method of tuning nonuniform structures is proposed on the basis of the above theory. The internal reflection coefficient of each cell is obtained from analyzing the normalized voltage distribution. A numerical simulation of tuning process according to the coupled cavity chain theory has been done and the result shows each cell is in right phase advance after tuning. The method will be used in the tuning of a disk-loaded traveling-wave structure being developed at the Accelerator Laboratory, Tsinghua University. (authors)

A modified symplectic PRK scheme for seismic wave modeling

Science.gov (United States)

Liu, Shaolin; Yang, Dinghui; Ma, Jian

2017-02-01

A new scheme for the temporal discretization of the seismic wave equation is constructed based on symplectic geometric theory and a modified strategy. The ordinary differential equation in terms of time, which is obtained after spatial discretization via the spectral-element method, is transformed into a Hamiltonian system. A symplectic partitioned Runge-Kutta (PRK) scheme is used to solve the Hamiltonian system. A term related to the multiplication of the spatial discretization operator with the seismic wave velocity vector is added into the symplectic PRK scheme to create a modified symplectic PRK scheme. The symplectic coefficients of the new scheme are determined via Taylor series expansion. The positive coefficients of the scheme indicate that its long-term computational capability is more powerful than that of conventional symplectic schemes. An exhaustive theoretical analysis reveals that the new scheme is highly stable and has low numerical dispersion. The results of three numerical experiments demonstrate the high efficiency of this method for seismic wave modeling.

Development of Extended Ray-tracing method including diffraction, polarization and wave decay effects

Science.gov (United States)

Yanagihara, Kota; Kubo, Shin; Dodin, Ilya; Nakamura, Hiroaki; Tsujimura, Toru

2017-10-01

Geometrical Optics Ray-tracing is a reasonable numerical analytic approach for describing the Electron Cyclotron resonance Wave (ECW) in slowly varying spatially inhomogeneous plasma. It is well known that the result with this conventional method is adequate in most cases. However, in the case of Helical fusion plasma which has complicated magnetic structure, strong magnetic shear with a large scale length of density can cause a mode coupling of waves outside the last closed flux surface, and complicated absorption structure requires a strong focused wave for ECH. Since conventional Ray Equations to describe ECW do not have any terms to describe the diffraction, polarization and wave decay effects, we can not describe accurately a mode coupling of waves, strong focus waves, behavior of waves in inhomogeneous absorption region and so on. For fundamental solution of these problems, we consider the extension of the Ray-tracing method. Specific process is planned as follows. First, calculate the reference ray by conventional method, and define the local ray-base coordinate system along the reference ray. Then, calculate the evolution of the distributions of amplitude and phase on ray-base coordinate step by step. The progress of our extended method will be presented.

The Nodal Polynomial Expansion method to solve the multigroup diffusion equations

International Nuclear Information System (INIS)

Ribeiro, R.D.M.

1983-03-01

The methodology of the solutions of the multigroup diffusion equations and uses the Nodal Polynomial Expansion Method is covered. The EPON code was developed based upon the above mentioned method for stationary state, rectangular geometry, one-dimensional or two-dimensional and for one or two energy groups. Then, one can study some effects such as the influence of the baffle on the thermal flux by calculating the flux and power distribution in nuclear reactors. Furthermore, a comparative study with other programs which use Finite Difference (CITATION and PDQ5) and Finite Element (CHD and FEMB) Methods was undertaken. As a result, the coherence, feasibility, speed and accuracy of the methodology used were demonstrated. (Author) [pt

Analysis of a plane stress wave by the moving least squares method

Directory of Open Access Journals (Sweden)

Wojciech Dornowski

2014-08-01

Full Text Available A meshless method based on the moving least squares approximation is applied to stress wave propagation analysis. Two kinds of node meshes, the randomly generated mesh and the regular mesh are used. The nearest neighbours’ problem is developed from a triangulation that satisfies minimum edges length conditions. It is found that this method of neighbours’ choice significantly improves the solution accuracy. The reflection of stress waves from the free edge is modelled using fictitious nodes (outside the plate. The comparison with the finite difference results also demonstrated the accuracy of the proposed approach.[b]Keywords[/b]: civil engineering, meshless method, moving least squares method, elastic waves

Optical wave microphone measurement during laser ablation of Si

Energy Technology Data Exchange (ETDEWEB)

Mitsugi, Fumiaki, E-mail: mitsugi@cs.kumamoto-u.ac.jp [Graduate School of Science and Technology, Kumamoto University, 2-39-1 Kurokami, Kumamoto, 860-8555 (Japan); Ide, Ryota; Ikegami, Tomoaki [Graduate School of Science and Technology, Kumamoto University, 2-39-1 Kurokami, Kumamoto, 860-8555 (Japan); Nakamiya, Toshiyuki; Sonoda, Yoshito [Graduate School of Industrial Engineering, Tokai University, 9-1-1 Toroku, Kumamoto, 862-8652 (Japan)

2012-10-30

Pulsed laser irradiation is used for surface treatment of a solid and ablation for particle formation in gas, liquid or supercritical phase media. When a pulsed laser is used to irradiate a solid, spatial refractive index variations (including photothermal expansion, shockwaves and particles) occur, which vary depending on the energy density of the pulsed laser. We focused on this phenomenon and applied an unique method for detection of refractive index variation using an optical wave microphone based on Fraunhofer diffraction. In this research, we analyzed the waveforms and frequencies of refractive index variations caused by pulsed laser irradiation of silicon in air and measured with an optical wave microphone.

Systems and methods for wave energy conversion

Science.gov (United States)

MacDonald, Daniel G.; Cantara, Justin; Nathan, Craig; Lopes, Amy M.; Green, Brandon E.

2017-02-28

Systems for wave energy conversion that have components that can survive the harsh marine environment and that can be attached to fixed structures, such as a pier, and having the ability to naturally adjust for tidal height and methods for their use are presented.

Conservation Laws and Traveling Wave Solutions of a Generalized Nonlinear ZK-BBM Equation

Directory of Open Access Journals (Sweden)

Khadijo Rashid Adem

2014-01-01

Full Text Available We study a generalized two-dimensional nonlinear Zakharov-Kuznetsov-Benjamin-Bona-Mahony (ZK-BBM equation, which is in fact Benjamin-Bona-Mahony equation formulated in the ZK sense. Conservation laws for this equation are constructed by using the new conservation theorem due to Ibragimov and the multiplier method. Furthermore, traveling wave solutions are obtained by employing the (G'/G-expansion method.

Generalization of the linear algebraic method to three dimensions

International Nuclear Information System (INIS)

Lynch, D.L.; Schneider, B.I.

1991-01-01

We present a numerical method for the solution of the Lippmann-Schwinger equation for electron-molecule collisions. By performing a three-dimensional numerical quadrature, this approach avoids both a basis-set representation of the wave function and a partial-wave expansion of the scattering potential. The resulting linear equations, analogous in form to the one-dimensional linear algebraic method, are solved with the direct iteration-variation method. Several numerical examples are presented. The prospect for using this numerical quadrature scheme for electron-polyatomic molecules is discussed

A fast method for linear waves based on geometrical optics

NARCIS (Netherlands)

Stolk, C.C.

2009-01-01

We develop a fast method for solving the one-dimensional wave equation based on geometrical optics. From geometrical optics (e.g., Fourier integral operator theory or WKB approximation) it is known that high-frequency waves split into forward and backward propagating parts, each propagating with the

Angular correlation methods

International Nuclear Information System (INIS)

Ferguson, A.J.

1974-01-01

An outline of the theory of angular correlations is presented, and the difference between the modern density matrix method and the traditional wave function method is stressed. Comments are offered on particular angular correlation theoretical techniques. A brief discussion is given of recent studies of gamma ray angular correlations of reaction products recoiling with high velocity into vacuum. Two methods for optimization to obtain the most accurate expansion coefficients of the correlation are discussed. (1 figure, 53 references) (U.S.)

The resonance expansion for the Green's function of the Schroedinger and wave equations

International Nuclear Information System (INIS)

Albeverio, S.; Aix-Marseille-2 Univ., 13 - Marseille; Hoeegh-Krohn, R.; Oslo Univ.

1984-01-01

We give a survey of some recent mathematical work on resonances, in particular on perturbation series, low energy expansions and on resonances for point interactions. Expansions of the kernels of esup(-it)√sup(H+) and esup(-itH) in terms of resonances are also given (where Hsub(+) is the positive part of the Hamiltonian). (orig.)

Nolinear stability analysis of nuclear reactors : expansion methods for stability domains

International Nuclear Information System (INIS)

Yang, Chae Yong

1992-02-01

Two constructive methods for estimating asymptotic stability domains of nonlinear reactor models are developed in this study: an improved Chang and Thorp's method based on expansion of a Lyapunov function and a new method based on expansion of any positive definite function. The methods are established on the concept of stability definitions of Lyapunov itself. The first method provides a sequence of stability regions that eventually approaches the exact stability domain, but requires many expansions in order to obtain the entire stability region because the starting Lyapunov function usually corresponds to a small stability region and because most dynamic systems are stiff. The second method (new method) requires only a positive definite function and thus it is easy to come up with a starting region. From a large starting region, the entire stability region is estimated effectively after sufficient iterations. It is particularly useful for stiff systems. The methods are applied to several nonlinear reactor models known in the literature: one-temperature feedback model, two-temperature feedback model, and xenon dynamics model, and the results are compared. A reactor feedback model for a pressurized water reactor (PWR) considering fuel and moderator temperature effects is developed and the nonlinear stability regions are estimated for the various values of design parameters by using the new method. The steady-state properties of the nonlinear reactor system are analyzed via bifurcation theory. The analysis of nonlinear phenomena is carried out for the various forms of reactivity feedback coefficients that are both temperature- (or power-) independent and dependent. If one of two temperature coefficients is positive, unstable limit cycles or multiplicity of the steady-state solutions appear when the other temperature coefficient exceeds a certain critical value. As an example, even though the fuel temperature coefficient is negative, if the moderator temperature

Traveling Wave Resonance and Simplified Analysis Method for Long-Span Symmetrical Cable-Stayed Bridges under Seismic Traveling Wave Excitation

Directory of Open Access Journals (Sweden)

Zhong-ye Tian

2014-01-01

Full Text Available The seismic responses of a long-span cable-stayed bridge under uniform excitation and traveling wave excitation in the longitudinal direction are, respectively, computed. The numerical results show that the bridge’s peak seismic responses vary significantly as the apparent wave velocity decreases. Therefore, the traveling wave effect must be considered in the seismic design of long-span bridges. The bridge’s peak seismic responses do not vary monotonously with the apparent wave velocity due to the traveling wave resonance. A new traveling wave excitation method that can simplify the multisupport excitation process into a two-support excitation process is developed.

Application of the Asymptotic Taylor Expansion Method to Bistable Potentials

Directory of Open Access Journals (Sweden)

Okan Ozer

2013-01-01

Full Text Available A recent method called asymptotic Taylor expansion (ATEM is applied to determine the analytical expression for eigenfunctions and numerical results for eigenvalues of the Schrödinger equation for the bistable potentials. Optimal truncation of the Taylor series gives a best possible analytical expression for eigenfunctions and numerical results for eigenvalues. It is shown that the results are obtained by a simple algorithm constructed for a computer system using symbolic or numerical calculation. It is observed that ATEM produces excellent results consistent with the existing literature.

Nonlocal symmetry and explicit solutions from the CRE method of the Boussinesq equation

Science.gov (United States)

Zhao, Zhonglong; Han, Bo

2018-04-01

In this paper, we analyze the integrability of the Boussinesq equation by using the truncated Painlevé expansion and the CRE method. Based on the truncated Painlevé expansion, the nonlocal symmetry and Bäcklund transformation of this equation are obtained. A prolonged system is introduced to localize the nonlocal symmetry to the local Lie point symmetry. It is proved that the Boussinesq equation is CRE solvable. The two-solitary-wave fusion solutions, single soliton solutions and soliton-cnoidal wave solutions are presented by means of the Bäcklund transformations.

A comparison of high-order polynomial and wave-based methods for Helmholtz problems

Science.gov (United States)

Lieu, Alice; Gabard, Gwénaël; Bériot, Hadrien

2016-09-01

The application of computational modelling to wave propagation problems is hindered by the dispersion error introduced by the discretisation. Two common strategies to address this issue are to use high-order polynomial shape functions (e.g. hp-FEM), or to use physics-based, or Trefftz, methods where the shape functions are local solutions of the problem (typically plane waves). Both strategies have been actively developed over the past decades and both have demonstrated their benefits compared to conventional finite-element methods, but they have yet to be compared. In this paper a high-order polynomial method (p-FEM with Lobatto polynomials) and the wave-based discontinuous Galerkin method are compared for two-dimensional Helmholtz problems. A number of different benchmark problems are used to perform a detailed and systematic assessment of the relative merits of these two methods in terms of interpolation properties, performance and conditioning. It is generally assumed that a wave-based method naturally provides better accuracy compared to polynomial methods since the plane waves or Bessel functions used in these methods are exact solutions of the Helmholtz equation. Results indicate that this expectation does not necessarily translate into a clear benefit, and that the differences in performance, accuracy and conditioning are more nuanced than generally assumed. The high-order polynomial method can in fact deliver comparable, and in some cases superior, performance compared to the wave-based DGM. In addition to benchmarking the intrinsic computational performance of these methods, a number of practical issues associated with realistic applications are also discussed.

Symmetry in the polarization expansion for intermolecular forces

International Nuclear Information System (INIS)

Chipman, D.M.; Hirschfelder, J.O.

1980-01-01

In the usual polarization expansion for intermolecular forces, exchange effects that determine the separations of energy levels within the manifold of interacting states are ignored. Previous low order calculations on simple physical systems have indicated that these exchange terms can be described reasonably well by an appropriate ad hoc symmetrization of the polarization wave function (the SYM-P method). But theoretical considerations suggest that the SYM-P method should be good for only one of the interacting states and not for the others in the manifold. Here this long standing apparent conflict between theoretical expectations and actual results is explained by consideration of a simple model system in which the relevant equations can be solved exactly. It is concluded that while the SYM-P method is potentially exact for only one of the interacting states, it may provide good approximations to the other states of the manifold in the case of large separations of the interacting subsystems

Modal Ring Method for the Scattering of Electromagnetic Waves

Science.gov (United States)

Baumeister, Kenneth J.; Kreider, Kevin L.

1993-01-01

The modal ring method for electromagnetic scattering from perfectly electric conducting (PEC) symmetrical bodies is presented. The scattering body is represented by a line of finite elements (triangular) on its outer surface. The infinite computational region surrounding the body is represented analytically by an eigenfunction expansion. The modal ring method effectively reduces the two dimensional scattering problem to a one-dimensional problem similar to the method of moments. The modal element method is capable of handling very high frequency scattering because it has a highly banded solution matrix.

Damage evaluation by a guided wave-hidden Markov model based method

Science.gov (United States)

Mei, Hanfei; Yuan, Shenfang; Qiu, Lei; Zhang, Jinjin

2016-02-01

Guided wave based structural health monitoring has shown great potential in aerospace applications. However, one of the key challenges of practical engineering applications is the accurate interpretation of the guided wave signals under time-varying environmental and operational conditions. This paper presents a guided wave-hidden Markov model based method to improve the damage evaluation reliability of real aircraft structures under time-varying conditions. In the proposed approach, an HMM based unweighted moving average trend estimation method, which can capture the trend of damage propagation from the posterior probability obtained by HMM modeling is used to achieve a probabilistic evaluation of the structural damage. To validate the developed method, experiments are performed on a hole-edge crack specimen under fatigue loading condition and a real aircraft wing spar under changing structural boundary conditions. Experimental results show the advantage of the proposed method.

Science with the space-based interferometer eLISA. III: probing the expansion of the universe using gravitational wave standard sirens

Energy Technology Data Exchange (ETDEWEB)

Tamanini, Nicola; Caprini, Chiara [Institut de Physique Théorique, CEA-Saclay, CNRS UMR 3681, Université Paris-Saclay, F-91191 Gif-sur-Yvette (France); Barausse, Enrico [Sorbonne Universités, UPMC Université Paris 6, UMR 7095, Institut d' Astrophysique de Paris, 98 bis Bd Arago, 75014 Paris (France); Sesana, Alberto [School of Physics and Astronomy, The University of Birmingham, Edgbaston, Birmingham B15 2TT (United Kingdom); Klein, Antoine [Department of Physics and Astronomy, The University of Mississippi, University, MS 38677 (United States); Petiteau, Antoine, E-mail: nicola.tamanini@cea.fr, E-mail: chiara.caprini@cea.fr, E-mail: barausse@iap.fr, E-mail: asesana@star.sr.bham.ac.uk, E-mail: aklein@physics.montana.edu, E-mail: antoine.petiteau@apc.univ-paris7.fr [APC, Université Paris Diderot, Observatoire de Paris, Sorbonne Paris Cité, 10 rue Alice Domon et Léonie Duquet, 75205 Paris Cedex 13 (France)

2016-04-01

We investigate the capability of various configurations of the space interferometer eLISA to probe the late-time background expansion of the universe using gravitational wave standard sirens. We simulate catalogues of standard sirens composed by massive black hole binaries whose gravitational radiation is detectable by eLISA, and which are likely to produce an electromagnetic counterpart observable by future surveys. The main issue for the identification of a counterpart resides in the capability of obtaining an accurate enough sky localisation with eLISA. This seriously challenges the capability of four-link (2 arm) configurations to successfully constrain the cosmological parameters. Conversely, six-link (3 arm) configurations have the potential to provide a test of the expansion of the universe up to z ∼ 8 which is complementary to other cosmological probes based on electromagnetic observations only. In particular, in the most favourable scenarios, they can provide a significant constraint on H{sub 0} at the level of 0.5%. Furthermore, (Ω{sub M}, Ω{sub Λ}) can be constrained to a level competitive with present SNIa results. On the other hand, the lack of massive black hole binary standard sirens at low redshift allows to constrain dark energy only at the level of few percent.

An adaptive Bayesian inversion for upper mantle structure using surface waves and scattered body waves

Science.gov (United States)

Eilon, Zachary; Fischer, Karen M.; Dalton, Colleen A.

2018-04-01

We present a methodology for 1-D imaging of upper mantle structure using a Bayesian approach that incorporates a novel combination of seismic data types and an adaptive parameterisation based on piecewise discontinuous splines. Our inversion algorithm lays the groundwork for improved seismic velocity models of the lithosphere and asthenosphere by harnessing the recent expansion of large seismic arrays and computational power alongside sophisticated data analysis. Careful processing of P- and S-wave arrivals isolates converted phases generated at velocity gradients between the mid-crust and 300 km depth. This data is allied with ambient noise and earthquake Rayleigh wave phase velocities to obtain detailed VS and VP velocity models. Synthetic tests demonstrate that converted phases are necessary to accurately constrain velocity gradients, and S-p phases are particularly important for resolving mantle structure, while surface waves are necessary for capturing absolute velocities. We apply the method to several stations in the northwest and north-central United States, finding that the imaged structure improves upon existing models by sharpening the vertical resolution of absolute velocity profiles, offering robust uncertainty estimates, and revealing mid-lithospheric velocity gradients indicative of thermochemical cratonic layering. This flexible method holds promise for increasingly detailed understanding of the upper mantle.

The extended hyperbolic function method and exact solutions of the long-short wave resonance equations

International Nuclear Information System (INIS)

Shang Yadong

2008-01-01

The extended hyperbolic functions method for nonlinear wave equations is presented. Based on this method, we obtain a multiple exact explicit solutions for the nonlinear evolution equations which describe the resonance interaction between the long wave and the short wave. The solutions obtained in this paper include (a) the solitary wave solutions of bell-type for S and L, (b) the solitary wave solutions of kink-type for S and bell-type for L, (c) the solitary wave solutions of a compound of the bell-type and the kink-type for S and L, (d) the singular travelling wave solutions, (e) periodic travelling wave solutions of triangle function types, and solitary wave solutions of rational function types. The variety of structure to the exact solutions of the long-short wave equation is illustrated. The methods presented here can also be used to obtain exact solutions of nonlinear wave equations in n dimensions

In situ measurement of plasma and shock wave properties inside laser-drilled metal holes

Science.gov (United States)

Brajdic, Mihael; Hermans, Martin; Horn, Alexander; Kelbassa, Ingomar

2008-10-01

High-speed imaging of shock wave and plasma dynamics is a commonly used diagnostic method for monitoring processes during laser material treatment. It is used for processes such as laser ablation, cutting, keyhole welding and drilling. Diagnosis of laser drilling is typically adopted above the material surface because lateral process monitoring with optical diagnostic methods inside the laser-drilled hole is not possible due to the hole walls. A novel method is presented to investigate plasma and shock wave properties during the laser drilling inside a confined environment such as a laser-drilled hole. With a novel sample preparation and the use of high-speed imaging combined with spectroscopy, a time and spatial resolved monitoring of plasma and shock wave dynamics is realized. Optical emission of plasma and shock waves during drilling of stainless steel with ns-pulsed laser radiation is monitored and analysed. Spatial distributions and velocities of shock waves and of plasma are determined inside the holes. Spectroscopy is accomplished during the expansion of the plasma inside the drilled hole allowing for the determination of electron densities.

Full wave simulation of waves in ECRIS plasmas based on the finite element method

Energy Technology Data Exchange (ETDEWEB)

Torrisi, G. [INFN - Laboratori Nazionali del Sud, via S. Sofia 62, 95123, Catania, Italy and Università Mediterranea di Reggio Calabria, Dipartimento di Ingegneria dell' Informazione, delle Infrastrutture e dell' Energia Sostenibile (DIIES), Via Graziella, I (Italy); Mascali, D.; Neri, L.; Castro, G.; Patti, G.; Celona, L.; Gammino, S.; Ciavola, G. [INFN - Laboratori Nazionali del Sud, via S. Sofia 62, 95123, Catania (Italy); Di Donato, L. [Università degli Studi di Catania, Dipartimento di Ingegneria Elettrica Elettronica ed Informatica (DIEEI), Viale Andrea Doria 6, 95125 Catania (Italy); Sorbello, G. [INFN - Laboratori Nazionali del Sud, via S. Sofia 62, 95123, Catania, Italy and Università degli Studi di Catania, Dipartimento di Ingegneria Elettrica Elettronica ed Informatica (DIEEI), Viale Andrea Doria 6, 95125 Catania (Italy); Isernia, T. [Università Mediterranea di Reggio Calabria, Dipartimento di Ingegneria dell' Informazione, delle Infrastrutture e dell' Energia Sostenibile (DIIES), Via Graziella, I-89100 Reggio Calabria (Italy)

2014-02-12

This paper describes the modeling and the full wave numerical simulation of electromagnetic waves propagation and absorption in an anisotropic magnetized plasma filling the resonant cavity of an electron cyclotron resonance ion source (ECRIS). The model assumes inhomogeneous, dispersive and tensorial constitutive relations. Maxwell's equations are solved by the finite element method (FEM), using the COMSOL Multiphysics{sup ®} suite. All the relevant details have been considered in the model, including the non uniform external magnetostatic field used for plasma confinement, the local electron density profile resulting in the full-3D non uniform magnetized plasma complex dielectric tensor. The more accurate plasma simulations clearly show the importance of cavity effect on wave propagation and the effects of a resonant surface. These studies are the pillars for an improved ECRIS plasma modeling, that is mandatory to optimize the ion source output (beam intensity distribution and charge state, especially). Any new project concerning the advanced ECRIS design will take benefit by an adequate modeling of self-consistent wave absorption simulations.

Application of the Most Likely Extreme Response Method for Wave Energy Converters: Preprint

Energy Technology Data Exchange (ETDEWEB)

Quon, Eliot; Platt, Andrew; Yu, Yi-Hsiang; Lawson, Michael

2016-07-01

Extreme loads are often a key cost driver for wave energy converters (WECs). As an alternative to exhaustive Monte Carlo or long-term simulations, the most likely extreme response (MLER) method allows mid- and high-fidelity simulations to be used more efficiently in evaluating WEC response to events at the edges of the design envelope, and is therefore applicable to system design analysis. The study discussed in this paper applies the MLER method to investigate the maximum heave, pitch, and surge force of a point absorber WEC. Most likely extreme waves were obtained from a set of wave statistics data based on spectral analysis and the response amplitude operators (RAOs) of the floating body; the RAOs were computed from a simple radiation-and-diffraction-theory-based numerical model. A weakly nonlinear numerical method and a computational fluid dynamics (CFD) method were then applied to compute the short-term response to the MLER wave. Effects of nonlinear wave and floating body interaction on the WEC under the anticipated 100-year waves were examined by comparing the results from the linearly superimposed RAOs, the weakly nonlinear model, and CFD simulations. Overall, the MLER method was successfully applied. In particular, when coupled to a high-fidelity CFD analysis, the nonlinear fluid dynamics can be readily captured.

Hydrodynamic analysis and shape optimization for vertical axisymmetric wave energy converters

Science.gov (United States)

Zhang, Wan-chao; Liu, Heng-xu; Zhang, Liang; Zhang, Xue-wei

2016-12-01

The absorber is known to be vertical axisymmetric for a single-point wave energy converter (WEC). The shape of the wetted surface usually has a great influence on the absorber's hydrodynamic characteristics which are closely linked with the wave power conversion ability. For complex wetted surface, the hydrodynamic coefficients have been predicted traditionally by hydrodynamic software based on the BEM. However, for a systematic study of various parameters and geometries, they are too multifarious to generate so many models and data grids. This paper examines a semi-analytical method of decomposing the complex axisymmetric boundary into several ring-shaped and stepped surfaces based on the boundary discretization method (BDM) which overcomes the previous difficulties. In such case, by using the linear wave theory based on eigenfunction expansion matching method, the expressions of velocity potential in each domain, the added mass, radiation damping and wave excitation forces of the oscillating absorbers are obtained. The good astringency of the hydrodynamic coefficients and wave forces are obtained for various geometries when the discrete number reaches a certain value. The captured wave power for a same given draught and displacement for various geometries are calculated and compared. Numerical results show that the geometrical shape has great effect on the wave conversion performance of the absorber. For absorbers with the same outer radius and draught or displacement, the cylindrical type shows fantastic wave energy conversion ability at some given frequencies, while in the random sea wave, the parabolic and conical ones have better stabilization and applicability in wave power conversion.

Application of wavelet scaling function expansion continuous-energy resonance calculation method to MOX fuel problem

International Nuclear Information System (INIS)

Yang, W.; Wu, H.; Cao, L.

2012-01-01

More and more MOX fuels are used in all over the world in the past several decades. Compared with UO 2 fuel, it contains some new features. For example, the neutron spectrum is harder and more resonance interference effects within the resonance energy range are introduced because of more resonant nuclides contained in the MOX fuel. In this paper, the wavelets scaling function expansion method is applied to study the resonance behavior of plutonium isotopes within MOX fuel. Wavelets scaling function expansion continuous-energy self-shielding method is developed recently. It has been validated and verified by comparison to Monte Carlo calculations. In this method, the continuous-energy cross-sections are utilized within resonance energy, which means that it's capable to solve problems with serious resonance interference effects without iteration calculations. Therefore, this method adapts to treat the MOX fuel resonance calculation problem natively. Furthermore, plutonium isotopes have fierce oscillations of total cross-section within thermal energy range, especially for 240 Pu and 242 Pu. To take thermal resonance effect of plutonium isotopes into consideration the wavelet scaling function expansion continuous-energy resonance calculation code WAVERESON is enhanced by applying the free gas scattering kernel to obtain the continuous-energy scattering source within thermal energy range (2.1 eV to 4.0 eV) contrasting against the resonance energy range in which the elastic scattering kernel is utilized. Finally, all of the calculation results of WAVERESON are compared with MCNP calculation. (authors)

Chiral metamaterials characterisation using the wave propagation retrieval method

DEFF Research Database (Denmark)

Andryieuski, Andrei; Lavrinenko, Andrei; Malureanu, Radu

2010-01-01

In this presentation we extend the wave propagation method for the retrieval of the effective properties to the case of chiral metamaterials with circularly polarised eigenwaves. The method is unambiguous, simple and provides bulk effective parameters. Advantages and constraints are discussed...

The development of efficient numerical time-domain modeling methods for geophysical wave propagation

Science.gov (United States)

Zhu, Lieyuan

This Ph.D. dissertation focuses on the numerical simulation of geophysical wave propagation in the time domain including elastic waves in solid media, the acoustic waves in fluid media, and the electromagnetic waves in dielectric media. This thesis shows that a linear system model can describe accurately the physical processes of those geophysical waves' propagation and can be used as a sound basis for modeling geophysical wave propagation phenomena. The generalized stability condition for numerical modeling of wave propagation is therefore discussed in the context of linear system theory. The efficiency of a series of different numerical algorithms in the time-domain for modeling geophysical wave propagation are discussed and compared. These algorithms include the finite-difference time-domain method, pseudospectral time domain method, alternating directional implicit (ADI) finite-difference time domain method. The advantages and disadvantages of these numerical methods are discussed and the specific stability condition for each modeling scheme is carefully derived in the context of the linear system theory. Based on the review and discussion of these existing approaches, the split step, ADI pseudospectral time domain (SS-ADI-PSTD) method is developed and tested for several cases. Moreover, the state-of-the-art stretched-coordinate perfect matched layer (SCPML) has also been implemented in SS-ADI-PSTD algorithm as the absorbing boundary condition for truncating the computational domain and absorbing the artificial reflection from the domain boundaries. After algorithmic development, a few case studies serve as the real-world examples to verify the capacities of the numerical algorithms and understand the capabilities and limitations of geophysical methods for detection of subsurface contamination. The first case is a study using ground penetrating radar (GPR) amplitude variation with offset (AVO) for subsurface non-aqueous-liquid (NAPL) contamination. The

An Unconditionally Stable Method for Solving the Acoustic Wave Equation

Directory of Open Access Journals (Sweden)

Zhi-Kai Fu

2015-01-01

Full Text Available An unconditionally stable method for solving the time-domain acoustic wave equation using Associated Hermit orthogonal functions is proposed. The second-order time derivatives in acoustic wave equation are expanded by these orthogonal basis functions. By applying Galerkin temporal testing procedure, the time variable can be eliminated from the calculations. The restriction of Courant-Friedrichs-Levy (CFL condition in selecting time step for analyzing thin layer can be avoided. Numerical results show the accuracy and the efficiency of the proposed method.

A Novel Grey Wave Method for Predicting Total Chinese Trade Volume

Directory of Open Access Journals (Sweden)

Kedong Yin

2017-12-01

Full Text Available The total trade volume of a country is an important way of appraising its international trade situation. A prediction based on trade volume will help enterprises arrange production efficiently and promote the sustainability of the international trade. Because the total Chinese trade volume fluctuates over time, this paper proposes a Grey wave forecasting model with a Hodrick–Prescott filter (HP filter to forecast it. This novel model first parses time series into long-term trend and short-term cycle. Second, the model uses a general GM (1,1 to predict the trend term and the Grey wave forecasting model to predict the cycle term. Empirical analysis shows that the improved Grey wave prediction method provides a much more accurate forecast than the basic Grey wave prediction method, achieving better prediction results than autoregressive moving average model (ARMA.

Refinements to the method of epicentral location based on surface waves from ambient seismic noise: introducing Love waves

Science.gov (United States)

Levshin, Anatoli L.; Barmin, Mikhail P.; Moschetti, Morgan P.; Mendoza, Carlos; Ritzwoller, Michael H.

2012-01-01

The purpose of this study is to develop and test a modification to a previous method of regional seismic event location based on Empirical Green’s Functions (EGFs) produced from ambient seismic noise. Elastic EGFs between pairs of seismic stations are determined by cross-correlating long ambient noise time-series recorded at the two stations. The EGFs principally contain Rayleigh- and Love-wave energy on the vertical and transverse components, respectively, and we utilize these signals between about 5 and 12 s period. The previous method, based exclusively on Rayleigh waves, may yield biased epicentral locations for certain event types with hypocentral depths between 2 and 5 km. Here we present theoretical arguments that show how Love waves can be introduced to reduce or potentially eliminate the bias. We also present applications of Rayleigh- and Love-wave EGFs to locate 10 reference events in the western United States. The separate Rayleigh and Love epicentral locations and the joint locations using a combination of the two waves agree to within 1 km distance, on average, but confidence ellipses are smallest when both types of waves are used.

On a generalized oscillator system: interbasis expansions

Energy Technology Data Exchange (ETDEWEB)

Kibler, M [Lyon-1 Univ., 69 - Villeurbanne (France). Inst. de Physique Nucleaire; Mardoyan, L G; Pogosyan, G S [Joint Inst. for Nuclear Research, Dubna (Russian Federation). Lab. of Theoretical Physics

1997-12-31

This article deals with a nonrelativistic quantum mechanical study of a dynamical system which generalizes the isotropic harmonic oscillator system in three dimensions. The Schroedinger equation for this generalized oscillator system is separable in spherical, cylindrical, and spheroidal (prolate and oblate) coordinates. The quantum mechanical spectrum of this system is worked out in some details. The problem of interbasis expansions of the wave functions is completely solved. The coefficients for the expansion of the cylindrical basis in terms of the spherical basis, and vice-versa, are found to be analytic continuations (to real values of their arguments) of Clebsch-Gordan coefficients for the group SU(2). The interbasis expansion coefficients for the prolate and oblate spheroidal bases in terms of the spherical or the cylindrical bases are shown to satisfy three-term recursion relations. Finally, a connection between the generalized oscillator system (projected on the z-line) and the Morse system (in one dimension) are discussed. 41 refs.,.

On a generalized oscillator system: interbasis expansions

International Nuclear Information System (INIS)

Kibler, M.; Mardoyan, L.G.; Pogosyan, G.S.

1996-01-01

This article deals with a nonrelativistic quantum mechanical study of a dynamical system which generalizes the isotropic harmonic oscillator system in three dimensions. The Schroedinger equation for this generalized oscillator system is separable in spherical, cylindrical, and spheroidal (prolate and oblate) coordinates. The quantum mechanical spectrum of this system is worked out in some details. The problem of interbasis expansions of the wave functions is completely solved. The coefficients for the expansion of the cylindrical basis in terms of the spherical basis, and vice-versa, are found to be analytic continuations (to real values of their arguments) of Clebsch-Gordan coefficients for the group SU(2). The interbasis expansion coefficients for the prolate and oblate spheroidal bases in terms of the spherical or the cylindrical bases are shown to satisfy three-term recursion relations. Finally, a connection between the generalized oscillator system (projected on the z-line) and the Morse system (in one dimension) are discussed. 41 refs.,

Modeling and analysis of waves in a heat conducting thermo-elastic plate of elliptical shape

Directory of Open Access Journals (Sweden)

R. Selvamani

Full Text Available Wave propagation in heat conducting thermo elastic plate of elliptical cross-section is studied using the Fourier expansion collocation method based on Suhubi's generalized theory. The equations of motion based on two-dimensional theory of elasticity is applied under the plane strain assumption of generalized thermo elastic plate of elliptical cross-sections composed of homogeneous isotropic material. The frequency equations are obtained by using the boundary conditions along outer and inner surface of elliptical cross-sectional plate using Fourier expansion collocation method. The computed non-dimensional frequency, velocity and quality factor are plotted in dispersion curves for longitudinal and flexural (symmetric and antisymmetric modes of vibrations.

The Peano-series solution for modeling shear horizontal waves in piezoelectric plates

Directory of Open Access Journals (Sweden)

Ben Ghozlen M.H.

2012-06-01

Full Text Available The shear horizontal (SH wave devices have been widely used in electroacoustic. To improve their performance, the phase velocity dispersion and the electromechanical coupling coefficient of the Lamb wave should be calculated exactly in the design. Therefore, this work is to analyze exactly the Lamb waves polarized in the SH direction in homogeneous plate pie.zoelectric material (PZT-5H. An alternative method is proposed to solve the wave equation in such a structure without using the standard method based on the electromechanical partial waves. This method is based on an analytical solution, the matricant explicitly expressed under the Peano series expansion form. Two types of configuration have been addressed, namely the open circuited and the short circuited. Results confirm that the SH wave provides a number of attractive properties for use in sensing and signal processing applications. It has been found that the phase velocity remains nearly constant for all values of h/λ (h is the plate thickness, λ the acoustic wavelength. Secondly the SH0 wave mode can provide very high electromechanical coupling. Graphical representations of electrical and mechanical amounts function of depth are made, they are in agreement with the continuity rules. The developed Peano technique is in agreement with the classical approach, and can be suitable with cylindrical geometry.

An alternative solver for the nodal expansion method equations - 106

International Nuclear Information System (INIS)

Carvalho da Silva, F.; Carlos Marques Alvim, A.; Senra Martinez, A.

2010-01-01

An automated procedure for nuclear reactor core design is accomplished by using a quick and accurate 3D nodal code, aiming at solving the diffusion equation, which describes the spatial neutron distribution in the reactor. This paper deals with an alternative solver for nodal expansion method (NEM), with only two inner iterations (mesh sweeps) per outer iteration, thus having the potential to reduce the time required to calculate the power distribution in nuclear reactors, but with accuracy similar to the ones found in conventional NEM. The proposed solver was implemented into a computational system which, besides solving the diffusion equation, also solves the burnup equations governing the gradual changes in material compositions of the core due to fuel depletion. Results confirm the effectiveness of the method for practical purposes. (authors)

Multiple travelling wave solutions of nonlinear evolution equations using a unified algebraic method

International Nuclear Information System (INIS)

Fan Engui

2002-01-01

A new direct and unified algebraic method for constructing multiple travelling wave solutions of general nonlinear evolution equations is presented and implemented in a computer algebraic system. Compared with most of the existing tanh methods, the Jacobi elliptic function method or other sophisticated methods, the proposed method not only gives new and more general solutions, but also provides a guideline to classify the various types of the travelling wave solutions according to the values of some parameters. The solutions obtained in this paper include (a) kink-shaped and bell-shaped soliton solutions, (b) rational solutions, (c) triangular periodic solutions and (d) Jacobi and Weierstrass doubly periodic wave solutions. Among them, the Jacobi elliptic periodic wave solutions exactly degenerate to the soliton solutions at a certain limit condition. The efficiency of the method can be demonstrated on a large variety of nonlinear evolution equations such as those considered in this paper, KdV-MKdV, Ito's fifth MKdV, Hirota, Nizhnik-Novikov-Veselov, Broer-Kaup, generalized coupled Hirota-Satsuma, coupled Schroedinger-KdV, (2+1)-dimensional dispersive long wave, (2+1)-dimensional Davey-Stewartson equations. In addition, as an illustrative sample, the properties of the soliton solutions and Jacobi doubly periodic solutions for the Hirota equation are shown by some figures. The links among our proposed method, the tanh method, extended tanh method and the Jacobi elliptic function method are clarified generally. (author)

Overdetermined shooting methods for computing standing water waves with spectral accuracy

International Nuclear Information System (INIS)

Wilkening, Jon; Yu Jia

2012-01-01

A high-performance shooting algorithm is developed to compute time-periodic solutions of the free-surface Euler equations with spectral accuracy in double and quadruple precision. The method is used to study resonance and its effect on standing water waves. We identify new nucleation mechanisms in which isolated large-amplitude solutions, and closed loops of such solutions, suddenly exist for depths below a critical threshold. We also study degenerate and secondary bifurcations related to Wilton's ripples in the traveling case, and explore the breakdown of self-similarity at the crests of extreme standing waves. In shallow water, we find that standing waves take the form of counter-propagating solitary waves that repeatedly collide quasi-elastically. In deep water with surface tension, we find that standing waves resemble counter-propagating depression waves. We also discuss the existence and non-uniqueness of solutions, and smooth versus erratic dependence of Fourier modes on wave amplitude and fluid depth. In the numerical method, robustness is achieved by posing the problem as an overdetermined nonlinear system and using either adjoint-based minimization techniques or a quadratically convergent trust-region method to minimize the objective function. Efficiency is achieved in the trust-region approach by parallelizing the Jacobian computation, so the setup cost of computing the Dirichlet-to-Neumann operator in the variational equation is not repeated for each column. Updates of the Jacobian are also delayed until the previous Jacobian ceases to be useful. Accuracy is maintained using spectral collocation with optional mesh refinement in space, a high-order Runge–Kutta or spectral deferred correction method in time and quadruple precision for improved navigation of delicate regions of parameter space as well as validation of double-precision results. Implementation issues for transferring much of the computation to a graphic processing units are briefly

Investigation on the reliability of expansion joint for piping with probabilistic method

International Nuclear Information System (INIS)

Ishii, Y.; Kambe, M.

1980-01-01

The reduction of the plant size is necessitated as one of the major targets in LMFBR design. Usually, piping work system is extensively used to absorb thermal expansion between two components anywhere. Besides above, expansion joint for piping seems to be attractive lately for the same object. This paper describes the significance of expansion joint with multiple boundaries, breakdown probability of expansion joint assembly and partly the bellows by introducing several hypothetical conditions in connection with piping. Also, an importance of in-service inspection (ISI) for expansion joint was discussed using a comparative table and probabilities on reliability from partly broken to full penetration. In conclusion, the expansion joint with ISI should be manufactured with excellent reliability in order to cope with piping work system; several conditions of the practical application for piping systems are suggested. (author)

Investigation on the reliability of expansion joint for piping with probabilistic method

Energy Technology Data Exchange (ETDEWEB)

Ishii, Y; Kambe, M

1980-02-01

The reduction of the plant size is necessitated as one of the major targets in LMFBR design. Usually, piping work system is extensively used to absorb thermal expansion between two components anywhere. Besides above, expansion joint for piping seems to be attractive lately for the same object. This paper describes the significance of expansion joint with multiple boundaries, breakdown probability of expansion joint assembly and partly the bellows by introducing several hypothetical conditions in connection with piping. Also, an importance of in-service inspection (ISI) for expansion joint was discussed using a comparative table and probabilities on reliability from partly broken to full penetration. In conclusion, the expansion joint with ISI should be manufactured with excellent reliability in order to cope with piping work system; several conditions of the practical application for piping systems are suggested. (author)

Investigation on the reliability of expansion joint for piping with probabilistic method

International Nuclear Information System (INIS)

Ishii, Yoichiro; Kambe, Mitsuru.

1979-11-01

The reduction of the plant size if necessitated as one of the major target in LMFBR design. Usually, piping work system is extensively used to absorb thermal expansion between two components anywhere. Besides above, expansion joint for piping seems to be attractive lately for the same object. This paper describes about the significance of expansion joint with multiple boundaries, breakdown probability of expansion joint assembly and partly the bellows by introducing several hypothetical conditions in connection with piping. Also, an importance of inservice inspection (ISI) for expansion joint was discussed using by comparative table and probabilities on reliability from partly broken to full penetration. In the conclusion, the expansion joint with ISI should be manufactured with excellent reliability in order to cope with piping work system, and several conditions of the practical application for piping systems are suggested. (author)

Methods of abdominal wall expansion for repair of incisional herniae: a systematic review.

Science.gov (United States)

Alam, N N; Narang, S K; Pathak, S; Daniels, I R; Smart, N J

2016-04-01

To systematically review the available literature regarding methods for abdominal wall expansion and compare the outcome of primary fascial closure rates. A systematic search of Pubmed and Embase databases was conducted using the search terms "Abdominal wall hernia", "ventral hernia", "midline hernia", "Botulinum toxin", "botox", "dysport", "progressive preoperative pneumoperitoneum", and "tissue expanders". Study quality was assessed using the Methodological Index for Non-Randomised Studies. 21 of the 105 studies identified met the inclusion criteria. Progressive preoperative pneumoperitoneum (PPP) was performed in 269 patients across 15 studies with primary fascial closure being achieved in 226 (84%). 16 patients had a recurrence (7.2%) and the complication rate was 12% with 2 reported mortalities. There were 4 studies with 14 patients in total undergoing abdominal wall expansion using tissue expanders with a fascial closure rate of 92.9% (n = 13). A recurrence rate of 10.0% (n = 1) was reported with 1 complication and no mortalities. Follow up ranged from 3 to 36 months across the studies. There were 2 studies reporting the use of botulinum toxin with 29 patients in total. A primary fascial closure rate of 100% (n = 29) was demonstrated although a combination of techniques including component separation and Rives-Stoppa repair were used. There were no reported complications related to the use of Botulinum Toxin. However, the short-term follow up in many cases and the lack of routine radiological assessment for recurrence suggests that the recurrence rate has been underestimated. PPP, tissue expanders and Botulinum toxin are safe and feasible methods for abdominal wall expansion prior to incisional hernia repair. In combination with existing techniques for repair, these methods may help provide the crucial extra tissue mobility required to achieve primary closure.

Coronary wave energy: a novel predictor of functional recovery after myocardial infarction.

Science.gov (United States)

De Silva, Kalpa; Foster, Paul; Guilcher, Antoine; Bandara, Asela; Jogiya, Roy; Lockie, Tim; Chowiencyzk, Phil; Nagel, Eike; Marber, Michael; Redwood, Simon; Plein, Sven; Perera, Divaka

2013-04-01

Revascularization after acute coronary syndromes provides prognostic benefit, provided that the subtended myocardium is viable. The microcirculation and contractility of the subtended myocardium affect propagation of coronary flow, which can be characterized by wave intensity analysis. The study objective was to determine in acute coronary syndromes whether early wave intensity analysis-derived microcirculatory (backward) expansion wave energy predicts late viability, defined by functional recovery. Thirty-one patients (58±11 years) were enrolled after non-ST elevation myocardial infarction. Regional left ventricular function and late-gadolinium enhancement were assessed by cardiac magnetic resonance imaging, before and 3 months after revascularization. The backward-traveling (microcirculatory) expansion wave was derived from wave intensity analysis of phasic coronary pressure and velocity in the infarct-related artery, whereas mean values were used to calculate hyperemic microvascular resistance. Twelve-hour troponin T, left ventricular ejection fraction, and percentage late-gadolinium enhancement mass were 1.35±1.21 µg/L, 56±11%, and 8.4±6.0%, respectively. The infarct-related artery backward-traveling (microcirculatory) expansion wave was inversely correlated with late-gadolinium enhancement infarct mass (r=-0.81; Pwave threshold of 2.8 W m(-2) s(-2)×10(5) predicted functional recovery with sensitivity and specificity of 0.91 and 0.82 (AUC 0.88). Hyperemic microvascular resistance correlated with late-gadolinium enhancement mass (r=0.48; P=0.03) but not left ventricular recovery (r=-0.34; P=0.07). The microcirculation-derived backward expansion wave is a new index that correlates with the magnitude and location of infarction, which may allow for the prediction of functional myocardial recovery. Coronary wave intensity analysis may facilitate myocardial viability assessment during cardiac catheterization.

TSOS and TSOS-FK hybrid methods for modelling the propagation of seismic waves

Science.gov (United States)

Ma, Jian; Yang, Dinghui; Tong, Ping; Ma, Xiao

2018-05-01

We develop a new time-space optimized symplectic (TSOS) method for numerically solving elastic wave equations in heterogeneous isotropic media. We use the phase-preserving symplectic partitioned Runge-Kutta method to evaluate the time derivatives and optimized explicit finite-difference (FD) schemes to discretize the space derivatives. We introduce the averaged medium scheme into the TSOS method to further increase its capability of dealing with heterogeneous media and match the boundary-modified scheme for implementing free-surface boundary conditions and the auxiliary differential equation complex frequency-shifted perfectly matched layer (ADE CFS-PML) non-reflecting boundaries with the TSOS method. A comparison of the TSOS method with analytical solutions and standard FD schemes indicates that the waveform generated by the TSOS method is more similar to the analytic solution and has a smaller error than other FD methods, which illustrates the efficiency and accuracy of the TSOS method. Subsequently, we focus on the calculation of synthetic seismograms for teleseismic P- or S-waves entering and propagating in the local heterogeneous region of interest. To improve the computational efficiency, we successfully combine the TSOS method with the frequency-wavenumber (FK) method and apply the ADE CFS-PML to absorb the scattered waves caused by the regional heterogeneity. The TSOS-FK hybrid method is benchmarked against semi-analytical solutions provided by the FK method for a 1-D layered model. Several numerical experiments, including a vertical cross-section of the Chinese capital area crustal model, illustrate that the TSOS-FK hybrid method works well for modelling waves propagating in complex heterogeneous media and remains stable for long-time computation. These numerical examples also show that the TSOS-FK method can tackle the converted and scattered waves of the teleseismic plane waves caused by local heterogeneity. Thus, the TSOS and TSOS-FK methods proposed in

A Full-wave Model for Wave Propagation and Dissipation in the Inner Magnetosphere Using the Finite Element Method

International Nuclear Information System (INIS)

Valeo, Ernest; Johnson, Jay R.; Kim, Eun-Hwa; Phillips, Cynthia

2012-01-01

A wide variety of plasma waves play an important role in the energization and loss of particles in the inner magnetosphere. Our ability to understand and model wave-particle interactions in this region requires improved knowledge of the spatial distribution and properties of these waves as well as improved understanding of how the waves depend on changes in solar wind forcing and/or geomagnetic activity. To this end, we have developed a two-dimensional, finite element code that solves the full wave equations in global magnetospheric geometry. The code describes three-dimensional wave structure including mode conversion when ULF, EMIC, and whistler waves are launched in a two-dimensional axisymmetric background plasma with general magnetic field topology. We illustrate the capabilities of the code by examining the role of plasmaspheric plumes on magnetosonic wave propagation; mode conversion at the ion-ion and Alfven resonances resulting from external, solar wind compressions; and wave structure and mode conversion of electromagnetic ion cyclotron waves launched in the equatorial magnetosphere, which propagate along the magnetic field lines toward the ionosphere. We also discuss advantages of the finite element method for resolving resonant structures, and how the model may be adapted to include nonlocal kinetic effects.

Two modified symplectic partitioned Runge-Kutta methods for solving the elastic wave equation

Science.gov (United States)

Su, Bo; Tuo, Xianguo; Xu, Ling

2017-08-01

Based on a modified strategy, two modified symplectic partitioned Runge-Kutta (PRK) methods are proposed for the temporal discretization of the elastic wave equation. The two symplectic schemes are similar in form but are different in nature. After the spatial discretization of the elastic wave equation, the ordinary Hamiltonian formulation for the elastic wave equation is presented. The PRK scheme is then applied for time integration. An additional term associated with spatial discretization is inserted into the different stages of the PRK scheme. Theoretical analyses are conducted to evaluate the numerical dispersion and stability of the two novel PRK methods. A finite difference method is used to approximate the spatial derivatives since the two schemes are independent of the spatial discretization technique used. The numerical solutions computed by the two new schemes are compared with those computed by a conventional symplectic PRK. The numerical results, which verify the new method, are superior to those generated by traditional conventional methods in seismic wave modeling.

A method of acoustic wave registration and determination their generation region

International Nuclear Information System (INIS)

Kozin, I.D.; Marchenko, M.V.

1998-01-01

Here is presented a method of acoustic wave registration with using of a synchronous LF broadcasting system. This method of detection and determination of underground nuclear explosion location is based on a registration of ionospheric disturbances induced by acoustic waves at the region of LF sign al reflection. The measuring complex created in the institute of the Ionosphere /1/ allows to register amplitude-frequency characteristics of composite signal from synchronous broadcasting net

Scattering of Lamb waves in a composite plate

Science.gov (United States)

Bratton, Robert; Datta, Subhendu; Shah, Arvind

1991-01-01

A combined analytical and finite element technique is developed to gain a better understanding of the scattering of elastic waves by defects. This hybrid method is capable of predicting scattered displacements from arbitrary shaped defects as well as inclusions of different material. The continuity of traction and displacements at the boundaries of the two areas provided the necessary equations to find the nodal displacements and expansion coefficients. Results clearly illustrate the influence of increasing crack depth on the scattered signal.

Experiments of Long-range Inspection Method in Straight Pipes using Ultrasonic Guided Waves

International Nuclear Information System (INIS)

Eom, H. S.; Lim, S. H.; Kim, J. H.; Joo, Y.S.

2006-02-01

This report describes experimental results of a long-range inspection method of pipes using ultrasonic guided waves. In chapter 2, theory of guided wave was reviewed. In chapter 3, equipment and procedures which were used in the experiments were described. Detailed specifications of the specimens described in chapter 4. In chapter 5, we analyzed characteristics of guided wave signals according to shapes and sizes of defects and presents results of various signal processing methods

Application of Wave Distribution Function Method to the ERG/PWE Data

Science.gov (United States)

Ota, M.; Kasahara, Y.; Matsuda, S.; Kojima, H.; Matsuoka, A.; Hikishima, M.; Kasaba, Y.; Ozaki, M.; Yagitani, S.; Tsuchiya, F.; Kumamoto, A.

2017-12-01

The ERG (Arase) satellite was launched on 20 December 2016 to study acceleration and loss mechanisms of relativistic electrons in the Earth's magnetosphere. The Plasma Wave Experiment (PWE), which is one of the science instruments on board the ERG satellite, measures electric field and magnetic field. The PWE consists of three sub-systems; EFD (Electric Field Detector), OFA/WFC (Onboard Frequency Analyzer and Waveform Capture), and HFA (High Frequency Analyzer).The OFA/WFC measures electromagnetic field spectra and raw waveforms in the frequency range from few Hz to 20 kHz. The OFA produces three kind of data; OFA-SPEC (power spectrum), OFA-MATRIX (spectral matrix), and OFA-COMPLEX (complex spectrum). The OFA-MATRIX measures ensemble averaged complex cross-spectra of two electric field components, and of three magnetic field components. The OFA-COMPLEX measures instantaneous complex spectra of electric and magnetic fields. These data are produced every 8 seconds in the nominal mode, and it can be used for polarization analysis and wave propagation direction finding.In general, spectral matrix composed by cross-spectra of observed signals is used for direction finding, and many algorithms have been proposed. For example, Means method and SVD method can be applied on the assumption that the spectral matrix is consists of a single plane wave, while wave distribution function (WDF) method is applicable even to the data in which multiple numbers of plane waves are simultaneously included. In this presentation, we introduce the results when the WDF method is applied to the ERG/PWE data.

Comparative numerical solutions of stiff Ordinary differential equations using magnus series expansion method

Directory of Open Access Journals (Sweden)

SURE KÖME

2014-12-01

Full Text Available In this paper, we investigated the effect of Magnus Series Expansion Method on homogeneous stiff ordinary differential equations with different stiffness ratios. A Magnus type integrator is used to obtain numerical solutions of two different examples of stiff problems and exact and approximate results are tabulated. Furthermore, absolute error graphics are demonstrated in detail.

Thermal gravitational waves in accelerating universe

Directory of Open Access Journals (Sweden)

B Ghayour

2013-10-01

Full Text Available Gravitational waves are considered in thermal vacuum state. The amplitude and spectral energy density of gravitational waves are found enhanced in thermal vacuum state compared to its zero temperature counterpart. Therefore, the allowed amount of enhancement depends on the upper bound of WMAP-5 and WMAP-7 for the amplitude and spectral energy density of gravitational waves. The enhancement of amplitude and spectral energy density of the waves in thermal vacuum state is consistent with current accelerating phase of the universe. The enhancement feature of amplitude and spectral energy density of the waves is independent of the expansion model of the universe and hence the thermal effect accounts for it. Therefore, existence of thermal gravitational waves is not ruled out

Thermal expansion of CeCu5.8Ag0.2

International Nuclear Information System (INIS)

Kuechler, R.; Gegenwart, P.; Heuser, K.; Scheidt, E.-W.; Stewart, G.R.; Steglich, F.

2005-01-01

We present low-temperature thermal expansion measurements on the heavy fermion system CeCu 5.8 Ag 0.2 , which is located at an antiferromagnetic (AF) quantum critical point (QCP). At zero magnetic field, the volume expansion coefficient divided by temperature shows a logarithmic divergence upon cooling below 1K. This temperature dependence is incompatible with the predictions of the itinerant spin-density wave theory for an AF QCP. The application of magnetic fields leads to a cross-over to Landau Fermi liquid behavior as expected for a zero-field QCP

Convergence of the multiple scattering expansion in XAFS and XANES

International Nuclear Information System (INIS)

Rehr, J.J.

1992-01-01

The convergence of the multiple-scattering expansion of XAFS and XANES by explicit path-bypath calculations. The approach is based on the fast scattering matrix formalism of Rehr and Albers, together with an automated path finder and filters that exclude negligible paths. High-order scattering terms are found to be essential, especially at low energies. Several factors including the magnitude of curved wave scattering amplitudes, inelastic losses and multiple-scattering Debye-Waller factors control convergence of the expansion. The convergence is illustrated explicitly for the case of diatomic molecules

Stability analysis of CMFD acceleration for the wavelet expansion method of neutron transport equation

International Nuclear Information System (INIS)

Zheng Youqi; Wu Hongchun; Cao Liangzhi

2013-01-01

This paper describes the stability analysis for the coarse mesh finite difference (CMFD) acceleration used in the wavelet expansion method. The nonlinear CMFD acceleration scheme is transformed by linearization and the Fourier ansatz is introduced into the linearized formulae. The spectral radius is defined as the stability criterion, which is the least upper bound (LUB) of the largest eigenvalue of Fourier analysis matrix. The stability analysis considers the effect of mesh size (spectral length), coarse mesh division and scattering ratio. The results show that for the wavelet expansion method, the CMFD acceleration is conditionally stable. The small size of fine mesh brings stability and fast convergent. With the increase of the mesh size, the stability becomes worse. The scattering ratio does not impact the stability obviously. It makes the CMFD acceleration highly efficient in the strong scattering case. The results of Fourier analysis are verified by the numerical tests based on a homogeneous slab problem.

Tunnel Vision Prismatic Field Expansion: Challenges and Requirements.

Science.gov (United States)

Apfelbaum, Henry; Peli, Eli

2015-12-01

No prismatic solution for peripheral field loss (PFL) has gained widespread acceptance. Field extended by prisms has a corresponding optical scotoma at the prism apices. True expansion can be achieved when each eye is given a different view (through visual confusion). We analyze the effects of apical scotomas and binocular visual confusion in different designs to identify constraints on any solution that is likely to meet acceptance. Calculated perimetry diagrams were compared to perimetry with PFL patients wearing InWave channel prisms and Trifield spectacles. Percept diagrams illustrate the binocular visual confusion. Channel prisms provide no benefit at primary gaze. Inconsequential extension was provided by InWave prisms, although accessible with moderate gaze shifts. Higher-power prisms provide greater extension, with greater paracentral scotoma loss, but require uncomfortable gaze shifts. Head turns, not eye scans, are needed to see regions lost to the apical scotomas. Trifield prisms provide field expansion at all gaze positions, but acceptance was limited by disturbing effects of central binocular visual confusion. Field expansion when at primary gaze (where most time is spent) is needed while still providing unobstructed central vision. Paracentral multiplexing prisms we are developing that superimpose shifted and see-through views may accomplish that. Use of the analyses and diagramming techniques presented here will be of value when considering prismatic aids for PFL, and could have prevented many unsuccessful designs and the improbable reports we cited from the literature. New designs must likely address the challenges identified here.

Tur\\'an type inequalities for regular Coulomb wave functions

OpenAIRE

Baricz, Árpád

2015-01-01

Tur\\'an, Mitrinovi\\'c-Adamovi\\'c and Wilker type inequalities are deduced for regular Coulomb wave functions. The proofs are based on a Mittag-Leffler expansion for the regular Coulomb wave function, which may be of independent interest. Moreover, some complete monotonicity results concerning the Coulomb zeta functions and some interlacing properties of the zeros of Coulomb wave functions are given.

Perfectly Matched Layer for the Wave Equation Finite Difference Time Domain Method

Science.gov (United States)

Miyazaki, Yutaka; Tsuchiya, Takao

2012-07-01

The perfectly matched layer (PML) is introduced into the wave equation finite difference time domain (WE-FDTD) method. The WE-FDTD method is a finite difference method in which the wave equation is directly discretized on the basis of the central differences. The required memory of the WE-FDTD method is less than that of the standard FDTD method because no particle velocity is stored in the memory. In this study, the WE-FDTD method is first combined with the standard FDTD method. Then, Berenger's PML is combined with the WE-FDTD method. Some numerical demonstrations are given for the two- and three-dimensional sound fields.

Runge-Kutta Integration of the Equal Width Wave Equation Using the Method of Lines

Directory of Open Access Journals (Sweden)

M. A. Banaja

2015-01-01

Full Text Available The equal width (EW equation governs nonlinear wave phenomena like waves in shallow water. Numerical solution of the (EW equation is obtained by using the method of lines (MOL based on Runge-Kutta integration. Using von Neumann stability analysis, the scheme is found to be unconditionally stable. Solitary wave motion and interaction of two solitary waves are studied using the proposed method. The three invariants of the motion are evaluated to determine the conservation properties of the generated scheme. Accuracy of the proposed method is discussed by computing the L2 and L∞ error norms. The results are found in good agreement with exact solution.

Measuring longitudinal wave speed in solids: two methods and a half

International Nuclear Information System (INIS)

Fazio, C; Guastella, I; Sperandeo-Mineo, R M; Tarantino, G

2006-01-01

Three methods to analyse longitudinal wave propagation in metallic rods are discussed. Two of these methods also prove to be useful for measuring the sound propagation speed. The experimental results, as well as some interpretative models built in the context of a workshop on mechanical waves at the Graduate School for Pre-Service Physics Teacher Education, Palermo University, are described. Some considerations about observed modifications in trainee teachers' attitudes to utilizing physics experiments to build pedagogical activities are discussed

Analysis and design of efficient planar leaky-wave antennas

NARCIS (Netherlands)

Ettore, M.

2008-01-01

This thesis deals with the effective design of planar leaky-wave antennas. The work describes a methodology based on the polar expansion of Green's function representations to address very different geometrical configurations which might appear to have little in common. In fact leaky waves with

Imaging of THz waves in 2D photonic crystal structures embedded in a slab waveguide

International Nuclear Information System (INIS)

Peier, P; Merbold, H; Feurer, T; Pahinin, V; Nelson, K A

2010-01-01

We present space- and time-resolved simulations and measurements of single-cycle terahertz (THz) waves propagating through two-dimensional (2D) photonic crystal structures embedded in a slab waveguide. Specifically, we use a plane wave expansion technique to calculate the band structure and a time-dependent finite-element method to simulate the temporal evolution of the THz waves. Experimentally, we measure the space-time evolution of the THz waves through a coherent time-resolved imaging method. Three different structures are laser machined in LiNbO 3 crystal slabs and analyzing the transmitted as well as the reflected THz waveforms allows determination of the bandgaps. Comparing the results with the calculated band diagrams and the time-dependent simulations shows that the experiments are consistent with 3D simulations, which include the slab waveguide geometry, the birefringence of the material, and a careful analysis of the excited modes within the band diagrams.

Piecewise parabolic method for simulating one-dimensional shear shock wave propagation in tissue-mimicking phantoms

Science.gov (United States)

Tripathi, B. B.; Espíndola, D.; Pinton, G. F.

2017-11-01

The recent discovery of shear shock wave generation and propagation in the porcine brain suggests that this new shock phenomenology may be responsible for a broad range of traumatic injuries. Blast-induced head movement can indirectly lead to shear wave generation in the brain, which could be a primary mechanism for injury. Shear shock waves amplify the local acceleration deep in the brain by up to a factor of 8.5, which may tear and damage neurons. Currently, there are numerical methods that can model compressional shock waves, such as comparatively well-studied blast waves, but there are no numerical full-wave solvers that can simulate nonlinear shear shock waves in soft solids. Unlike simplified representations, e.g., retarded time, full-wave representations describe fundamental physical behavior such as reflection and heterogeneities. Here we present a piecewise parabolic method-based solver for one-dimensional linearly polarized nonlinear shear wave in a homogeneous medium and with empirical frequency-dependent attenuation. This method has the advantage of being higher order and more directly extendable to multiple dimensions and heterogeneous media. The proposed numerical scheme is validated analytically and experimentally and compared to other shock capturing methods. A Riemann step-shock problem is used to characterize the numerical dissipation. This dissipation is then tuned to be negligible with respect to the physical attenuation by choosing an appropriate grid spacing. The numerical results are compared to ultrasound-based experiments that measure planar polarized shear shock wave propagation in a tissue-mimicking gelatin phantom. Good agreement is found between numerical results and experiment across a 40 mm propagation distance. We anticipate that the proposed method will be a starting point for the development of a two- and three-dimensional full-wave code for the propagation of nonlinear shear waves in heterogeneous media.

Ultrasound viscoelasticity assessment using an adaptive torsional shear wave propagation method

Energy Technology Data Exchange (ETDEWEB)

Ouared, Abderrahmane [Laboratory of Biorheology and Medical Ultrasonics, University of Montréal Hospital Research Center (CRCHUM), Montréal, Québec H2X 0A9, Canada and Institute of Biomedical Engineering, University of Montréal, Montréal, Québec H3T 1J4 (Canada); Kazemirad, Siavash; Montagnon, Emmanuel [Laboratory of Biorheology and Medical Ultrasonics, University of Montréal Hospital Research Center (CRCHUM), Montréal, Québec H2X 0A9 (Canada); Cloutier, Guy, E-mail: guy.cloutier@umontreal.ca [Laboratory of Biorheology and Medical Ultrasonics, University of Montréal Hospital Research Center (CRCHUM), Montréal, Québec H2X 0A9 (Canada); Department of Radiology, Radio-Oncology and Nuclear Medicine, University of Montréal, Montréal, Québec H3T 1J4 (Canada); Institute of Biomedical Engineering, University of Montréal, Montréal, Québec H3T 1J4 (Canada)

2016-04-15

Purpose: Different approaches have been used in dynamic elastography to assess mechanical properties of biological tissues. Most techniques are based on a simple inversion based on the measurement of the shear wave speed to assess elasticity, whereas some recent strategies use more elaborated analytical or finite element method (FEM) models. In this study, a new method is proposed for the quantification of both shear storage and loss moduli of confined lesions, in the context of breast imaging, using adaptive torsional shear waves (ATSWs) generated remotely with radiation pressure. Methods: A FEM model was developed to solve the inverse wave propagation problem and obtain viscoelastic properties of interrogated media. The inverse problem was formulated and solved in the frequency domain and its robustness to noise and geometric constraints was evaluated. The proposed model was validated in vitro with two independent rheology methods on several homogeneous and heterogeneous breast tissue-mimicking phantoms over a broad range of frequencies (up to 400 Hz). Results: Viscoelastic properties matched benchmark rheology methods with discrepancies of 8%–38% for the shear modulus G′ and 9%–67% for the loss modulus G″. The robustness study indicated good estimations of storage and loss moduli (maximum mean errors of 19% on G′ and 32% on G″) for signal-to-noise ratios between 19.5 and 8.5 dB. Larger errors were noticed in the case of biases in lesion dimension and position. Conclusions: The ATSW method revealed that it is possible to estimate the viscoelasticity of biological tissues with torsional shear waves when small biases in lesion geometry exist.

The principal component analysis method used with polynomial Chaos expansion to propagate uncertainties through critical transport problems

Energy Technology Data Exchange (ETDEWEB)

Rising, M. E.; Prinja, A. K. [Univ. of New Mexico, Dept. of Chemical and Nuclear Engineering, Albuquerque, NM 87131 (United States)

2012-07-01

A critical neutron transport problem with random material properties is introduced. The total cross section and the average neutron multiplicity are assumed to be uncertain, characterized by the mean and variance with a log-normal distribution. The average neutron multiplicity and the total cross section are assumed to be uncorrected and the material properties for differing materials are also assumed to be uncorrected. The principal component analysis method is used to decompose the covariance matrix into eigenvalues and eigenvectors and then 'realizations' of the material properties can be computed. A simple Monte Carlo brute force sampling of the decomposed covariance matrix is employed to obtain a benchmark result for each test problem. In order to save computational time and to characterize the moments and probability density function of the multiplication factor the polynomial chaos expansion method is employed along with the stochastic collocation method. A Gauss-Hermite quadrature set is convolved into a multidimensional tensor product quadrature set and is successfully used to compute the polynomial chaos expansion coefficients of the multiplication factor. Finally, for a particular critical fuel pin assembly the appropriate number of random variables and polynomial expansion order are investigated. (authors)

A parallel orbital-updating based plane-wave basis method for electronic structure calculations

International Nuclear Information System (INIS)

Pan, Yan; Dai, Xiaoying; Gironcoli, Stefano de; Gong, Xin-Gao; Rignanese, Gian-Marco; Zhou, Aihui

2017-01-01

Highlights: • Propose three parallel orbital-updating based plane-wave basis methods for electronic structure calculations. • These new methods can avoid the generating of large scale eigenvalue problems and then reduce the computational cost. • These new methods allow for two-level parallelization which is particularly interesting for large scale parallelization. • Numerical experiments show that these new methods are reliable and efficient for large scale calculations on modern supercomputers. - Abstract: Motivated by the recently proposed parallel orbital-updating approach in real space method , we propose a parallel orbital-updating based plane-wave basis method for electronic structure calculations, for solving the corresponding eigenvalue problems. In addition, we propose two new modified parallel orbital-updating methods. Compared to the traditional plane-wave methods, our methods allow for two-level parallelization, which is particularly interesting for large scale parallelization. Numerical experiments show that these new methods are more reliable and efficient for large scale calculations on modern supercomputers.