Prediction of Roll Damping Using Viscous Flow Solvers

2012 ◽  
Author(s):  
Manuel Manzke ◽  
Thomas Rung
Keyword(s):  
Domain Size ◽  
Building Blocks ◽  
Sensitivity Study ◽  
Convergence Criterion ◽  
Parameter Study ◽  
Floating Bodies ◽  
Time Step ◽  
Roll Damping ◽  
Damping Coefficients ◽  
Ranse Solver

This article illustrates the use of a RANSE solver coupled to a motion solver to predict the free roll decay and the associated damping coefficients of floating bodies. The necessary building blocks to perform such a prediction are described briefly. A sensitivity study for the convergence criterion, the time step, the domain size and the grid resolution is performed for a simple 2-dimensional barge. The results are compared to results from experiments. Furthermore a simulation for a free roll decay of a Navy Combatant is performed, considering the results of the parameter study for the Barge. Overall results indicate that the natural roll frequency can be well predicted, while the prediction of the roll damping coefficients is afflicted with some uncertainties.

Using CFD to Assess Low Frequency Damping

2012 ◽  
Author(s):  
Alessio Pistidda ◽  
Harald Ottens ◽  
Richard Zoontjes
Keyword(s):  
Standard Procedure ◽  
Low Frequency ◽  
Added Mass ◽  
Model Tests ◽  
Viscous Damping ◽  
Viscous Effects ◽  
Floating Bodies ◽  
Time Step ◽  
Damping Coefficients ◽  
Quadratic Component

During offshore installation operations, floating bodies are often moored using soft mooring which are designed to withstand the environmental forces. Large amplitude motions often occur due to excitation by slowly varying wind and wave drift forces. To analyze these motions the dynamic system has to be accurately described, which includes an estimation of the added mass and damping coefficients. In general, the added mass can be accurately calculated with traditional potential theory. However for the damping this method is not adequate because viscous effects play an important role. Generally these data are obtained using model tests. This paper validates the CFD methodology as an alternative to model tests to evaluate the viscous damping. The aim is to define a standard procedure to derive viscous damping coefficients for surge, sway and yaw motion of floating bodies. To estimate viscous damping in CFD, a 3D model of the launch and float-over barge H-851 was used. For this barge, model test data is available which could be compared with the results of the CFD analysis. For the simulations, the commercial package STAR-CCM+ with the implicit unsteady solver for Reynolds-Averaged Navier-Stokes (RANS) equations was used. The turbulence model implemented was the k-Omega-SST. Numerical errors have been assessed performing sensitivity analysis on time step and grid size. Damping has been investigated by performing decay simulations as in the model tests, taking the effect of coupling among all motions into account. The P-Q fitting method has been used to determine the linear and quadratic component of the damping. Numerical results are validated with those obtained from the towing tank. Results show that CFD is an adequate tool to estimate the low frequency damping in terms of equivalent damping. More investigations are required to determine the linear and quadratic component.

Fully Coupled Time Domain Modelling of 3D Floating Bodies and Mooring Systems in Regular and Irregular Sea States

2012 ◽  
Author(s):  
Kameswara S. Vepa ◽  
Diederik Van Nuffel ◽  
Wim Van Paepegem ◽  
Joris Degrieck
Keyword(s):  
Time Domain ◽  
Domain Size ◽  
Floating Body ◽  
Computational Time ◽  
Floating Bodies ◽  
Step Size ◽  
Coupled Models ◽  
Time Step ◽  
Time Step Size ◽  
Fully Coupled

Research on floating bodies like Wave Energy Converters (WECs) and Laser Imaging Detection And Ranging (LIDAR) systems has recently known a large growth. To study the minute details of the working model, it is important to study the effect of interactions between the waves, floating bodies and the mooring systems that are controlling the motion of the floating body. To achieve a more realistic numerical model in the time domain, a number of programs are linked together. The idea is to use the strength of each individual program for better results and also reduce the computational time. This paper provides a solution in the direction of using a fully coupled time domain coupling code that controls the data flow between a fluid solver, a structural solver, and a kinematic system simulator. Two- and three-dimensional fully coupled models are studied for calculation times and accuracy of results, and scaling is tested through parallelization on a large HPC cluster. The time step size of the whole model can be controlled by the user. Calculation times and memory requirements vary largely based on the factors like: domain size, SPH particle size, material model used for the floating body and the mooring system, complexity of the mechanical system inside the floating body. As a test case, a rigid body model is presented in this paper.

