Application of the characteristic CIP method to a shallow water model on the sphere

Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the efficiency of the characteristic constrained interpolation profile (CIP) method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.

Study on the Mechanism of Detonation to Quasi-detonation Transition

In this paper, the mechanism of detonation to quasi-detonation transition was discussed, a new physical model to simulate quasi-detonation was proposed, and one-dimensional theoretical and numerical simulation was conducted. This study firstly demonstrates that the quasi-detonation is of thermal choking. If the conditions of thermal choking are created by some disturbances, the supersonic flow is then unable to accept additional thermal energy, and the CJ detonation becomes the unstable quasi-detonation precipitately. The kinetic energy loss caused by this transition process is firstly considered in this new physical model. The numerical results are in good agreement with previous experimental observations qualitatively, which demonstrates that the quasi-detonation model is physically correct and the study are fundamentally important for detonation and supersonic combustion research.

Heat Transfer Problems Induced By Multi-shocks Interaction

Numerical simulations of the multi-shock interactions observable around hypersonic vehicles were carried out by solving Navier-Stokes equations with the AUSMPW scheme and the new type of the IV interaction created by two incident shock waves was investigated in detail. Numerical results show that the intersection point of the second incident shock with the bow shock plays important role on the flow pattern, peak pressures and heat fluxes. In the case of two incident shocks interacting with the bow shock at the same position, the much higher peak pressure and more severe heat transfer rate are induced than the classical IV interaction. The phenomenon is referred to as the multi-shock interaction and higher requirements will be imposed on thermal protection systems.

Thermally induced evolution of phase transformations in gas hydrate sediment

Thermally induced evolution of phase transformations is a basic physical-chemical process in the dissociation of gas hydrate in sediment (GHS). Heat transfer leads to the weakening of the bed soil and the simultaneous establishment of a time varying stress field accompanied by seepage of fluids and deformation of the soil. As a consequence, ground failure could occur causing engineering damage or/and environmental disaster. This paper presents a simplified analysis of the thermal process by assuming that thermal conduction can be decoupled from the flow and deformation process. It is further assumed that phase transformations take place instantaneously. Analytical and numerical results are given for several examples of simplified geometry. Experiments using Tetra-hydro-furan hydrate sediments were carried out in our laboratory to check the theory. By comparison, the theoretical, numerical and experimental results on the evolution of dissociation fronts and temperature in the sediment are found to be in good agreement.

Isoscalar Giant Monopole Resonance in Relativistic Continuum Random Phase Approximation

The isoscalar giant monopole resonance (ISGMR) in nuclei is studied in the framework of a fully consistent relativistic continuum random phase approximation (RCRPA). In this method the contribution of the continuum spectrum to nuclear excitations is treated exactly by the single particle Green's function technique. The negative energy states in the Dirac sea are also included in the single particle Green's function in the no-sea approximation. The single particle Green's function is calculated numerically by a proper product of the regular and irregular solutions of the Dirac equation. The strength distributions in the RCRPA calculations, the inverse energy-weighted sum rule m(-1) and the centroid energy of the ISGMR in Sn-120 and Pb-208 are analysed. Numerical results of the RCRPA are checked with the constrained relativistic mean field model and relativistic random phase approximation with a discretized spectrum in the continuum. Good agreement between them is achieved.

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.

Molecular dynamics and control following excitation with an ultra-short intense laser pulse

Vibronic excitations of the tri-atomic molecule OClO (A(2)A(2)(nu(1), nu(2), nu(3)) <-- (XB1)-B-2 (0, 0, 0)) with weak and strong ultra-short laser fields are studied within full quantum wavepacket dynamics in hyperspherical coordinates. Different dynamics is observed following excitation with laser pulses of different intensities. With a strong laser pulse, many vibrational states are excited and a spatially more localised wavepacket arises. The numerical results show that the population of different vibrational states of the wavepacket on the excited potential energy surface is altered by the intensity of the laser pulse. The numerical results also suggest a related effect on the phase of the wavepacket. These interesting phenomena can be understood by an analysis of the corresponding results for two model diatomic molecules. The possible physical mechanisms of control of chemical processes using strong laser fields are discussed. (C) 2004 Elsevier B.V. All rights reserved.

Studies on coarsening of microdomains in binary polymer mixtures by Lattice Boltzmann methods

An improved free energy approach Lattice Boltzmann model(LBM) is proposed by introducing a forcing term instead of the pressure tensor. This model can reach the proper thermodynamic equilibrium after enough simulation time. On the basis of this model, the phase separation in binary polymer mixtures is studied by applying a Flory-Huggins-type free energy. The numerical results show good agreement with the analytic coexistence curve. This model can also be used to study the coarsening of microdomains in binary polymer mixtures at the early and intermediate stages.

Wave generation by the falling rock in the two-dimensional wave tank

Wave generation by the falling rock in the two-dimensional wave tank is experimentally and numerically studied, where the numerical model utilizes the boundary element method to solve the fully nonlinear potential flow theory. The wave profiles at different times are measured in the laboratory, which are also used to test the numerical model. Comparisons show that the experimental and numerical results are in good agreement, and the numerical model can be used to simulate the wave generation due to the submarine rock falling. Further numerical tests on the influences of the rock size, density, initial position and the falling angle on the wave elevation of the generated waves are performed, respectively. The results show that the size and density of the rock have strong effects on the maximum elevation of the generated wave, while the effects of the initial position and the falling angle of the rock are also significant. When the size or the density of the rock increases, the maximum elevation of the generated wave increases. The same effect on the generated wave would be produced if the initial position of the rock becomes closer to the surface, or the falling angle between the falling route and the vertical direction turns larger. In addition, the present numerical tests reveal that the submarine rock falling provides a new generation method for the breaking wave in the wave tank.

