Study of Polarized Wave with a Hydrodynamic Model and Fourier Spectral Method

1 Department of Fundamental and Applied Science, Universiti Teknologi PETRONAS, 31750 Tronoh, Perak, Malaysia 2 UTM Centre for Industrial and Applied Mathematics, Department of Mathematical Sciences, Faculty of Science, Universiti Teknologi Malaysia (UTM), 81310 Johor Bahru, Malaysia 3 Department of Mathematics, Faculty of Science, University of Kordofan, North Kordofan State, El-Obeid 51111, Sudan


Introduction
Via the massive energy released by the faults, energy of the earthquake is usually measured in magnitude or Rachel scale.Released energy in less than one percent of is said transmitted to form Tsunamis [1].Therefore, we assume that the Tsunami and seismic waves should exist together in this paper during the derivation of polarized model.The notion is strong ground motion by the wave field polarization that generates Tsunami.Reflection being transmitted or polarization gives strong earth motion.Seismic waves due to strong polarized earthquakes propagating in seabed layers will reduce soil stiffness and increase the energy dissipation into the soil [2,3].
Lately, recording of polarization of the Earth's teleseismic wave fields and strong ground motion is imaged by applying a vector stacking method to three component broadband seismograms [4].Polarization effect of the P waves and S waves was also studied analytically for the dispersion and attenuation effects in fluid saturated medium [5][6][7][8].
Besides, wave field polarization may both amplify the displacement and increase the attenuation.It is found that multiple wave fields impinging each other give amplification effect while wave field separation gives the negative damping effects [9].By considering the polarization of seismic wave to the generation of water displacement, we analyze the effect of strong polarized seabed motion towards the generation of Tsunamis.
There are experimental approaches for exploring nonlinear wave's motion [15][16][17][18].Experimental results by Hwang et al. [19] indicated that the Tsunamis run-ups are directly proportional to the incline beach plane and the run-ups are higher than the run-ups height proposed by Synolakis [13], and Hwang et al. [19] had conducted a study in water tank with 1/20 beach plane.Chang and Hwung [17] and Hsiao et al. [18] implemented water tank experiments for wave propagation in a water tank (300 m × 5 m × 5.2 m) with incline beach planes 1/40 and 1/60 at Tainan Hydraulics Laboratory of National Cheng Kung University.Besides, the wave flume (2 m × 2 m × 200 m) is as well implemented by Chang and Hwung [17] with the purpose of reducing the scale effect.
There are approaches to modeling run-up on beach: the nonlinear shallow water equations (nondispersive system) and Boussinesq equations to include the dispersion in weakly nonlinear waves.More sophisticated methods such as Navier-Stokes equations are used to simulate wave run-up, that is, UNESCO TSUNAMI model [20].The tsunami model TsunAWI by Alfred Wegener Institute simulates all stages of a tsunami from the origin and the propagation in the ocean to the arrival at the coast and the inundation on land.Solved nonlinear shallow water equations on an unstructured finite element grid allows changing the resolution seamlessly between a coarse grid in the deep ocean and a fine representation of coastal structures.Besides, it provides fast computation [21].
Here, we use the dimensionless equation of hydrodynamics (1) which is derived in accordance with Boussinesq equations by Thierry and Michel [22]:

Polarization
However, a simple calculation based on considering the tsunami has a wave profile that reads [23]  0 = 2sech 2 () , given for the -axis direction propagation only which is not right to explain the waves propagates in the - directions or run-ups at the incline beach plane.Figures 1(a) and 1(b) illustrate the difference between constant amplitude and variable amplitude for nonlinear waves.The conventional wave (3) that is used in modelling is shown in Figure 1(a).Figure 1(b) illustrates the schematic model of this research.The nonlinear wave is variable because of the reflection at beach plane.Incidence of seismic plane wave is not left behind whereas both the nonlinear waves and seismic waves make a polarization angle ] during the Tsunami run-ups.
Complying with the Snell laws, the variation of density will lead to the seismic wave diffraction during incidence.Here, the diffraction is due to polarization and consequently the purposed objective is to derive the polarized nonlinear waves and seismic waves as initial waves in hydrodynamics simulation.In next session, the derivation for polarization is introduced before the application of polarization in the nonlinear waves and the seismic waves.Then, numerical simulation and discussion will be shown.
In Figure 1,  is the incident wave coefficient with  = (, ) and  is the coefficients for wave height variation  and distance travelled .Initial wave height (, ) with minimal variation  reads  (, )  ≅  (, 0)  =  () . ( The wave height depends on the distance travelled or  = ().Essentially, we are going to find the polarization angle V from total variation angle   .In accordance with Figure 1, shock dynamics relations with the angle (, ) between the -axis and the ray were derived: In this section, we will show the derivation for polarization along with the angle between the beach plane and -axis (, ) and polarization angle ] according to Han and Yin [24].The total variation angle   is For the wave orientation, it is appropriate to write the variation as Respectively,  and  are the variation in horizontal and vertical direction.Hence, relations (5) become Subsequently, the polarization angle ] will be extracted from the (8a) and (8b), given that   =  + ].
By inserting  = () into (8a), the relations read Equations ( 9a) and (9b) can be written in characteristic form.By introducing (9a) ±  × (9b) with  being the coefficient of incident waves, one gets the relations of variation in characteristic form: Equations (10a) and (10b) will limit minimum beach plane gradient to  = 1 such that The characteristic curves of incident waves,  and Riemann invariants  and  are Consequently, the total variation angle is From (10a) and (10b), it is noted that The value for ] > 1 will be used for the numerical simulation later.By solving (10a) and (10b), the total variation due to polarization reads Next, the derivation of initial waves with polarization will be explained.

