Importance of self-consistency in relativistic continuum random-phase approximation calculations

A fully consistent relativistic continuum random phase approximation (RCRPA) is constructed, where the contribution of the continuum spectrum to nuclear excitations is treated exactly by the single-particle Green's function technique. The full consistency of the calculations is achieved that the same effective Lagrangian is adopted for the ground state and the excited states. The negative energy states in the Dirac sea are also included in the single-particle Green's function in the no-sea approximation. The currents from the vector meson and photon exchanges and the Coulomb interaction in RCRPA are treated exactly. The spin-orbit interaction is included naturally in the relativistic frame. Numerical results of the RCRPA are checked with the constrained relativistic mean-field theory. We study the effects of the inconsistency, particularly the currents and Coulomb interaction in various collective multipole excitations.

Relativistic Continuum Random Phase Approximation and Applications II. Applications

The fully consistent relativistic continuum random phase approximation (RCRPA) has been constructed in the momentum representation in the first part of this paper. In this part we describe the numerical details for solving the Bethe-Salpeter equation. The numerical results are checked by the inverse energy weighted sum rules in the isoscalar giant monopole resonance, which are obtained from the constraint relativistic mean field theory and also calculated with the integration of the RCRPA strengths. Good agreement between the misachieved. We study the effects of the self-consistency violation, particularly the currents and Coulomb interaction to various collective multipole excitations. Using the fully consistent RCRPA method, we investigate the properties of isoscalar and isovector collective multipole excitations for some stable and exotic from light to heavy nuclei. The properties of the resonances, such as the centroid energies and strength distributions are compared with the experimental data as well as with results calculated in other models.

Relativistic Continuum Random Phase Approximation and Applications I. Formalism

A fully consistent relativistic continuum random phase approximation (RCRPA) is constructed in terms of the Green's function technique. In this method the contribution of the continuum spectrum to nuclear excitations is treated exactly by the single particle Green's function, which includes also the negative states in the Dirac sea in the nose aapproximation. The theoretical formalism of RCRPA and numerical details are presented. The single particle Green's function is calculated numerically by a proper product of regular and irregular solutions of the Dirac equation. The numerical details and the formalism of RCRPA in the momentum representation are presented.

Digital simulation in a thin layer spectroelectrochemical cell with minigrid platinum electrode by microregion approximation explicit finite difference method

The microregion approximation explicit finite difference method is used to simulate cyclic voltammetry of an electrochemical reversible system in a three-dimensional thin layer cell with minigrid platinum electrode. The simulated CV curve and potential scan-absorbance curve were in very good accordance with the experimental results, which differed from those at a plate electrode. The influences of sweep rate, thickness of the thin layer, and mesh size on the peak current and peak separation were also studied by numerical analysis, which give some instruction for choosing experimental conditions or designing a thin layer cell. The critical ratio (1.33) of the diffusion path inside the mesh hole and across the thin layer was also obtained. If the ratio is greater than 1.33 by means of reducing the thickness of a thin layer, the electrochemical property will be far away from the thin layer property.

Applications of improved 2-D numerical wave tanks in wave breaking

Song and Banner (2002, henceforth referred to as SB02) used a numerical wave tank (developed by Drimer and Agnon, and further refined by Segre, henceforth referred to as DAS) to study the wave breaking in the deep water, and proposed a dimensionless breaking threshold that based on the behaviour of the wave energy modulation and focusing during the evolution of the wave group. In this paper, two modified DAS models are used to further test the SB02's results, the first one (referred to MDAS1) corrected many integral calculation errors appeared in the DAS code, and the second one (referred to MDAS2) replaced the linear boundary element approximation of DAS into the cubic element on the free surface. Researches show that the results of MDAS1 are the same with those of DAS for the simulations of deep water wave breaking, but, the different values of the wavemaker amplitude, the breaking time and the maximum local average energy growth rate delta(max) for the marginal breaking cases are founded by MDAS2 and MDAS1. However, MDAS2 still satisfies the SB02' s breaking threshold. Furthermore, MDAS1 is utilized to study the marginal breaking case in the intermediate water depth when wave passes over a submerged slope, where the slope is given by 1 : 500, 1 : 300, 1 : 150 or 1 : 100. It is found that the maximum local energy density U increases significantly if the slope becomes steeper, and the delta(max) decreases weakly and increases intensively for the marginal recurrence case and marginal breaking case respectively. SB02's breaking threshold is still valid for the wave passing over a submerged slope gentler than 1 : 100 in the intermediate water depth.

Numerical Analysis of a Dual Polarization Mode-Locked Laser with a Quarter Wave Plate

