scholarly journals A Simple and Unified Linear Solver for Free-Surface and Pressurized Mixed Flows in Hydraulic Systems

Water ◽  
2019 ◽  
Vol 11 (10) ◽  
pp. 1979 ◽  
Author(s):  
Dechao Hu ◽  
Songping Li ◽  
Shiming Yao ◽  
Zhongwu Jin
Keyword(s):  
Free Surface ◽  
Linear System ◽  
Cross Sections ◽  
Sensitivity Study ◽  
Slot Width ◽  
Width Ratio ◽  
Hydraulic Systems ◽  
Time Step ◽  
Linear Solver ◽  
Mixed Flows

A semi–implicit numerical model with a linear solver is proposed for the free-surface and pressurized mixed flows in hydraulic systems. It solves the two flow regimes within a unified formulation, and is much simpler than existing similar models for mixed flows. Using a local linearization and an Eulerian–Lagrangian method, the new model only needs to solve a tridiagonal linear system (arising from velocity-pressure coupling) and is free of iterations. The model is tested using various types of mixed flows, where the simulation results agree with analytical solutions, experiment data and the results reported by former researchers. Sensitivity studies of grid scales and time steps are both performed, where a common grid scale provides grid-independent results and a common time step provides time-step-independent results. Moreover, the model is revealed to achieve stable and accurate simulations at large time steps for which the CFL is greater than 1. In simulations of a challenging case (mixed flows characterized by frequent flow-regime conversions and a closed pipe with wide-top cross-sections), an artificial slot (A-slot) technique is proposed to cope with possible instabilities related to the discontinuous main-diagonal coefficients of the linear system. In this test, a slot-width sensitivity study is also performed, and the suitable slot-width ratio (ε) for the linear solver is suggested to be 0.05–0.1.

Scalability studies and large grid computations for surface combatant using CFDShip-Iowa

2011 ◽  
Vol 25 (4) ◽  
pp. 466-487 ◽  
Author(s):  
Shanti Bhushan ◽  
Pablo Carrica ◽  
Jianming Yang ◽  
Frederick Stern
Keyword(s):  
Free Surface ◽  
Message Passing Interface ◽  
Free Surface Flows ◽  
Processing Unit ◽  
Mean Flow ◽  
Time Step ◽  
Vortical Structures ◽  
Cpu Time ◽  
Surface Combatant ◽  
Straight Ahead

Scalability studies and computations using the largest grids to date for free-surface flows are performed using message-passing interface (MPI)-based CFDShip-Iowa toolbox curvilinear (V4) and Cartesian (V6) grid solvers on Navy high-performance computing systems. Both solvers show good strong scalability up to 2048 processors, with V6 showing somewhat better performance than V4. V6 also outperforms V4 in terms of the memory requirements and central processing unit (CPU) time per time-step per grid point. The explicit solvers show better scalability than the implicit solvers, but the latter allows larger time-step sizes, resulting in a lower total CPU time. The multi-grid HYPRE solver shows better scalability than the portable, extensible toolkit for scientific computation solver. The main scalability bottleneck is identified to be the pressure Poisson solver. The memory bandwidth test suggests that further scalability improvements could be obtained by using hybrid MPI/open multi-processing (OpenMP) parallelization. V4-detached eddy simulation (DES) on a 300 M grid for the surface combatant model DTMB 5415 in the straight-ahead condition provides a plausible description of the vortical structures and mean flow patterns observed in the experiments. However, the vortex strengths are over predicted and the turbulence is not resolved. V4-DESs on up to 250 M grids for DTMB 5415 at 20° static drift angle significantly improve the forces and moment predictions compared to the coarse grid unsteady Reynolds averaged Navier–Stokes, due to the improved resolved turbulence predictions. The simulations provide detailed resolution of the free-surface and breaking pattern and vortical and turbulent structures, which will guide planned experiments. V6 simulations on up to 276 M grids for DTMB 5415 in the straight-ahead condition predict diffused vortical structures due to poor wall-layer predictions. This could be due to the limitations of the wall-function implementation for the immersed boundary method.

NUMERICAL SOLUTIONS OF 2D AND 3D SLAMMING PROBLEMS

2021 ◽  
Vol 153 (A2) ◽  
Author(s):  
Q Yang ◽  
W Qiu
Keyword(s):  
Free Surface ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Type Equation ◽  
The Body ◽  
Navier Stokes ◽  
Time Step ◽  
Navier Stokes Equations ◽  
Highly Nonlinear ◽  
2D And 3D

