windVelocity = (0.0, 0.0)
depth = 1.0
inflowHeightMean = 1.0
inflowVelocityMean = (0.0, 0.0)
period = 1.33
waveheight = 0.062
amplitude = waveheight / 2.0
wavelength = 2.77
waveDir = numpy.array([1, 0, 0])
g = numpy.array([0, -9.81, 0])
waves = WT.RandomWaves(
Tp=period, # Peak period
Hs=waveheight, # Height
d=depth, # Depth
fp=1. / period, #peak Frequency
bandFactor=2.0, #fmin=fp/Bandfactor, fmax = Bandfactor * fp
N=101, #No of frequencies for signal reconstruction
mwl=inflowHeightMean, # Sea water level
waveDir=waveDir, # waveDirection
g=g, # Gravity vector, defines the vertical
gamma=1.0, #Pierson Moskowitz spectum for gamma=1.0
spec_fun=WT.JONSWAP)
# Discretization -- input options
genMesh = True
movingDomain = False
applyRedistancing = True
useOldPETSc = False
useSuperlu = False
timeDiscretization = 'be' #'vbdf'#'be','flcbdf'
spaceOrder = 1
useHex = False
0.00001078,
0.00000151,
0.00000023,
0.00000007
],
Bcoeff=[
0.05395079, 0.00357780, 0.00020506, 0.00000719, -0.00000016,
-0.00000005, 0.00000000, 0.00000000
])
elif opts.wave_type == 'single-peaked':
waves = WT.RandomWaves(
Tp=period, # Peak period
Hs=waveheight, # Height
mwl=inflowHeightMean, # Sea water level
depth=depth, # Depth
#fp = 1./period, #peak Frequency
bandFactor=2.0, #fmin=fp/Bandfactor, fmax = Bandfactor * fp
N=101, #No of frequencies for signal reconstruction
#mwl = inflowHeightMean, # Sea water level
waveDir=waveDir, # waveDirection
g=g,
spectName="JONSWAP") # Gravity vector, defines the vertical
#gamma=3.3,
#spec_fun = JONSWAP)
elif opts.wave_type == 'double-peaked':
waves = WT.DoublePeakedRandomWaves( #Tp = period, # Peak period
Hs=waveheight, # Height
d=depth, # Depth
fp=1. / period, #peak Frequency
bandFactor=2.0, #fmin=fp/Bandfactor, fmax = Bandfactor * fp
N=101, #No of frequencies for signal reconstruction
mwl=inflowHeightMean, # Sea water level
#wave generator
windVelocity = (0.0,0.0)
inflowHeightMean = 1.0
inflowVelocityMean = (0.0,0.0)
period = 1.94
omega = 2.0*math.pi/period
waveheight = 0.1
amplitude = waveheight/ 2.0
wavelength = 5.0
k = 2.0*math.pi/wavelength
waves = WaveTools.RandomWaves(Tp = 1.94,
Hs = 0.1,
d = 1.0,
fp = 1.0/1.94,
bandFactor = 2.0,
N = 101,
mwl = 1.0,
g = 9.8)#shouldn't mwl = d always?
# Discretization -- input options
genMesh=True
movingDomain=False
applyRedistancing=True
useOldPETSc=False
useSuperlu=False
timeDiscretization='be'#'be','vbdf','flcbdf'
spaceOrder = 1
useHex = False
useRBLES = 0.0
useMetrics = 1.0