Initial Waves
In this modelling, we assume that the nonlinear waves and the seismic plane waves will share single polarization angle V.
For the cases of multiple polarization angle ], it will be studied in further research.
For deriving the polarized nonlinear waves, the amplitude in (3) is replaced with / or tan( + ]) as illustrated in Figure 1(b).The initial waves read and the polarized seismic waves are where  is the phase velocity of the wave.We are assuming that the initial seismic waves are plane waves with the amplitude depending on the polarization or /.Both ( 17) and ( 18) are used as initial waves in this analysis in conjunction with the hydrodynamics equation (1).Fourier spectral method and Runge-Kutta 4 methods are used to solve the identified model.

Fourier Spectral and Runge-Kutta 4 Methods
Hydrodynamics equation ( 1) is rewritten as    Spatial part of hydrodynamics equation will be discretized by using Fourier spectral method.The Fourier space gives 20) or, in standard form of   =  + (, ), where  and  are the linear and nonlinear operators, respectively.Then, we use Runge-Kutta 4 for time stepping.Initial waves from ( 17) and ( 18) will be used to analyze the effect of polarization.

Results and Discussions
In this section, the configurations of hydrodynamics models as shown in Figures 1(a) and 1(b) will be shown.The initial waves of (3) will produce the initial waves with dimensionless amplitude 2. However, the dispersion of the waves is still significant in both Figures 2(a) and 2(b).There is no beach plane or run-ups in Figures 2(a For these flow cases, we should see a beach plane when  = 0 as shown in Figures 3(a Consequently, dispersion is reduced when the size of the domain increases.By comparing Figures 3(a) and 3(b), the dispersive waves are more significant in Figure 3(a).Nevertheless, the initial wave comes with the displacement height around 1.2 units which is higher than the initial wave as shown in Figure 3   However, the reflected waves are still visible.In Figures 2(a) and 2(b), the reflected waves are shown to diminish the incoming waves when the domain increases.In contrast, the incoming wave run-ups on beach plane without the polarization effect are greater and higher in displacement when the domain increases as shown in Figures 3(a The nonpolarized initial waves were shown previously.Here, the consequence of polarized initial waves in hydrodynamics will be shown.Figures 4(a  The polarization effect in Tsunami wave is referred to the situation for the near-shore Tsunami.From the source of the earthquake, the Indian Ocean Tsunami is known as near-shore Tsunami.If the polarization effect is ignored, the inundation in Sumatra in 2004 is less disastrous as shown in Figure 3(b).However, the inundation is deeper and comes in earlier once the polarization effect is considered as shown in Figure 4(b).

Conclusion
The development of hydrodynamics model with the existence of seismic waves is a promising task.Type of setup especially the structure of the bottom geography is important.The nonlinear wave in the domain without an inclined beach plane is responsive to the reflected waves.However, an incline beach plane will show no reflected waves when the polarization in waves is considered.
The size of the domain denotes the roles of polarization.The Tsunami wave propagates in the smaller domain exclusive of polarization effect is longer in size, but the inundation is less profound once compared with the Tsunami wave embedded with polarization.With the polarization effects, the tsunami wave is catastrophic: the height of the displacement is greater than before and the inundation comes earlier and deeper for which we take for granted the hazard is being brought forward by the high displacement.

Figure 1 :
Figure 1: The nonlinear waves.(a) The nonlinear wave propagates with amplitude 2 in -direction.(b)The nonlinear waves propagate with variation amplitude /, initiated by reflection at incline beach plane.There is an incidence of seismic plane waves,  from the seabed whereas the seismic plane waves suffer diffraction initially.

Figure 2 :
Figure 2: Time evolution for the hydrodynamics equation together with initial wave  = 2 sech 2 ().Time runs from 0 at the bottom of the figure to 100 at the top.

Figure 2
shows the initial waves generated for two different domains: (a) [0, 32pi], (b) [0, 128pi].Time runs from 0 at the bottom of the Figure to 100 at the top and the initial wave's propagation began from the right of both Figures 2(a) and 2(b).From the Figures 2(a) and 2(b), stability of the initial waves is shown to be greatly influenced by the size of the domain; the reflected waves tend to diminish the waves when the size of the domain increases as shown in Figure 2(b).
Figures3(a) and 3(b) are plotted to show the run-ups by initial waves(22).For these flow cases, we should see a beach plane when  = 0 as shown in Figures3(a) and 3(b).Figures 3(a) and 3(b) demonstrate the time evolution of the wave's run-ups without polarization by letting ] = 0.The angle between beach plane and -axis is  = 0.1.Consequently, dispersion is reduced when the size of the domain increases.By comparing Figures3(a) and 3(b), the dispersive waves are more significant in Figure3(a).Nevertheless, the initial wave comes with the displacement height around 1.2 units which is higher than the initial wave as shown in Figure3(b) with the displacement height around 0.33 in dimensionless.When the length of the wave is evaluated, the wave shown in Figure3(b) is found longer.Apparently, the weak dispersion has led to the long waves in Figure3(b).
Figures3(a) and 3(b) are plotted to show the run-ups by initial waves(22).For these flow cases, we should see a beach plane when  = 0 as shown in Figures3(a) and 3(b).Figures 3(a) and 3(b) demonstrate the time evolution of the wave's run-ups without polarization by letting ] = 0.The angle between beach plane and -axis is  = 0.1.Consequently, dispersion is reduced when the size of the domain increases.By comparing Figures3(a) and 3(b), the dispersive waves are more significant in Figure3(a).Nevertheless, the initial wave comes with the displacement height around 1.2 units which is higher than the initial wave as shown in Figure3(b) with the displacement height around 0.33 in dimensionless.When the length of the wave is evaluated, the wave shown in Figure3(b) is found longer.Apparently, the weak dispersion has led to the long waves in Figure3(b).
Figures3(a) and 3(b) are plotted to show the run-ups by initial waves(22).For these flow cases, we should see a beach plane when  = 0 as shown in Figures3(a) and 3(b).Figures 3(a) and 3(b) demonstrate the time evolution of the wave's run-ups without polarization by letting ] = 0.The angle between beach plane and -axis is  = 0.1.Consequently, dispersion is reduced when the size of the domain increases.By comparing Figures3(a) and 3(b), the dispersive waves are more significant in Figure3(a).Nevertheless, the initial wave comes with the displacement height around 1.2 units which is higher than the initial wave as shown in Figure3(b) with the displacement height around 0.33 in dimensionless.When the length of the wave is evaluated, the wave shown in Figure3(b) is found longer.Apparently, the weak dispersion has led to the long waves in Figure3(b).
) and 4(b) illustrate the polarized initial waves with  = 0.1 and ] = 1.Similarly, the initial waves are simulated in two domains: (a) [0, 32pi] and (b) [0, 16pi].Apparently, the Figures4(a) and 4(b) show dissimilar illustration as Figures2 and 3, respectively.The initial waves when  = 0 show extraordinary initial displacement with displacement height around 7 in dimensionless as shown in Figure4(a).Nevertheless, the displacement height is around 6 units as shown in Figure4(b).These high displacements have totally ignored the existence of reflected waves.By comparing the inundation time of Figure3(a) with Figure4(a), the paces of the initial waves are quicker.This indicates that the velocity of the waves is higher.The initial waves arrive at the shore or  = 100 when  = 40 units.When the domain reduces, the arrival time of the waves at the shore is reduced as well.

Figure 4 (
b) shows the inundation time at  = 50.
in which the dispersion term   represents shock formation produced by convective term.  represents the steepening of the wave while   represents the depressiveness of the wave.The last phase of a tsunami is its run-ups and inundation.With linear speed  0 = √ℎ, dimensionless variables are