|
[Sponsors] |
Choosing a stretching factor in DIVEMesh - wave number |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
September 6, 2022, 04:13 |
Choosing a stretching factor in DIVEMesh - wave number
|
#1 |
New Member
Felix Daddo
Join Date: Aug 2022
Location: Sydney, Australia
Posts: 25
Rep Power: 3 |
High there,
I've been reading through the paper "Reduction of the wave propagation error of a sigma grid based numerical tank using a vertical spacing based on the constant truncation error" by Pakozdi et al. (2021.https://doi.org/10.1016/j.oceaneng.2021.109741) and am trying to work out the way to compute the optimum grid spacing profile so that I can compare this profile with that computed by the grid stretching function in DIVEMesh and estimate the ideal stretching factor for my simulation. I am confused in the paper what the "wave number" is and how this property is found and used to normalise a number of variables. I am trying to understand after kdz (normalised grid spacing) is computed, it can be converted back to z (not normalised) to compared to the DIVEMesh profile. I am also not sure what value to set k (wave number) to when trying to numerically solve the equation for kdz. I hope this makes sense and there is a simple way to wrap my head around this so that I can find the optimum grid spacing. Cheers, Felix |
|
September 6, 2022, 04:29 |
|
#2 |
Senior Member
Arun Kamath
Join Date: Nov 2014
Location: Trondheim, Norway
Posts: 265
Rep Power: 13 |
Hi Felix,
I am not sure if I understood your question correctly. Wave number k= 2*pi/L. In the case of regular waves, this is the wavelength L of the desired input wave. In the case of irregular waves, 5.4, 5.5 and 5.6 in the paper you mention offers some guidelines for the basis of L for the stretching. Then, to use your custom stretching in Divemesh, you need to then create a file with one column going 0 to 1 based on the optimal stretching. This file is to named z-spacing.dat and the stretching for the z-direction using this file is activated by B 103 10 Hope that helps
__________________
Arun X years with REEF3D Last edited by kamath; September 6, 2022 at 04:32. Reason: added info for irregular waves |
|
September 7, 2022, 03:24 |
|
#3 |
New Member
Felix Daddo
Join Date: Aug 2022
Location: Sydney, Australia
Posts: 25
Rep Power: 3 |
Thanks Kamath,
That makes a lot of sense and clear it all up! Cheers Felix |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
use the message in macro DEFINE_PROFILE with parallel processor | alireza_T | Fluent UDF and Scheme Programming | 3 | May 11, 2022 02:08 |
[Other] Can't Shake Erros: patch type 'patch' not constraint type 'empty' | BrendaEM | OpenFOAM Meshing & Mesh Conversion | 12 | April 3, 2022 18:32 |
GenerateVolumeMesh Error - Surface Wrapper Self Interacting (?) | AndreP | STAR-CCM+ | 10 | August 2, 2018 07:48 |
Compressor Simulation using rhoPimpleDyMFoam | Jetfire | OpenFOAM Running, Solving & CFD | 107 | December 9, 2014 13:38 |
DecomposePar unequal number of shared faces | maka | OpenFOAM Pre-Processing | 6 | August 12, 2010 09:01 |