Three-dimensional numerical modeling of thermohaline and wind-driven circulations in the Persian Gulf

The Persian Gulf circulation is investigated with respect to the relevant forcing mechanism including wind stress and thermohaline surface fluxes by using a three-dimensional numerical hydrodynamic model. The model results show a correlation between the strength of the bottom layer outflow of the Persian Gulf and that of the Indian Ocean Surface Water (IOSW) inflow into the Gulf. The inflow of IOSW into the Gulf attain maximum values in May–June in conjunction with peak bottom outflow through the Hormuz Strait. The results of sensitivity experiment indicate that circulation is dominated by thermohaline flows at almost all parts of the Gulf. The heat fluxes play an essential role on the general circulation of the Persian Gulf. In spring and summer, the wind stress generates southeast-flowing surface currents of magnitude about 5 cm/s along the Saudi Arabia and Iranian coasts on the northern Gulf. In winter and autumn, due to weak static stability, the wind produces mesoscale eddies in most parts of the Gulf. In winter and spring the wind stress acts to reinforce the thermohaline circulation of deep outflow. Conversely, in summer and autumn the wind forcing acts in opposition to the thermohaline forcing and causes a bottom inflow from Oman Sea into the Gulf.     

Numerical modeling of hydrodynamic and hydrothermal characteristics in subtropical alpine lake

A three-dimensional, time-dependent hydrodynamic and hydrothermal model was performed and applied to the subtropical alpine Yuan-Yang Lake (YYL) in northeastern region of Taiwan. The model was driven with discharge inflow, heat, and wind stress to simulate the hydrodynamic and hydrothermal in the lake. The model was validated with measured water surface elevation, current, and temperature in 2008. The overall model simulation results are in quantitative agreement with the available field data. The validated model was then used to investigate wind-driven current, mean circulation, and residence time in the YYL. The modeling results reveal that the velocity field along the wind axis present the variations over depth with return current where the velocity at the surface layer is along the wind direction while it is opposite near 1 m below water surface. The simulated mean current indicates that the surface currents flow towards the southwest direction and form a clock-wise rotation. The calculated residence time is strongly dependent on the inflows and wind effects. Regression analysis of model results reveals that an exponential regression equation can be employed to correlate the residence time to change of discharge input. The residence time without wind stress is higher than that with wind effect, indicating that wind plays an important role in lake mixing. The calculated residence time is approximately 2-2.5 days under low inflow with wind effect.     

Studying the role of temperature and salinity anomalies in the formation of world ocean meridional circulation modes

A numerical model of the ocean circulation is considered for using in a coupled oceanice-atmosphere climatic model. Preliminary experiments conducted on a numerical grid with a coarse spatial resolution over a period of one thousand years reproduce the World ocean dynamics under climatic surface forcing. By introducing anomalies into the surface sources, the sensitivity of the global ocean circulation to thermohaline variability in the northern regions of the Atlantic Ocean is investigated.     

Study of storm surge due to Typhoon Linda (1997) in the Gulf of Thailand using a three dimensional ocean model

Numerical integrations using the three dimensional ocean model based on the princeton ocean model (POM) were applied for the study of both sea level elevation and ocean circulation patterns forced by the wind fields during typhoons that moved over the Gulf of Thailand (GoT). The simulation concerned a case of Typhoon Linda which occurred during November 1-4, 1997. Typhoon Linda was one of the worst storms that passed the Gulf of Thailand and hit the southern coastal provinces of Thailand on November 3, 1997. It caused flooding and a strong wind covering large areas of agriculture and fisheries, which destroyed households, utilities and even human lives. The model is the time-dependent, primitive equation, Cartesian coordinates in a horizontal and sigma coordinate in the vertical. The model grid has 37 × 97 orthogonal curvilinear grid points in the horizontal, with variable spacing from 2 km near the head of the GoT to 55 km at the eastern boundary, with 10 sigma levels in the vertical conforming to a realistic bottom topography. Open boundary conditions are determined by using radiation conditions, and the sea surface elevation is prescribed from the archiving, validation and interpretation of satellite oceanographic data (AVISO). The initial condition is determined from the spin up phase of the first model run, which was executed by using wind stress calculated from climatological monthly mean wind, restoring-type surface heat and salt and climatological monthly mean freshwater flux. The model was run in spin up phase until an ocean model reached an equilibrium state under the applied force. A spatially variable wind field taken from the European Centre for Medium-Range Weather Forecasts (ECMWF) is used to compute the wind stress directly from the velocity fluctuations. Comparison of tendency between the sea surface elevations from model and the observed significant wave heights of moored buoys in the Gulf of Thailand under Seawatch project is investigated. The model predicts the sea level elevation up to 68.5 cm at the Cha-Am area located in the north of where the typhoon strands to the shore. Results of sea level elevation show that there is an area of peak set-up in the upper gulf, particularly in the western coast, and the effects of the storm surge are small at the lower gulf. During the entire period of this study, the surge in the gulf was induced by the northeasterly wind blowing over it.     

