A non-equilibrium sediment transport model for rill erosion

Sediment transport in rill flows exhibits the characteristics of non-equilibrium transport, and the sediment transport rate of rill flow gradually recovers along the flow direction by erosion. By employing the concept of partial equilibrium sediment transport from open channel hydraulics, a dynamic model of rill erosion on hillslopes was developed. In the model, a parameter, called the restoration coefficient of sediment transport capacity, was used to express the recovery process of sediment transport rate, which was analysed by dimensional analysis and determined from laboratory experimental data. The values of soil loss simulated by the model were in agreement with observed values. The model results showed that the length and gradient of the hillslope and rainfall intensity had different influences on rill erosion. Copyright (c) 2006 John Wiley & Sons, Ltd.

Kinetics and transport model for the chemical vapor epitaxy of GexSi1-x

A numerical model that combines mass transport and surface kinetics was applied, for the first time, to the chemical vapor epitaxy of GexSi1-x. The temperature, velocity and concentration fields were calculated from the conservation equations for energy, momentum and species coupled with the boundary conditions on the growth surface which were determined by surface kinetics. The deposition rates of Si and Ge were assumed to be limited, respectively, by surface kinetics and mass transport. A theoretical relation between the initial conditions and the Ge composition in the solid was established. The calculated growth rate as well as the Ge composition in the solid and its dependence on growth temperature agree well with experimental data.

Circumstantial Evidence for a Soft Nuclear Symmetry Energy at Suprasaturation Densities

Within an isospin- and momentum-dependent hadronic transport model, it is shown that the recent FOPI data on the pi(-)/pi(+) ratio in central heavy-ion collisions at SIS/GSI energies [Willy Reisdorf , Nucl. Phys. A 781, 459 (2007)] provide circumstantial evidence suggesting a rather soft nuclear symmetry energy E-sym(rho) at rho >= 2 rho(0) compared to the Akmal-Pandharipande-Ravenhall prediction. Some astrophysical implications and the need for further experimental confirmations are discussed.

Collective Transverse Flow in Intermediate Energy Heavy-Ion Collisions

Within the hadronic transport model IBUU04, we study the density-dependent symmetry energy by using the neutron-proton differential flow from the Sn-132+Sn-124 reactions at beam energies of 200, 400, 600 and 800MeV per nucleon. The strong effect of the symmetry energy is shown at the incident beam energy of 400 MeV/A. The small medium-effect of the neutron-proton differential flow is also found. We also study the neutron-proton differential flows with impact parameters of 3, 5, 7 fm. It is found that in semi-central collisions the sensitivity of the neutron-proton differential flow to the symmetry energy is larger.

Entrance-channel dependence of isospin effects on double n/p ratio in HIC

Within the hadronic transport model IBUU04, we investigate the effect of density-dependent symmetry energy on double neutron/proton (n/p) ratio of free nucleons in heavy ion collisions by taking four isotopic Sn+Sn reaction systems. Especially the entrance-channel asymmetry and impact-parameter dependence of the effect of symmetry energy are discussed. It is found that in both central and semi-central collisions the sensitivity of the double n/p ratio to the density-dependent symmetry energy is more pronounced in neutron-richer systems. Our results also indicate clearly that the effect of symmetry energy is stronger in central collisions than that in semi-central collisions.

Effect of entrance-channel asymmetry on the isospin dependence of nucleon emission in heavy ion collisions

Using the isospin- and momentum-dependent hadronic transport model 1BUU04, we have investigated the influence of the entrance-channel isospin asymmetry on the sensitivity of the pre-equilibrium neutron/proton ratio to symmetry energy in central heavy-ion collisions induced by high-energy radioactive beams. Our analysis and discussion are based on the dynamical simulations of the three isotopic reaction Systems Sn-132+Sn-124, Sn-124+Sn-112 and Sn-112+(112)Su which are of the same total proton number but, different isospin asymmetry. We find that, the kinetic-energy distributions of the pre-equilibrium neutron/proton ratio are quite sensitive to the density-dependence of symmetry energy at incident beam energy E/A = 400 MeV, and the sensitivity increases as the isospin asymmetry of the reaction system increases.

Double neutron-proton differential transverse flow as a probe for the high density behavior of the nuclear symmetry energy

The double neutron-proton differential transverse flow taken from two reaction systems using different isotopes of the same element is studied at incident beam energies of 400 and 800 MeV/nucleon within the framework of an isospin- and momentum-dependent hadronic transport model IBUU04. The double differential flow is found to retain about the same sensitivity to the density dependence of the nuclear symmetry energy as the single differential flow in the more neutron-rich reaction. Because the double differential flow reduces significantly both the systematic errors and the influence of the Coulomb force, it is thus more effective probe for the high-density behavior of the nuclear symmetry energy.

Beam energy dependence of the pi(-)/pi(+) ratio

Using the isopin- and momentum-dependant hadronic transport model IBUU04, the effect of symmetry energy on the pi(-)/pi(+) ratio are studied. Our investigations are based on the calculations of the Sn-132+Sn-124 semi-central collisions at beam energies of 400/ A MwV, 600/ A MeV and 800MeV. It is found that both the transverse momentum and kinetic energy distributions of the pi(-)/pi(+) ratio are rather sensitive to the symmetry energy, especially around the Colomb peaks. The position of the coulomb peak is shown to be nearly indrpendant of beam energy. The sesitivity of the pi(-)/pi(+) ratio to the symmetry ebergy decreases as the beam energy increases from 600/ A MeV to 800/ A MeV.

