High-dimensional sensitivity analysis and uncertainty quantification for multiscale gas dynamics, spanning the continuum to rarefied regimes, require computationally efficient and mathematically consistent gradient evaluation. This paper develops a discrete adjoint method for the unified gas-kinetic scheme (UGKS) based on a dual-consistent formulation. The adjoint system is derived directly from the discrete microscopic velocity-distribution equation coupled with the macroscopic-moment compatibility conditions. To resolve the stiff cross-scale coupling, we propose an asymptotic-preserving (AP) adjoint formulation constructed via macroscopic-moment projection and microscopic lifting. Under this framework, the AP adjoint formulation eliminates the stiff collision coupling and removes the collision-time step restriction in the continuum regime. Numerically, a memory-efficient residual-evaluation algorithm that mirrors the forward UGKS cell-vertex data structure is implemented to bypass the memory bottleneck in velocity space. Furthermore, a macroscopic--microscopic predictor--corrector implicit marching scheme is designed to accelerate convergence without solving a globally coupled system. The accuracy, consistency, and robustness of the proposed AP-adjoint scheme are rigorously verified against an independent linearized UGKS solver across a wide range of Knudsen numbers, including lid-driven cavity heat conduction, microchannel thermal creep flow, and hypersonic flow past a circular cylinder.
A fast-converging and asymptotic-preserving adjoint shape optimization method is proposed for drag reduction of multiscale gas flows in vibrating micro-electro-mechanical systems. The convergence of the Boltzmann kinetic equation is accelerated by macroscopic synthetic equations, whose constitutive relations integrate...
In this work, a globally stiffly accurate Implicit-Explicit (IMEX) Runge-Kutta scheme is developed and implemented in the GBS code [Ricci et al., Plasma Phys. Control. Fusion, 2012], for two-fluid plasma turbulence simulations. The stiffest phenomena, governed by shear Alfv\'en waves and parallel diffusion, are treated...
Micol Bassanini, Simone Deparis, Paolo Ricci et al.· 0 citations
A smooth artificial viscosity treatment based on smooth approximations of non-differentiable pointwise operations is introduced based on smooth approximations of non-differentiable pointwise operations to ensure robust and differentiable nonlinear solves in the presence of shocks.
Julian Andrej, Jean-Sylvain Camier, V. Dobrev et al.· 0 citations
Pressure anisotropy and off-diagonal pressure stresses are ubiquitous and play important roles in collisionless/weakly collisional plasmas. The Chew-Goldberger-Low (CGL) MHD model is often used; however, it can lose hyperbolicity when the pressure anisotropy or plasma beta becomes large, making it hard to develop appro...
K. Akutagawa, Shinsuke Imada, M. Shoda· 0 citations
This work introduces a hyperbolic relaxation system designed for the simulation of kinetic equations in the low-Mach number limit. The methodology is built upon a micro-macro decomposition of the scaled BGK model, which reformulates the kinetic distribution function into a coupled system consisting of a macroscopic equ...
G. Dimarco, Axel Klar, Theresa Kofler et al.· 0 citations
We present a new semi-implicit structure-preserving (SP) finite volume discretization for a unified first-order hyperbolic model of continuum mechanics. This framework provides a common mathematical description for fluids and solids and incorporates transport, viscous effects, heat conduction, and elastic deformations...
M. Dumbser, W. Boscheri, M. Tavelli et al.· 1 citation
We use cookies to run the site and, with your consent, for analytics and to show ads.
See our Cookie Policy.