Package Bees for Honey Production and Pollination

PDF pages: 1

On the modelling of isothermal gas flows at the microscale

This paper makes two new propositions regarding the modelling of rarefied (non-equilibrium) isothermal gas flows at the microscale. The first is a new test case for benchmarking high-order, or extended, hydrodynamic models for these flows. This standing time-varying shear-wave problem does not require boundary conditions to be specified at a solid surface, so is useful for assessing whether fluid models can capture rarefaction effects in the bulk flow. We assess a number of different proposed extended hydrodynamic models, and we find the R13 equations perform the best in this case. Our second proposition is a simple technique for introducing non-equilibrium effects caused by the presence of solid surfaces into the computational fluid dynamics framework. By combining a new model for slip boundary conditions with a near-wall scaling of the Navier--Stokes constitutive relations, we obtain a model that is much more accurate at higher Knudsen numbers than the conventional second-order slip model. We show that this provides good results for combined Couette/Poiseuille flow, and that the model can predict the stress/strain-rate inversion that is evident from molecular simulations. The model's generality to non-planar geometries is demonstrated by examining low-speed flow around a micro-sphere. It shows a marked improvement over conventional predictions of the drag on the sphere, although there are some questions regarding its stability at the highest Knudsen numbers

A wall-function approach to incorporating Knudsen-layer effects in gas micro flow simulations

For gas flows in microfluidic configurations, the Knudsen layer close to the wall can comprise a substantial part of the entire flow field and has a major effect on quantities such as the mass flow rate through micro devices. The Knudsen layer itself is characterized by a highly nonlinear relationship between the viscous stress and the strain rate of the gas, so even if the Navier-Stokes equations can be used to describe the core gas flow they are certainly inappropriate for the Knudsen layer itself. In this paper we propose a "wall-function" model for the stress/strain rate relations in the Knudsen layer. The constitutive structure of the Knudsen layer has been derived from results from kinetic theory for isothermal shear flow over a planar surface. We investigate the ability of this simplified model to predict Knudsen-layer effects in a variety of configurations. We further propose a semi-empirical Knudsen-number correction to this wall function, based on high-accuracy DSMC results, to extend the predictive capabilities of the model to greater degrees of rarefaction

Capturing the Knudsen layer in continuum-fluid models of non-equilibrium gas flows

In hypersonic aerodynamics and microflow device design, the momentum and energy fluxes to solid surfaces are often of critical importance. However, these depend on the characteristics of the Knudsen layer - the region of local non-equilibrium existing up to one or two molecular mean free paths from the wall in any gas flow near a surface. While the Knudsen layer has been investigated extensively using kinetic theory, the ability to capture it within a continuum-fluid formulation (in conjunction with slip boundary conditions) suitable for current computational fluid dynamics toolboxes would offer distinct and practical computational advantages

The Galaxy-Halo/Subhalo connection: Mass relations and implications for some satellite occupational distributions

We infer the local stellar-to-halo/subhalo mass relations (MRs) for central and satellite galaxies separately. We constraint this relations by using several combinations of observational data, consisting of the total galaxy stellar mass function, its decomposition into centrals and satellites, and the projected two-point correlation functions measured in different stellar mass bins. The differences among the resulting MRs are within the model-fit uncertainties (which are very small, smaller than the intrinsic scatter between galaxy and halo mass), no matter what combination of data are used. We also constrain the values of the intrinsic scatter around the central-halo (CH) and satellite-subhalo (SS) MRs assuming them to be constant: sigma_c= 0.168+-0.051 dex and sigma_s = 0.172+-0.057 dex, respectively. The CH and SS MRs are actually different, in particular when we take the subhalo mass at the present-day epoch instead of at their accretion time. We calculate several population statistics. We find that the central galaxy is not on average within the mass distribution of the most-massive satellite, even for cluster-sized halos, i.e., centrals are not a mere realization of the high-end of the satellite mass function; however for > 3x10^13 msun halos, ~15% of centrals could be. We also find that the probabilities of MW-sized halos of having N Magellanic-Clouds (MCs)-sized satellites agree well with observational measures; for a halo mass of 2x10^12 msun, the probability to have 2 MCs is 5.4%, but if we exclude those systems with satellites larger than the MCs, then the probability decreases to <2.2%.Comment: 18 pages, 7 figures. Accepted for publication in ApJ. Some references were adde

Coupled continuum hydrodynamics and molecular dynamics method for multiscale simulation

We present a new hybrid methodology for carrying out multiscale simulations of flow problems lying between continuum hydrodynamics and molecular dynamics, where macro/micro lengthscale separation exists only in one direction. Our multiscale method consists of an iterative technique that couples mass and momentum flux between macro and micro domains, and is tested on a converging/diverging nanochannel case containing flow of a simple Lennard-Jones liquid. Comparisons agree well with a full MD simulation of the same test case

The usefulness of higher-order constitutive relations for describing the Knudsen layer

The Knudsen layer is an important rarefaction phenomenon in gas flows in and around microdevices. Its accurate and efficient modeling is of critical importance in the design of such systems and in predicting their performance. In this paper we investigate the potential that higher-order continuum equations may have to model the Knudsen layer, and compare their predictions to high-accuracy DSMC (direct simulation Monte Carlo) data, as well as a standard result from kinetic theory. We find that, for a benchmark case, the most common higher-order continuum equation sets (Grad's 13 moment, Burnett, and super-Burnett equations) cannot capture the Knudsen layer. Variants of these equation families have, however, been proposed and some of them can qualitatively describe the Knudsen layer structure. To make quantitative comparisons, we obtain additional boundary conditions (needed for unique solutions to the higher-order equations) from kinetic theory. However, we find the quantitative agreement with kinetic theory and DSMC data is only slight

A rapidly moving spot on jupiter's north temperate belt

Photographic observations of rapidly moving spot in Jupiter atmospher

Time-step coupling for hybrid simulations of multiscale flows

A new method is presented for the exploitation of time-scale separation in hybrid continuum-molecular models of multiscale flows. Our method is a generalisation of existing approaches, and is evaluated in terms of computational efficiency and physical/numerical error. Comparison with existing schemes demonstrates comparable, or much improved, physical accuracy, at comparable, or far greater, efficiency (in terms of the number of time-step operations required to cover the same physical time). A leapfrog coupling is proposed between the ‘macro’ and ‘micro’ components of the hybrid model and demonstrates potential for improved numerical accuracy over a standard simultaneous approach. A general algorithm for a coupled time step is presented. Three test cases are considered where the degree of time-scale separation naturally varies during the course of the simulation. First, the step response of a second-order system composed of two linearly-coupled ODEs. Second, a micro-jet actuator combining a kinetic treatment in a small flow region where rarefaction is important with a simple ODE enforcing mass conservation in a much larger spatial region. Finally, the transient start-up flow of a journal bearing with a cylindrical rarefied gas layer. Our new time-stepping method consistently demonstrates as good as or better performance than existing schemes. This superior overall performance is due to an adaptability inherent in the method, which allows the most-desirable aspects of existing schemes to be applied only in the appropriate conditions