Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Study of turbulent heat transfer of aviation kerosene flows in a curved pipe at supercritical pressure

The now and heat transfer characteristics of China No. 3 aviation kerosene in a heated curved tube under supercritical pressure are numerically investigated by a finite volume method. A two-layer turbulence model, consisting of the RNG k-epsilon two-equation model and the Wolfstein one-equation model, is used for the simulation of turbulence. A 10-species kerosene surrogate model and the NIST Supertrapp software are applied to obtain the thermophysical and transport properties of the kerosene at various temperature under a supercritical pressure of 4 MPa. The large variation of thermophysical properties of the kerosene at the supercritical pressure make the flow and heat transfer more complicated, especially under the effects of buoyancy and centrifugal force. The centrifugal force enhances the heat transfer, but also increases the friction factors. The rise of the velocity caused by the variation of the density does not enhance the effects of the centrifugal force when the curvature ratios are less than 0.05. On the contrary, the variation of the density increases the effects of the buoyancy. (C) 2010 Elsevier Ltd. All rights reserved.

Characteristic rainfall for warning of debris flows

A characteristic rainfall is introduced to overcome the difficulties encountered in determining a critical rainfall value for triggering debris flow. The characteristic value is defined as the rainfall at which debris-flow occurrence probability shows a rapid increase, and can be used as a warning rainfall threshold for debris flows. Investigation of recorded debris flows and 24-hour rainfall data at Jiangjia basin, Yunnan Province, in southwestern China, demonstrates the existence of such a characteristic rainfall. It was found that the characteristic rainfall corresponds to the daily rainfall of 90% cumulative probability by analyzing the basin's daily rainfall histogram. The result provides a simple and useful method for estimating a debris-flow warning rainfall threshold from the daily rainfall distribution. It was applied to estimate the debris-flow warning rainfall threshold for the Subaohe basin, a watershed in the 2008 Wenchuan earthquake zone with many physical characteristics similar to those of the Jiangjia basin.

Probability of Sediment Incipient Motion Under Complex Flows

Presented in this paper is a mathematical model to calculate the probability of the sediment incipient motion, in which the effects of the fluctuating pressure and the seepage are considered. The instantaneous bed shear velocity and the pressure gradient on the bed downstream of the backward-facing step flow are obtained according to the PIV measurements. It is found that the instantaneous pressure gradient on the bed obeys normal distribution. The probability of the sediment incipient motion on the bed downstream of the backward-facing step flow is given by the mathematical model. The predicted results agree well with the experiment in the region downstream of the reattachment point while a large discrepancy between the theory and experiment is seen in the region near the reattachment point. The possible reasons for this discrepancy are discussed.

A finite-element algorithm for modeling variably saturated flows

A general numerical algorithm in the context of finite element scheme is developed to solve Richards’ equation, in which a mass-conservative, modified head based scheme (MHB) is proposed to approximate the governing equation, and mass-lumping techniques are used to keep the numerical simulation stable. The MHB scheme is compared with the modified Picard iteration scheme (MPI) in a ponding infiltration example. Although the MHB scheme is a little inferior to the MPI scheme in respect of mass balance, it is superior in convergence character and simplicity. Fully implicit, explicit and geometric average conductivity methods are performed and compared, the first one is superior in simulation accuracy and can use large time-step size, but the others are superior in iteration efficiency. The algorithm works well over a wide variety of problems, such as infiltration fronts, steady-state and transient water tables, and transient seepage faces, as demonstrated by its performance against published experimental data. The algorithm is presented in sufficient detail to facilitate its implementation.

Triton-He-3 relative and differential flows as probes of the nuclear symmetry energy at supra-saturation densities

Using a transport model coupled with a phase-space coalescence afterburner, we study the triton-He-3 (t-He-3) ratio with both relative and differential transverse flows in semicentral Sn-132 + Sn-124 reactions at a beam energy of 400 MeV/nucleon. The neutron-proton ratios with relative and differential flows are also discussed as a reference. We find that similar to the neutron-proton pairs, the t-He-3 pairs also carry interesting information regarding the density dependence of the nuclear symmetry energy. Moreover, the nuclear symmetry energy affects more strongly the t-He-3 relative and differential flows than the pi(-)/pi(+) ratio in the same reaction. The t-He-3 relative flow can be used as a particularly powerful probe of the high-density behavior of the nuclear symmetry energy.

TRITON-(HE)-H-3 RELATIVE AND DIFFERENTIAL FLOWS AND THE HIGH DENSITY BEHAVIOR OF NUCLEAR SYMMETRY ENERGY

Using a transport model coupled with a phase-space coalescence after-burner we study the triton-He-3 relative and differential transverse flows in semi-central Sn-132 + Sn-124 reactions at a beam energy of 400 MeV/nucleon. We find that the triton-He-3 pairs carry interesting information about the density dependence of the nuclear symmetry energy. The t-He-3 relative flow can be used as a particularly powerful probe of the high-density behavior of the nuclear symmetry energy.

Pressure drops of single and two-phase flows through T-type microchannel mixers

In this paper, preliminary experimental results are presented on pressure drop characteristics of single and two-phase flows through two T-type rectangular microchannel mixers with hydraulic diameters of 528 and 333 mum, respectively. It is shown that both N-2 and water single-phase laminar flows in microchannels, with consideration of experimental uncertainties, are consistent with classic theory, if additional effects, such as entrance effects that will interfere with the interpretation of experimental results, are eliminated by carefully designing the experiments. The obtained pressure drop data of N-2-water two-phase flow in micromixers are analyzed and compared with existing flow pattern-independent models. It is found that the Lockhart-Martinelli method generally underpredicts the frictional pressure drop. Thereafter, a modified correlation of C value in the Chisholm's equation based on linear regression of experimental data is proposed to provide a better prediction of the two-phase frictional pressure drop. Also among the homogeneous flow models investigated, the viscosity correlation of McAdams indicates the best performance in correlating the frictional pressure drop data (mean deviations within +/-20% for two micromixers both). Finally it is suggested that systematic studies are still required to accurately predict two-phase frictional performance in microchannels. (C) 2004 Elsevier B.V. All rights reserved.

Spatial and temporal flows of China's forest resources: Development of a framework for evaluating resource efficiency

Science & Technology Basic Work Program of China: Scientific Survey of the Middle-lower Reaches of Lantsang River and the Great Shangri-La Region [2008FY110300]; National Basic Research Program of China (973 Program): Ecosystem Services and Ecological Safety of the Major Terrestrial Ecosystems of China [2009CB421106]; National Natural Science Foundation of China [30670374]; EU ; European Commission, DG Research [003874]

Strongly nonlinear long interfacial waves in uniform basic flows and their dispersion relationship

In this paper, long interfacial waves of finite amplitude in uniform basic flows are considered with the assumption that the aspect ratio between wavelength and water depth is small. A new model is derived using the velocities at arbitrary distances from the still water level as the velocity variables instead of the commonly used depth-averaged velocities. This significantly improves the dispersion properties and makes them applicable to a wider range of water depths. Since its derivation requires no assumption on wave amplitude, the model thus can be used to describe waves with arbitrary amplitude.

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.