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