Large Vortex Research Articles

О неэффективности нагнетательного способа проветривания при решении задачи нормализации пылевой обстановки в тупиковом забое калийного рудника

When potash ore is mined, a large amount of the harmful impurities, including dust, are released into the mine environment. Although potash dust is not toxic, at constant high concentrations it adversely affects the miner health. Hygroscopic properties of rocks in the potash mines do not allow the use of water to control dust. Ventilation approaches to the normalization of the dust environment also do not give positive results. Forced method of the ventilation of the dead-end mine workings used in all the mines in Russia is distinguished by the long-range air jet from the ventilation stack. This allows to create intensive vortex formation in the bottomhole part of the mine working for avoiding the occurrence of stagnant zones. However, during forced ventilation, a large vortex occurs. Results of studying the dynamics of dusty air streams in the space of a dead-end tunneling and cleaning working are presented in the article. In the mining face, there are two main sources of dust formation: the working face, where the rock mass is destroyed, and the bunker-loader, where the broken ore is poured. Both the first and second sources are identical in terms of dust formation power. The vortex formed during the forced ventilation method returns to the combine operator cab through the ore pouring area into the bunker. As a result, the composition of the atmosphere in the bottom hole significantly deteriorates. The greater the amount of air supplied during the forced ventilation method, the greater the dustiness of the bottomhole part of the mine. Conducted studies allow to conclude that the injection method of the ventilation of dead-end combine workings of the potash mines is ineffective to normalize the dust environment.

Numerical Study on the Shear Stress Characteristics of Open-Channel Flow over Rough Beds

Bed shear stress is an important measure of benthic habitats since it is related to many ecological processes. In this study, we focused on the fluctuating characteristics of shear stress in rough-bed open-channel flows. The roughness element method was adopted to mimic natural rough beds and the Improved Delayed Detached Eddy Simulation (IDDES) model was used to obtain comprehensive information about shear stress near the rough bed. Three arrangement patterns of the roughness elements were simulated to compare their effects on flow structure and shear stresses. The arrangements of the roughness elements altered the Reynold stress and turbulent kinetic energy characteristics, due to the variance of blockage in lateral directions that led to flow detachment and changes in the flow directions. Quadrant analysis revealed the spatial variations of the instantaneous shear stress burst events at different locations in the wake. By using spectrum analysis, the accumulation of shear-stress energy from small to large vortex scales was estimated, which revealed that the instantaneous effect of the shear stress was significantly stronger than the effect of the time-averaged shear stress, especially on small scales. The results of this study suggest the significance of the fluctuation part of shear stress in further studies on ecological processes.

Vortex-induced vibration characteristics of two open girders: A comparison of experimental and numerical investigation

Open girder sections are vulnerable to vortex-induced vibrations (VIV) because of their bluff aerodynamic characteristics. A comparative investigation was conducted by wind tunnel tests and numerical simulation. Firstly, based on sectional model wind tunnel tests, the VIV performance of semi-open girder and separated edge-boxes open girder was studied with consideration of the influence of the equivalent mass, wind attack angle and damping ratio. The aerodynamic parameters Scruton number ( S c) and Strouhal number ( S t) of two girders were analyzed contrastively. In addition, VIV amplitudes corresponding to the real bridge girders were calculated by the linear and nonlinear theories. Then, flow fields around two semi-open girder sections were investigated on the basis of computational fluid dynamics. The results show that vertical and torsional VIVs of two open girders are observed, at +3°and +5° wind attack angles. There were two vertical VIV regions and the maximum amplitude in the second vertical VIV region is significantly larger than that in the first one. For the semi-open girder and separated edged-box open girder, the vertical VIV amplitudes at +5° wind attack angle is 115% and 75% larger than at +3° wind attack angle, respectively. The damping ratio could obviously mitigate the VIV of two girders and VIV amplitudes are decreased linearly with the S C number increases. There are vortexes above the girder deck at the opening position of the semi-open girder and the separated edge-boxes open girder. The oblique web and wind faring may break the large vortex into several smaller vortexes at the opening position, thus optimizing the vortex vibration performance of the girder. A thicker shear layer is formed at separated open-girder during the development of airflow towards backward of bridge deck, due to the effect of railings and the blunt boxes on the lower surface.

Numerical Investigation of 3D Flow Properties around Finite Emergent Vegetation by Using the Two-Phase Volume of Fluid (VOF) Modeling Technique