Numerical modelling of residual flow and salinity in the Río de la Plata

The Río de la Plata discharges into the Atlantic Ocean. The particular characteristics of the study area, the variable width and shallowness of the river, the high fluvial discharges and the dynamic processes involving interactions between river discharges, tidal currents and wind, generate complex velocity and salinity fields. We applied the hydrodynamic model RMA-10 to examine the effects of various forcing (tides, flow discharge and winds) on residual currents and salinity fields in the Río de la Plata, focusing on the outer zone of the river. The RMA-10 code, developed by Ian King, is a multiparameter finite element model representing estuarine flow in three dimensions. In this study the model has been applied in a depth-averaged-baroclinic mode and a series of observed data is used for model calibration and verification. The model result shows that it is able to simulate velocity and the salinity fields with a reasonable accuracy. The analysis of residual currents in the river, when forced by freshwater discharge and astronomical tide, shows that the flow discharge takes place mainly over the shallower areas of the river and that the saline water is advected up-river through the deeper channels. The numerical simulations show that the winds from the South-West and North-East quadrants have a great influence over the salinity and velocity fields.     

Dynamic modelling of dissolved oxygen in the creeks of Sagar island,Hooghly–Matla estuarine system,West Bengal,India

Hooghly–Matla estuarine ecosystem is one of the largest estuarine ecosystems of the world. Sagar island is the largest delta in this estuarine complex. This island is criss-crossed by small and large creeks with mangrove vegetation and all are connected to the principal estuarine water. Decomposition of mangrove litter in soil is major source of inorganic nutrient to phytoplankton of the adjacent creeks. Deforestation of mangrove affects the primary production, which in turn reduces the availability of dissolved oxygen for the organisms residing in the estuary. Considering the importance of dissolved oxygen in various aspects of aquatic life, a dynamic model of dissolved oxygen at Sagar island of Hooghly–Matla estuarine complex with the help of single dimension differential equation is proposed in the present paper. Different physical, chemical and biological factors such as solar irradiance, temperature, salinity of water, particulate organic matter, re-aeration, wind velocity, phytoplankton and zooplankton, which control the fluctuation of dissolved oxygen, are included in the present model. Most of the parameter values are collected directly from the field surveys. The parameter values which are not able to collect from the field, obtained from literatures are calibrated. To make the model realistic it is properly validated with observed data and to know the statistical significance, chi square goodness fit test is performed. Field surveys are performed over two years. During calibration and validation, two sets of data (first year and second year data) are used. Chi-square values are 5.97 and 6.17 for first and second sets of data respectively (p < 0.05). Sensitivity analysis reveals that optimal light intensity is the most sensitive parameter for dissolved oxygen dynamics. Results also show that wind velocity, solar irradiation, salinity of water and temperature are important factors for controlling the dynamics of dissolved oxygen. Macrophytes have very little contribution to oxygen production in the creeks of Sagar island. Model reveals that low dissolved oxygen in the creek water is one of the causes of decline in fish population of the estuary.     

Three‐Dimensional Surface Water Waves Governed by the Forced Benney–Luke Equation

This paper studies the propagation of three‐dimensional surface waves in water with an ambient current over a varying bathymetry. When the ambient flow is near the critical speed, under the shallow water assumptions, a forced Benney–Luke (fBL) equation is derived from the Euler equations. An asymptotic approximation of the water's reaction force over the varying bathymetry is derived in terms of topographic stress. Numerical simulations of the fBL equation over a trough are compared to those using a forced Kadomtsev–Petviashvilli equation. For larger variations in the bathymetry that upstream‐radiating three‐dimensional solitons are observed, which are different from the upstream‐radiating solitons simulated by the forced Kadomtsev–Petviashvilli equation. In this case, we show the fBL equation is a singular perturbation of the forced Kadomtsev–Petviashvilli equation which explains the significant differences between the two flows.     

