Shallow Water Model On Cubed-Sphere By Multi-Moment Finite Volume Method

A global numerical model for shallow water flows on the cubed-sphere grid is proposed in this paper. The model is constructed by using the constrained interpolation profile/multi-moment finite volume method (CIP/MM FVM). Two kinds of moments, i.e. the point value (PV) and the volume-integrated average (VIA) are defined and independently updated in the present model by different numerical formulations. The Lax-Friedrichs upwind splitting is used to update the PV moment in terms of a derivative Riemann problem, and a finite volume formulation derived by integrating the governing equations over each mesh element is used to predict the VIA moment. The cubed-sphere grid is applied to get around the polar singularity and to obtain uniform grid spacing for a spherical geometry. Highly localized reconstruction in CIP/MM FVM is well suited for the cubed-sphere grid, especially in dealing with the discontinuity in the coordinates between different patches. The mass conservation is completely achieved over the whole globe. The numerical model has been verified by Williamson's standard test set for shallow water equation model on sphere. The results reveal that the present model is competitive to most existing ones. (C) 2008 Elsevier Inc. All rights reserved.

A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.

A micromechanical model of influence of particle fracture and particle cluster on mechanical properties of metal matrix composites

A general analytical model for a composite with an isotropic matrix and two populations of spherical inclusions is proposed. The method is based on the second order moment of stress for evaluating the homogenised effective stress in the matrix and on the secant moduli concept for the plastic deformation. With Webull's statistical law for the strength of SiCp particles, the model can quantitatively predict the influence of particle fracture on the mechanical properties of PMMCs. Application of the proposed model to the particle cluster shows that the particle cluster has neglected influence on the strain and stress curves of the composite. (C) 1998 Elsevier Science B.V.

Numerical simulation of g-jitter effects on half floating zone convection under microgravity environment

The g-jitter effects on the thermocapillary convection in liquid bridge of floating half zone were studied by numerical simulation for unsteady and axi-symmetric model in the cylindrical coordinate system. The g-jitter field was given by a steady microgravity field in addition to an oscillatory low-gravity field, and the effects on the flow field, temperature distribution and free surface deformation were analyzed numerically.

G-jitter effects on half floating-zone convection in intermediate-frequency range

The g-jitter influence on thermocapillary convection and critical Marangoni number in a liquid bridge of half-floating rone was discussed in the low frequency range of 0.4 to 1.5 Hz in a previous paper. This paper extended the experiments to the intermediate frequency range of 2 to 18 Hz, which htrs often been recorded as vibration environment of spacecrafts. The experiment was completed on the deck of a vibration machine, which gave a periodical applied acceleration to simulate the effects of g-jitter. The experimental results in the intermediate frequency range are different from that in the low frequency range. The velocity field and the shape of the free surface have periodical fluctuations in response to g-jitter. The amplitude of the periodical varying part of the temperature response decreases obviously with increasing frequency of g-jitter and vanishes almost when the frequency of g-jitter is high enough. The critical Marangoni number is defined to describe the transition from a periodical convection in response to g-jitter to an oscillatory convection due to internal instability, and will increase with increasing g-jitter frequency. According to the spectral analysis, it can be found that the oscillatory part of temperature is a superposition of two harmonic waves if the Marangoni number is larger than a critical value.

Effects of g-jitter on the critical marangoni number

A half floating zone is fixed on a vibrational deck, which supports a periodical applied acceleration to simulate the effect of g-jitter. This paper deals with the effects of g-jitter on the fluid fields and the critical Marangoni number, which describes the transition from a forced oscillation of thermocapillary convection into an instability oscillatory convection in a liquid bridge of half floating zone with top rod heated. The responses of g-jitter field on the temperature profiles and flow pattern in the liquid bridge were obtained experimentally. The results indicated that the critical Marangoni number decreases with the increasing of g-jitter effect and is slightly smaller for higher frequency of g-jitter with fixed strength of applied gravity.

Quantum and variational monte-carlo interaction potentials for li2 (x1-sigma-g+)

Optimized trial functions are used in quantum Monte Carlo and variational Monte Carlo calculations of the Li₂(X ¹Σ⁺_g) potential curve. The trial functions used are a product of a Slater determinant of molecular orbitals multiplied by correlation functions of electron—nuclear and electron—electron separation. The parameters of the determinant and correlation functions are optimized simultaneously by reducing the deviations of the local energy E_L (E_L  Ψ⁻¹_THΨ_T, where Ψ_T denotes a trial function) over a fixed sample. At the equilibrium separation, the variational Monte Carlo and quantum Monte Carlo methods recover 68% and 98% of the correlation energy, respectively. At other points on the curves, these methods yield similar accuracies.

源于G.I.Taylor工作的现代经典物理学

一位科学家的工作提供了一门完整课程的素材，涉及从流体动力学稳定性、湍流到流体电动力学、微生物的运动．

Starting Nozzle Flow Simulation Using K-G Two-Equation Turbulence Model

