Nose micro-blowing for asymmetric vortices control on blunt-nose slender body at high angle of attack

The asymmetric vortices over blunt-nose slender body at high angles of attack result in random side force. In this paper, a nose micro-blowing technology is used to control the asymmetric flow. Pressure measurement and particle image velocimetry(PIV) experiments are conducted in a low-speed wind tunnel to research effects of jet flow rate on asymmetric vortices over blunt-nose slender body. The angle of attack of the model is fixed at 50?and the Reynolds number for the experiments is 1.6×10~5 based on diameter of aftbody. A blow hole(5 mm in diameter) on the nose is processed at circumferential angle θ_b= 90?and meridian angle γ_b= 20?with jet momentum ratio C_μ ranging from 5.30×10~(-7) to1.19×10~(-4). Tests are made under two kinds of perturbations. One is called single perturbation with only blow hole and the other is called combined perturbation consists of blow hole and additional granules set on nose. The results show that whether the model has the single perturbation or the combined one,the sectional side force of x/D = 3 varies in the same direction with the increasement of C_μ and remains stable when C_μ is greater than 3.29×10~(-6). But the stable force values are different according to various perturbations. The fact proves that the size and direction of the side force of blunt-nose slender body can be controlled by the nose micro-blowing.     

Computation of field structure and aerodynamic characteristics of delta wings at high angles of attack

A numerical investigation of the structure of the vortical flowfield over delta wings at high angles of attack in longitudinal and with small sideslip angle is presented. Three-dimensional Navier-Stokes numerical simulations were carried out to predict the complex leeward-side flowfield characteristics that are dominated by the effect of the breakdown of the leading-edge vortices. The methods that analyze the flowfield structure quantitatively were given by using flowfield data from the computational results. In the region before the vortex breakdown, the vortex axes are approximated as being straight line. As the angle of attack increases, the vortex axes are closer to the root chord, and farther away from the wing surface. Along the vortex axes, as the adverse pressure gradients occur, the axial velocity decreases, that is, A is negativee, so the vortex is unstable, and it is possible to breakdown. The occurrence of the breakdown results in the instability of lateral motion for a delta wing, and the lateral moment diverges after a small perturbation occurs at high angles of attack. However, after a critical angle of attack is reached the vortices breakdown completely at the wing apex, and the instability resulting from the vortex breakdown disappears.     

Numerical simulations of gas-particle flow behavior created by low-level rotary-winged aircraft flight over particle bed

The aerodynamics of gas-particle suspensions is simulated as an Euler-Euler two-fluid model in a revolving rotor over a particle bed. The interactions of collisions between the blade and particles and particle-particle interactions are modeled using the kinetic theory of granular flow(KTGF). The gas turbulence induced by the rotation of the rotor is modeled using the kg-εg model. The flow field of a revolving rotor is simulated using the multiple reference frame(MRF) method. The distributions of velocities, volume fractions, and gas pressure are predicted while the aircraft hovers at different altitudes.The gas pressure decreases from the hub to the tip of the blade, and it is higher at the pressure side than that at the suction side of the rotor. The turbulent kinetic energy of the gas increases toward the blade tip. The volume fraction of particles decreases as the hovering altitude increases. The simulated pressure coefficient is compared with that in experimental measurements.     

Experimental study of the boundary layer over an airfoil in plunging motion

This is an experimental study on the boundary layer over an airfoil under steady and unsteady conditions.It specifically deals with the effect of plunging oscillation on the laminar/turbulent characteristics of the boundary layer.The wind tunnel measurements involved surfacemounted hot-film sensors and boundary-layer rake.The experiments were conducted at Reynolds numbers of 0.42×10 6 to 0.84 × 10 6 and the reduced frequency was varied from 0.01 to 0.11.The results of the quasi-wall-shear stress as well as the boundary layer velocity profiles provided important information about the state of the boundary layer over the suction surface of the airfoil in both static and dynamic cases.For the static tests,boundary layer transition occurred through a laminar separation bubble.By increasing the angle of attack,disturbances and the transition location moved toward the leading edge.For the dynamic tests,earlier transition occurred with increasing rather than decreasing effective angle of attack.The mean angle of attack and the oscillating parameters significantly affected the state of the boundary layer.By increasing the reduced frequency,the boundary layer transition was promoted to the upstroke portion of the equivalent angle of attack,but the quasi skin friction coefficient was decreased.     