This study predicts how the Free Surface Level (FSL) variations around finite length vegetation affect flow structure by using a numerical simulation. The volume of fluid (VOF) technique with the Reynolds stress model (RSM) was used for the simulation. Multizone Hexahedral meshing was adopted to accurately track the free surface level with minimum numerical diffusion at the water–air interface. After the validation, finite length emergent vegetation patches were selected based on the aspect ratio (AR = vegetation width-length ratio) under constant subcritical flow conditions for an inland tsunami flow. The results showed that the generation of large vortices was predominated in wider vegetation patches (AR > 1) due to the increase and decrease in the FSL at the front and back of the vegetation compared to longer vegetation patches (AR ≤ 1), as this offered more resistance against the approaching flow. The wider vegetation patches (AR > 1) are favorable in terms of generating a large area of low velocity compared to the longer vegetation patch (AR < 1) directly downstream of the vegetation patch. On the other hand, it has a negative impact on the adjacent downstream gap region, where a 14.3–34.9% increase in velocity was observed. The longer vegetation patches (AR < 1) generate optimal conditions within the vegetation region due to great velocity reduction. Moreover, in all the AR vegetation cases, the water turbulent intensity was maximum in the vegetation region compared to the adjacent gap region and air turbulent intensity above the FSL, suggesting strong air entrainment over this region. The results of this study are important in constructing vegetation layouts based on the AR of the vegetation for tsunami mitigation.

A Study on Swirling Jets Flow Controlled by Two Tandem Bodies

This study presents the effects of using two centerbodies arranged in tandem on the flow and velocity properties of swirling jets. The centerbodies are installed concentrically downstream of the jets to modify the wake. Smoke flow visualization was employed to illustrate the flow behaviors, while the particle image velocimetry technique was utilized to investigate the velocity fields, turbulent properties, vorticities, and Lagrangian time and length scales. The toroidal recirculation and vortex shedding are found near the field using one centerbody. The recirculation zone is evolved from the edge of the centerbody, while the vortex shedding is formed by the interaction between the central and the annular jets. Since two centerbodies are installed, two four-way saddles are located near the rim of the upstream centerbody to associate two pairs of vortices. A large momentum of the annular jet dominates the flow field; therefore, the central jet forms dual-rotating vortices at downstream centerbody. The turbulent intensity and vorticity along the central axis in the case of two centerbodies are extremely lower than those in the case of one centerbody. Large turbulent intensity and vorticity are located around the shear layers of jets. This promotes a better mixing efficiency of the swirling jets.

Simultaneous velocity–density measurements of downslope density clouds

We present results from laboratory experiments of the downslope propagation of a finite volume released gravity current by means of combined PIV/PLIF measurements. The experimental data were used to estimate the global characteristics of the current, such as the propagation speed, the lateral surface and the buoyancy, revealing that thermal theory is a robust model that can predict such quantities properly, especially the models of Beghin et al. (1981), Maxworthy (2010) and the theoretical model of Dai (2013a). The collected data allowed also for the determination of the entrainment rates based on the turbulent fluxes instead using the variation of the volume flux, and the regions of entrainment/detrainment have been identified. The results support previous observations of the mechanism responsible for the buoyancy loss of the cloud and its consequent deceleration down the slope, as a large scale recirculation vortex at the back of the cloud. A very complex interior structure of the cloud is observed, with a large variety of turbulent processes taking place, such as the large-scale convectively unstable recirculation vortex at the scale of the current itself, the small-scale convectively unstable motions inside the head of the cloud and close to the bottom boundary and shear (Kelvin–Helmholtz) instabilities at the boundary between the current and the ambient water. The synoptic velocity and density measurements allowed also to test existent parametrizations of turbulent fluxes, that have been quantified in 2D fields with a high spatial resolution. Results confirm that parametrization laws based on the assumption of a constant turbulent diffusivity or mixing length do not apply for buoyancy clouds due to their high spatial heterogeneity. Hence, the parametrization of the turbulent diffusivities in such flows should be based on scalar quantities that avoid the problem of spatial heterogeneity and takes into account the different sources of turbulence production, e.g. using energetic considerations that compare the terms having a definite exchange of energy and acts as a source or a sink.

Existence of stationary vortex sheets for the 2D incompressible Euler equation

