scholarly journals On the computation of the barotropic mode of a free-surface world ocean model

1995 ◽  
Vol 13 (6) ◽  
pp. 675-688 ◽  
Author(s):  
E. Deleersnijder ◽  
J.-M. Campin

Abstract. The free-surface formulation of the equations of our world ocean model is briefly described. The barotropic mode equations are solved according to the split-explicit method, using different time steps for the external and internal modes. Because the numerical algorithm is implemented on the B-grid, a spurious, free-surface, two-grid interval mode may develop. This mode must be filtered out. The properties of two filters are theoretically investigated and their actual performance is tested in a series of numerical experiments. It is seen that one of these filters may severely perturb the local mass conservation, rendering it impossible to enforce the impermeability of the surface or the bottom of the ocean. The dynamics of the external mode is also examined, by studying the depth-integrated momentum equations. The depth-integral of the pressure force due to the slope of the ocean surface is approximately balanced by the depth-integral of the force ensuing from the horizontal variations of the density. The depth-integral of the Coriolis force is an order of magnitude smaller, except in the Southern Ocean. Two variational principles are resorted to for computing the fictitious ocean surface elevation corresponding to the approximate equilibrium between the dominant forces of the barotropic momentum equations.

2017 ◽  
Author(s):  
Ralph Timmermann ◽  
Sebastian Goeller

Abstract. A Regional Antarctic and Global Ocean (RAnGO) model has been developed to study the interaction between the world ocean and the Antarctic ice sheet. The coupled model is based on a global implementation of the Finite Element Sea-ice Ocean Model (FESOM) with a mesh refinement in the Southern Ocean, particularly in its marginal seas and in the sub-ice shelf cavities. The cryosphere is represented by a regional setup of the ice flow model RIMBAY comprising the Filchner-Ronne Ice Shelf and the grounded ice in its catchment area up to the ice divides. At the base of the RIMBAY ice shelf, melt rates from FESOM's ice-shelf component are supplied. RIMBAY returns ice thickness and the position of the grounding line. The ocean model uses a pre-computed mesh to allow for an easy adjustment of the model domain to a varying cavity geometry. RAnGO simulations with a 20th-century climate forcing yield realistic basal melt rates and a quasi-stable grounding line position close to the presently observed state. In a centennial-scale warm-water-inflow scenario, the model suggests a substantial thinning of the ice shelf and a local retreat of the grounding line. The potentially negative feedback from ice-shelf thinning through a rising in-situ freezing temperature is more than outweighed by the increasing water column thickness in the deepest parts of the cavity. Compared to a control simulation with fixed ice-shelf geometry, the coupled model thus yields a slightly stronger increase of ice-shelf basal melt rates.


2009 ◽  
Vol 625 ◽  
pp. 435-443 ◽  
Author(s):  
MARK A. KELMANSON

A novel pseudo-three-timescale asymptotic procedure is developed and implemented for obtaining accurate approximations to solutions of an evolution equation arising in thin-film free-surface viscous flow. The new procedure, which employs strained fast and slow timescales, requires considerably fewer calculations than its standard three-timescale counterpart employing fast, slow and slower timescales and may readily be applied to other evolution equations of fluid mechanics possessing wave-like solutions exhibiting exponential decay in amplitude and variations in phase over disparate timescales. The new method is validated on the evolution of free-surface waves on a thin, viscous film coating the exterior of a horizontal rotating cylinder and is shown to yield accurate solutions up to non-dimensional times exceeding by an order of magnitude those of previous related studies. Results of the new method applied to this test problem are demonstrated to be in excellent agreement, over large timescales, with those of corroborative spectrally accurate numerical integrations.


2021 ◽  
Author(s):  
Nicholas Keville-Reynolds Kevlahan ◽  
Florian Lemarié