Side force control on slender body by self-excited oscillation flag

Strong asymmetrical vortices appear on the leeward of slender body at high angles of attack,which has very unfavorable effect on the stability and control of the aircraft.A method is developed to control the side force of slender body at high angles of attack,and is verified in wind tunnel.A thin-film triangular self-excited oscillation flag is fixed at the tip of the slender body model whose semi-apex angle is 10°.Side force is approximately linearly proportional to roll-setting angle of self-excited oscillation flag at high angles of attack,and the slop of fitting straight line obtained by the least square method is-0.158.The linear relationship between side force and roll-setting angle provides convenience for developing side force control law of slender body at high angles of attack.Experimental data shows that the side force coefficients vary linearly with roll-setting angles when a specific plastic self-excited oscillation flag is used as the control flag.The range of side force coefficient and roll-setting angle are,respectively,-3.2to 3.0 and-20° to 20°.The device is simple,effective,and is of great potential in engineering application.     

Effect of inclined ribs on heat transfer coefficient in stationary square channel

The main objective of this research is to study the effect of rib arrangement on the distributions of the local heat transfer coefficient in a stationary channel. In this study, the ribs with square cross section were used to place on two side walls for study. The rib height-to-hydraulic diameter ratio(e/Dh) and the rib pitchto-height(p/e) ratio were fixed at 0.133 and 10, respectively. Three different types of rib arrangement for inclined ribs, V-shaped ribs and inverted V-shaped ribs were investigated. The rib angle of attack(α)was varied from 30?to 90?for inclined ribs and 45?and 60?for both V-shaped and inverted V-shaped ribs, and compared at constant Reynolds number Re =30000. Thermal Liquid Crystal sheet was applied for evaluating the heat transfer distributions. The results showed that the average Nusselt number on surface with rib inclined angle at 60?, 45?, and 60?V-shaped ribs was improved up to about 20%, 25% and 30%higher than case of angle 90?and the rib inclined angle at 60?V-shaped ribs provided the highest Nusselt number covering largest area when compared to the other cases.     

Similarity solutions to viscous flow and heat transfer of nanofluid over nonlinearly stretching sheet

The boundary-layer flow and heat transfer in a viscous fluid containing metallic nanoparticles over a nonlinear stretching sheet are analyzed. The stretching velocity is assumed to vary as a power function of the distance from the origin. The governing partial differential equation and auxiliary conditions are reduced to coupled nonlinear ordinary differential equations with the appropriate corresponding auxiliary conditions. The resulting nonlinear ordinary differential equations (ODEs) are solved numerically. The effects of various relevant parameters, namely, the Eckert number Ec, the solid volume fraction of the nanoparticles φ, and the nonlinear stretching parameter n are discussed. The comparison with published results is also presented. Different types of nanoparticles are studied. It is shown that the behavior of the fluid flow changes with the change of the nanoparticles type.     

Effect of disturbances at inlet on hypersonic boundary layer transition on a blunt cone at small angle of attack

To investigate the effect of different disturbances in the upstream, we present numerical simulation of transition for a hypersonic boundary layer on a 5-degree half-angle blunt cone in a freestream with Mach number 6 at 1-degree angle of attack. Evolution of small disturbances is simulated to compare with the linear stability theory （LST）, indicating that LST can provide a good prediction on the growth rate of the disturbance. The effect of different disturbances on transition is investigated. Transition onset distributions along the azimuthal direction are obtained with two groups of disturbances of different frequencies. It shows that transition onset is relevant to frequencies and amplitudes of the disturbances at the inlet, and is decided by the amplitudes of most unstable waves at the inlet. According to the characteristics of environmental disturbances in most wind tunnels, we explain why transition occurs leeside-forward and windside-aft over a circular cone at an angle of attack. Moreover, the indentation phenomenon in the transition curve on the leeward is also revealed.     