Sensitivity of hypersaline Arabian Gulf to seawater desalination plants

Due to the hot, arid nature of its bordering lands, seawater in the Arabian Gulf can have significant evaporation rates leading to hypersaline conditions. If additional desalination plants were to operate along its coast, then the extraction of desalinated water and returned brine waste stream into the Gulf would increase the salinity. This paper uses a tidally and cross-sectionally averaged mathematical model that reveals multiplicative dependence of the salinity on factors associated with river flow, evaporation rates and each of the desalination plants. Present-day desalinated water production rates are in the linear regime, but hypersalinity has exponential sensitivity to the position and volumetric rate of desalinated water extraction.     

Numerical processing of thin-film thermometer data for determining transient heat fluxes

The purpose of this study is the evaluation of transient heat fluxes using thin-film thermometer data. A thin-film temperature sensor operated in transient mode enables very accurate surface temperature measurements (e.g., in turbine blades or in hypersonic wind tunnels), which are obtained because the gauge is non-intrusive and has a high frequency response; however, the complexity of data processing is the major drawback for the popular use of the thin-film sensor.     

Seasonal variation of temperature and salinity fields and water masses in the Bien Dong (South China) sea

Seasonal and bimonthly analyses of the CTD and STD data in the Bien Dong Sea have been performed using a variational inverse method (VIM). The results, materialized as gridded data sets (horizontal resolution at 1/4 of a degree at 26 depth levels from 5m to 4500m), show the basic features of the hydrology of the Bien Dong Sea and its trend of seasonal variabilities affecting the properties of the water masses. The seasonally reversing monsoon wind and the variation of net air-sea fluxes play an important role in determining the transformation and spreading of the water masses. The present work is conceived as a basic support to more advanced studies, including diagnostic calculations and assimilation of hydrological data in 3D primitive equation models for the Bien Dong Sea.     

An analogy between vortical and thermal fields subject to the large pressure gradient

Local thermal disturbances generated by a small confined mechanical forcing on the curved surface are under consideration within the framework of the interacting boundary layer in the limit of large Reynolds number. The variations of the self-induced pressure and temperature are found to be explicitly related, and an analytic solution of a linear perturbation analysis is presented. The heat flux associated with the self-induced pressure may attain values comparable with the heat flux existing in the oncoming boundary layer. The analogy between the normal-to-wall vorticity and spanwise derivative of temperature is rigorously derived to supersede locally the well-known analogy by Reynolds. Heavily increased thermal loads can cause severe damage to the engine in operation.     

含正三角形孔口热电材料的断裂力学分析

采用复变函数方法研究了含正三角形孔口热电材料中受到无穷远处均匀电流密度和能量流作用下的断裂问题.在电绝缘和热绝缘边界条件下,获得了温度场和应力场的解析表达式,分析了三角形尺寸、能量流载荷对热电材料性能的影响.结果表明,能量流载荷和三角形尺寸的变化对环向能量流、环向应力、环向热流有较明显的影响.     

Retracted: Unsteady flow and heat transfer on a rotating disk in a porous medium

This article has been retracted. See retraction notice DOI: 10.1002/mma.850 . An unsteady flow and heat transfer in a porous medium of a viscous incompressible fluid over a rotating disk in an otherwise ambient fluid are studied. The unsteadiness in the flow field is caused by the angular velocity of the disk which varies with time. The new self‐similar solution of the Navier–Stokes and energy equations is obtained numerically. The solution obtained here is not only the solution of the Navier–Stokes equations, but also of the boundary layer equations. Also, for a simple scaling factor, it represents the solution of the flow and heat transfer in the forward stagnation‐point region of a rotating sphere or over a rotating cone. The asymptotic behaviour of the solution for a large porosity or for a large independent variable is also examined. The surface shear stresses in the radial and tangential directions and the surface heat transfer increase as the acceleration parameter increases. Also, the surface shear stress in the radial direction and the surface heat transfer decrease with increasing porosity, but the surface shear stress in the tangential direction increases. Copyright © 2006 John Wiley & Sons, Ltd.     

