On an integral equation arising in the transport of radiation through a slab involving internal reflection

In an earlier contribution to this journal [M.M.R. Williams, Eur. Phys. J. B 47, 291 (2005)], we derived an integral equation for the transmission of radiation through a slab of finite thickness which incorporated internal reflection at the surfaces. Here we generalise the problem to the case when there is a source on each face and the reflection coefficients are different at each face. We also discuss numerical and analytic solutions of the equation discussed in [M.M.R. Williams, Eur. Phys. J. B 47, 291 (2005)] when the reflection is governed by the Fresnel conditions. We obtain numerical and graphical results for the reflection and transmission coefficients, the scalar intensity and current and the emergent angular distributions at each face. The incident source is either a mono-directional beam or a smoothly varying distribution which goes from isotropic to a normal beam. Of particular interest is the philosophy of the numerical solution and whether a direct numerical approach is more effective than one involving a more elegant analytical solution using replication and the Hilbert problem. We also develop the solution of this problem using diffusion theory and compare the results with the exact transport solution. An erratum to this article is available at .     

Nonlinear Flow Numerical Simulation of an Ultra-Low Permeability Reservoir

A nonlinear flow mathematical model is established and the grid equation is deduced. A nonlinear flow reservoir numerical simulation program is compiled. The permeability loss coefficient is used to describe the permeability loss. A pilot calculation is made on the basis of actual field data, which reflects the reservoir development characteristics. The numerical simulation program based on nonlinear flow can anticipate the dynamic characteristics of the ultra-low permeability reservoir exploitation more exactly.     

Coherent Structures in Transition of a Flat-Plate Boundary Layer at Ma = 0.7

Direct numerical simulation （DNS） of a spatially evolving fiat-plate boundary layer transition process at free stream Mach number 0. 7 is performed. Tollmien-Schlichting （T-S） waves are added on the inlet boundary as the disturbances before transition. Typical coherent structures in the transition process are investigated based on the second invariant of velocity gradient tensor. The instantaneous shear stress and the mean velocity profile in the transition region are studied. In our view, the fact that the peak value of shear stress in the stress concentration area increases and exceeds a threshold value during the later stage of the trallsition process plays an important role in the laminar breakdown process.     

Turbulence in diffusion replicator equation

Dynamical behaviors in the diffusion replicator equation of three species are numerically studied. We point out the significant role of the heteroclinic cycle in the equation, and analyze the details of the turbulent solution that appeared in this system. Firstly, the bifurcation diagram for a certain parameter setting is drawn. Then it is shown that the turbulence appears with the supercritical Hopf bifurcation of a stationary uniform solution and it disappears under a subcritical-type bifurcation. Secondly, the statistical property of the turbulence near the supercritical Hopf onset point is analyzed precisely. Further, the correlation lengths and correlation times obey some characteristic scaling laws.     

Directed transient long-range transport in a slowly driven Hamiltonian system of interacting particles

We study the Hamiltonian dynamics of a one-dimensional chain of linearly coupled particles in a spatially periodic potential which is subjected to a time-periodic mono-frequency external field. The average over time and space of the related force vanishes and hence, the system is effectively without bias which excludes any ratchet effect. We pay special attention to the escape of the entire chain when initially all of its units are distributed in a potential well. Moreover for an escaping chain we explore the possibility of the successive generation of a directed flow based on large accelerations. We find that for adiabatic slope-modulations due to the ac-field transient long-range transport dynamics arises whose direction is governed by the initial phase of the modulation. Most strikingly, that for the driven many particle Hamiltonian system directed collective motion is observed provides evidence for the existence of families of transporting invariant tori confining orbits in ballistic channels in the high-dimensional phase spaces.     

Convective heat transport in a fluid layer of infinite Prandtl number: upper bounds for the case of rigid lower boundary and stress-free upper boundary