AbstractWe construct a new type of planar Euler flows with localized vorticity. Let $\kappa _i\not =0$ , $i=1,\ldots , m$ , be m arbitrarily fixed constants. For any given nondegenerate critical point $\mathbf {x}_0=(x_{0,1},\ldots ,x_{0,m})$ of the Kirchhoff–Routh function defined on $\Omega ^m$ corresponding to $(\kappa _1,\ldots , \kappa _m)$ , we construct a family of stationary planar flows with vortex sheets that have large vorticity amplitude and concentrate on curves perturbed from small circles centered near $x_{0,i}$ , $i=1,\ldots ,m$ . The proof is accomplished via the implicit function theorem with suitable choice of function spaces.

On the nonlinear thermal stress, thermal creep, and thermal edge flows in triangular cavities

Here, rarefied thermally driven flow is investigated in two-dimensional equilateral triangular cavities with different uniform wall temperatures. We used three different solvers, i.e., the direct simulation Monte Carlo solver, discrete unified gas kinetic scheme solver, and continuum set of equations of a slow non-isothermal flow solver. Two main cases were considered; in the first case, the cavity's base is considered hot, and the other sides were set cold. In the second case, the right half of the bottom wall was regarded as a diffuse reflector with high temperature, while the left half of the bottom border was set as a specular reflector. The adjacent side walls were set cold with diffuse reflector boundary conditions. The imposed temperature difference/wall boundary condition induces various vortices in the geometry. In case 1, we observe that principal vortices appearing in the triangle are due to nonlinear thermal stress effects, and the thermal creep effects cause other smaller, confined ones. In case 2, a thermal edge flow is set up from the specular wall on the way to the diffusive hot wall, creating a large vortex in the geometry. As the Knudsen number decreases, another small vortex appears near the left cold border.

Wind Pressure Characteristics Based on the Rise–Span Ratio of Spherical Domes with Openings on the Roof

Wind loads are a primary concern in dome roof structures with openings such as retractable dome roofs. This is because the openings can cause damage to the cladding owing to high internal pressure. In this study, the wind pressure characteristics of a dome with an opening that varied based on the opening, rise–span ratio, and height span were examined by comparing the results from wind tunnel tests with those from previous studies. The negative pressure dominated the internal pressure of the roof in all regions and was not significantly affected by changes in the rise–span and height–span ratios. The reattachment distance of the windward region increased as the rise–span ratio increased, increasing the negative net pressure and decreasing the positive net pressure owing to a relatively large vortex. The roof inclination angle of the leeward region decreased as the rise–span ratio decreased, resulting in a decrease in the negative net pressure and an increase in the positive net pressure owing to a relatively small vortex. Based on the experimental results, a peak net pressure coefficient for cladding design was proposed for an open dome roof with a rise–span ratio of 0.05.

Hotspots of the stokes rotating circulation in a large marginal sea

Marginal seas, surrounded by continents with dense populations, are vulnerable and have a quick response to climate change effects. The seas typically have alternatively rotating layered circulations to regulate regional heat and biogeochemical transports. The circulations are composed of dynamically active hotspots and governed by the couplings between unique extrinsic inflow and intrinsic dynamic response. Ambiguities about the circulations’ structure, composition, and physics still exist, and these ambiguities have led to poor numerical simulation of the marginal sea in global models. The South China Sea is an outstanding example of a marginal sea that has this typical rotating circulation. Our study demonstrates that the rotating circulation is structured by energetic hotspots with large vorticity arising from unique dynamics in the marginal sea and is identifiable by the constraints of Stokes Theorem. These hotspots contribute most of the vorticity and most of energy needed to form and maintain the rotating circulation pattern. Our findings provide new insights on the distinguishing features of the rotating circulation and the dominant physics with the objectives of advancing our knowledge and improving modeling of marginal seas.

A global climatology of polar lows investigated for local differences and wind-shear environments

Abstract. Polar lows are intense mesoscale cyclones developing in marine polar air masses. This study presents a new global climatology of polar lows based on the ERA5 reanalysis for the years 1979–2020. Criteria for the detection of polar lows are derived based on a comparison of five polar-low archives with cyclones derived by a mesoscale tracking algorithm. The characteristics associated with polar lows are considered by the following criteria: (i) intense cyclone (large relative vorticity), (ii) mesoscale (small vortex diameter), and (iii) development in the marine polar air masses (a combination of low potential static stability and low potential temperature at the tropopause). Polar lows develop in all marine areas adjacent to sea ice or cold landmasses, mainly in the winter half year. The length and intensity of the season are regionally dependent. The highest density appears in the Nordic Seas. For all ocean sub-basins, forward-shear polar lows are the most common, whereas weak-shear polar lows and those propagating towards warmer environments are second and third most frequent, depending on the area. Reverse-shear polar lows and those propagating towards colder environments are rather seldom, especially in the Southern Ocean. Generally, polar lows share many characteristics across ocean basins and wind-shear categories. The most remarkable difference is that forward-shear polar lows often occur in a stronger vertical wind shear, whereas reverse-shear polar lows feature lower static stability. Hence, the contribution to a fast baroclinic growth rate is slightly different for the shear categories.

