Vibration decay along deepwater moorings with complex geometry and configuration, as a tool for truncation

This paper is concerned with the difficulties in model testing deepwater structures at reasonable scales. An overview of recent research efforts to tackle this challenge is given first, introducing the concept of line truncation. Passive truncation has traditionally been the preferred method by industry; however, these techniques tend to suffer in capturing accurately line dynamic response and so reproducing peak tensions. In an attempt to improve credibility of model test data the proposed truncation procedure sets up the truncated model, based on line dynamic response rather than quasi-static system stiffness. Vibration decay of transverse elastic waves due to fluid drag forces is assessed and it is found that below a certain length criterion, the transverse vibrational characteristics for each line are inertia driven, hence with respect to these motions the truncated model can assume a linear damper whose coefficient depends on the local line properties and vibration frequency. Initially a simplified taut string model is assumed for which the line is submerged in still water, one end fixed at the bottom the other assumed to follow the vessel response, which can be harmonic or random. A dimensional analysis, supported by exact benchmark numerical solutions, has shown that it is possible to produce a general guideline for the truncation length criterion, which is suitable for any kind of line with any top motion. The focus of this paper is to extend this work to a more complex line configuration of a conventional deepwater mooring line and so enhance the generality of the truncation guideline. The paper will close with an example case study of a spread mooring system, applying this method to create an equivalent numerical model at a reduced depth that replicates exactly the static and dynamic characteristics of the full depth system. Copyright © 2012 by ASME.

地震波传播模拟的边界元―体积元积分方程法及其近似求解

To improve the efficiency of boundary-volume integral equation technique, this paper is involved in the approximate solutions of boundary-volume integral equation technique. Firstly, based on different interpretations of the self-interaction and extrapolation operators of the resulting boundary integral equation matrix, two different hybrid BEM+Born series modeling schemes are formulated and validated through comparisons with the full-waveform BE numerical solutions for wave propagation simulation in a semicircular alluvial valley and a complex fault model respectively. Numerical experiments indicate that both the BEM+Born series modeling schemes are suitable for complex geological structures and significantly improve computational efficiency especially for the cases of high frequencies and multisource seismic survey. Then boundary-volume integral equation technique is illuminated in detail and verified by modeling wave propagation in complex media. Furthermore, the first-order and second-order Born approximate solutions for the volume-scattering waves are studied and quantified by numerical simulation in different random medium models. Finally, preconditioning generalized minimal residual method is applied to solve boundary-volume integral equation and compared with Gaussian elimination method. Numerical experiments indicate this method makes the calculations more efficient.

特大地震激发地球自由振荡的数值模拟研究

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.

各向异性与非线性介质中弹性波场传播特征研究

The topic of this study is about the propagation features of elastic waves in the anisotropic and nonlinear media by numerical methods with high accuracy and stability. The main achievements of this paper are as followings: Firstly, basing on the third order elastic energy formula, principle of energy conservation and circumvolved matrix method, we firstly reported the equations of non-linear elastic waves with two dimensions and three components in VTI media. Secondly, several conclusions about some numerical methods have been obtained in this paper. Namely, the minimum suitable sample stepth in space is about 1/8-1/12 of the main wavelength in order to distinctly reduce the numerical dispersion resulted from the numerical mehtod, at the same time, the higher order conventional finite difference (CFD) schemes will give little contribution to avoid the numerical solutions error accumulating with time. To get the similar accuracy with the fourth order center finite difference method, the half truncation length of SFFT should be no less than 7. The FDFCT method can present with the numerical solutions without obvious dispersion when the paprameters of FCT is suitable (we think they should be in the scope from 0.0001 to 0.07). Fortunately, the NADM method not only can reported us with the higher order accuracy solutions (higher than that of the fourth order finite difference method and lower than that of the sixth order finite difference method), but also can distinctly reduce the numerical dispersion. Thirdly, basing on the numerial and theoretical analysis, we reported such nonlinear response accumulating with time as waveform aberration, harmonic generation and resonant peak shift shown by the propagation of one- and two-dimensional non-linear elasticwaves in this paper. And then, we drew the conclusion that these nonlinear responses are controlled by the product between nonlinear strength (SN) and the amplitude of the source. At last, the modified FDFCT numerical method presented by this paper is used to model the two-dimensional non-linear elastic waves propagating in VTI media. Subsequently, the wavelet analysis and polarization are adopted to investigate and understand the numerical results. And then, we found the following principles (attention: the nonlinear strength presented by this paper is weak, the thickness of the -nonlinear media is thin (200m), the initial energy of the source is weak and the anisotropy of the media is weak too): The non-linear response shown by the elastic waves in VTI media is anisotropic too; The instantaneous main frequency sections of seismic records resulted from the media with a non-linear layer have about 1/4 to 1/2 changes of the initial main frequency of source with that resulted from the media without non-linear layer; The responses shown by the elasic waves about the anisotropy and nonlinearity have obvious mutual reformation, namely, the non-linear response will be stronger in some directions because of the anisotropy and the anisotropic strength shown by the elastic waves will be stronger when the media is nonlinear.