We present the theory of the multi--solutions of the variational problem for the upper bounds on the convective heat transport in a heated from below horizontal fluid layer with rigid lower boundary and stress-free upper boundary. A sequence of upper bounds on the convective heat transport is obtained. The highest bound is between the bounds for the case of a fluid layer with two rigid boundaries and for the case of a fluid layer with two stress-free boundaries. As an additional result of the presented theory we obtain small corrections of the boundary layer thicknesses of the optimum fields for the case of fluid layer with two rigid boundaries. These corrections lead to systematically lower upper bounds on the convective heat transport in comparison to the bounds obtained in [5]. Received 29 September 1999     

A computationally efficient scheme for the non-linear diffusion equation

This Letter proposes a new numerical scheme for integrating the non-linear diffusion equation. It is shown that it is linearly stable. Some tests are presented comparing this scheme to a popular decentered version of the linearized Crank-Nicholson scheme, showing that, although this scheme is slightly less accurate in treating the highly resolved waves, (i) the new scheme better treats highly non-linear systems, (ii) better handles the short waves, (iii) for a given test bed turns out to be three to four times more computationally cheap, and (iv) is easier in implementation.     

Multi-Bifurcation Effect of Blood Flow by Lattice Boltzmann Method

The multi-bifurcation effect of blood flow is investigated by lattice Boltzmann method at Re=200 with six different bifurcation angles α, which are 22.5°, 25°, 28°, 30°, 33°, 35°, respectively. The velocities and ratios of average velocity at various bifurcations are discussed. It is indicated that the maximum velocity at the section near the first divider increases and shifts towards the walls of branch with the increase of α. At the first bifurcation, the average horizontal velocities increase with the increase of α. The average horizontal velocities of outer branches at the secondary bifurcation decrease at 22.5°≤ α≤30° and increase at 30°≤α≤ 35°, whereas those of inner branches at the secondary bifurcation have the opposite variation, as the same as the above variations of the ratios of average horizontal velocities at various bifurcations. The ratios of average vertical velocities of branch at first bifurcation to that of outer branches at the secondary bifurcation increase at 22.5°≤α≤30° and decrease at 30°≤ α ≤ 35°, whereas the ratios of average vertical velocities of branch at first bifurcation to that of inner branches at the secondary bifurcation always decrease.     

Two Phases of Coherent Structure Motions in Turbulent Boundary Layer

Two phases of coherent structure motion are acquired after obtaining conditional phase-averaged waveforms for longitudinal velocity of coherent structures in turbulent boundary layer based on Harr wavelet transfer. The correspondences of the two phases to the two processes （i.e. ejection and sweep） during a burst are determined.     

Piecewise parabolic method on a local stencil for magnetized supersonic turbulence simulation

Stable, accurate, divergence-free simulation of magnetized supersonic turbulence is a severe test of numerical MHD schemes and has been surprisingly difficult to achieve due to the range of flow conditions present. Here we present a new, higher order-accurate, low dissipation numerical method which requires no additional dissipation or local “fixes” for stable execution. We describe PPML, a local stencil variant of the popular PPM algorithm for solving the equations of compressible ideal magnetohydrodynamics. The principal difference between PPML and PPM is that cell interface states are evolved rather that reconstructed at every timestep, resulting in a compact stencil. Interface states are evolved using Riemann invariants containing all transverse derivative information. The conservation laws are updated in an unsplit fashion, making the scheme fully multidimensional. Divergence-free evolution of the magnetic field is maintained using the higher order-accurate constrained transport technique of Gardiner and Stone. The accuracy and stability of the scheme is documented against a bank of standard test problems drawn from the literature. The method is applied to numerical simulation of supersonic MHD turbulence, which is important for many problems in astrophysics, including star formation in dark molecular clouds. PPML accurately reproduces in three-dimensions a transition to turbulence in highly compressible isothermal gas in a molecular cloud model. The low dissipation and wide spectral bandwidth of this method make it an ideal candidate for direct turbulence simulations.     

Direct Numerical Simulation of the Rayleigh-Taylor Instability with the Spectral Element Method