Numerical Simulation of Turbulent Flow in Eccentric Co-Rotating Heat Transfer

Heat transfer engineering is significant in many applications, especially in buoyancy natural convection in concentric and eccentric cavities. The biggest practical challenges, in this context, are capturing the self-natural flow, estimating the mixing performance, and determining what parameters affect the temperature distribution in the cavity. In this paper, we focus on the improvement of a mathematical model, in order to enhance the accuracy of the solution, by investigating a new source term in the SST k−ω turbulence model based on the finite volume technique. The commercial numerical simulation software ANSYS Fluent 2021R1 is implemented to validate the accuracy. A concentric cavity was chosen for validation, the obtained temperature profiles at θ=0°, θ=30°, θ=60°, θ=90°, θ=120°, θ=150°, and θ=180° were compared with previous experimental data. We applied this model to four eccentric rotating scenarios, including inner counterclockwise rotation, outer counterclockwise rotation, inner–outer clockwise rotation, and inner clockwise–outer counterclockwise rotation. The numerical simulation results reveal that the new source term in the momentum equation can produce superior results in the concentric test-case. The proposed mathematical model can describe the heat transfer under the eccentric co-rotation scenario well. Furthermore, the results for eccentric cases confirm that the rotational direction affects the mixing temperature by generating a large vortex in the cavity, which increases the temperature mixing performance.

Reducing labyrinth seal leakage with application of axially arranged baffles with increasing height

Abstract For turbomachines, reducing seal leakage is beneficial for increasing the efficiency and reducing the energy consumption. This study proposes the design of the axially arranged baffles with increasing height to reduce the leakage of the labyrinth seal. Due to the existence of the axially arranged baffles with increasing height, the leakage fluid is guided to flow into the seal cavities. In each seal cavity, one large vortex is induced, and additional small vortexes are induced in the regions between adjacent baffles. The leakage flow resistance and kinetic energy dissipation of the leakage fluid are enhanced. To demonstrate the effectiveness of the new design, numerical comparisons are conducted between the labyrinth seals with and without axially arranged baffles with increasing height. Results show that the new seal can reduce the leakage flow rate by 37.7-38.8% compared with the labyrinth seal.

Potential vorticity fronts and the late-time evolution of large-scale quasi-geostrophic flows

The late-time behaviour of freely evolving quasi-geostrophic flows with initial characteristic length scale $L_0$ larger than or equal to the deformation radius $L_D$ , $L_0/L_D \geq 1$ , is studied. At late time the flows are dominated by large multi-level vortices consisting of ascending terraces of well-mixed potential vorticity (PV), i.e. PV staircases. We examine how the number of mixed PV levels depends on the initial conditions, in particular $L_0/L_D$ . For sufficiently large values of $L_0/L_D \approx 5$ , a complete staircase with regular steps forms, but as $L_0/L_D$ decreases, the staircase becomes more irregular, with fewer mixed levels and the appearance of a large step centred on zero PV, corresponding to large regions of near-zero PV separating the multi-level vortices. This occurs because weak PV features in the initial field with scales smaller than $L_D$ undergo filamentation and are coarse-grained away or homogenised. For all values of $L_0/L_D$ considered, inverse cascades of potential energy commence at sufficiently late times. The onset of these cascades, even when the flow is initialised well within the ‘asymptotic model’ (AM) regime, suggests that the AM regime is not self-consistent: when potential vorticity fronts are well-resolved, frontal dynamics eventually drive ongoing flow evolution.

Numerical simulation of unsteady flow characteristics in a centrifugal pump during runaway process caused by power failure