MHD natural convection in porous media-filled enclosures

In this work, the magnetohydrodynamics （MHD） natural convection heat transfer problem inside a porous medium filled with inclined rectangular enclosures is investigated numerically. The boundary conditions selected on the enclosure are two adiabatic and two isothermal walls. The governing equations, continuity, and Forchheimer extension of the Darcy law and energy are transformed into dimensionless forms by using a set of suitable variables, and then solved by using a finite difference scheme. The governing parameters are the magnetic influence number, the Darcy Rayleigh number, the inclination angle, and the aspect ratio of the enclosure. It is found that the magnetic influence number and the inclination angle have pronounced effects on the fluid flow and heat transfer in porous media-filled enclosures.     

Comparative assessment of SAS and DES turbulence modeling for massively separated flows

Numerical studies of the flow past a circular cylinder at Reynolds number 1.4 × 10~5 and NACA0021airfoil at the angle of attack 60?have been carried out by scale-adaptive simulation (SAS) and detached eddy simulation (DES), in comparison with the existing experimental data. The new version of the model developed by Egorov and Menter is assessed, and advantages and disadvantages of the SAS simulation are analyzed in detail to provide guidance for industrial application in the future. Moreover, the mechanism of the scale-adaptive characteristics in separated regions is discussed, which is obscure in previous analyses. It is concluded that: the mean flow properties satisfactorily agree with the experimental results for the SAS simulation, although the prediction of the second order turbulent statistics in the near wake region is just reasonable. The SAS model can produce a larger magnitude of the turbulent kinetic energy in the recirculation bubble, and, consequently, a smaller recirculation region and a more rapid recovery of the mean velocity outside the recirculation region than the DES approach with the same grid resolution. The vortex shedding is slightly less irregular with the SAS model than with the DES approach,probably due to the higher dissipation of the SAS simulation under the condition of the coarse mesh.     

On the Justification of Viscoelastic Elliptic Membrane Shell Equations

We consider a family of linearly viscoelastic shells with thickness \(2\varepsilon\), clamped along their entire lateral face, all having the same middle surface \(S=\boldsymbol{\theta}(\bar{\omega})\subset \mathbb{R}^{3}\), where \(\omega\subset\mathbb{R}^{2}\) is a bounded and connected open set with a Lipschitz-continuous boundary \(\gamma\). We make an essential geometrical assumption on the middle surface \(S\), which is satisfied if \(\gamma\) and \(\boldsymbol{\theta}\) are smooth enough and \(S\) is uniformly elliptic. We show that, if the applied body force density is \(O(1)\) with respect to \(\varepsilon\) and surface tractions density is \(O(\varepsilon)\), the solution of the scaled variational problem in curvilinear coordinates, \(\boldsymbol{u}( \varepsilon)\), defined over the fixed domain \(\varOmega=\omega\times (-1,1)\) for each \(t\in[0,T]\), converges to a limit \(\boldsymbol{u}\) with \(u_{\alpha}(\varepsilon)\rightarrow u_{\alpha}\) in \(W^{1,2}(0,T,H ^{1}(\varOmega))\) and \(u_{3}(\varepsilon)\rightarrow u_{3}\) in \(W^{1,2}(0,T,L^{2}(\varOmega))\) as \(\varepsilon\to0\). Moreover, we prove that this limit is independent of the transverse variable. Furthermore, the average \(\bar{\boldsymbol{u}}= \frac{1}{2}\int_{-1}^{1} \boldsymbol{u}dx_{3}\), which belongs to the space \(W^{1,2}(0,T, V_{M}( \omega))\), where

