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

We study the translational and rotational dynamics of neutrally buoyant finite-size spheroids in hydrodynamic turbulence by means of fully resolved numerical simulations. We examine axisymmetric shapes, from oblate to prolate, and the particle volume dependences. We show that the accelerations and rotations experienced by non-spherical inertial-scale particles result from volume filtered fluid forces and torques, similar to spherical particles. However, the particle orientations carry signatures of preferential alignments with the surrounding flow structures, which are reflected in distinct axial and lateral fluctuations for accelerations and rotation rates. The randomization of orientations does not occur even for particles with volume-equivalent diameter size in the inertial range, here up to 60 dissipative units (eta) at Taylor-scale Reynolds number Re-lambda = 120. Additionally, we demonstrate that the role of fluid boundary layers around the particles cannot be neglected in reaching a quantitative understanding of particle statistical dynamics, as they affect the intensities of the angular velocities and the relative importance of tumbling with respect to spinning rotations. This study brings to the fore the importance of inertial-scale flow structures in homogeneous and isotropic turbulence and their impacts on the transport of neutrally buoyant bodies with sizes in the inertial range.

For a fixed geometric configuration, hydrodynamic instabilities and bifurcation processes of laminar isothermal planar opposed jet flows with symmetric and slightly asymmetric inlet boundary conditions are investigated numerically using a high-resolution approach based on spectral element method. In current configuration, when inlet boundary conditions are symmetric, in the range of the Reynolds number considered (Re <= 200), multiple new symmetry-breaking bifurcations are observed and four new flow patterns are identified. Their hydrodynamic characteristics are analyzed, in particular their symmetries. In addition, the case that inlet boundary conditions are slightly asymmetric is investigated. It is found that bifurcation processes are extremely sensitive to this small symmetry-breaking imperfection and much different from those in the symmetric case. Furthermore, model equations are constructed by symmetry consideration to explain the numerical results based on hydrodynamic equations.