Abstract Centrifugal pump is widely used in agricultural irrigation, petrochemical, mining, metallurgy and other industries. A special centrifugal pump with elbow pipe is researched in this paper. The flow rate and rotating speed change rapidly, especially in elbow pipe, causing serious pressure fluctuation and complex vorticity in the centrifugal pump. A user-defined function to solve the angular momentum equation was calculated the pump rotational speed. Unsteady numerical simulation with SST turbulence model is carried out and Short time Fourier transform (STFT) and vorticity strength are used to analyse the pressure fluctuation and vortexes respectively. The results indicate that the rotational speed first decreases to zero and then reversely increases to runaway speed. Large area vorticity appeared at the elbow pipe when the high amplitude frequency of pressure fluctuation increased. The change law of main frequency of pressure fluctuation at the elbow pipe is consistent with the rotational speed change ruler, which is equal to blade passing frequency (BPF). The main source of pressure fluctuation was rotor stator interaction. This study’s results provide a meaningful reference about pump runaway process caused by power failure, leading to the prevention of dangerous situations, for safe operations of the pumping unit.

Vortex dynamics of tandem bare and spiral finned cylinders in cross-flow and their susceptibility to acoustic resonance excitation

The vortex dynamics of tandem bare and spiral finned cylinders in the cross flow are experimentally investigated at a spacing ratio of L/Deq=2. Three different fin pitch-to-root diameter ratios (p/Dr) are considered. For high p/Dr, the fins on the downstream cylinder inhibit the shear layer growth and flapping in the cylinders' gap. Shear layer flapping is only observed for finned cylinders with the lowest p/Dr. Moreover, the gap shear layers enclose the downstream cylinder, enabling the formation of large vortex cores. The flow characteristics of finned cylinders are found to be independent of the Reynolds number with higher values of the Strouhal number due to an enlarged vortex formation length. For tandem bare cylinders, acoustic resonance excitation is triggered by two sources: (1) shear layer flapping in the cylinders' gap and (2) vortex shedding in the wake of the cylinders. All the tandem finned cylinders were able to excite acoustic resonance by the latter source. However, only tandem finned cylinders with the lowest p/Dr were capable of exciting resonance by the former source. During this excitation, weaker shear layer flapping resulted in the formation of smaller vortex cores in the finned cylinders' gap. During the second resonance, single vortex pairs formed per cycle in the wake of the downstream finned cylinder, whereas two vortex pairs formed per cycle in the wake of the downstream bare cylinder. This work shows that the addition of fins changes the impinging flow mechanism and the flow topology, which cannot be captured using the equivalent diameter approach.

Mesoscopic Fluid-Particle Flow and Vortex Structural Transmission in a Submerged Entry Nozzle of Continuous Caster.

Understanding the essence of the flow oscillations within a submerged-entry nozzle (SEN) is essential to control flow patterns in the continuous casting mold and consequently increase the superficial quality of steel products. A numerical study of the mesoscopic fluid-particle flow in a bifurcated pool-type SEN under steady operating conditions is conducted using the lattice Boltzmann method (LBM) coupled with the large eddy simulation (LES) model. The accuracy of the model has been verified by comparing vortex structures and simulated velocities with published experimental values. The LBM modeling is also verified by comparing the “stair-step” jet patterns observed in the experiment. The geometrical parameters and operational conditions of physical experiments are reproduced in the simulations. By comparing the time-averaged velocities of Reynolds-averaged Navier–Stokes equations (RANS) with LBM models, transient mesoscopic fluid-particles and related vortex structures can be better reproduced within the SEN. The visualization of internal flow within the SEN is illustrated through the mass-less Discrete Phase Model (DPM) model. The trajectories show that the LBM–LES–DPM coupled model is good at predicting the transient vortical flow within the SEN. A large vortex is found inside the exit port and continuously changes in shape and size therein. The monitoring points and lines within the SEN are selected to illustrate the velocity variations and effective viscosity, which can reflect the oscillating characteristics even under stable operating conditions without changes at the exit from the SEN. Furthermore, the formation, development, diffusion, and dissipation of the vortex structures from the exit port of the SEN are also investigated using the Q criteria. The comparison of the power spectrum with high-frequency components along the exit port indicates that the flow oscillations must originate from within the SEN and are intensified in the exit port. The mesoscopic LBM model can replicate the fluid-particle flow and vortex structure transmission as well as their turbulence effects inside the SEN in detail.

Experimental and Numerical Simulation Study on the Flow Characteristics of the Draft Tube in Francis Turbine

