Source author record

Sha Liu

Sha Liu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Investigation of nonlinear squeeze-film damping involving rarefied gas effect in micro-electro-mechanical-systems

In this paper, the nonlinear squeeze-film damping (SFD) involving rarefied gas effect in the micro-electro-mechanical-systems (MEMS) is investigated. Considering the motion of structures (beam, cantilever, and membrane) in MEMS, the dynamic response of structure will be influenced largely by the squeeze-film damping. In the traditional model, a viscous damping assumption that damping force is linear with moving velocity is used. As the nonlinear damping phenomenon is observed for a micro-structure oscillating with a high-velocity, this assumption is invalid and will generates error result for predicting the response of micro-structure. In addition, due to the small size of device and the low pressure of encapsulation, the gas in MEMS usually is rarefied gas. Therefore, to correctly predict the damping force, the rarefied gas effect must be considered. To study the nonlinear SFD phenomenon involving the rarefied gas effect, a kinetic method, namely discrete unified gas kinetic scheme (DUGKS), is introduced. And based on DUGKS, two solving methods, a traditional decoupled method (Eulerian scheme) and a coupled framework (arbitrary Lagrangian-Eulerian scheme), are adopted. With these two methods, two basic motion forms, linear (perpendicular) and tilting motions of a rigid micro-beam, are studied with forced and free oscillations.

preprint2020arXiv

A multi-scale kinetic inviscid flux extracted from the gas-kinetic scheme for simulating incompressible and compressible flows

A Kinetic Inviscid Flux (KIF) is proposed for simulating incompressible and compressible flows. It is constructed based on the direct modeling of multi-scale flow behaviors, which is used in the Gas-Kinetic Scheme (GKS), the Unified Gas-Kinetic Scheme (UGKS), the Discrete Unified Gas-Kinetic Scheme (DUGKS), etc.. In KIF, the discontinuities (such as the shock wave) that can not be well resolved by mesh cells are mainly solved by the Kinetic Flux Vector Splitting (KFVS) method representing the free transport mechanism (or micro-scale mechanism), while in other flow regions that are smooth, the flow behavior is solved mainly by the central-scheme-like Totally Thermalized Transport (TTT). The weights of KFVS and TTT in KIF is automatically determined by those in the theory of direct modeling. Two ways of choosing the weights in KIF are proposed, which are actually the weights adopted in the UGKS and the DUGKS, respectively. By using the test cases of Sod shock tube, rarefaction wave, the boundary layer of flat plate, the cavity flow and hypersonic flow over circular cylinder, the validity and accuracy of the present method are examined. The KIF does not suffer from the carbuncle phenomenon, and does not introduce extra numerical viscosity in smooth regions. Especially, in the case of hypersonic cylinder, it gives a quite sharp and clear density and temperature contours. The KIF can be viewed as an inviscid-viscous splitting version of the GKS. By the doing this splitting, it is easy to be used in the traditional CFD frameworks. It can also be classified as a new type in the numerical schemes based on the kinetic theory that are represented by the works in Ref. \cite{SunA} and Ref. \cite{ohwada2018a}, except the weights are determined by the weights of direct modeling.

preprint2020arXiv

A multiscale discrete velocity method for model kinetic equations

In this paper, authors focus effort on improving the conventional discrete velocity method (DVM) into a multiscale scheme in finite volume framework for gas flow in all flow regimes. Unlike the typical multiscale kinetic methods unified gas-kinetic scheme (UGKS) and discrete unified gas-kinetic scheme (DUGKS), which concentrate on the evolution of the distribution function at the cell interface, in the present scheme the flux for macroscopic variables is split into the equilibrium part and the nonequilibrium part, and the nonequilibrium flux is calculated by integrating the discrete distribution function at the cell center, which overcomes the excess numerical dissipation of the conventional DVM in the continuum flow regime. Afterwards, the macroscopic variables are finally updated by simply integrating the discrete distribution function at the cell center, or by a blend of the increments based on the macroscopic and the microscopic systems, and the multiscale property is achieved. Several test cases, involving unsteady, steady, high speed, low speed gas flows in all flow regimes, have been performed, demonstrating the good performance of the multiscale DVM from free molecule to continuum Navier-Stokes solutions and the multiscale property of the scheme is proved.