$$V_{M}(\omega)=H^{1}_{0}(\omega)\times H^{1}_{0}(\omega)\times L ^{2}(\omega), $$

satisfies what we have identified as (scaled) two-dimensional equations of a viscoelastic membrane elliptic shell, which includes a long-term memory that takes into account previous deformations. We finally provide convergence results which justify those equations.

     

Uniaxial experimental study of the acoustic emission and deformation behavior of composite rock based on 3D digital image correlation (DIC)

In this paper, uniaxial compression tests were carried out on a series of composite rock specimens with different dip angles, which were made from two types of rock-like material with different strength. The acoustic emission technique was used to monitor the acoustic signal characteristics of composite rock specimens during the entire loading process. At the same time, an optical non-contact 3 D digital image correlation technique was used to study the evolution of axial strain field and the maximal strain field before and after the peak strength at different stress levels during the loading process. The effect of bedding plane inclination on the deformation and strength during uniaxial loading was analyzed. The methods of solving the elastic constants of hard and weak rock were described. The damage evolution process, deformation and failure mechanism, and failure mode during uniaxial loading were fully determined. The experimental results show that the θ = 0?–45?specimens had obvious plastic deformation during loading, and the brittleness of the θ = 60?–90?specimens gradually increased during the loading process. When the anisotropic angle θincreased from 0?to 90?, the peak strength, peak strain,and apparent elastic modulus all decreased initially and then increased. The failure mode of the composite rock specimen during uniaxial loading can be divided into three categories:tensile fracture across the discontinuities(θ = 0?–30?), slid-ing failure along the discontinuities(θ = 45?–75?), and tensile-split along the discontinuities(θ = 90?). The axial strain of the weak and hard rock layers in the composite rock specimen during the loading process was significantly different from that of the θ = 0?–45?specimens and was almost the same as that of the θ = 60?–90?specimens. As for the strain localization highlighted in the maximum principal strain field, the θ = 0?–30?specimens appeared in the rock matrix approximately parallel to the loading direction,while in the θ = 45?–90?specimens it appeared at the hard and weak rock layer interface.     

Long-time Dynamics of Resonant Weakly Nonlinear CGL Equations

Consider a weakly nonlinear CGL equation on the torus \(\mathbb {T}^d\):

$$\begin{aligned} u_t+i\Delta u=\epsilon [\mu (-1)^{m-1}\Delta ^{m} u+b|u|^{2p}u+ ic|u|^{2q}u]. \end{aligned}$$

(*)

Here \(u=u(t,x)\), \(x\in \mathbb {T}^d\), \(0<\epsilon <<1\), \(\mu \geqslant 0\), \(b,c\in \mathbb {R}\) and \(m,p,q\in \mathbb {N}\). Define \(I(u)=(I_{\mathbf {k}},\mathbf {k}\in \mathbb {Z}^d)\), where \(I_{\mathbf {k}}=v_{\mathbf {k}}\bar{v}_{\mathbf {k}}/2\) and \(v_{\mathbf {k}}\), \(\mathbf {k}\in \mathbb {Z}^d\), are the Fourier coefficients of the function \(u\) we give. Assume that the equation \((*)\) is well posed on time intervals of order \(\epsilon ^{-1}\) and its solutions have there a-priori bounds, independent of the small parameter. Let \(u(t,x)\) solve the equation \((*)\). If \(\epsilon \) is small enough, then for \(t\lesssim {\epsilon ^{-1}}\), the quantity \(I(u(t,x))\) can be well described by solutions of an effective equation:

$$\begin{aligned} u_t=\epsilon [\mu (-1)^{m-1}\Delta ^m u+ F(u)], \end{aligned}$$

where the term \(F(u)\) can be constructed through a kind of resonant averaging of the nonlinearity \(b|u|^{2p}+ ic|u|^{2q}u\).

     

Insight into Heterogeneity Effects in Methane Hydrate Dissociation via Pore-Scale Modeling

