Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHyperbolic Model Reduction for Kinetic Equations
288
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We make a brief historical review to the moment model reduction to the kinetic equations, particularly the Grads moment method for Boltzmann equation. The focus is on the hyperbolicity of the reduced model, which is essential to the existence of its classical solution as a Cauchy problem. The theory of the framework we developed in last years is then introduced, which may preserve the hyperbolic nature of the kinetic equations with high universality. Some lastest progress on the comparison between models with/without hyperbolicity is presented to validate the hyperbolic moment models for rarefied gases.
rate research
Read More
Entropy production during the process of thermal phase-separation of multiphase flows is investigated by means of a discrete Boltzmann kinetic model. The entropy production rate is found to increase during the spinodal decomposition stage and to decrease during the domain growth stage, attaining its maximum at the crossover between the two. Such behaviour provides a natural criterion to identify and discriminate between the two regimes. Furthermore, the effects of heat conductivity, viscosity and surface tension on the entropy production rate are investigated by systematically probing the interplay between non-equilibrium energy and momentum fluxes. It is found that the entropy production rate due to energy fluxes is an increasing function of the Prandtl number, while the momentum fluxes exhibit an opposite trend. On the other hand, both contributions show an increasing trend with surface tension. The present analysis inscribes within the general framework of non-equilibrium thermodynamics and consequently it is expected to be relevant to a broad class of soft-flowing systems far from mechanical and thermal equilibrium.
We expand on two previous developments in the modeling of discrete-time Langevin systems. One is the well-documented Gr{o}nbech-Jensen Farago (GJF) thermostat, which has been demonstrated to give robust and accurate configurational sampling of the phase space. Another is the recent discovery that also kinetics can be accurately sampled for the GJF method. Through a complete investigation of all possible finite difference approximations to the velocity, we arrive at two main conclusions:~1) It is not possible to define a so-called on-site velocity such that kinetic temperature will be correct and independent of the time step, and~2) there exists a set of infinitely many possibilities for defining a two-point (leap-frog) velocity that measures kinetic energy correctly for linear systems in addition to the correct configurational statistics obtained from the GJF algorithm. We give explicit expressions for the possible definitions, and we incorporate these into convenient and practical algorithmic forms of the normal Verlet-type algorithms along with a set of suggested criteria for selecting a useful definition of velocity.
We show how the basic idea of parabolic Jacobi relaxation can be modified to obtain a new class of hyperbolic relaxation schemes that are suitable for the solution of elliptic equations. Some of the analytic and numerical properties of hyperbolic relaxation are examined. We describe its implementation as a first order system in a pseudospectral evolution code, demonstrating that certain elliptic equations can be solved within a framework for hyperbolic evolution systems. Applications include various initial data problems in numerical general relativity. In particular we generate initial data for the evolution of a massless scalar field, a single neutron star, and binary neutron star systems.
We present modifications to the Athena++ framework to enable use of general equations of state (EOS). Part of our motivation for doing so is to model transient astrophysics phenomena, as these types of events are often not well approximated by an ideal gas. This necessitated changes to the Riemann solvers implemented in Athena++. We discuss the adjustments made to the HLLC, and HLLD solvers and EOS calls required for arbitrary EOS. We demonstrate the reliability of our code in a number of tests which utilize a relatively simple, but non-trivial EOS based on hydrogen ionization, appropriate for the transition from atomic to ionized hydrogen. Additionally, we perform tests using an electron-positron Helmholtz EOS, appropriate for regimes where nuclear statistical equilibrium is a good approximation. These new complex EOS tests overall show that our modifications to Athena++ accurately solve the Riemann problem with linear convergence and linear-wave tests with quadratic convergence. We provide our test solutions as a means to check the accuracy of other hydrodynamic codes. Our tests and additions to Athena++ will enable further research into (magneto)hydrodynamic problems where realistic treatments of the EOS are required.
We present an analytical derivation of the transport coefficients of a relativistic gas in (2+1) dimensions for both Chapman-Enskog (CE) asymptotics and Grads expansion methods. Moreover, we develop a systematic calibration method, connecting the relaxation time of relativistic kinetic theory to the transport parameters of the associated dissipative hydrodynamic equations. Comparison between the analytical results and numerical simulations, shows that the CE method correctly captures dissipative effects, while Grads method does not. The resulting calibration procedure based on the CE method opens the way to the quantitative kinetic description of dissipative relativistic fluid dynamics under fairly general conditions, namely flows with strongly non-linearities, in non-ideal geometries, across both ultra-relativistic and near-non-relativistic regimes.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