A novel method is proposed to simulate Rayleigh-Taylor instabilities using a specially-developed unsteady three-dimensional high-order spectral element method code. The numerical model used consists of Navier-Stokes equations and a transport-diffusive equation. The code is first validated with the results of linear stability perturbation theory. Then several characteristics of the Rayleigh-Taylor instabilities are studied using this three-dimensional unsteady code, including instantaneous turbulent structures and statistical turbulent mixing heights under different initial wave numbers. These results indicate that turbulent structures of Rayleigh-Taylor instabilities are strongly dependent on the initial conditions. The results also suggest that a high-order numerical method should provide the capability of simulating small scale fluctuations of Rayleigh-Taylor instabilities of turbulent flows.     

Upper bounds on the convective heat transport in a rotating fluid layer of infinite Prandtl number: case of large Taylor numbers

By means of the Howard-Busse method of the optimum theory of turbulence we obtain upper bounds on the convective heat transport in a horizontal fluid layer heated from below and rotating about a vertical axis. We consider the interval of large Taylor numbers where the intermediate layers of the optimum fields expand in the direction of the corresponding internal layers. We consider the 1 - α-solution of the arising variational problem for the cases of rigid-stress-free, stress-free, and rigid boundary conditions. For each kind of boundary condition we discuss four cases: two cases where the boundary layers are thinner than the Ekman layers of the optimum field and two cases where the boundary layers are thicker than the Ekman layers. In most cases we use an improved solution of the Euler-Lagrange equations of the variational problem for the intermediate layers of the optimum fields. This solution leads to corrections of the thicknesses of the boundary layers of the optimum fields and to lower upper bounds on the convective heat transport in comparison to the bounds obtained by Chan [J. Fluid Mech. 64, 477 (1974)] and Hunter and Riahi [J. Fluid Mech. 72, 433 (1975)]. Compared to the existing experimental data for the case of a fluid layer with rigid boundaries the corresponding upper bounds on the convective heat transport is less than two times larger than the experimental results, the corresponding upper bound on the convective heat transport, obtained by Hunter and Riahi is about 10% higher than the bound obtained in this article. When Rayleigh number and Taylor number are high enough the upper bound on the convective heat transport ceases to depend on the boundary conditions. Received 30 January 2001 and Received in final form 28 May 2001     

Scaling of global momentum transport in Taylor-Couette and pipe flow

We interpret measurements of the Reynolds number dependence of the torque in Taylor-Couette flow by Lewis and Swinney [Phys. Rev. E 59, 5457 (1999)] and of the pressure drop in pipe flow by Smits and Zagarola [Phys. Fluids 10, 1045 (1998)] within the scaling theory of Grossmann and Lohse [J. Fluid Mech. 407, 27 (2000)], developed in the context of thermal convection. The main idea is to split the energy dissipation into contributions from a boundary layer and the turbulent bulk. This ansatz can account for the observed scaling in both cases if it is assumed that the internal wind velocity introduced through the rotational or pressure forcing is related to the external (imposed) velocity U, by with and for the Taylor-Couette (U inner cylinder velocity) and pipe flow (U mean flow velocity) case, respectively. In contrast to the Rayleigh-Bénard case the scaling exponents cannot (yet) be derived from the dynamical equations. Received 9 September 2000     

Vortex Structures and Behavior of a Flow Past Two Rotating Circular Cylinders Arranged Side-by-Side