Slamming forces on 2D and 3D bodies have been computed based on a CIP method. The highly nonlinear water entry problem governed by the Navier-Stokes equations was solved by a CIP based finite difference method on a fixed Cartesian grid. In the computation, a compact upwind scheme was employed for the advection calculations and a pressure-based algorithm was applied to treat the multiple phases. The free surface and the body boundaries were captured using density functions. For the pressure calculation, a Poisson-type equation was solved at each time step by the conjugate gradient iterative method. Validation studies were carried out for 2D wedges with various deadrise angles ranging from 0 to 60 degrees at constant vertical velocity. In the cases of wedges with small deadrise angles, the compressibility of air between the bottom of the wedge and the free surface was modelled. Studies were also extended to 3D bodies, such as a sphere, a cylinder and a catamaran, entering calm water. Computed pressures, free surface elevations and hydrodynamic forces were compared with experimental data and the numerical solutions by other methods.

Computation of Wave Added Resistance by Control Surface Integration

2016 ◽  
Author(s):  
Zhiyuan Pan ◽  
Torgeir Vada ◽  
Kaijia Han
Keyword(s):  
Free Surface ◽  
Near Field ◽  
Momentum Conservation ◽  
Ship Hull ◽  
Control Surface ◽  
Floating Body ◽  
Added Resistance ◽  
Time Step ◽  
Rankine Source ◽  
Linearization Methods

A time domain Rankine source solver is extended to compute the wave added resistance of ships. The proposed approach applies the momentum conservation principle on the near field fluid volume enclosed by the wet surface of a floating body, the free surface and a control surface. The wave added resistance is then calculated by the integration over the control surface of the fluid velocities and free surface elevations. To be able to incorporate the proposed method with the Rankine source code, an interpolation scheme has been developed to compute the kinematics for the off-body points close to (or on) the free surface. Two Wigley ship models, a containership model S175 and a tanker model KVLCC2 are used to validate the present method. In general good agreement is found comparing with the model test data. The convergence behavior is examined for the proposed method including the selection of the time step and location of the control surface. Both Neumann-Kelvin and double body linearization methods are evaluated with the proposed method. It is found that the Neumann-Kelvin linearization can only be applied for slender ship hull, whereas double body method fits also for blunt ships. It is suggested to apply the proposed method with double body linearization to evaluate the wave added resistance of ships with a control surface close to the ship hull.

scholarly journals WAVETRISK-2.1: an adaptive dynamical core for ocean modelling

2021 ◽  
Author(s):  
Nicholas Keville-Reynolds Kevlahan ◽  
Florian Lemarié
Keyword(s):  
Free Surface ◽  
Shallow Water ◽  
Time Integration ◽  
Ocean Model ◽  
Mesh Adaptivity ◽  
Time Step ◽  
Vertical Coordinate ◽  
Slip Boundary ◽  
Modelling Framework ◽  
Atmosphere Model

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.

Non-Boussinesq gravity currents and surface waves generated by lock release in a circular-section channel: theoretical and experimental investigation

2019 ◽  
Vol 869 ◽  
pp. 610-633 ◽  
Author(s):  
L. Chiapponi ◽  
M. Ungarish ◽  
D. Petrolo ◽  
V. Di Federico ◽  
S. Longo
Keyword(s):  
Free Surface ◽  
Surface Waves ◽  
Cross Section ◽  
Cross Sections ◽  
Density Ratio ◽  
Circular Cross Section ◽  
Gravity Currents ◽  
Numerical Code ◽  
Long Distance ◽  
Front Speed

We present a combined theoretical and experimental study of lock-release inertial gravity currents (GCs) propagating in a horizontal channel of circular cross-section with open-top surface in the non-Boussinesq regime. A two-layer shallow-water (SW) model is developed for a generic shape of the cross-section with open top, and then implemented in a finite difference numerical code for the solution in a circular-cross-section channel of the type used in the experiments. The model predicts propagation with (almost) constant speed for a fairly long distance, accompanied by a depression of the ambient free open-top surface behind the front of the current. Sixteen experiments were conducted with a density ratio $r=0.587{-}0.939$ in full-depth and part-depth release conditions, measuring the front speed and the free-surface time series at four cross-sections. The channel was a circular tube 409 cm long, with a radius of 9.5 cm; the lengths of the locks were 52 and 103.5 cm. Density contrast was obtained by adding sodium chloride and dipotassium phosphate to fresh water. The theoretical values of the front speed and of the depression overestimate the experimental values, but they predict correctly their trend for varying parameters and provide reliable insights into the underlying mechanisms. In particular, we demonstrate that the circular cross-section increases the speed of propagation as compared to the standard rectangular cross-section case (for the same initial height and density ratio). The discrepancies between the SW predictions and the present experiments are of the same order of magnitude as those of previously published results for simpler systems (Boussinesq, rectangular). In addition to the depression, which is a wave bound to, and following the front of, the GC, the system also displays two kinds of free-surface waves, namely the initial bump (its amplitude is of the same order as the depression) and some short-length and low-amplitude waves in the tail of the bump. These free waves propagate with a celerity well predicted by the ‘fast’ eigenvalues of the mathematical model. Comparison is provided with the celerity of a solitary wave. It is expected that discrepancies between theory and experiments can be partly attributed to the presence of these waves. The reported insights and SW prediction method can be applied to a variety of cross-sections of practical interest (triangles, trapezoids, etc.).