scholarly journals Transient Two-Way Molecular-Continuum Coupling with OpenFOAM and MaMiCo: A Sensitivity Study

Computation ◽  
2021 ◽  
Vol 9 (12) ◽  
pp. 128
Author(s):  
Helene Wittenberg ◽  
Philipp Neumann
Keyword(s):  
Domain Size ◽  
Sensitivity Study ◽  
Computational Effort ◽  
Lessons Learned ◽  
Computational Domain ◽  
Time Step ◽  
Coupling Algorithm ◽  
Long Term Behavior

Molecular-continuum methods, as considered in this work, decompose the computational domain into continuum and molecular dynamics (MD) sub-domains. Compared to plain MD simulations, they greatly reduce computational effort. However, the quality of a fully two-way coupled simulation result strongly depends on a variety of system-specific parameters, and the corresponding sensitivity is only rarely addressed in the literature. Using a state-flux molecular-continuum coupling algorithm, we investigated the influences of various parameters, such as the size of the overlapping region, the coupling time step and the quality of ensemble-based sampling of flow velocities, in a Couette flow scenario. In particular, we considered a big setup in terms of domain size and number of time steps, which allowed us to investigate the long-term behavior of the coupling algorithm close to the incompressible regime. While mostly good agreement was reached on short time scales, it was the long-term behavior which differed even with slightly differently parametrized simulations. We demonstrated our findings by measuring the error in velocity, and we summarize our main observations with a few lessons learned.

scholarly journals Bilge Keel Induced Roll Damping of an FPSO With Sponsons

2016 ◽  
Author(s):  
Babak Ommani ◽  
Nuno Fonseca ◽  
Trygve Kristiansen ◽  
Christopher Hutchison ◽  
Hanne Bakksjø
Keyword(s):  
Free Surface ◽  
Two Dimensions ◽  
The Body ◽  
Proximity Effects ◽  
Roll Damping ◽  
Damping Coefficients ◽  
Free Decay ◽  
Strip Theory ◽  
Decay Tests ◽  
Bilge Keel

The bilge keel induced roll damping of an FPSO with sponsons is investigated numerically and experimentally. The influence of the bilge keel size, on the roll damping is studied. Free decay tests of a three-dimensional ship model, for three different bilge keel sizes are used to determine roll damping coefficients. The dependency of the quadratic roll damping coefficient to the bilge keel height and the vertical location of the rotation center is studied using CFD. A Navier-Stokes solver based on the Finite Volume Method is adopted for solving the laminar flow of incompressible water around a section of the FPSO undergoing forced roll oscillations in two-dimensions. The free-surface condition is linearized by neglecting the nonlinear free-surface terms and the influence of viscous stresses in the free surface zone, while the body-boundary condition is exact. An averaged center of rotation is estimated by comparing the results of the numerical calculations and the free decay tests. The obtained two-dimensional damping coefficients are extrapolated to 3D by use of strip theory argumentations and compared with the experimental results. It is shown that this simplified approach can be used for evaluating the bilge keel induced roll damping with efficiency, considering unconventional ship shapes and free-surface proximity effects.

Experimental Investigation of Hydrodynamic Coefficients of a Wave-Piercing Tumblehome Hull Form

2007 ◽  
Author(s):  
Christopher C. Bassler ◽  
Jason B. Carneal ◽  
Paisan Atsavapranee
Keyword(s):  
Experimental Investigation ◽  
System Modeling ◽  
Flow Simulation ◽  
Roll Damping ◽  
Hull Form ◽  
Damping Coefficients ◽  
Calm Water ◽  
Forced Roll ◽  
Future Work ◽  
Damping Model