Single and double pi(-)/pi(+) ratios in heavy-ion reactions as probes of the high-density behavior of the nuclear symmetry energy

Based on the isospin- and momentum-dependent hadronic transport model IBUU04, effects of the nuclear symmetry energy on the single and double pi(-)/pi(+) ratios in central reactions of Sn-132+Sn-124 and Sn-112+Sn-112 at a beam energy of 400 MeV/nucleon are studied. It is found that around the Coulomb peak of the single pi(-)/pi(+) ratio the double pi(-)/pi(+) ratio taken from the two isotopic reactions retains about the same sensitivity to the density dependence of nuclear symmetry energy. Because the double pi(-)/pi(+) ratio can significantly reduce the systematic errors, it is thus a more effective probe for the high-density behavior of the nuclear symmetry energy.

System size and beam energy effects on probing the high-density behavior of nuclear symmetry energy with pion ratio

Based on the isospin-and momentum-dependent hadronic transport model IBUU04, we have investigated the pi(-)/pi(+) ratio in the following three reactions: Ca-48+Ca-48, Sn-124 +Sn-124 and Au-197+Au-197 with nearly the same isospin asymmetry but different masses, at the bombarding energies from 0.25 to 0.6 A GeV. It is shown that the sensitivity of probing the E-sym (rho) with pi(-)/pi(+) increases with increasing the system size or decreasing the beam energy, showing a correlation to the degree of isospin fractionation. Therefore, with a given isospin asymmetry, heavier system at energies near the pion threshold is preferential to study the behavior Of nuclear symmetry energy at supra-saturation densities.

The effects of sediment transport on grain-size distributions

Changes in statistics (mean, sorting, and skewness) describing grain-size distributions have long been used to speculate on the direction of sediment transport. We present a simple model whereby the distributions of sediment in transport are related to their source by a sediment transfer function which defines the relative probability that a grain within each particular class interval will be eroded and transported. A variety of empirically derived transfer functions exhibit negatively skewed distributions (on a phi scale). Thus, when a sediment is being eroded, the probability of any grain going into transport increases with diminishing grain size throughout more than half of its size range. This causes the sediment in transport to be finer and more negatively skewed than its source, whereas the remaining sediment (a lag) must become relatively coarser and more positively skewed. Flume experiments show that the distributions of transfer functions change from having a highly negative skewness to being nearly symmetrical (although still negatively skewed) as the energy of the transporting process increases. We call the two extremes low-energy and high-energy transfer functions , respectively. In an expanded sediment-transport model, successive deposits in the direction of transport are related by a combination of two transfer functions. If energy is decreasing and the transfer functions have low-energy distributions, successive deposits will become finer and more negatively skewed. If, however, energy is decreasing, but the initial transfer function has a high-energy distribution, successive deposits will become coarser and more positively skewed. The variance of the distributions of lags, sediment in transport, and successive deposits in the down-current direction must eventually decrease (i.e., the sediments will become better sorted). We demonstrate that it is possible for variance first to increase, but suggest that, in reality, an increasing variance in the direction of transport will seldom be observed, particularly when grain-size distributions are described in phi units. This model describing changes in sediment distributions was tested in a variety of environments where the transport direction was known. The results indicate that the model has real-world validity and can provide a method to predict the directions of sediment transport

Diagnostic experiments for transport mechanisms of suspended sediment discharged from the Yellow River in the Bohai Sea

Five diagnostic experiments with a 3D baroclinic hydrodynamic and sediment transport model ECOMSED in couple with the third generation wave model SWAN and the Grant-Madsen bottom boundary layer model driven by the monthly sediment load of the Yellow River, were conducted to separately diagnose effects of different hydrodynamic factors on transport of suspended sediment discharged from the Yellow River in the Bohai Sea. Both transport and spatio-temporal distribution of suspended sediment concentration in the Bohai Sea were numerially simulated. It could be concluded that suspended sediment discharged from the Yellow River cannot be delivered in long distance under the condition of tidal current. Almost all of sediments from the Yellow River are deposited outside the delta under the condition of wind-driven current, and only very small of them are transported faraway. On the basis of wind forcing, sediments from the Yellow River are mainly transported north-northwestward, and others which are first delivered to the Laizhou Bay are continuously moved northward. An obvious 3D structure characteristic of sediment transport is produced in the wind-driven and tide-induced residual circulation condition. Transport patterns at all layers are generally consistent with circulation structure, but there is apparent deviation between the depth-averaged sediment flux and the circulation structure. The phase of temporal variation of sediment concentration is consistent with that of the bottom shear stress, both of which are proved to have a ten-day cycle in wave and current condition.

Large-Eddy-Simulation of 3-Dimensional Rayleigh-Taylor Instability in Incompressible Fluids

The 3-dimensiqnal incompressible Rayleigh-Taylor instability is numerically studied through the large-eddy-simulation (LES) approach based on the passive scalar transport model. Both the instantaneous velocity and the passive scalar fields excited by sinu