The starting process of two-dimensional nozzle flows has been simulated with Euler, laminar and k - g two-equation turbulence Navier-Stokes equations. The flow solver is based on a combination of LUSGS subiteration implicit method and five spatial discretized schemes, which are Roe, HLLE, MHLLE upwind schemes and AUSM+, AUSMPW schemes. In the paper, special attention is for the flow differences of the nozzle starting process obtained from different governing equations and different schemes. Two nozzle flows, previously investigated experimentally and numerically by other researchers, are chosen as our examples. The calculated results indicate the carbuncle phenomenon and unphysical oscillations appear more or less near a wall or behind strong shock wave except using HLLE scheme, and these unphysical phenomena become more seriously with the increase of Mach number. Comparing the turbulence calculation, inviscid solution cannot simulate the wall flow separation and the laminar solution shows some different flow characteristics in the regions of flow separation and near wall.

Spatio-Temporal Scanning and Statistical Test of the Accelerating Moment Release (AMR) Model Using Australian Earthquake Data

The Accelerating Moment Release (AMR) preceding earthquakes with magnitude above 5 in Australia that occurred during the last 20 years was analyzed to test the Critical Point Hypothesis. Twelve earthquakes in the catalog were chosen based on a criterion for the number of nearby events. Results show that seven sequences with numerous events recorded leading up to the main earthquake exhibited accelerating moment release. Two occurred near in time and space to other earthquakes preceded by AM R. The remaining three sequences had very few events in the catalog so the lack of AMR detected in the analysis may be related to catalog incompleteness. Spatio-temporal scanning of AMR parameters shows that 80% of the areas in which AMR occurred experienced large events. In areas of similar background seismicity with no large events, 10 out of 12 cases exhibit no AMR, and two others are false alarms where AMR was observed but no large event followed. The relationship between AMR and Load-Unload Response Ratio (LURR) was studied. Both methods predict similar critical region sizes, however, the critical point time using AMR is slightly earlier than the time of the critical point LURR anomaly.

液氮温区G-M制冷机的数值分析及新型结构制冷机的研制

G-M制冷机是回热式的小型低温制冷机，它利用绝热放气膨胀原理（又称为西蒙膨胀法）获得低温，具有振动小、运行稳定、寿命长、操作方便等特点，在对效率、重量、尺寸等没有太高要求的场合应用非常广泛。在八十年代末，G-M制冷机已经突破了液氦温度，非常适合在液氦温区为超导器件或电子元件提供冷量，在低温真空泵、MRI超导磁体冷却系统、SQUID再冷凝系统等方面有良好的应用前景，推广其应用已属当务之急。本文围绕着探索氦温区G-M制冷机工作机理、提高低温蓄冷器性能及新型结构G-M制冷机的研制等方面，进行了系统地理论分析和纳归纳，以及初步地实验研究：一. 首次比较完善地建立了液氦温区G-M型制冷机整机数值模拟方法，数值模拟方法给出了制冷机内工质参数瞬态分布及动态变化曲线，为分析制冷机独特的循环特性提供了直观的依据，为探讨运行及结构参数对制冷性能的影响机理提供了强有力的工具。模型中考虑了制冷机中的阻力、实际进排气角、物性变化及蓄冷器内空隙率的存在等多种实际因素：解决了一、二级耦合及部件交界处物性变化不连续对计算影响等问题；采取网格非均匀化、牛顿迭代法以及负反馈原理等措施，提高了计算精度和收敛性。在微机上成功地模拟了液氦温区G-M制冷机的工作过程，程序采用模块化编程，具有一定通用性。二. 运用液氦温区整套机数值模拟方法，计算了制冷机内工质参数（氦流温度、压力、流速等）周期性变化和空间分布，采用一级蓄冷器与低温蓄冷器工作特性对比法，从整机内工质参数动态变化规律分析的角度，验证了液氦温区G-M制冷机工质氦的循环主要分为两部分：常规循环、低温循环。并详细地讨论了运行参数（频率、工作压力）及一级制冷机结构参数对整机内工数动态变化的影响规律以及制冷机性能的影响机理。三. 首次建立了较完善的液氦温区多层混合填料型低温蓄冷器的数值模拟方法，运用数值模拟方法，首次详细地研究了常用填料的不同组合、一定组合下填料比例以及蓄冷器结构对制冷机性能的影响机理，提出了不同填料的最佳节组合及一定组合下填料最佳比例的确定和低温蓄冷器结构和填料优化的原则，为合理有效地设计高性能液氦温区低温蓄冷器提供了依据。四. 提出并设计、加工了一种新型结构的液氦温区双级G-M型制冷机，该机结构在国内外属于首创。其具有以下主要结构特点：一、二级分别独立驱动；一、二级之间通过热桥连接；二级蓄冷器外置于汽缸等。同时建立了新型结构液氦温区双级G-M型制冷机实验系统，为今后整机性能和内部动态过程的研究奠定了基础。五. 新型结构双级G-M型制冷机二级单机动转频率为0.6Hz时，制冷温度为13.6K，且在20K时有4.4W的制冷量；制冷机已达到液氦温区，f = 0.4Hz时，最低温度为4.6K；f = 1Hz，制冷温度为10K时，制冷量大于6W，上述结果目前均未见有文献报道。在新型G-M制冷机上，初步进行了低温蓄冷器性能测试实验及运行参数对制冷机性能影响的实验研究。