A systematic series of calm-water forced roll model tests were performed over a range of forward speeds using an advanced tumblehome hull form (DTMB model #5613-1) to examine the mechanisms of roll damping. This experimental investigation is part of an ongoing effort to advance the capability to assess seakeeping, maneuvering, and dynamic stability characteristics of an advanced surface combatant. The experiment was performed to provide data for development and validation of a semi-empirical roll damping model for use in validation of ship motion and viscous flow simulation codes, as well as to provide a basis for future work with additional experiments, contributing to the development of an improved analytical roll damping model. Two hull configurations were tested: barehull with skeg, and bare hull with skeg and bilge keels. Measurements of forces and moments were obtained over a range of forward speeds, roll frequencies, and roll amplitudes. Stereo particle-image velocimetry (SPIV) measurments were also taken for both zero and forward speeds. Test data was used to calculate added mass/inertia and damping coefficients. Two different system modeling techniques were used. The first method modeled the system as an equivalent linearly-damped second-order harmonic oscillator with the time-varying total stiffness coefficient considered linear. The second technique used equivalent linear damping, including higher-order Fourier components, and a non-linear stiffness formulation. Results are shown, including plots of added inertia and damping coefficients as functions of roll frequency, roll amplitude, and forward speed and SPIV measurements. Trends from the experimental data are compared to results from traditional component roll damping formulations for conventional hull from geometries and differences are discussed.

The Meshless Dynamic Relaxation Techniques for Simulating Atomic Structures of Materials

2002 ◽  
Author(s):  
Li Pan ◽  
Don R. Metzger ◽  
Marek Niewczas
Keyword(s):  
Potential Energy ◽  
Boundary Conditions ◽  
Internal Force ◽  
Crystalline Lattice ◽  
Structural Defects ◽  
Convergence Criterion ◽  
Relaxation Techniques ◽  
Dynamic Relaxation ◽  
Time Step ◽  
Eam Potential

Traditionally, Molecular Dynamics combined with pair potential functions or the Embedded Atom Method (EAM) is applied to simulate the motion of atoms. When a defect is generated in the crystalline lattice, the equilibrium of atoms around it is destroyed. The atoms move to find a new place where the potential energy in the system is minimum, which could result in a change of the local atomic structure. The present paper introduces new Dynamic Relaxation algorithm, which is based on explicit Finite Element Analysis, and pair or EAM potential function, to find equilibrium positions of the block of atoms containing different structural defects. The internal force and stiffness at the atoms (nodes) are obtained by the first and second derivatives of the potential energy functions. The convergence criterion is based on the Euclidean norm of internal force being close to zero when the potential energy is minimum. The damping ratio affects the solution path so that different damping ratios could lead to different minimum potential energy and equilibrium shapes. The numerical responses and results by applying free boundary conditions and certain periodic boundary conditions are presented. The choice of scaled mass of atoms, proper time step and damping appropriate for the efficient and stable simulation is studied.

Parameter Study of a Thermal Analysis of a Bead-on-Plate Weld

2018 ◽  
Author(s):  
Spyridon A. Alexandratos ◽  
Lei Shi ◽  
Noel P. O’Dowd
Keyword(s):  
Grain Size ◽  
Steel Plate ◽  
Martensitic Steel ◽  
Sensitivity Study ◽  
Parameter Study ◽  
Stainless Steel 316L ◽  
Material Microstructure ◽  
Simulation Techniques ◽  
Subsequent Effect ◽  
Thermal Histories

Failures in engineering components operating at high temperature often initiate in welded joints, particularly in the heat-affected zone (HAZ) adjacent to a weld. This is due to the inhomogenous microstructure of the weld and adjacent material due to the different thermal histories experienced during the weld cycle. It is therefore important to accurately predict the temperature distributions arising during welding and the subsequent effect on material microstructure. The NET TG1 bead-on-plate weld geometry is examined in this work. This geometry is a single weld bead laid on the surface of an AISI 316L austenite steel plate. Experimental data from the TG1 study are available to validate different weld simulation techniques. Here, a sensitivity study to the thermal properties is carried out and the influence on the HAZ temperatures and grain size is examined. The study shows that the conductivity and the specific heat capacity significantly affect the temperature prediction in the HAZ with a similar influence on predicted grain size following welding. Results are presented for a stainless steel (316L) and a martensitic steel (P91) plate.

A computational method for earthquake cycles within anisotropic media

2019 ◽  
Vol 219 (2) ◽  
pp. 816-833 ◽  
Author(s):  
Maricela Best Mckay ◽  
Brittany A Erickson ◽  
Jeremy E Kozdon
Keyword(s):  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Elastic Solid ◽  
Elliptic Partial Differential Equation ◽  
Computational Method ◽  
Parameter Study ◽  
Time Step ◽  
Time Stepping ◽  
Static Elasticity ◽  
Earthquake Cycles

SUMMARY We present a numerical method for the simulation of earthquake cycles on a 1-D fault interface embedded in a 2-D homogeneous, anisotropic elastic solid. The fault is governed by an experimentally motivated friction law known as rate-and-state friction which furnishes a set of ordinary differential equations which couple the interface to the surrounding volume. Time enters the problem through the evolution of the ordinary differential equations along the fault and provides boundary conditions for the volume, which is governed by quasi-static elasticity. We develop a time-stepping method which accounts for the interface/volume coupling and requires solving an elliptic partial differential equation for the volume response at each time step. The 2-D volume is discretized with a second-order accurate finite difference method satisfying the summation-by-parts property, with boundary and fault interface conditions enforced weakly. This framework leads to a provably stable semi-discretization. To mimic slow tectonic loading, the remote side-boundaries are displaced at a slow rate, which eventually leads to earthquake nucleation at the fault. Time stepping is based on an adaptive, fourth-order Runge–Kutta method and captures the highly varying timescales present. The method is verified with convergence tests for both the orthotropic and fully anisotropic cases. An initial parameter study reveals regions of parameter space where the systems experience a bifurcation from period one to period two behaviour. Additionally, we find that anisotropy influences the recurrence interval between earthquakes, as well as the emergence of aseismic transients and the nucleation zone size and depth of earthquakes.

scholarly journals WRF Hindcasts of Cold Front Passages over the ARM Eastern North Atlantic Site: A Sensitivity Study

2018 ◽  
Vol 146 (8) ◽  
pp. 2417-2432 ◽  
Author(s):  
Fayçal Lamraoui ◽  
James F. Booth ◽  
Catherine M. Naud
Keyword(s):  
Case Studies ◽  
North Atlantic ◽  
Cold Front ◽  
Domain Boundary ◽  
Wrf Model ◽  
Domain Size ◽  
Sensitivity Study ◽  
Spectral Nudging ◽  
Eastern North Atlantic ◽  
The Impact

Abstract The present study explores the ability of the Weather Research and Forecasting (WRF) Model to accurately reproduce the passage of extratropical cold fronts at the DOE ARM eastern North Atlantic (ENA) observation site on the Azores. An analysis of three case studies is performed in which the impact of the WRF domain size, position of the model boundary relative to the ENA site, grid spacing, and spectral nudging conditions are explored. The results from these case studies indicate that model biases in the timing and duration of cold front passages change with the distance between the model domain boundary and the ENA site. For these three cases, if the western model boundary is farther than 1500 km from the site, the front becomes too meridional and fails to reach the site, making 1000 or 1500 km the optimal distances. In contrast, integrations with small distances (e.g., 500 km) between the site and domain boundaries have inadequate spatial spinup (i.e., the domain is too small for the model to properly stabilize). For all three cases, regardless of domain size, the model has biases in its upper-level circulation that impact the position and timing of the front. However, this issue is most serious for 4000-km2 domains and larger. For these domains, prolonged spectral nudging can correct cold front biases. As such, this analysis provides a framework to optimize the WRF Model configuration necessary for a realistic hindcast of a cold front passage at a fixed location centered in a domain as large as computationally possible.

Nonlinear-wave effects on fixed and floating bodies

1982 ◽  
Vol 120 ◽  
pp. 267-281 ◽  
Author(s):  
Michael De St Q. Isaacson
Keyword(s):  
Initial Conditions ◽  
Integral Equation Method ◽  
Ocean Waves ◽  
Floating Bodies ◽  
Time Step ◽  
Fluid Forces ◽  
Time Stepping ◽  
Transient Problem ◽  
Wave Effects ◽  
Incident Waves

A numerical method for calculating the interaction of steep (nonlinear) ocean waves with large fixed or floating structures of arbitrary shape is described. The interaction is treated as a transient problem with known initial conditions corresponding to still water in the vicinity of the structure and a prescribed incident waveform approaching it. The development of the flow, together with the associated fluid forces and structural motions, are obtained by a time-stepping procedure in which the flow at each time step is calculated by an integral-equation method based on Green's theorem. A few results are presented for two reference situations and these serve to illustrate the effects of nonlinearities in the incident waves.

Sign in / Sign up

Export Citation Format

Share Document