preprint2020arXiv

Simplified Unified Wave-Particle Method with Quantified Model-Competition Mechanism for Numerical Calculation of Multi-Scale Flows

A Quantified Model-Competition (QMC) mechanism for multi-scale flows is extracted from the integral (analytical) solution of the Boltzmann-BGK model equation. In the QMC mechanism, the weight of the rarefied model and the weight of the continuum (aerodynamic/hydrodynamic) model are quantified. Then, a Simplified Unified Wave-Particle method (SUWP) is constructed based the on the QMC mechanism. In the SUWP, the stochastic particle method and the continuum Navier-Stokes method are combined together. Their weights are determined by the QMC mechanism quantitatively in every discrete cells of the computational domain. The validity and accuracy of the present numerical method are examined using a series of test cases including the high non-equilibrium shock wave structure case, the unsteady Sod shock-tube case with a wide range of Kn number, the hypersonic flow around the circular cylinder from the free-molecular regime to the near continuum regime, and the viscous boundary layer case. In the construction process of the present method, an anti-dissipation effect in the continuum mechanism is also discussed.

preprint2015arXiv

1D momentum-conserving systems: the conundrum of anomalous versus normal heat transport

Transport and diffusion of heat in one dimensional (1D) nonlinear systems which {\it conserve momentum} is typically thought to proceed anomalously. Notable exceptions, however, exist of which the rotator model is a prominent case. Therefore, the quest arises to identify the origin of manifest anomalous transport in those low dimensional systems. Here, we develop the theory for both, momentum/heat diffusion and its corresponding momentum/heat transport. We demonstrate that the second temporal derivative of the mean squared deviation of the momentum spread is proportional to the equilibrium correlation of the total momentum flux. This result in turn relates, via the integrated momentum flux correlation, to an effective viscosity, or equivalently, to the underlying momentum diffusivity. We put forward the intriguing hypothesis that a fluid-like momentum dynamics with a {\it finite viscosity} causes {\it normal} heat transport; its corollary being that superdiffusive momentum diffusion with an intrinsic {\it diverging viscosity} in turn yields {\it anomalous} heat transport. This very hypothesis is corroborated over wide extended time scales by use of precise molecular dynamics simulations. The numerical validation of the hypothesis involves three distinct archetype classes of nonlinear pair-interaction potentials: (i) a globally bounded pair interaction (the noted rotator model), (ii) unbounded interactions acting at large distances (the Fermi-Pasta-Ulam $β$ model, or the rotator model amended with harmonic pair interactions) and (iii), a pair interaction potential being unbounded at short distances while displaying an asymptotic free part (Lennard-Jones model).

preprint2015arXiv

Renormalized phonons in nonlinear lattices: A variational approach

We propose a variational approach to study renormalized phonons in momentum conserving nonlinear lattices with either symmetric or asymmetric potentials. To investigate the influence of pressure to phonon properties, we derive an inequality which provides both the lower and upper bound of the Gibbs free energy as the associated variational principle. This inequality is a direct extension to the Gibbs-Bogoliubov inequality. Taking the symmetry effect into account, the reference system for the variational approach is chosen to be harmonic with an asymmetric quadratic potential which contains variational parameters. We demonstrate the power of this approach by applying it to one dimensional nonlinear lattices with a symmetric or asymmetric Fermi-Pasta- Ulam type potential. For a system with a symmetric potential and zero pressure, we recover existing results. For other systems which beyond the scope of existing theories, including those having the symmetric potential and pressure, and those having the asymmetric potential with or without pressure, we also obtain accurate sound velocity.

preprint2014arXiv

Triggering waves in nonlinear lattices: Quest for anharmonic phonons and corresponding mean free paths

