A hybrid coupled ocean-atmosphere model and ENSO prediction study

A hybrid coupled ocean-atmosphere model is designed, which consists of a global AGCM and a simple anomaly ocean model in the tropical Pacific. Retroactive experimental predictions initiated in each season from 1979 to 1994 are performed. Analyses indicate that: (1) The overall predictive capability of this model for SSTA over the central-eastern tropical Pacific can reach one year, and the error is not larger than 0.8 degrees C. (2) The prediction skill depends greatly on the season when forecasts start. However, the phenomenon of SPB (spring prediction barrier) is not found in the model. (3) The ensemble forecast method can effectively improve prediction results. A new initialization scheme is discussed.

Submesoscale activity over the shelf of the northern South China Sea in summer: simulation with an embedded model

We applied a primitive equation ocean model to simulate submesoscale activities and processes over the shelf of the northern South China Sea (NSCS) with a one-way nesting technology for downscaling. The temperature and density fields showed that submesoscale activities were ubiquitous in the NSCS shelf. The vertical velocity was considerably enhanced in submesoscale processes and could reach an average of 58 m per day in the subsurface. At this point, the mixed layer depth also was deepened along the front, and the surface kinetic energy also increased with the intense vertical movement induced by submesoscale activity. Thus, submesoscale stirring/mixing is important for tracers, such as temperature, salinity, nutrients, dissolved organic, and inorganic carbon. This result may have implication for climate and biogeochemical investigations.

Dust-Cloud Structures behind a Shock Wave Moving Over a Deposited Layer of Fine Particles

The present paper investigates dispersed-phase flow structures of a dust cloud induced by a normal shock wave moving at a constant speed over a flat surface deposited with fine particles. In the shock-fitted coordinates, the general equations of dusty-gas

A coupling model for terrestrial processes in arid areas and its application

In this paper, the importance of investigation on terrestrical processes in arid areas for mankind's living environment protection and local economy development as well as its present state of the art are elucidated. A coupling model, which evaluates heat, mass, momentum and radiative fluxes in the SPAC system, is developed for simulating microclimate over plant and bare soil. Especially, it is focussed on the details of turbulence transfer. For illustration, numerical simulation of the water-heat exchange processes at Shapotou Observatory, GAS, Ninxia Province are conducted, and the computational results show that the laws of land-surface processes are rather typical in the arid areas.

Shock-Induced Near-Wall Two-Phase Flow Structure Over a Micron-Sized Particles Bed

The present paper studies numerical modelling of near-wall two-phase flows induced by a normal shock wave moving at a constant speed, over a micronsized particles bed. In this two-fluid model, the possibility of particle trajectory intersection is considered and a full Lagrangian formulation of the dispersed phase is introduced. The finiteness of the Reynolds and Mach numbers of the flow around a particle as well as the fineness of the particle sizes are taken into account in describing the interactions between the carrier- and dispersed- phases. For the small mass-loading ratio case, the numerical simulation of flow structure of the two phases is implemented and the profiles of the particle number density are obtained under the constant-flux condition on the wall. The effects of the shock Mach number and the particle size and material density on particle entrainment motion are discussed in detail.The obtained results indicate that interphase non-equilibrium in the velocity and temperature is a common feature for this type of flows and a local particle accumulation zone may form near the envelope of the particle trajectory family.

Investigation of Vaporized Kerosene Injection and Combustion in a Supersonic Model Combustor

Injection and combustion of vaporized kerosene was experimentally investigated in a Mach 2.5 model combustor at various fuel temperatures and injection pressures. A unique kerosene heating and delivery system, which can prepare heated kerosene up to 820 K at a pressure of 5.5 MPa with negligible fuel coking, was developed. A three-species surrogate was employed to simulate the thermophysical properties of kerosene. The calculated thermophysical properties of surrogate provided insight into the fuel flow control in experiments. Kerosene jet structures at various preheat temperatures injecting into both quiescent environment and a Mach 2.5 crossflow were characterized. It was shown that the use ofvaporized kerosene injection holds the potential of enhancing fuel-air mixing and promoting overall burning. Supersonic combustion tests further confirmed the preceding conjecture by comparing the combustor performances of supercritical kerosene with those of liquid kerosene and effervescent atomization with hydrogen barbotage. Under the similar flow conditions and overall kerosene equivalence ratios, experimental results illustrated that the combustion efficiency of supercritical kerosene increased approximately 10-15% over that of liquid kerosene, which was comparable to that of effervescent atomization.

A Numerical Study for Turbulent Flow and Thermal Influence Over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier-Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052-2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349-377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369-398)'s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243-258)'s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.

Chemical kinetics and self-ignition in a model supersonic hydrogen-air combustor