A large number (1253) of high-quality streaming potential coefficient (\(C_\mathrm{sp})\) measurements have been carried out on Berea, Boise, Fontainebleau, and Lochaline sandstones (the latter two including both detrital and authigenic overgrowth forms), as a function of pore fluid salinity (\(C_\mathrm{f})\) and rock microstructure. All samples were saturated with fully equilibrated aqueous solutions of NaCl (10\(^{-5}\) and 4.5 mol/dm\(^{3})\) upon which accurate measurements of their electrical conductivity and pH were taken. These \(C_\mathrm{sp}\) measurements represent about a fivefold increase in streaming potential data available in the literature, are consistent with the pre-existing 266 measurements, and have lower experimental uncertainties. The \(C_\mathrm{sp}\) measurements follow a pH-sensitive power law behaviour with respect to \(C_\mathrm{f}\) at medium salinities (\(C_\mathrm{sp} =-\,1.44\times 10^{-9} C_\mathrm{f}^{-\,1.127} \), units: V/Pa and mol/dm\(^{3})\) and show the effect of rock microstructure on the low salinity \(C_\mathrm{sp}\) clearly, producing a smaller decrease in \(C_\mathrm{sp}\) per decade reduction in \(C_\mathrm{f}\) for samples with (i) lower porosity, (ii) larger cementation exponents, (iii) smaller grain sizes (and hence pore and pore throat sizes), and (iv) larger surface conduction. The \(C_\mathrm{sp}\) measurements include 313 made at \(C_\mathrm{f} > 1\) mol/dm\(^{3}\), which confirm the limiting high salinity \(C_\mathrm{sp}\) behaviour noted by Vinogradov et al., which has been ascribed to the attainment of maximum charge density in the electrical double layer occurring when the Debye length approximates to the size of the hydrated metal ion. The zeta potential (\(\zeta \)) was calculated from each \(C_\mathrm{sp}\) measurement. It was found that \(\zeta \) is highly sensitive to pH but not sensitive to rock microstructure. It exhibits a pH-dependent logarithmic behaviour with respect to \(C_\mathrm{f}\) at low to medium salinities (\(\zeta =0.01133 \log _{10} \left( {C_\mathrm{f} } \right) +0.003505\), units: V and mol/dm\(^{3})\) and a limiting zeta potential (zeta potential offset) at high salinities of \({\zeta }_\mathrm{o} = -\,17.36\pm 5.11\) mV in the pH range 6–8, which is also pH dependent. The sensitivity of both \(C_\mathrm{sp}\) and \(\zeta \) to pH and of \(C_\mathrm{sp}\) to rock microstructure indicates that \(C_\mathrm{sp}\) and \(\zeta \) measurements can only be interpreted together with accurate and equilibrated measurements of pore fluid conductivity and pH and supporting microstructural and surface conduction measurements for each sample.     

Bifurcations of a creeping air–water flow in a conical container

This numerical study describes the eddy emergence and transformations in a slow steady axisymmetric air–water flow, driven by a rotating top disk in a vertical conical container. As water height \(H_{\mathrm{w}}\) and cone half-angle \(\beta \) vary, numerous flow metamorphoses occur. They are investigated for \(\beta =30^{\circ }, 45^{\circ }\), and \(60^{\circ }\). For small \(H_{\mathrm{w}}\), the air flow is multi-cellular with clockwise meridional circulation near the disk. The air flow becomes one cellular as \(H_{\mathrm{w}}\) exceeds a threshold depending on \(\beta \). For all \(\beta \), the water flow has an unbounded number of eddies whose size and strength diminish as the cone apex is approached. As the water level becomes close to the disk, the outmost water eddy with clockwise meridional circulation expands, reaches the interface, and induces a thin layer with anticlockwise circulation in the air. Then this layer expands and occupies the entire air domain. The physical reasons for the flow transformations are provided. The results are of fundamental interest and can be relevant for aerial bioreactors.     

Characterization for elastic constants of fused deposition modelling-fabricated materials based on the virtual fields method and digital image correlation