The flow characteristics of the draft tube of a Francis turbine have a significant influence on turbine stability. Numerical simulations were performed for a Francis turbine under three different output conditions of 20%, 100%, and 120% at the rated and maximum heads. Laser Doppler velocimetry (LDV) tests were conducted to test the flow characteristics of the draft tube of the Francis turbine model. The flow characteristics in the draft tube, the mechanism of the flow characteristics change, and the effect of the opening on the vortex rope were analyzed. The results showed that the large and invisible vortex in the conical cross-section at the inlet of the draft tube gradually changed to a tangible vortex rope as the guide vane opening (GVO) increased. The pressure and velocity are significantly influenced by the GVO, and the flow characteristics in the draft tube improve as the GVO increases. Simultaneously, the influence range of the vortex rope increased as the head increased.

Hydrodynamic modeling of non-swirling and swirling gas-particle two-phase turbulent flow using large eddy simulation

Flow structures of non-swirling and swirling particle-laden turbulence are modeled using the large eddy simulation. Four-way coupling to reveal the interaction mechanism of gas to particle, particle to gas and particle-particle collisions is adopted by a new proposed subgrid scale (SGS) particle kinetic energy-granular temperature model. Evolution of particle vortices, particle coherent structures and particle dynamics are predicted numerically, and comparisons are also conducted. Predictions are in good agreement with experimental measurements. For non-swirling flows, the typical coherent structures, and larger vortices at downward region of gas flow are dominantly. However, particle transport characteristics are significantly different, distinctly coherent structures and more vortices have not been generated in additions to fewer ones at far distance from entrance. Meanwhile, fewer coherent structures and vortices that dispersed at corner recirculation are observed for swirling flow. Compared to gas flow, the number of particle vortices is less and Reynolds stress and kinetic energy transport behaviors of particle flow are relatively independent due to inertia. Gas fluctuations of shear stresses are approximately twice greater than those of particle in non-swirling flows. A novelty finding is that turbulent fluctuations and instantaneous gas vortices of non-swirling flow are stronger than swirling case.

Контроль течений на поверхности прибрежных морских акваторий дистанционными методами

Object and purpose of research. The compare analysis carried out on the vorticies structure optical satellite pictures observation of the sea surface and the velocities field obtained by the ground HF Doppler radar. The purpose of this research is estimation of the possibility on fast monitoring feasibility assessment for the changes in littoral waters dynamics through joint application of satellite surveillance, radar probing of sea currents, wind data and numerical simulation of water circulation in the area under investigation. Materials and methods. The Sea surface observation data obtained by the satellites Sentinel 1 and 2 and by the coastal Doppler HF radar Sea Sonde in the South-east Black Sea area are analyzing. The additional data on the satellite observation of this area and on the wind condition is used from LANCE processing system [9, 10] and reanalysis NCEP/NCAR [11, 12]. The comparative analysis is carried out on the large vertical structures (LVS) establishing in the research coastal area. The numerical modeling of the water mass circulation in this area is carried out by the software package COHERENS (Belgian Institute of Natural Science) to investigate the role of the wind conditions at LVS appearance in the research area. Main results. The LVS in the vicinity of Sea coast line could be generated, as satellites observations demonstrated, by the vortices arrived along the coast line from the east-south direction in accordance with the Black Sea general circulation flow. Such vortices could be induced by the bottom roughness, rivers discharges and another local sources. Another LVS source could be connected with the wind impact in some affordable direction as demonstrated by the numerical modeling of the coastal water circulation. The continuous currents measurements at the Sea surface by Doppler HF coastal Radar supply the valuable data for further details on the space –time variation of LVS parameters and operational forecast of the hydrology regime in the aquatic areas. The operative numerical modelling of the wind impact to the surface currents structures formation provides the efficient estimation of the local large vorticity structure appearance in the research area. Conclusion. Satellite radar and optical observation provide the large scale control on the Sea surface dynamics conditions, but that is limited by relatively long interval between satellite tracks over the region of interest and do not acquire optical data at the cloud cover. The additional data could be operative obtain from the open internet resource in satellite optical, radar and meteorological data (LANCE, NCEP/NCAR), to get the more details development of the Sea surface dynamics. Sub-satellite Sea currents measurements by HF Doppler radar are useful for the reliable interpretation of satellite pictures and for the proper forecast of the aquatic area dynamics. Operative numerical modelling of the Sea coastal water circulation in accordance with variable meteo data combining with dynamic processes remote measurements could be efficient for the forecast of the environment parameters variation due the water mass transfer. In general the presented results demonstrates the combined application possibilities of the satellite observations, coastal Doppler radar measurements, internet data and operative numerical modelling for operational oceanographic and ecological monitoring of Sea coastal areas with heavy pollution load.