Self-ignition tests of a model scramjet combustor were conducted by using parallel sonic injection of gaseous hydrogen from the base of a blade-like strut into a supersonic vitiated airstream. The range of stagnation pressure and temperature studied varied from 1.0 to 4.5 MPa and from 1300 to 2200 K, respectively. Experimental results show that the self-ignition limit, in terms of either global or local quantities of pressure and temperature, exhibits a nonmonotonic behavior resembling the classical homogeneous explosion limit of the hydrogen-oxygen system. Specifically, for a given temperature, increasing pressure from a low value can render a nonignitable mixture to first become ignitable, then nonignitable again, This correspondence shows that, despite the globally supersonic nonpremixed configuration studied herein, ignition is strongly influenced by the intricate chemical reaction mechanism and thereby exhibits the homogeneous explosion character. Consequently, self-ignition criteria based on a global reaction rate approximating the complex chemistry are inadequate. An auxiliary computational study on counterflow ignition was also conducted to systematically investigate the contamination effects of vitiated air. Results indicate that the net contamination effects for the present experimental data are expected to be substantially smaller than contributions from the individual contamination species because of the counterbalancing influences of the H2O-inhibition and NO-promotion reactions in effecting ignition.

Combustion and Ignition of Thermally Cracked Kerosene in Supersonic Model Combustors

A series of experiments were conducted to characterize the self-ignition and combustion of thermally cracked kerosene in both a Mach 2.5 model combustor with a combustor entrance height of 51 mm and a Mach 3.0 model combustor with an entrance height of 70 mm. A unique kerosene heating and delivery system was developed, which can prepare heated kerosene up to 950 K at a pressure of 5.5 MPa with negligible fuel coking. The extent of China no. 3 kerosene conversion under supercritical conditions was measured using a specially designed system. The compositions of gaseous products as a result of thermal cracking were analyzed using gas chromatography. The mass flow rates of cracked kerosene were also calibrated and measured using sonic nozzles. With the injection of thermally cracked kerosene, the ability to achieve enhanced combustion performance was demonstrated under a variety of airflow and fuel conditions. Furthermore, self-ignition tests of cracked kerosene in a Mach 2.5 model combustor over a range of fuel injection conditions and with the help of different amounts of pilot hydrogen were conducted and discussed.

Turbulent Modeling of ABL over Wavy Water Surface

A new model accounting for both turbulence and sea state effects has been proposed in this paper to describe momentum, heat and moisture exchanges through air-sea interface. While long wave components mainly change air flow profile, short wave components

Study Of The Damage-Induced Anisotropy Of Quasi-Brittle Materials Using The Component Assembling Model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

High Aspect Ratio (L/D) Riser Viv Prediction Using Wake Oscillator Model

A two-dimensional (2-D) vortex-induced vibration (VIV) prediction model for high aspect ratio (LID) riser subjected to uniform and sheared flow is studied in this paper. The nonlinear structure equations are considered. The near wake dynamics describing the fluctuating nature of vortex shedding is modeled using classical van der Pol equation. A new approach was applied to calibrate the empirical parameters in the wake oscillator model. Compared the predicted results with the experimental data and computational fluid dynamic (CFD) results. Good agreements are observed. It can be concluded that the present model can be used as simple computational tool in predicting some aspects of VIV of long flexible structures. (C) 2008 Elsevier Ltd. All rights reserved.

A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid

A numerical model for shallow-water equations has been built and tested on the Yin-Yang overset spherical grid. A high-order multimoment finite-volume method is used for the spatial discretization in which two kinds of so-called moments of the physical field [i.e., the volume integrated average ( VIA) and the point value (PV)] are treated as the model variables and updated separately in time. In the present model, the PV is computed by the semi-implicit semi-Lagrangian formulation, whereas the VIA is predicted in time via a flux-based finite-volume method and is numerically conserved on each component grid. The concept of including an extra moment (i.e., the volume-integrated value) to enforce the numerical conservativeness provides a general methodology and applies to the existing semi-implicit semi-Lagrangian formulations. Based on both VIA and PV, the high-order interpolation reconstruction can only be done over a single grid cell, which then minimizes the overlapping zone between the Yin and Yang components and effectively reduces the numerical errors introduced in the interpolation required to communicate the data between the two components. The present model completely gets around the singularity and grid convergence in the polar regions of the conventional longitude-latitude grid. Being an issue demanding further investigation, the high-order interpolation across the overlapping region of the Yin-Yang grid in the current model does not rigorously guarantee the numerical conservativeness. Nevertheless, these numerical tests show that the global conservation error in the present model is negligibly small. The model has competitive accuracy and efficiency.

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.

Quantum Monte Carlo study of the CO interaction with a dimer model surface for Cr(110)

The chemisorption of CO on a Cr( 110) surface is investigated using the quantum Monte Carlo method in the diffusion Monte Carlo (DMC) variant and a model Cr2CO cluster. The present results are consistent with the earlier ab initio HF study with this model that showed the tilted/ near-parallel orientation as energetically favoured over the perpendicular arrangement. The DMC energy difference between the two orientations is larger (1.9 eV) than that computed in the previous study. The distribution and reorganization of electrons during CO adsorption on the model surface are analysed using the topological electron localization function method that yields electron populations, charge transfer and clear insight on the chemical bonding that occurs with CO adsorption and dissociation on the model surface.