Far-field analysis of axially symmetric three-dimensional directional cloaks.

Axisymmetric radiating and scattering structures whose rotational invariance is broken by non-axisymmetric excitations present an important class of problems in electromagnetics. For such problems, a cylindrical wave decomposition formalism can be used to efficiently obtain numerical solutions to the full-wave frequency-domain problem. Often, the far-field, or Fraunhofer region is of particular interest in scattering cross-section and radiation pattern calculations; yet, it is usually impractical to compute full-wave solutions for this region. Here, we propose a generalization of the Stratton-Chu far-field integral adapted for 2.5D formalism. The integration over a closed, axially symmetric surface is analytically reduced to a line integral on a meridional plane. We benchmark this computational technique by comparing it with analytical Mie solutions for a plasmonic nanoparticle, and apply it to the design of a three-dimensional polarization-insensitive cloak.

A defect equation approach for the coupling of subdomains in domain decomposition methods

A defect equation for the coupling of nonlinear subproblems defined in nonoverlapped subdomains arise in domain decomposition methods is presented. Numerical solutions of defect equations by means of quasi-Newton methods are considered.

Performance evaluation of a distributed algorithm for an inverse heat conduction problem

Numerical solutions of realistic 2-D and 3-D inverse problems may require a very large amount of computation. A two-level concept on parallelism is often used to solve such problems. The primary level uses the problem partitioning concept which is a decomposition based on the mathematical/physical problem. The secondary level utilizes the widely used data partitioning concept. A theoretical performance model is built based on the two-level parallelism. The observed performance results obtained from a network of general purpose Sun Sparc stations are compared with the theoretical values. Restrictions of the theoretical model are also discussed.

Comparison of encapsulant curing with convention and microwave systems

Comparison of the performance of a conventional convection oven system with a dual-section microwave system for curing thermosetting polymer encapsulant materials has been performed numerically. A numerical model capable of analysing both the convection and microwave cure processes has been developed and is breifly outliines. The model is used to analyse the curing of a commercially available encapsulant material using both systems. Results obtained from numerical solutions are presented, confirming that the VFM system enables the cure process to be carried out far more rapidly than with the convection oven system. This capability stems from the fundamental heating processes involved, namely that microwave processing enables the heating rate to be varied independently of the material temperature. Variations in cure times, curing rates, maximum temperatures and residual stresses between the processes are fully discussed.

High-energy cutoff in the spectrum of strong-field nonsequential double ionization

Electron energy distributions of singly and doubly ionized helium in an intense 390 nm laser field have been measured at two intensities (0.8 PW/cm(2) and 1.1 PW/cm(2), where PW equivalent to 10(15) W/cm(2)). Numerical solutions of the full-dimensional time-dependent helium Schrodinger equation show excellent agreement with the experimental measurements. The high-energy portion of the two-electron energy distributions reveals an unexpected 5U(p) cutoff for the double ionization (DI) process and leads to a proposed model for DI below the quasiclassical threshold.

Surface plasmons in layered structures with semiconductor and metallic films