Fused deposition modelling (FDM), a widely used rapid prototyping process, is a promising technique in manufacturing engineering. In this work, a method for characterizing elastic constants of FDM-fabricated materials is proposed. First of all, according to the manufacturing process of FDM, orthotropic constitutive model is used to describe the mechanical behavior. Then the virtual fields method (VFM) is applied to characterize all the mechanical parameters \((Q_{11}\), \(Q_{22}\), \(Q_{12}\), \(Q_{66})\) using the full-field strain, which is measured by digital image correlation (DIC). Since the principal axis of the FDM-fabricated structure is sometimes unknown due to the complexity of the manufacturing process, a disk in diametrical compression is used as the load configuration so that the loading angle can be changed conveniently. To verify the feasibility of the proposed method, finite element method (FEM) simulation is conducted to obtain the strain field of the disk. The simulation results show that higher accuracy can be achieved when the loading angle is close to \(30^{\circ }\). Finally, a disk fabricated by FDM was used for the experiment. By rotating the disk, several tests with different loading angles were conducted. To determine the position of the principal axis in each test, two groups of parameters \((Q_{11}\), \(Q_{22}\), \(Q_{12}\), \(Q_{66})\) are calculated by two different groups of virtual fields. Then the corresponding loading angle can be determined by minimizing the deviation between two groups of the parameters. After that, the four constants \((Q_{11}\), \(Q_{22}\), \(Q_{12}\), \(Q_{66})\) were determined from the test with an angle of \(27^{\circ }\).     

Influence of hydrodynamic slip on convective transport in flow past a circular cylinder

The presence of a finite tangential velocity on a hydrodynamically slipping surface is known to reduce vorticity production in bluff body flows substantially while at the same time enhancing its convection downstream and into the wake. Here, we investigate the effect of hydrodynamic slippage on the convective heat transfer (scalar transport) from a heated isothermal circular cylinder placed in a uniform cross-flow of an incompressible fluid through analytical and simulation techniques. At low Reynolds (\({\textit{Re}}\ll 1\)) and high Péclet (\({\textit{Pe}}\gg 1\)) numbers, our theoretical analysis based on Oseen and thermal boundary layer equations allows for an explicit determination of the dependence of the thermal transport on the non-dimensional slip length \(l_s\). In this case, the surface-averaged Nusselt number, Nu transitions gradually between the asymptotic limits of \(Nu \sim {\textit{Pe}}^{1/3}\) and \(Nu \sim {\textit{Pe}}^{1/2}\) for no-slip (\(l_s \rightarrow 0\)) and shear-free (\(l_s \rightarrow \infty \)) boundaries, respectively. Boundary layer analysis also shows that the scaling \(Nu \sim {\textit{Pe}}^{1/2}\) holds for a shear-free cylinder surface in the asymptotic limit of \({\textit{Re}}\gg 1\) so that the corresponding heat transfer rate becomes independent of the fluid viscosity. At finite \({\textit{Re}}\), results from our two-dimensional simulations confirm the scaling \(Nu \sim {\textit{Pe}}^{1/2}\) for a shear-free boundary over the range \(0.1 \le {\textit{Re}}\le 10^3\) and \(0.1\le {\textit{Pr}}\le 10\). A gradual transition from the lower asymptotic limit corresponding to a no-slip surface, to the upper limit for a shear-free boundary, with \(l_s\), is observed in both the maximum slip velocity and the Nu. The local time-averaged Nusselt number \(Nu_{\theta }\) for a shear-free surface exceeds the one for a no-slip surface all along the cylinder boundary except over the downstream portion where unsteady separation and flow reversal lead to an appreciable rise in the local heat transfer rates, especially at high \({\textit{Re}}\) and Pr. At a Reynolds number of \(10^3\), the formation of secondary recirculating eddy pairs results in appearance of additional local maxima in \(Nu_{\theta }\) at locations that are in close proximity to the mean secondary stagnation points. As a consequence, Nu exhibits a non-monotonic variation with \(l_s\) increasing initially from its lowermost value for a no-slip surface and then decreasing before rising gradually toward the upper asymptotic limit for a shear-free cylinder. A non-monotonic dependence of the spanwise-averaged Nu on \(l_s\) is observed in three dimensions as well with the three-dimensional wake instabilities that appear at sufficiently low \(l_s\), strongly influencing the convective thermal transport from the cylinder. The analogy between heat transfer and single-component mass transfer implies that our results can directly be applied to determine the dependency of convective mass transfer of a single solute on hydrodynamic slip length in similar configurations through straightforward replacement of Nu and \({\textit{Pr}}\) with Sherwood and Schmidt numbers, respectively.     