Abstract. This paper introduces WAVETRISK-2.1 (i.e. WAVETRISK-OCEAN), an incompressible version of the atmosphere model wavetrisk-1.x with free-surface. This new model is built on the same wavelet-based dynamically adaptive core as wavetrisk, which itself uses DYNANICO's mimetic vector-invariant multilayer rotating shallow water formulation. Both codes use a Lagrangian vertical coordinate with conservative remapping. The ocean variant solves the incompressible multilayer shallow water equations with inhomogeneous density layers. Time integration uses barotropic--baroclinic mode splitting via an semi-implicit free surface formulation, which is about 34–44 times faster than an unsplit explicit time-stepping. The barotropic and baroclinic estimates of the free surface are reconciled at each time step using layer dilation. No slip boundary conditions at coastlines are approximated using volume penalization. The vertical eddy viscosity and diffusivity coefficients are computed from a closure model based on turbulent kinetic energy (TKE). Results are presented for a standard set of ocean model test cases adapted to the sphere (seamount, upwelling and baroclinic turbulence). An innovative feature of wavetrisk-ocean is that it could be coupled easily to the wavetrisk atmosphere model, thus providing a first building block toward an integrated Earth-system model using a consistent modelling framework with dynamic mesh adaptivity and mimetic properties.


2019 ◽  
Vol 5 (3) ◽  
pp. eaav5335 ◽  
Author(s):  
Yan Kuai ◽  
Junxue Chen ◽  
Xi Tang ◽  
Yifeng Xiang ◽  
Fengya Lu ◽  
...  

Surface plasmon resonance microscopy (SPRM) with single-direction illumination is a powerful platform for biomedical imaging because of its wide-field, label-free, and high-surface-sensitivity imaging capabilities. However, two disadvantages prevent wider use of SPRM. The first is its poor spatial resolution that can be as large as several micrometers. The second is that SPRM requires use of metal films as sample substrates; this introduces working wavelength limitations. In addition, cell culture growth on metal films is not as universally available as growth on dielectric substrates. Here we show that use of azimuthal rotation illumination allows SPRM spatial resolution to be enhanced by up to an order of magnitude. The metal film can also be replaced by a dielectric multilayer and then a different label-free surface-sensitive photonic microscopy is developed, which has more choices in terms of the working wavelength, polarization, and imaging section, and will bring opportunities for applications in biology.


2020 ◽  
Vol 2020 ◽  
pp. 1-18
Author(s):  
Puyang Gao

In this paper, we develop a new computational framework to investigate the sloshing free surface flow of Newtonian and non-Newtonian fluids in the rectangular tanks. We simulate the flow via a two-phase model and employ the fixed unstructured mesh in the computation to avoid the mesh distortion and reconstruction. As for the solution of Navier–Stokes equation, we utilize the SUPG finite element method based on the splitting scheme. The same order interpolation functions are then used for velocity and pressure. Moreover, the moving interface is captured via the concise level set method. We take advantage of the implicit discontinuous Galerkin method to handle the solution of level set and its reinitialization equations. A mass correction technique is also added to ensure the mass conservation property. The dam break-free surface flow is simulated firstly to demonstrate the validity of our mathematical model. In addition, the sloshing Newtonian fluid in the tank with flat and rough bottoms is considered to illustrate the feasibility and robustness of our computational scheme. Finally, the development of free surface for non-Newtonian fluid is also studied in the two tanks, and the influence of power-law index on the sloshing fluid flow is analyzed.


Ocean Science ◽  
2017 ◽  
Vol 13 (5) ◽  
pp. 765-776 ◽  
Author(s):  
Ralph Timmermann ◽  
Sebastian Goeller