Dynamical behaviors and frequency characteristics of an active mode-locked laser with a quarter wave plate (QWP) are numerically studied by using a set pf vectorial laser equation. Like a polarization self-modulated laser, a frequency shift of half the cavity mode spacing exists between the eigen-modes in the two neutral axes of QWP. Within the active medium, the symmetric gain and cavity structure maintain the pulse's circular polarization with left-hand and right-hand in turn for each round trip. Once the left-hand or right-hand circularly polarized pulse passes through QWP, its polarization is linear and the polarized direction is in one of the directions of i45o with respect to the neutral axes of QWP. The output components in the directions of i45" from the mirror close to QWP are all linearly polarized with a period of twice the round-trip time.

Numerical Investigation of Richtmyer-Meshkov Instability Driven by Cylindrical Shocks

In this paper, a numerical method with high order accuracy and high resolution was developed to simulate the Richtmyer-Meshkov(RM) instability driven by cylindrical shock waves. Compressible Euler equations in cylindrical coordinate were adopted for the cylindrical geometry and a third order accurate group control scheme was adopted to discretize the equations. Moreover, an adaptive grid technique was developed to refine the grid near the moving interface to improve the resolution of numerical solutions. The results of simulation exhibited the evolution process of RM instability, and the effect of Atwood number was studied. The larger the absolute value of Atwood number, the larger the perturbation amplitude. The nonlinear effect manifests more evidently in cylindrical geometry. The shock reflected from the pole center accelerates the interface for the second time, considerably complicating the interface evolution process, and such phenomena of reshock and secondary shock were studied.

Two-Point Closure Strategy in the Mapping Closure Approximation Approach

A two-point closure strategy in mapping closure approximation (MCA) approach is developed for the evolution of the probability density function (PDF) of a scalar advected by stochastic velocity fields. The MCA approach is based on multipoint statistics. We formulate a MCA modeled system using the one-point PDFs and two-point correlations. The MCA models can describe both the evolution of the PDF shape and the rate at which the PDF evolves.

Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma=6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at free-stream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibilty effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely to occur.

Structural Failure Analysis and Numerical Simulation of Micro-Accelerometers under Impulsive Loading

Micromachined accelerometer is a kind of inertial MEMS devices, which usually operate under intensive impact loading. The reliability of micromachined accelerometers is one of the most important performance indices for their design, manufacture and commer

Numerical Analysis of Rainfall Infiltration in the Slope with a Fracture

With the finite volume method, a 2D numerical model for seepage in unsaturated soil has been established to study the rainfall infiltration in the fractured slope.The result shows that more rain may infiltrate into the slope due to existing fracture and then the pore pressure rises correspondingly. Very probably, it is one of the crucial factors accounting for slope failure. Furthermore a preliminary study has been conducted to investigate the influence of various fracture and rainfall factors such as the depth, width and location of a crack, surface condition, rainfall intensity and duration. Pore pressure and water volumetric content during the transient seepage are carefully examined to reveal the intrinsic mechanism.

Intensity Distribution Design for Laser-Induced Thermal Loading Based on Numerical Simulation

To accomplish laser-induced thermal loading simulation tests for pistons,the Gaussian beam was modulated into multi-circular beam with specific intensity distribution.A reverse method was proposed to design the intensity distribution for the laser-induced thermal loading based on finite element(FE) analysis.Firstly,the FE model is improved by alternating parameters of boundary conditions and thermal-physical properties of piston material in a reasonable range,therefore it can simulate the experimental resul...

Numerical simulation of an adaptive optics system with laser propagation in the atmosphere

A comprehensive model of laser propagation in the atmosphere with a complete adaptive optics (AO) system for phase compensation is presented, and a corresponding computer program is compiled. A direct wave-front gradient control method is used to reconstruct the wave-front phase. With the long-exposure Strehl ratio as the evaluation parameter, a numerical simulation of an AO system in a stationary state with the atmospheric propagation of a laser beam was conducted. It was found that for certain conditions the phase screen that describes turbulence in the atmosphere might not be isotropic. Numerical experiments show that the computational results in imaging of lenses by means of the fast Fourier transform (FFT) method agree well with those computed by means of an integration method. However, the computer time required for the FFT method is 1 order of magnitude less than that of the integration method. Phase tailoring of the calculated phase is presented as a means to solve the problem that variance of the calculated residual phase does not correspond to the correction effectiveness of an AO system. It is found for the first time to our knowledge that for a constant delay time of an AO system, when the lateral wind speed exceeds a threshold, the compensation effectiveness of an AO system is better than that of complete phase conjugation. This finding indicates that the better compensation capability of an AO system does not mean better correction effectiveness. (C) 2000 Optical Society of America.