A semi-analytical theory of coastal upwelling and jet

The general forms of the conservation of momentum, temperature and potential vorticity of coastal ocean are obtained in the x-z plane for the nonlinear ocean circulation of Boussinesq fluid, and a elliptic type partial differential equations of second order are derived. Solution of the partial differential equations are obtained under the conditions that the fluid moves along the topography. The numerical results show that there exist both upwelling and downwelling along coastline that mainly depends on the large scale ocean condition. Numerically results of the upwelling (downwelling), coastal jet and temperature front zone are favorable to the observations.

SUSPENDED SEDIMENT TRANSPORT IN THE OFFSHORE NEAR YANGTZE ESTUARY

Based on the Estuarine, Coastal and Ocean Modeling System with Sediments (ECOMSED) model, a 3-D hydrodynamic-transport numerical model was established for the offshore area near the Yangtze Estuary in the East China Sea. The hydrodynamic module was driven by tide and wind. Sediment module included sediment resuspension, transport and deposition of cohesive and non-cohesive sediment. The settling of cohesive sediment in the water column was modeled as a function of aggregation (flocculation) and deposition. The numerical results were compared with observation data for August, 2006. It shows that the sediment concentration reduces gradually from the seashore to the offshore area. Numerical results of concentration time series in the observation stations show two peaks and two valleys, according with the observation data. It is mainly affected by tidal current. The suspended sediment concentration is related to the tidal current during a tidal cycle, and the maximum concentration appears 1 h-4 h after the current maximum velocity has reached.

Algorithm Studies on How to Obtain a Conditional Nonlinear Optimal Perturbation (CNOP)

The conditional nonlinear optimal perturbation (CNOP), which is a nonlinear generalization of the linear singular vector (LSV), is applied in important problems of atmospheric and oceanic sciences, including ENSO predictability, targeted observations, and ensemble forecast. In this study, we investigate the computational cost of obtaining the CNOP by several methods. Differences and similarities, in terms of the computational error and cost in obtaining the CNOP, are compared among the sequential quadratic programming (SQP) algorithm, the limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm, and the spectral projected gradients (SPG2) algorithm. A theoretical grassland ecosystem model and the classical Lorenz model are used as examples. Numerical results demonstrate that the computational error is acceptable with all three algorithms. The computational cost to obtain the CNOP is reduced by using the SQP algorithm. The experimental results also reveal that the L-BFGS algorithm is the most effective algorithm among the three optimization algorithms for obtaining the CNOP. The numerical results suggest a new approach and algorithm for obtaining the CNOP for a large-scale optimization problem.

基于各向异性非线性扩散的图像消抖算法

对于图像抖动产生偏移,提出了一种基于各向异性非线性扩散以及抖动估计的抖动消除算法。这种各向异性非线性扩散的模型由两项组成,即扩散项以及强制项。基本思想就是对于边界点以及图像内部的点分别进行处理。利用Newton-Raphson算法最小化抖动误差,估计出抖动偏移量。实验结果表明本文的抖动消除技术比其他方法的消除性能好,恢复效果接近于理想图像且性能稳定。

复杂条件地层滤波预测的理论与关键技术研究

The theory researches of prediction about stratigraphic filtering in complex condition are carried out, and three key techniques are put forward in this dissertation. Theoretical aspects: The prediction equations for both slant incidence in horizontally layered medium and that in laterally variant velocity medium are expressed appropriately. Solving the equations, the linear prediction operator of overlaid layers, then corresponding reflection/transmission operators, can be obtained. The properties of linear prediction operator are elucidated followed by putting forward the event model for generalized Goupillaud layers. Key technique 1: Spectral factorization is introduced to solve the prediction equations in complex condition and numerical results are illustrated. Key technique 2: So-called large-step wavefield extrapolation of one-way wave under laterally variant velocity circumstance is studied. Based on Lie algebraic integral and structure preserving algorithm, large-step wavefield depth extrapolation scheme is set forth. In this method, the complex phase of wavefield extrapolation operator’s symbol is expressed as a linear combination of wavenumbers with the coefficients of this linear combination in the form of the integral of interval velocity and its derivatives over depth. The exponential transform of the complex phase is implemented through phase shifting, BCH splitting and orthogonal polynomial expansion. The results of numerical test show that large-step scheme takes on a great number of advantages as low accumulating error, cheapness, well adaptability to laterally variant velocity, small dispersive, etc. Key technique 3: Utilizing large-step wavefield extrapolation scheme and based on the idea of local harmonic decomposition, the technique generating angle gathers for 2D case is generalized to 3D case so as to solve the problems generating and storing 3D prestack angle gathers. Shot domain parallel scheme is adopted by which main duty for servant-nodes is to compute trigonometric expansion coefficients, while that for host-node is to reclaim them with which object-oriented angle gathers yield. In theoretical research, many efforts have been made in probing into the traits of uncertainties within macro-dynamic procedures.