scholarly journals Experimental study of the wake behind a surface-piercing cylinder for a clean and contaminated free surface

2000 ◽  
Vol 402 ◽  
pp. 109-136 ◽  
Author(s):  
AMY WARNCKE LANG ◽  
MORTEZA GHARIB
Keyword(s):  
Surface Tension ◽  
Free Surface ◽  
Turbulent Flows ◽  
Cross Sections ◽  
Surface Deformation ◽  
Surface Condition ◽  
Free Surface Flows ◽  
Shear Stresses ◽  
Surface Contaminants ◽  
A New Technique

This experimental investigation into the nature of free-surface flows was to study the effects of surfactants on the wake of a surface-piercing cylinder. A better understanding of the process of vorticity generation and conversion at a free surface due to the absence or presence of surfactants has been gained. Surfactants, or surface contaminants, have the tendency to reduce the surface tension proportionally to the respective concentration at the free surface. Thus when surfactant concentration varies across a free surface, surface tension gradients occur and this results in shear stresses, thus altering the boundary condition at the free surface. A low Reynolds number wake behind a surface-piercing cylinder was chosen as the field of study, using digital particle image velocimetry (DPIV) to map the velocity and vorticity field for three orthogonal cross-sections of the flow. Reynolds numbers ranged from 350 to 460 and the Froude number was kept below 1.0. In addition, a new technique was used to simultaneously map the free surface deformation. Shadowgraph imaging of the free surface was also used to gain a better understanding of the flow. It was found that, depending on the surface condition, the connection of the shedding vortex filaments in the wake of the cylinder was greatly altered with the propensity for surface tension gradients to redirect the vorticity near the free surface to that of the surface-parallel component. This result has an impact on the understanding of turbulent flows in the vicinity of a free surface with varying surface conditions.

Numerical Simulation of Nonlinear Water Wave Propagation Over Rippled Bed

2007 ◽  
Author(s):  
Gerasimos A. Kolokythas ◽  
Athanassios A. Dimas
Keyword(s):  
Free Surface ◽  
Boundary Conditions ◽  
Water Waves ◽  
Stokes Equations ◽  
Surface Flow ◽  
Surface Boundary ◽  
Computational Domain ◽  
Time Step ◽  
Outflow Boundary ◽  
Rippled Bed

In the present study, numerical simulations of the free-surface flow, developing by the propagation of nonlinear water waves over a rippled bottom, are performed assuming that the corresponding flow is two-dimensional, incompressible and viscous. The simulations are based on the numerical solution of the Navier-Stokes equations subject to the fully-nonlinear free-surface boundary conditions and the suitable bottom, inflow and outflow boundary conditions. The equations are properly transformed so that the computational domain becomes time-independent. For the spatial discretization, a hybrid scheme with finite-differences and Chebyshev polynomials is applied, while a fractional time-step scheme is used for the temporal discretization. A wave absorption zone is placed at the outflow region in order to efficiently minimize reflection of waves by the outflow boundary. The numerical model is validated by comparison to the analytical solution for the laminar, oscillatory, current flow which develops a uniform boundary layer over a horizontal bottom. For the propagation of finite-amplitude waves over a rigid rippled bed, the case with wavelength to water depth ratio λ/d0 = 6 and wave height to wavelength ratio H0/λ = 0.05 is considered. The ripples have parabolic shape, while their dimensions — length and height — are chosen accordingly to fit laboratory and field data. Results indicate that the wall shear stress over the ripples and the form drag forces on the ripples increase with increasing ripple height, while the corresponding friction force is insensitive to this increase. Therefore, the percentage of friction in the total drag force decreases with increasing ripple height.

Towards Development of Uncertainty Library for Nuclear Reactor Core Simulation

2018 ◽  
Author(s):  
Hany S. Abdel-Khalik ◽  
Dongli Huang ◽  
Ondrej Chvala ◽  
G. Ivan Maldonado
Keyword(s):  
Cross Sections ◽  
Degrees Of Freedom ◽  
Nuclear Reactor ◽  
Reactor Core ◽  
Sensitivity Study ◽  
Reactor Physics ◽  
Rigorous Approach ◽  
Uncertainty Space ◽  
Core Simulation