Heat and mass transfer in MHD free convection from a moving permeable vertical surface by a perturbation technique

Numerical results are presented for heat and mass transfer effect on hydromagnetic flow of a moving permeable vertical surface. An analysis is performed to study the momentum, heat and mass transfer characteristics of MHD natural convection flow over a moving permeable surface. The surface is maintained at linear temperature and concentration variations. The non-linear coupled boundary layer equations were transformed and the resulting ordinary differential equations were solved by perturbation technique [Aziz A, Na TY. Perturbation methods in heat transfer. Berlin: Springer-Verlag; 1984. p. 1–184; Kennet Cramer R, Shih-I Pai. Magneto fluid dynamics for engineers and applied physicists 1973;166–7]. The solution is found to be dependent on several governing parameter, including the magnetic field strength parameter, Prandtl number, Schmidt number, buoyancy ratio and suction/blowing parameter, a parametric study of all the governing parameters is carried out and representative results are illustrated to reveal a typical tendency of the solutions. Numerical results for the dimensionless velocity profiles, the temperature profiles, the concentration profiles, the local friction coefficient and the local Nusselt number are presented for various combinations of parameters.     

Three-dimensional exact thermo-elastic analysis of multilayered two-dimensional quasicrystal nanoplates

A three-dimensional thermo-elastic analytical solution for two-dimensional quasicrystal simply supported nanoplates subjected to a temperature change on their top surface is presented. The nonlocal theory and pseudo-Stroh formalism are used to obtain the exact solution for a homogeneous two-dimensional decagonal quasicrystal nanoplate with its thickness direction as a quasi-periodic direction. The propagator matrix method is introduced to deal with the corresponding multilayered nanoplates. Comprehensive numerical results show that nonlocal parameters, stress-temperature coefficients, stacking sequences have great influence on the stress, displacement components and heat fluxes of the nanoplates. In addition, the stacking sequences also influence the temperature and heat fluxes of the nanoplate. The exact thermo-elastic solution should be of interest to the design of the two-dimensional quasicrystal homogeneous and multilayered plates. The mechanical behaviors of the nanoplates in numerical results can also serve as benchmarks to verify various thin-plate theories or other numerical methods.     

表面粗糙度变化对近地层大气的影响——(Ⅰ)—风速和剪应力变异

本文系统地分析了任意分布的表面粗糙度变化,对三维和二维大气边界层近地风速和剪应力的影响,提出了一个四层结构的线性理论.由该理论得到的解与已有的结果符合得很好.     

Oscillatory reaction-diffusion equations with temporally varying parameters

Periodic wave trains are the generic one-dimensional solution form for reaction-diffusion equations with a limit cycle in the kinetics. Such systems are widely used as models for oscillatory phenomena in chemistry, ecology, and cell biology. In this paper, we study the way in which periodic wave solutions of such systems are modified by periodic forcing of kinetic parameters. Such forcing will occur in many ecological applications due to seasonal variations. We study temporal forcing in detail for systems of two reaction diffusion equations close to a supercritical Hopf bifurcation in the kinetics, with equal diffusion coefficients. In this case, the kinetics can be approximated by the Hopf normal form, giving reaction-diffusion equations of λ-ω type. Numerical simulations show that a temporal variation in the kinetic parameters causes the wave train amplitude to oscillate in time, whereas in the absence of any temporal forcing, this wave train amplitude is constant. Exploiting the mathematical simplicity of the λ-ω form, we derive analytically an approximation to the amplitude of the wave train oscillations with small forcing. We show that the amplitude of these oscillations depends crucially on the period of forcing.     

Model order reduction for thermomechanical problems including radiation

Thermal field problems including heat exchange by radiation lead to nonlinear system equations with a high number of inputs and outputs as radiation heat fluxes correspond to the fourth power of the temperature and thermal loads are distributed over the whole surface. In an alternative approach presented here, radiation is defined as a part of the load vector. Thus, the system matrices are constant. Furthermore, loads changing synchronously during operation are grouped into one column of the input matrix and load vector snapshots are used to consider the radiation heat fluxes. Hence, the Krylov Subspace Method can be applied to significantly reduce the system dimension and the computation times allowing transient thermal parameter studies. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)