The complete spectrum of eigenwaves including surface plasmon polaritons (SPP), dynamic (bulk) and complex waves in the layered structures containing semiconductor and metallic films has been explored. The effects of loss, geometry and the parameters of dielectric layers on the eigenmode spectrum and, particularly, on the SPP modes have been analysed using both the asymptotic and rigorous numerical solutions of the full-wave dispersion equation. The field and Poynting vector distributions have been examined to identify the modes and elucidate their properties. It has been shown that losses and dispersion of permittivity qualitatively alter the spectral content and the eigenwave properties. The SPP counter-directional power fluxes in the film and surrounding dielectrics have been attributed to vortices of power flow, which are responsible for the distinctive features of SPP modes. It has been demonstrated for the first time that the maximal attainable slow-wave factor of the SPP modes guided by thin Au films at optical frequencies is capped not by losses but the frequency dispersion of the actual Au permittivity. © 2009 EDP Sciences.

Control of ionization and dissociation of H2+ by elliptically polarized ultra-short VUV laser pulses

Resonance-enhanced multiphoton ionization of H₂ ⁺ exposed to elliptically polarized VUV laser pulses is investigated. Differential cross sections for nuclei and electron are obtained using numerical solutions of the time-dependent Schrödinger equation. In this work in progress, we explore the dependence of the dissociative ionization observables with the polarization of the light.

Thermally driven surface winds in the tropics

Thesis (Ph. D.)--University of Washington, 1998

Grid structure impact in sparse point representation of derivatives

In the Sparse Point Representation (SPR) method the principle is to retain the function data indicated by significant interpolatory wavelet coefficients, which are defined as interpolation errors by means of an interpolating subdivision scheme. Typically, a SPR grid is coarse in smooth regions, and refined close to irregularities. Furthermore, the computation of partial derivatives of a function from the information of its SPR content is performed in two steps. The first one is a refinement procedure to extend the SPR by the inclusion of new interpolated point values in a security zone. Then, for points in the refined grid, such derivatives are approximated by uniform finite differences, using a step size proportional to each point local scale. If required neighboring stencils are not present in the grid, the corresponding missing point values are approximated from coarser scales using the interpolating subdivision scheme. Using the cubic interpolation subdivision scheme, we demonstrate that such adaptive finite differences can be formulated in terms of a collocation scheme based on the wavelet expansion associated to the SPR. For this purpose, we prove some results concerning the local behavior of such wavelet reconstruction operators, which stand for SPR grids having appropriate structures. This statement implies that the adaptive finite difference scheme and the one using the step size of the finest level produce the same result at SPR grid points. Consequently, in addition to the refinement strategy, our analysis indicates that some care must be taken concerning the grid structure, in order to keep the truncation error under a certain accuracy limit. Illustrating results are presented for 2D Maxwell's equation numerical solutions.

Reinforcement Learning Solution for Unit Commitment Problem Considering Minimum Start Up and Shut Down Times

Unit Commitment Problem (UCP) in power system refers to the problem of determining the on/ off status of generating units that minimize the operating cost during a given time horizon. Since various system and generation constraints are to be satisfied while finding the optimum schedule, UCP turns to be a constrained optimization problem in power system scheduling. Numerical solutions developed are limited for small systems and heuristic methodologies find difficulty in handling stochastic cost functions associated with practical systems. This paper models Unit Commitment as a multi stage decision making task and an efficient Reinforcement Learning solution is formulated considering minimum up time /down time constraints. The correctness and efficiency of the developed solutions are verified for standard test systems

Reinforcement Learning solution for Unit Commitment Problem through pursuit method

Unit commitment is an optimization task in electric power generation control sector. It involves scheduling the ON/OFF status of the generating units to meet the load demand with minimum generation cost satisfying the different constraints existing in the system. Numerical solutions developed are limited for small systems and heuristic methodologies find difficulty in handling stochastic cost functions associated with practical systems. This paper models Unit Commitment as a multi stage decision task and Reinforcement Learning solution is formulated through one efficient exploration strategy: Pursuit method. The correctness and efficiency of the developed solutions are verified for standard test systems