We present a study on the vortex structures and behavior of a flow past two rotating circular cylinders arranged side-by-side at a range of absolute rotational speeds （｜a｜≤2） for two different gap spacings g＊= 1.5 and 0.7 at Reynolds numbers Re = 160 and 200. The results show that the flow becomes stabilized and finally steady beyond the critical rotational speed as ｜a｜ increases, regardless of the variation in Re and g＊. The value of critical rotational speed increases with increasing Re. The wake patterns change in the unsteady regimes for g＊ =1.5 and 0.7. With increasing ｜a｜, the time-averaged drag coemcient CD decreases and the lift coemcient CL increases, respectively. CD at Re = 160 and g＊ = 0.7 decreases rapidly, resulting in the smallest value at the same ｜a｜ for 1 ≤ ｜a｜≤ 2. CD augments with increasing g＊ at the same ｜a｜. For g＊=1.5, CD has a little disparity between the cases of Re= 160 and 200. For the flow past two still cylinders, CL is inversely proportional to g＊ of two cylinders for a fixed ｜a｜, and is not dependent on Re.     

Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma = 6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at freestream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibility effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely tO OCCUr.     

Transport in a ratchet with spatial disorder and time-delayed feedback

A Brownian motor with Gaussian short-range correlated spatial disorder and time-delayed feedback is investigated. The effects of disorder intensity, correlation strength and delay time on the transport properties of an overdamped periodic ratchet are discussed for different driving force. For small driving force, the disorder intensity can induce a peak in the drift motion and a linear increasing function in diffusion motion. For large driving force, the disorder intensity can suppress the drift motion but enhance the diffusion motion. For both small and large driving forces, the correlation strength of the spatial disorder can enhance the drift motion but suppress the diffusion motion. While the delay time can reduce the drift motion to a small value and enhance the diffusion motion to a large value. The drift motion increases as the driving force increases. However, the diffusion motion is either decreases or only increases slightly when the driving force increases.     

Enskog-Landau kinetic equation for multicomponent mixture. Analytical calculation of transport coefficients

The Enskog-Landau kinetic equation is considered to describe non-equilibrium processes of a mixture of charged hard spheres. This equation has been obtained in our previous papers by means of the non-equilibrium statistical operator method. The normal solution of this kinetic equation found in the first approximation using the standard Chapman-Enskog method is given. On the basis of the found solution the flows and transport coefficients have been calculated. All transport coefficients for multicomponent mixture of spherical Coulomb particles are presented analytically for the first time. Numerical calculations of thermal conductivity and thermal diffusion coefficient are performed for some specific mixtures of noble gases of high density. We compare the calculations with those ones for point-like neutral and charged particles. Received 10 June 1999 and Received in final form 15 October 1999     

Stochastic resonance in a bistable system with time-delayed feedback and non-Gaussian noise

The phenomenon of stochastic resonance in a bistable system with time-delayed feedback driven by non-Gaussian noise is investigated. Combining the small time delay approximation, the path-integral approach and the unified colored noise approximation, a general approximate Fokker–Planck equation of a stochastic system is obtained. The effects of the parameter q indicating the departure from the Gaussian noise, the delay time τ  , and the correlation time τ₀${\tau }_0$ of the non-Gaussian noise on the quasi-steady-state probability distribution function (SPD) and the signal-to-noise ratio (SNR) are discussed. It is found that the number of peaks in SPD and the reentrant transition between one peak and two peaks and then to one peak again in the curve of SNR depends on the parameter q, the delay time τ  , and the noise correlation time τ₀${\tau }_0$.     

Results for a fractional diffusion equation with a nonlocal term in spherical symmetry

We obtain time dependent solutions for a fractional diffusion equation containing a nonlocal term by considering the spherical symmetry and using the Green function approach. The nonlocal term incorporated in the diffusion equation may also be related to the spatial and time fractional derivative and introduces different regimes of spreading of the solution with the time evolution. In addition, a rich class of anomalous diffusion processes may be described from the results obtained here.     

Distribution function of fusion reaction products and entropy evolution

One of the outcomes of nuclear reactions is that reaction products have at birth distribution functions far from Maxwellian. What role do those distribution functions play in the evolution of the entropy of the system? The purpose of this work is to show the effect of the distribution functions of reactant and reaction products on the entropy of a system undergoing DD nuclear fusion reactions. This analysis is conducted with the help of the H-theorem, in the framework of kinetic theory. It will be found that at the onset of this reaction, generalized system entropy decreases markedly.