Importance Sampling Variance Reduction for the Fokker-Planck Rarefied Gas Particle Method

Models and methods that are able to accurately and efficiently predict the flows of low-speed rarefied gases are in high demand, due to the increasing ability to manufacture devices at micro and nano scales. One such model and method is a Fokker-Planck approximation to the Boltzmann equation, which can be solved numerically by a stochastic particle method. The stochastic nature of this method leads to noisy estimates of the thermodynamic quantities one wishes to sample when the signal is small in comparison to the thermal velocity of the gas. Recently, Gorji et al have proposed a method which is able to greatly reduce the variance of the estimators, by creating a correlated stochastic process which acts as a control variate for the noisy estimates. However, there are potential difficulties involved when the geometry of the problem is complex, as the method requires the density to be solved for independently. Importance sampling is a variance reduction technique that has already been shown to successfully reduce the noise in direct simulation Monte Carlo calculations. In this paper we propose an importance sampling method for the Fokker-Planck stochastic particle scheme. The method requires minimal change to the original algorithm, and dramatically reduces the variance of the estimates. We test the importance sampling scheme on a homogeneous relaxation, planar Couette flow and a lid-driven-cavity flow, and find that our method is able to greatly reduce the noise of estimated quantities. Significantly, we find that as the characteristic speed of the flow decreases, the variance of the noisy estimators becomes independent of the characteristic speed

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

Sexual Assault Center Advocates and Psychotherapists: An Exploratory Study of Interventions

In the following qualitative study, the researcher attempted to identify and compare the self-identified approaches taken by Sexual Assault Center (SAC) advocates (n = 4) and outpatient psychotherapists (n = 4) in working with adult survivors of sexual violence. Using grounded theory methodology, data from semi-structured interviews were coded using an open coding process. From this, six themes emerged from the interviews with SAC advocates and ten themes emerged from the interviews with psychotherapists, suggesting approaches used when providing crisis intervention and counseling to adult survivors of sexual violence. These themes and their corresponding subthemes are discussed in this report. The study demonstrated that when providing crisis intervention and counseling to adult survivors of sexual violence, SAC advocates and psychotherapists take many of the same approaches, but that the work of SAC advocates is structured to be more short-term and to meet the immediate needs of survivors of sexual violence, whereas the work of psychotherapists tends to be more long-term and with the use of specific therapeutic treatment interventions and assessment tools. Several implications can be drawn from this study. This study informs direct practice of and referrals made by social workers, service choices made by survivors, policy surrounding the funding of SACs, and future research related to crisis intervention and counseling approaches used with survivors of sexual violence

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

Dynamics of liquid nano-threads : fluctuation-driven instability and rupture

The instability and rupture of nanoscale liquid threads is shown to strongly depend on thermal fluctuations. These fluctuations are naturally occurring within molecular dynamics (MD) simulations and can be incorporated via fluctuating hydrodynamics into a stochastic lubrication equation (SLE). A simple and robust numerical scheme is developed for the SLE that is validated against MD for both the initial (linear) instability and the nonlinear rupture process. Particular attention is paid to the rupture process and its statistics, where the `double-cone’ profile reported by Moseler & Landmann [Science, 2000, 289(5482): 1165-1169] is observed, as well as other distinct profile forms depending on the flow conditions. Comparison to the Eggers’ similarity solution [Physical Review Letters, 2002, 89(8): 084502], a power law of the minimum thread radius against time to rupture, shows agreement only at low surface tension; indicating that surface tension cannot generally be neglected when considering rupture dynamics

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

Sexual Assault Center Advocates and Psychotherapists: An Exploratory Study of Interventions

In the following qualitative study, the researcher attempted to identify and compare the self- identified approaches taken by Sexual Assault Center (SAC) advocates (n = 4) and outpatient psychotherapists (n = 4) in working with adult survivors of sexual violence. Using grounded theory methodology, data from semi-structured interviews were coded using an open coding process. From this, six themes emerged from the interviews with SAC advocates and ten themes emerged from the interviews with psychotherapists, suggesting approaches used when providing crisis intervention and counseling to adult survivors of sexual violence. These themes and their corresponding subthemes are discussed in this report. The study demonstrated that when providing crisis intervention and counseling to adult survivors of sexual violence, SAC advocates and psychotherapists take many of the same approaches, but that the work of SAC advocates is structured to be more short-term and to meet the immediate needs of survivors of sexual violence, whereas the work of psychotherapists tends to be more long-term and with the use of specific therapeutic treatment interventions and assessment tools. Several implications can be drawn from this study. This study informs direct practice of and referrals made by social workers, service choices made by survivors, policy surrounding the funding of SACs, and future research related to crisis intervention and counseling approaches used with survivors of sexual violence

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