11,082 research outputs found
Numerical simulation of liquid sloshing in a partially filled container with inclusion of compressibility effects
A numerical scheme of study is developed to model compressible two-fluid flows simulating liquid sloshing in a partially filled tank. For a two-fluid system separated by an interface as in the case of sloshing, not only a Mach-uniform scheme is required, but also an effective way to eliminate unphysical numerical oscillations near the interface. By introducing a preconditioner, the governing equations expressed in terms of primitive variables are solved for both fluids (i.e. water, air, gas etc.) in a unified manner. In order to keep the interface sharp and to eliminate unphysical numerical oscillations in unsteady fluid flows, the non-conservative implicit Split Coefficient Matrix Method (SCMM) is modified to construct a flux difference splitting scheme in the dual time formulation. The proposed numerical model is evaluated by comparisons between numerical results and measured data for sloshing in an 80% filled rectangular tank excited at resonance frequency. Through similar comparisons, the investigation is further extended by examining sloshing flows excited by forced sway motions in two different rectangular tanks with 20% and 83% filling ratios. These examples demonstrate that the proposed method is suitable to capture induced free surface waves and to evaluate sloshing pressure loads acting on the tank walls and ceiling
MFC: An open-source high-order multi-component, multi-phase, and multi-scale compressible flow solver
MFC is an open-source tool for solving multi-component, multi-phase, and bubbly compressible flows. It is capable of efficiently solving a wide range of flows, including droplet atomization, shock–bubble interaction, and bubble dynamics. We present the 5- and 6-equation thermodynamically-consistent diffuse-interface models we use to handle such flows, which are coupled to high-order interface-capturing methods, HLL-type Riemann solvers, and TVD time-integration schemes that are capable of simulating unsteady flows with strong shocks. The numerical methods are implemented in a flexible, modular framework that is amenable to future development. The methods we employ are validated via comparisons to experimental results for shock–bubble, shock–droplet, and shock–water-cylinder interaction problems and verified to be free of spurious oscillations for material-interface advection and gas–liquid Riemann problems. For smooth solutions, such as the advection of an isentropic vortex, the methods are verified to be high-order accurate. Illustrative examples involving shock–bubble-vessel-wall and acoustic–bubble-net interactions are used to demonstrate the full capabilities of MFC
A Moving Boundary Flux Stabilization Method for Cartesian Cut-Cell Grids using Directional Operator Splitting
An explicit moving boundary method for the numerical solution of
time-dependent hyperbolic conservation laws on grids produced by the
intersection of complex geometries with a regular Cartesian grid is presented.
As it employs directional operator splitting, implementation of the scheme is
rather straightforward. Extending the method for static walls from Klein et
al., Phil. Trans. Roy. Soc., A367, no. 1907, 4559-4575 (2009), the scheme
calculates fluxes needed for a conservative update of the near-wall cut-cells
as linear combinations of standard fluxes from a one-dimensional extended
stencil. Here the standard fluxes are those obtained without regard to the
small sub-cell problem, and the linear combination weights involve detailed
information regarding the cut-cell geometry. This linear combination of
standard fluxes stabilizes the updates such that the time-step yielding
marginal stability for arbitrarily small cut-cells is of the same order as that
for regular cells. Moreover, it renders the approach compatible with a wide
range of existing numerical flux-approximation methods. The scheme is extended
here to time dependent rigid boundaries by reformulating the linear combination
weights of the stabilizing flux stencil to account for the time dependence of
cut-cell volume and interface area fractions. The two-dimensional tests
discussed include advection in a channel oriented at an oblique angle to the
Cartesian computational mesh, cylinders with circular and triangular
cross-section passing through a stationary shock wave, a piston moving through
an open-ended shock tube, and the flow around an oscillating NACA 0012 aerofoil
profile.Comment: 30 pages, 27 figures, 3 table
- …