Uncertainty quantification is an indispensable analysis for nuclear reactor simulation as it provides a rigorous approach by which the credibility of the predictions can be assessed. Focusing on propagation of multi-group cross-sections, the major challenge lies in the enormous size of the uncertainty space. Earlier work has explored the use of the physics-guided coverage mapping (PCM) methodology to assess the quality of the assumptions typically employed to reduce the size of the uncertainty space. A reduced order modeling (ROM) approach has been further developed to identify the active degrees of freedom (DOFs) of the uncertainty space, comprising all the cross-section few-group parameters required in core-wide simulation. In the current work, a sensitivity study, based on the PCM and ROM results, is applied to identify a suitable compressed representation of the uncertainty space to render feasible the quantification and prioritization of the various sources of uncertainties. While the proposed developments are general to any reactor physics computational sequence, the proposed approach is customized to the TRITON-NESTLE computational sequence, simulating the BWR lattice model and the core model, which will serve as a demonstrative tool for the implementation of the algorithms.

Three-Dimensional Large Amplitude Body Motions in Waves

2007 ◽  
Author(s):  
Xinshu Zhang ◽  
Robert F. Beck
Keyword(s):  
Free Surface ◽  
Three Dimensional ◽  
Surface Boundary ◽  
Forward Speed ◽  
Time Step ◽  
Surface Boundary Conditions ◽  
Submerged Body ◽  
Calm Water ◽  
Wigley Hull ◽  
Constant Strength

Three-dimensional, time-domain, wave-body interactions are studied in this paper for cases with and without forward speed. In the present approach, an exact body boundary condition and linearized free surface boundary conditions are used. By distributing desingularized sources above the calm water surface and using constant-strength panels on the exact submerged body surface, the boundary integral equations are solved numerically at each time step. Once the fluid velocities on the free surface are computed, the free surface elevation and potential are updated by integrating the free surface boundary conditions. After each time step, the body surface and free surface are regrided due to the instantaneous changing submerged body geometry. The desingularized method applied on the free surface produces non-singular kernels in the integral equations by moving the fundamental singularities a small distance outside of the fluid domain. Constant strength panels are used for bodies with any arbitrary shape. Extensive results are presented to validate the efficiency of the present method. These results include the added mass and damping computations for a hemisphere. The calm water wave resistance for a submerged spheroid and a Wigley hull are also presented. All the computations with forward speed are started from rest and proceed until a steady state is reached. Finally, the time-domain forced motion results for a modified Wigley hull with forward speed are shown and compared with the experiments for both linear computations and body-exact computations.

Development of a Time Domain Boundary Element Method for the Seakeeping Behavior of a Cruise Ship Under Combined Strong Wind and Extreme Irregular Beam Seas

2014 ◽  
Author(s):  
G. D. Gkikas ◽  
F. van Walree
Keyword(s):  
Free Surface ◽  
Time Domain ◽  
Model Tests ◽  
The Body ◽  
Computational Method ◽  
Green Functions ◽  
Strong Wind ◽  
Cruise Ship ◽  
Time Step ◽  
Beam Seas

A computational method for the seakeeping behavior of a cruise ship at zero speed and under severe wind and oblique wave loads is presented. The proposed methodology is a time-domain panel method where the transient Green functions used for the estimation and implementation of the free surface effects on the vessel’s motions are estimated assuming constant low lateral speed, instead of the common practice zero speed influence functions. For the evaluation of the overall hydrodynamic forces, the so called “blended approach” is followed in the sense that the induced hydrodynamic pressures due to the scattering and radiation phenomena are calculated over the linearized position of the body, ignoring any displacements with respect to its mean position, while the hydrostatic and non-linear Froude-Krylov forces are considered at the actual body location and taking into account the free surface elevation at each time step. For the validation of the proposed methodology, heave and roll motions, the drift velocity as well as lateral accelerations of the vessel were investigated for two cases of severe beam seas combined with a constant strong wind load and the results were compared against experimental model tests. The model tests were performed to investigate the vessel’s behavior under extreme weather conditions. The low lateral speed Green functions were estimated for a speed similar to the one that the vessel was expected to drift, an estimation based on the model tests, as well as for the case where the input speed corresponded to the half of the expected speed. Good agreement was presented for both cases, showing that accurate and computationally efficient numerical simulations of the vessel’s motions under severe wind and wave excitations can be obtained by using low lateral speed transient Green functions.

Sign in / Sign up

Export Citation Format

Share Document