No abstract uploaded!

In this paper, we propose a rigorous and accurate non-local (in the oversampled region) upscaling framework based on some recently developed multiscale methods [10]. Our proposed method consists of identifying multi-continua parameters via local basis functions and constructing non-local (in the oversampled region) transfer and effective properties. To achieve this, we significantly modify our recent work proposed within Generalized Multiscale Finite Element Method (GMsFEM) in [10] and derive appropriate local problems in oversampled regions once we identify important modes representing each continuum. We use piecewise constant functions in each fracture network and in the matrix to write an upscaled equation. Thus, the resulting upscaled equation is of minimal size and the unknowns are average pressures in the fractures and the matrix. Note that the use of non-local upscaled model for porous

The bifurcations of penetrative Rayleigh-Benard convection in cylindrical containers are studied by the linear stability analysis (LSA) combined with the direct numerical simulation (DNS) method. The working fluid is cold water near 4 degrees C, where the Prandtl number Pr is 11.57, and the aspect ratio (radius/height) of the cylinder ranges from 0.66 to 2. It is found that the critical Rayleigh number increases with the increase in the density inversion parameter (m). The relationship between the normalized critical Rayleigh number (Ra-c((m))/Ra-c(0)) and (m) is formulated, which is in good agreement with the stability results within a large range of (m). The aspect ratio has a minor effect on Ra-c((m))/Ra-c(0). The bifurcation processes based on the axisymmetric solutions are also investigated. The results show that the onset of axisymmetric convection occurs through a trans-critical bifurcation due to the top-bottom symmetry breaking of the present system. Moreover, two kinds of qualitatively different steady axisymmetric solutions are identified.

We introduce a unified particle framework which integrates the phase-field method with multi-material simulation to allow modeling of both liquids and solids, as well as phase transitions between them. A simple elastoplastic model is used to capture the behavior of various kinds of solids, including deformable bodies, granular materials, and cohesive soils. States of matter or phases, particularly liquids and solids, are modeled using the non-conservative Allen-Cahn equation. In contrast, materials—made of different substances—are advected by the conservative Cahn-Hilliard equation. The distributions of phases and materials are represented by a phase variable and a concentration variable, respectively, allowing us to represent commonly observed fluid-solid interactions. Our multi-phase, multi-material system is governed by a unified Helmholtz free energy density. This framework provides the first method in computer graphics capable of modeling a continuous interface between phases. It is versatile and can be readily used in many scenarios that are challenging to simulate. Examples are provided to demonstrate the capabilities and effectiveness of this approach.

In this numerical study on Rayleigh-Benard convection, we seek to improve the heat transfer by passive means. To this end we introduce a single tilted conductive barrier centered in an aspect ratio one cell, breaking the symmetry of the geometry and to channel the ascending hot and descending cold plumes. We study the global and local heat transfer and the flow organization for Rayleigh numbers 10(5) <= Ra <= 10(9) for a fixed Prandtl number of Pr = 4.3. We find that the global heat transfer can be enhanced up to 18%, and locally around 800%. The averaged Reynolds number is always decreased when a barrier is introduced, even for those cases where the global heat transfer is increased. We map the entire parameter space spanned by the orientation and the size of a single barrier for Ra = 10(8).