Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSome improvements on Moment-of-Fluid method in 3D rectangular hexahedrons
59
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The moment-of-fluid method (MOF) is an extension of the volume-of-fluid method with piecewise linear interface construction (VOF-PLIC). In MOF reconstruction, the optimized normal vector is determined from the reference centroid and the volume fraction by iteration. The state-of-art work by citet{milcent_moment--fluid_2020} proposed an analytic gradient of the objective function, which greatly reduces the computational cost. In this study, we further accelerate the MOF reconstruction algorithm by using Gauss-Newton iteration instead of Broyden-Fletcher-Goldfarb-Shanno (BFGS) iteration. We also propose an improved initial guess for MOF reconstruction, which improves the efficiency and the robustness of the MOF reconstruction algorithm. Our implementation of the code and test cases are available on our Github repository.
rate research
Read More
The moment-of-fluid (MOF) method is an extension of the volume-of-fluid method with piecewise linear interface construction (VOF-PLIC). By minimizing the least square error of the centroid of the cutting polyhedron, the MOF method reconstructs the linear interface without using any neighboring information. Traditional MOF involves iteration while finding the optimized linear reconstruction. Here, we propose an alternative approach based on a machine learning algorithm: Decision Tree algorithm. A training data set is generated from a list of random cuts of a unit cube by plane. The Decision Tree algorithm extracts the input-output relationship from the training data, so that the resulting function determines the normal vector of the reconstruction plane directly, without any iteration. The present method is tested on a range of popular interface advection test problems. Numerical results show that our approach is much faster than the iteration-based MOF method while provides compatible accuracy with the conventional MOF method.
Simulating inhomogeneous flows with different characteristic scales in different coordinate directions using the collide-and-stream based lattice Boltzmann methods (LBM) can be accomplished efficiently using rectangular lattice grids. We develop and investigate a new rectangular central moment LBM based on non-orthogonal moment basis (referred to as RC-LBM). The equilibria to which the central moments relax under collision in this approach are obtained from matching with those corresponding to the continuous Maxwell distribution. A Chapman-Enskog analysis is performed to derive the correction terms to the second order moment equilibria involving the grid aspect ratio and velocity gradients that restores the isotropy of the viscous stress tensor and eliminates the non-Galilean invariant cubic velocity terms of the resulting hydrodynamical equations. A special case of this rectangular formulation involving the raw moments (referred to as the RNR-LBM) is also constructed. The resulting schemes represent a considerable simplification, especially for the transformation matrices and isotropy corrections, and improvement over the existing MRT-LB schemes on rectangular lattice grids that use orthogonal moment basis. Numerical validation study of both the RC-LBM and RNR-LBM for a variety of benchmark flow problems are performed that show good accuracy at various grid aspect ratios. The ability of our proposed schemes to simulate flows using relatively lower grid aspect ratios than considered in prior rectangular LB approaches is demonstrated. Furthermore, simulations reveal the superior stability characteristics of the RC-LBM over RNR-LBM in handling shear flows at lower viscosities and/or higher characteristic velocities. In addition, computational advantages of using our rectangular LB formulation in lieu of that based on the square lattice is shown.
We studied the role of electron physics in 3D two-fluid 10-moment simulation of the Ganymedes magnetosphere. The model captures non-ideal physics like the Hall effect, the electron inertia, and anisotropic, non-gyrotropic pressure effects. A series of analyses were carried out: 1) The resulting magnetic field topology and electron and ion convection patterns were investigated. The magnetic fields were shown to agree reasonably well with in-situ measurements by the Galileo satellite. 2) The physics of collisionless magnetic reconnection were carefully examined in terms of the current sheet formation and decomposition of generalized Ohms law. The importance of pressure anisotropy and non-gyrotropy in supporting the reconnection electric field is confirmed. 3) We compared surface brightness morphology, represented by surface electron and ion pressure contours, with oxygen emission observed by the Hubble Space Telescope (HST). The correlation between the observed emission morphology and spatial variability in electron/ion pressure was demonstrated. Potential extension to multi-ion species in the context of Ganymede and other magnetospheric systems is also discussed.
The computation of flows with large density contrasts is notoriously difficult. To alleviate the difficulty we consider a discretization of the Navier-Stokes equation that advects mass and momentum in a consistent manner. Incompressible flow with capillary forces is modeled and the discretization is performed on a staggered grid of Marker and Cell type. The Volume-of-Fluid method is used to track the interface and a Height-Function method is used to compute surface tension. The advection of the volume fraction is performed using either the Lagrangian-Explicit / CIAM (Calcul dInterface Affine par Morceaux) method or the Weymouth and Yue (WY) Eulerian-Implicit method. The WY method conserves fluid mass to machine accuracy provided incompressibility is satisfied. To improve the stability of these methods momentum fluxes are advected in a manner consistent with the volume-fraction fluxes, that is a discontinuity of the momentum is advected at the same speed as a discontinuity of the density. To find the density on the staggered cells on which the velocity is centered, an auxiliary reconstruction of the density is performed. The method is tested for a droplet without surface tension in uniform flow, for a droplet suddenly accelerated in a carrying gas at rest at very large density ratio without viscosity or surface tension, for the Kelvin-Helmholtz instability, for a 3mm-diameter falling raindrop and for an atomizing flow in air-water conditions.
We introduce a numerical method and python package, https://github.com/andillio/CHiMES, that simulates quantum systems initially well approximated by mean field theory using a second order extension of the classical field approach. We call this the field moment expansion method. In this way, we can accurately approximate the evolution of first and second field moments beyond where the mean field theory breaks down. This allows us to estimate the quantum breaktime of a classical approximation without any calculations external to the theory. We investigate the accuracy of the field moment expansion using a number of well studied quantum test problems. Interacting Bosonic systems similar to scalar field dark matter are chosen as test problems. We find that successful application of this method depends on two conditions: the quantum system must initially be well described by the classical theory, and that the growth of the higher order moments be hierarchical.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