Abstract. The Regional Antarctic ice and Global Ocean (RAnGO) model has been developed to study the interaction between the world ocean and the Antarctic ice sheet. The coupled model is based on a global implementation of the Finite Element Sea-ice Ocean Model (FESOM) with a mesh refinement in the Southern Ocean, particularly in its marginal seas and in the sub-ice-shelf cavities. The cryosphere is represented by a regional setup of the ice flow model RIMBAY comprising the Filchner–Ronne Ice Shelf and the grounded ice in its catchment area up to the ice divides. At the base of the RIMBAY ice shelf, melt rates from FESOM's ice-shelf component are supplied. RIMBAY returns ice thickness and the position of the grounding line. The ocean model uses a pre-computed mesh to allow for an easy adjustment of the model domain to a varying cavity geometry. RAnGO simulations with a 20th-century climate forcing yield realistic basal melt rates and a quasi-stable grounding line position close to the presently observed state. In a centennial-scale warm-water-inflow scenario, the model suggests a substantial thinning of the ice shelf and a local retreat of the grounding line. The potentially negative feedback from ice-shelf thinning through a rising in situ freezing temperature is more than outweighed by the increasing water column thickness in the deepest parts of the cavity. Compared to a control simulation with fixed ice-shelf geometry, the coupled model thus yields a slightly stronger increase in ice-shelf basal melt rates.


2019 ◽  
Vol 141 (11) ◽  
Author(s):  
Jieun Song ◽  
Seung Jin Song

AbstractAn integrated analytical model to predict non-axisymmetric flow fields and rotordynamic forces in a shrouded centrifugal compressor has been newly developed and validated. The model is composed of coupled, conservation law-based, bulk-flow submodels, and the model takes into account the flow coupling among the blades, labyrinth seals, and shroud cavity. Thus, the model predicts the entire flow field in the shrouded compressor when given compressor geometry, operating conditions, and eccentricity. When compared against the experimental data from part 1, the new model accurately predicts the evolution of the pressure perturbations along the shroud and labyrinth seal cavities as well as the corresponding rotordynamic stiffness coefficients. For the test compressor, the cross-coupled stiffness rotordynamic excitation is positive; the contribution of the shroud is the highest; the contribution of the seals is less than but on the same order of magnitude as that of the shroud; and contribution of impeller blades is insignificant. The new model also enables insight into the physical mechanism for pressure perturbation development. The labyrinth seal pressure distribution becomes non-axisymmetric to satisfy mass conservation in the seal cavity, and this non-axisymmetry, in turn, serves as the influential boundary condition for the pressure distribution in the shroud cavity. Therefore, for accurate flow and rotordynamic force predictions, it is important to model the flow coupling among the components (e.g., impeller, shroud, labyrinth seal, etc.), which determines the non-axisymmetric boundary conditions for the components.


2019 ◽  
Vol 49 (9) ◽  
pp. 2365-2376 ◽  
Author(s):  
A. J. George Nurser ◽  
Stephen M. Griffies

AbstractWe detail the physical means whereby boundary transfers of freshwater and salt induce diffusive fluxes of salinity. Our considerations focus on the kinematic balance between the diffusive fluxes of salt and freshwater, with this balance imposed by mass conservation for an element of seawater. The flux balance leads to a specific balanced form for the diffusive salt flux immediately below the ocean surface and, in the Boussinesq approximation, to a specific form for the salinity flux. This balanced form should be used in specifying the surface boundary condition for the salinity equation and the contribution of freshwater to the buoyancy budget.


2018 ◽  
Vol 33 (6) ◽  
pp. 325-331
Author(s):  
Ilya A. Chernov ◽  
Nikolay G. Iakovlev

Abstract In the present paper we consider the first results of modelling the World Ocean biogeochemistry system within the framework of the Earth system model: a global atmosphere-ocean-ice-land-biogeochemistry model. It is based on the INMCM climate model (version INMCM39) coupled with the pelagic ecosystem model BFM. The horizontal resolution was relatively low: 2∘ × 2.5∘ for the ‘longitude’ and ‘latitude’ in transformed coordinates with the North Pole moved to land, 33 non-equidistant σ-horizons, 1 hour time step. We have taken into account 54 main rivers worldwide with run–off supplied by the atmosphere submodel. The setup includes nine plankton groups, 60 tracers in total. Some components sink with variable speed. We discuss challenges of coupling the BFM with the σ-coordinate ocean model. The presented results prove that the model output is realistic in comparison with the observed data, the numerical efficiency is high enough, and the coupled model may serve as a basis for further simulations of the long-term climate change.


Sign in / Sign up

Export Citation Format

Share Document