Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userImproved phase-field models of melting and dissolution in multi-component flows
108
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We develop and analyse the first second-order phase-field model to combine melting and dissolution in multi-component flows. This provides a simple and accurate way to simulate challenging phase-change problems in existing codes. Phase-field models simplify computation by describing separate regions using a smoothed phase field. The phase field eliminates the need for complicated discretisations that track the moving phase boundary. However standard phase-field models are only first-order accurate. They often incur an error proportional to the thickness of the diffuse interface. We eliminate this dominant error by developing a general framework for asymptotic analysis of diffuse-interface methods in arbitrary geometries. With this framework we can consistently unify previous second-order phase-field models of melting and dissolution and the volume-penalty method for fluid-solid interaction. We finally validate second-order convergence of our model in two comprehensive benchmark problems using the open-source spectral code Dedalus.
rate research
Read More
In this paper, a high order quasi-conservative discontinuous Galerkin (DG) method using the non-oscillatory kinetic flux is proposed for the 5-equation model of compressible multi-component flows with Mie-Gruneisen equation of state. The method mainly consists of three steps: firstly, the DG method with the non-oscillatory kinetic flux is used to solve the conservative equations of the model; secondly, inspired by Abgralls idea, we derive a DG scheme for the volume fraction equation which can avoid the unphysical oscillations near the material interfaces; finally, a multi-resolution WENO limiter and a maximum-principle-satisfying limiter are employed to ensure oscillation-free near the discontinuities, and preserve the physical bounds for the volume fraction, respectively. Numerical tests show that the method can achieve high order for smooth solutions and keep non-oscillatory at discontinuities. Moreover, the velocity and pressure are oscillation-free at the interface and the volume fraction can stay in the interval [0,1].
We present a particle method for estimating the curvature of interfaces in volume-of-fluid simulations of multiphase flows. The method is well suited for under-resolved interfaces, and it is shown to be more accurate than the parabolic fitting that is employed in such cases. The curvature is computed from the equilibrium positions of particles constrained to circular arcs and attracted to the interface. The proposed particle method is combined with the method of height functions at higher resolutions, and it is shown to outperform the current combinations of height functions and parabolic fitting. The algorithm is conceptually simple and straightforward to implement on new and existing software frameworks for multiphase flow simulations thus enhancing their capabilities in challenging flow problems. We evaluate the proposed hybrid method on a number of two- and three-dimensional benchmark flow problems and illustrate its capabilities on simulations of flows involving bubble coalescence and turbulent multiphase flows.
A computational fluid dynamics (CFD) simulation framework for predicting complex flows is developed on the Tensor Processing Unit (TPU) platform. The TPU architecture is featured with accelerated performance of dense matrix multiplication, large high bandwidth memory, and a fast inter-chip interconnect, which makes it attractive for high-performance scientific computing. The CFD framework solves the variable-density Navier-Stokes equation using a Low-Mach approximation, and the governing equations are discretized by a finite difference method on a collocated structured mesh. It uses the graph-based TensorFlow as the programming paradigm. The accuracy and performance of this framework is studied both numerically and analytically, specifically focusing on effects of TPU-native single precision floating point arithmetic on solution accuracy. The algorithm and implementation are validated with canonical 2D and 3D Taylor Green vortex simulations. To demonstrate the capability for simulating turbulent flows, simulations are conducted for two configurations, namely the decaying homogeneous isotropic turbulence and a turbulent planar jet. Both simulations show good statistical agreement with reference solutions. The performance analysis shows a linear weak scaling and a super-linear strong scaling up to a full TPU v3 pod with 2048 cores.
We introduce in this work the normalizing field flows (NFF) for learning random fields from scattered measurements. More precisely, we construct a bijective transformation (a normalizing flow characterizing by neural networks) between a Gaussian random field with the Karhunen-Lo`eve (KL) expansion structure and the target stochastic field, where the KL expansion coefficients and the invertible networks are trained by maximizing the sum of the log-likelihood on scattered measurements. This NFF model can be used to solve data-driven forward, inverse, and mixed forward/inverse stochastic partial differential equations in a unified framework. We demonstrate the capability of the proposed NFF model for learning Non Gaussian processes and different types of stochastic partial differential equations.
Mineral precipitation and dissolution processes in a porous medium can alter the structure of the medium at the scale of pores. Such changes make numerical simulations a challenging task as the geometry of the pores changes in time in an apriori unknown manner. To deal with such aspects, we here adopt a two-scale phase-field model, and propose a robust scheme for the numerical approximation of the solution. The scheme takes into account both the scale separation in the model, as well as the non-linear character of the model. After proving the convergence of the scheme, an adaptive two-scale strategy is incorporated, which improves the efficiency of the simulations. Numerical tests are presented, showing the efficiency and accuracy of the scheme in the presence of anisotropies and heterogeneities.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