Further Dense Properties of the Space of Circle Diffeomorphisms with a Liouville Rotation Number

In continuation of Matsumoto’s paper (Nonlinearity 25:1495–1511, 2012) we show that various subspaces are \(C^{\infty }\)-dense in the space of orientation-preserving \(C^{\infty }\)-diffeomorphisms of the circle with rotation number \(\alpha \), where \(\alpha \in {\mathbb {S}}^1\) is any prescribed Liouville number. In particular, for every odometer \({\mathcal {O}}\) of product type we prove the denseness of the subspace of diffeomorphisms which are orbit-equivalent to \({\mathcal {O}}\).     

Approximation of Flat <Emphasis Type="Italic">W</Emphasis><Superscript>2,2</Superscript> Isometric Immersions by Smooth Ones

Let \({S\subset\mathbb{R}^2}\) be a bounded Lipschitz domain and denote by \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\) the set of mappings \({u\in W^{2,2}(S;\mathbb{R}^3)}\) which satisfy \({(\nabla u)^T(\nabla u) = Id}\) almost everywhere. Under an additional regularity condition on the boundary \({\partial S}\) (which is satisfied if \({\partial S}\) is piecewise continuously differentiable), we prove that the strong W ^2,2 closure of \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)\cap C^{\infty}(\overline{S};\mathbb{R}^3)}\) agrees with \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\).     

Displacement and Dissolution Characteristics of CO$$_{}$$/Brine System in Unconsolidated Porous Media

This study investigated the dynamic displacement and dissolution of \(\hbox {CO}_{2}\) in porous media at 313 K and 6/8 MPa. Gaseous (\(\hbox {gCO}_{2}\)) at 6 MPa and supercritical \(\hbox {CO}_{2 }(\hbox {scCO}_{2}) \) at 8 MPa were injected downward into a glass bead pack at different flow rates, following upwards brine injection. The processes occurring during \(\hbox {CO}_{2}\) drainage and brine imbibition were visualized using magnetic resonance imaging. The drainage flow fronts were strongly influenced by the flow rates, resulting in different gas distributions. However, brine imbibition proceeded as a vertical compacted front due to the strong effect of gravity. Additionally, the effects of flow rate on distribution and saturation were analyzed. Then, the front movement of \(\hbox {CO}_{2}\) dissolution was visualized along different paths after imbibition. The determined \(\hbox {CO}_{2}\) concentrations implied that little \(\hbox {scCO}_{2}\) dissolved in brine after imbibition. The dissolution rate was from \(10^{-8}\) to \(10^{-9}\, \hbox {kg}\, \hbox {m}^{-3} \, \hbox {s}^{-1}\) and from \(10^{-6}\) to \(10^{-8}\, \hbox {kg}\, \hbox {m}^{-3} \, \hbox {s}^{-1}\) for \(\hbox {gCO}_{2}\) at 6 MPa and \(\hbox {scCO}_{2 }\) at 8 MPa, respectively. The total time for the \(\hbox {scCO}_{2}\) dissolution was short, indicating fast mass transfer between the \(\hbox {CO}_{2}\) and brine. Injection of \(\hbox {CO}_{2}\) under supercritical conditions resulted in a quick establishment of a steady state with high storage safety.