Guided by a stylized experiment we develop a self-consistent anharmonic phonon concept for nonlinear lattices which allows for explicit "visualization." The idea uses a small external driving force which excites the front particles in a nonlinear lattice slab and subsequently one monitors the excited wave evolution using molecular dynamics simulations. This allows for a simultaneous, direct determination of the existence of the phonon mean free path with its corresponding anharmonic phonon wavenumber as a function of temperature. The concept for the mean free path is very distinct from known prior approaches: the latter evaluate the mean free path only indirectly, via using both, a scale for the phonon relaxation time and yet another one for the phonon velocity. Notably, the concept here is neither limited to small lattice nonlinearities nor to small frequencies. The scheme is tested for three strongly nonlinear lattices of timely current interest which either exhibit normal or anomalous heat transport.

preprint2013arXiv

Anomalous Heat Diffusion

Consider anomalous energy spread in solid phases, i.e., $MSD= \int (x -{< x >}_E)^2 ρ_E(x,t)dx \propto t^β$, as induced by a small initial excess energy perturbation distribution $ρ_{E}(x,t=0)$ away from equilibrium. The associated total thermal equilibrium heat flux autocorrelation function $C_{JJ}(t)$ is shown to obey rigorously the intriguing relation, $d^2 MSD/dt^2 = 2C_{JJ}(t)/(k_BT^2c)$, where $c$ is the specific volumetric heat capacity. Its integral assumes a time-local Helfand-moment relation; i.e. $ dMSD/dt|_{t=t_s} = 2/(k_BT^2c)\int_0^{t_s} C_{JJ}(s)ds$, where the chosen cut-off time $t_s$ is determined by the maximal signal velocity for heat transfer. Given the premise that the averaged nonequilibrium heat flux is governed by an anomalous heat conductivity, energy diffusion scaling necessarily determines a corresponding anomalous thermal conductivity scaling behavior.

preprint2013arXiv

Classical Heat Transport in Anharmonic Molecular Junctions: Exact Solutions

We study full counting statistics for classical heat transport through anharmonic/nonlinear molecular junctions formed by interacting oscillators. Analytical result of the steady state heat flux for an overdamped anharmonic junction with arbitrary temperature bias is obtained. It is found that the thermal conductance can be expressed in terms of temperature dependent effective force constant. The role of anharmonicity is identified. We also give the general formula for the second cumulant of heat in steady state, as well as the average geometric heat flux when two temperatures are modulated adiabatically. We present an anharmonic example for which all cumulants for heat can be obtained exactly. For a bounded single oscillator model with mass we found that the cumulants are independent of the nonlinear potential.

preprint2012arXiv

Anomalous Heat Conduction and Anomalous Diffusion in Low Dimensional Nanoscale Systems

Thermal transport is an important energy transfer process in nature. Phonon is the major energy carrier for heat in semiconductor and dielectric materials. In analogy to Ohm's law for electrical conductivity, Fourier's law is a fundamental rule of heat transfer in solids. It states that the thermal conductivity is independent of sample scale and geometry. Although Fourier's law has received great success in describing macroscopic thermal transport in the past two hundreds years, its validity in low dimensional systems is still an open question. Here we give a brief review of the recent developments in experimental, theoretical and numerical studies of heat transport in low dimensional systems, include lattice models, nanowires, nanotubes and graphenes. We will demonstrate that the phonon transports in low dimensional systems super-diffusively, which leads to a size dependent thermal conductivity. In other words, Fourier's law is breakdown in low dimensional structures.

preprint2012arXiv

Geometric Heat Flux for Classical Thermal Transport in Interacting Open Systems

We study classical heat conduction in a dissipative open system composed of interacting oscillators. By exactly solving a twisted Fokker-Planck equation which describes the full counting statistics of heat flux flowing through the system, we identify the geometric-phase-like effect and examine its impact on the classical heat transport. Particularly, we find that the nonlinear interaction as well as the closely related temperature-dependence of system-parameters are crucial in manifesting the geometric-phase contribution of heat flux. Finally, we propose an electronic experiment based on RC circuits to verify our theoretical predictions.