Experiments on exactly computing non-linear energy transfer rate in MASNUM-WAM*
|
|
- Ginger Caldwell
- 6 years ago
- Views:
Transcription
1 Chinese Journal of Oceanology and Limnology Experiments on exactly computing non-linear energy transfer rate in MASNUM-WAM* JIANG Xingjie ( 江兴杰 ) 1, 2, **,WANG Daolong ( 王道龙 ) 2, GAO Dalu ( 高大鲁 ) 1, 2, ZHANG Tingting ( 张婷婷 ) 2 1 Ocean University of China, Qingdao 26671, China 2 First Institute of Oceanography, State Oceanic Administration, Qingdao 26661, China Received Jan. 8, 215; accepted in principle Mar. 27, 215; accepted for publication Jun. 5, 215 Chinese Society for Oceanology and Limnology, Science Press, and Springer-Verlag Berlin Heidelberg 216 Abstract The Webb-Resio-Tracy (WRT) method for exact computation of the non-linear energy transfer rate was implemented in MASNUM-WAM, which is a third-generation wave model solving the discrete spectral balance equation. In this paper, we describe the transformation of the spectral space in the original WRT method. Four numerical procedures were developed in which the acceleration techniques in the original WRT method, such as geometric scaling, pre-calculating, and grid-searching, are all reorganized. A series of numerical experiments including two simulations based on real data were performed. The availability of such implementation in both serial and parallel versions of the wave model was proved, and a comparison of computation times showed that some of the developed procedures provided good efficacy. With exact computation of non-linear energy transfer, MASNUM-WAM now can be used to perform numerical experiments for research purposes, which augurs well for further developments of the model. Keyword: nonlinear energy transfer; the WRT method; geometric scaling; MASNUM-WAM 1 INTRODUCTION It is nowadays widely accepted that weak-resonance non-linear energy transfer between sea waves plays a very important role in energy dispersion. The theoretical framework of the wave-wave nonlinear interaction was developed in the 196s (Hasselmann, 1962, 1963a). A 6-fold Boltzmann integral was formulated to describe these non-linear interactions and many researchers spent much effort developing various schemes to compute the integral on computers (Hasselmann, 1963b; Webb, 1978; Masuda, 198). Termed exact computations or full Boltzmann integral (FBI), these schemes considered all resonance configurations. Because of their complexity and time demands, these exact methods are not suitable for inclusion in operational spectral wave models. Various approximate methods have been developed to improve the computational efficiency (Barnett, 1968; Ewing, 1971; Hasselmann and Hasselmann, 1981). Finally, the discrete interaction approximation (DIA) (Hasselmann and Hasselmann, 1985a; Hasselmann et al., 1985) was widely accepted and triggered the development of the third-generation wave models, like WAM (The Wamdi Group, 1988), the WaveWatch (Tolman, 1991) and the SWAN model (Booij et al., 1999). The DIA is much more efficient than exact computations because it considers only one pair of resonance configurations. The DIA preserves a few but important characteristics of the full solution, such as the slow downshifting of the peak frequency and the shape stabilization during wave growth. Nevertheless, the DIA also shows some deficiencies, such as allocating too much energy from the spectral region near the spectral peak to higher frequencies (van Vledder et al., 2). The known deficiencies of the DIA in a prediction model will not influence bulk wave parameters, such as the significant wave height, the zero-crossing period, and the mean wave direction, but it does distort the shape of energy spectra, which are vital in research on waves. * Supported by the State Oceanic Administration Marine Science Foundation for Young Scientists (Nos , ) ** Corresponding author: jiangxj@fio.org.cn
2 CHIN. J. OCEANOL. LIMNOL., Therefore, it is preferable to include an exact computation or FBI in a wave model for purposes of research. There are some third-generation wave models that have included exact computations, such as WaveWatchIII (Tolman, 22), SWAN (Booij et al., 24), CREST (Ardhuin et al., 21), and PROWAM (Monbaliu et al., 1999). The performance of exact computations has been well understood (Hasselmann and Hasselmann, 1985c; van Vledder, 26), but the effects of such methods in a certain wave model still need to be studied individually. In this paper, we describe the implementation of the Webb-Resio-Tracy (WRT) method (Webb, 1978; Tracy and Resio, 1982; Resio and Perrie, 1991) into MASNUM (the Key Laboratory of Marine Science and Numerical Modeling, First Institute of Oceanography, SOA, China) wave model (Yuan et al., 1992a, b; Yang et al., 25). MASNUM- WAM is a third-generation model which calculates energy spectrum balance equation using a characteristics inlaid scheme; the dissipation source function of this model is based on the wave broken statistic theory (Yuan et al., 1986); MASNUM-WAM adopts wave number spectral space which is coherent with the wave-wave nonlinear interaction theory (Hasselmann, 1962, 1963a). The WRT programs we used are free software distributed by Gerbrant Ph. van Vledder, ALKYON Hydraulic Consultancy & Research; the latest version is 5.4. In Section 2, we briefly introduce the theoretical background of the WRT method, the modifications we made and numerical procedures we developed. In Section 3, we describe a series of numerical experiments to illustrate the ease and efficiency of the implementation. Finally, in the last section, we discuss current and future work. 2 THEORETICAL BACKGROUND AND METHODS 2.1 Theoretical background The WRT method The WRT method is widely accepted as a mature and stable exact-computation scheme for non-linear energy transfer. Following Hasselmann (1962), the original 6-fold Boltzmann integral is n( k1 ) = G( k1, k2, k3, k4) P( k1, k2, k3, k4) t k1 + k2 ω1 + ω2 δ δ dk2dk3dk, (1) 4 k3 k 4 ω 3 ω 4 where n( k i ) is the action density and k i the ith interacting vector wave number corresponding to angular frequency ω i satisfying the dispersion relation ω 2 =gktanh(kd), (2) in which g is the acceleration of gravity and d local depth. The wave number and frequency resonance condition k1 + k2 = k3 + k4 ω1 + ω2 = ω3 + ω4, (3) are contained in the two delta-functions of the integral. The coupling coefficient G is a complicated function of wave numbers and frequencies, and P( k1, k2, k3, k4) is the density function Pk ( 1, k2, k3, k4) = nk ( 1) nk ( 2) nk ( 3) + nk ( 2) nk ( 3) nk ( 4) nk ( 1) nk ( 3) nk ( 4) nk nk nk. (4) ( 1) ( 2) ( 4) A key aspect in the WRT method is to consider the integration space for each ( k1, k3) combination n( k1 ) = 2 T( k1, k3) dk3 t, (5) where W T( k1, k3) = G P ( k1, k3, k4) ds n θ, (6) with W = ω1+ ω2 ω3 ω k1 + k2 k. The locations in 3 wave number space where W= trace out a closed curve or locus. Unit vector s is along the interaction locus, and n is normal to it. The function θ is defined as 1 when k1 k3 k1 k4 θ ( k1, k3, k4) = when k1 k3 > k1 k 4. (7) The above summarizes the basic theory of Webb (1978); additional work has been done by Tracy and Resio (1982) and Resio and Perrie (1991) Geometric scaling For convenience, we use a geometric scaling spectral space, in which k i+1 =λk i with λ>1. Such scaling provides a higher spectral resolution near the peak of the spectrum and less resolution in the highfrequency tail. An optimization algorithm was developed based on the geometric spacing spectral
3 JIANG et al.: WRT method in MASNUM-WAM grid (Tracy and Resio, 1982). In the geometrical scaling k-space, if k ' ' 1 k3 = λ k1 k3, the locus based on ( ' ' k1, k3 ) has be proven to scale linearly with λ for the locus based on ( k1, k3 ). Specifically, for each point along the original locus, a geometrically similar ' point exists in a scaling locus such as k2 = λk2. We easily deduce that k ' 4 = λk4, according to the wave number resonance condition Eq.3. For each combination of ( ' ' ' ' ) ( k ) 1, k2, k3, k4 = λ k1, k2, k3, k4, we have G( k ) ', k ', k ', k ' = λ 6 G( k, k, k, k ) ' W W 2 = λ n n d s' = λ d s. (8) Therefore, with some basic pairs of k 1 and k3, a -1 table of all possible values for G W / n d s can be 15 λ 2 j constructed easily by multiplying by a simple -based factor. Nevertheless, the geometric scaling can only be used when the dispersion relation ω 2 =gk is adopted. Because the coupling coefficient and locus solutions are all based on wave numbers, the Jacobian term, in which W=ω 1 +ω 2 ω 3 ω 4 is related to the angular frequency. Geometrical spacing is not allowed in both wave-number space and angular-frequency space in Eq.2 where a finite water depth is considered. Mismatching in the geometrical spacing may lead to incorrect results in Eq.8 and W= is the basic resonance condition that cannot be avoided by any scheme Pre-calculation and grid-searching Another important property of the WRT method is that all of the loci, interpolation coefficients, coupling coefficients, and Jacobian terms are only associated with the (wave number) spectral grid and water depth, but have nothing to do with the shape of the spectrum. For a given discrete spectral grid and water depth, all of those parameters mentioned above can be precomputed and stored in memory or in a database on disks, and then a search-and-read process could be performed instead of a repeated calculation if the current spectral grid and depth match. The grid-searching technique is based on some pre-calculated depths. New computations for the local depth are avoided, and the program tries to find a nearest pre-calculated depth record instead. The final results will then be scaled by a factor based on the local and nearest depths. More details can be found in Section Modifications Spectral space transformation In the originally distributed WRT subroutines, the angular-frequency spectral space (ω-space) is adopted, which is also adopted in WaveWatchIII and SWAN models. However, as we have seen in Section 2.1, the core integral processes are all based on the wave number spectral space (k-space). The action density terms can be transformed between ω-space and k-space following d k nω ( ωθ, ) = nk ( k, θ) dω dω nk ( k, θ) = nω ( ωθ, ), (9) d k with the dispersion relation Eq.2. For wave models like WaveWatchIII and SWAN, it may be convenient to interact with the outside main program with ω-space. However, once the ω-space is set in the main program and fixed during the whole computation, the k-space needs to be reinitialized according to Eq.9 every time when the WRT computation takes place at a geographic grid where a new depth d is met, which seems necessary but time consuming. Another disadvantage of frequent reinitialization is that the pre-calculation technique is not well exploited, because all of the pre-calculated parameters are dependent on the core k-space, and once the k-space is changed, those pre-calculated values will no longer match and hence become useless. Thus, we transformed the whole WRT scheme into the k-space including introducing a fixed wave number spectral grid that is independent of the original, transforming all spectral related variables into the new spectral space and simplifying some spectral grid checking processes, which is unnecessary now because the spectral space and related parameters will not be changed any more during the whole computation. However, the variables of the group speed, which depend on local depth, still need to be revalued for each grid. As mentioned above, wave spectra in MASNUM- WAM are also solved in a wave number space. After the transformation, the interaction between the main program and the WRT subroutines are more convenient now. We adopt the same spectral space
4 CHIN. J. OCEANOL. LIMNOL., used in the main program and the initialization is performed only once at the initializing part of the model. The pre-calculated records also can be accessed more efficiently now. More details are given in Section Parallel computation Parallel computing is a useful technique to improve computation speed, and MASNUM-WAM supports the MPI Parallelization (Wang et al., 27). In the MPI version of MASNUM-WAM, each CPU is assigned to a part of the modeling sea area and each CPU process has its own independent task, except data communicated on boundaries. When the original WRT method is parallelizing, each CPU manages its own memory but accesses files on the disk only according to filenames, such as error files, log files, and also pre-computed files. If two or more CPUs access the same file at the same time, the whole computation job stops. To solve this problem, we added CPU ids in all of the filenames stored on the disk allowing each CPU use its own files without interference. Those modifications made it possible to perform large-scale computation for real case modeling. 2.3 Numerical procedures design The original WRT subroutines can be separated into two kinds of functions, those that initialize the configurations and those that execute the main computation. They also provide some acceleration techniques such as geometric scaling, the precalculation, and grid-searching techniques. To produce better computational efficiency, these functions and techniques need to be well reorganized and optimized under the computation structure of the implemented model. Thus, we developed four numerical procedures which can implement the basic functions and use the acceleration techniques in the best possible way. Furthermore, the spectral space transformation mentioned above may give rise to additional efficiencies. Experiments described in Section 3 will verify efficacy and compare procedure efficiencies. The four numerical procedures are denoted NP1, NP2, NP3 and NP4, and outlined in Fig.1a d. NP1 is designed for infinite depth whereas the other three are for finite depth. Each of the procedures has an initial part and a main part, both of which will access the memory or the disks for data reading and writing. The initial parts are similar in that the basic configurations of the WRT scheme are set and the loci-related grids and coefficients are prepared for one or several depths. Because of a fixed spectral space has been adopted, the initial part needs to be executed only once and most of the spectral grid-related configurations need no revaluation. The main parts, however, are very different and complicated, having to be repeated many times for each time loop and topographical grid. NP1 (Fig.1a) is designed for deep water. Thus, only max-depth, which is set to m, is precalculated and related information is stored in memory. In the main part, all needed data for the main computation can be loaded from memory, which is much faster than looking up from disks. For infinite depth, the dispersion relationship ω 2 =gk is adopted; hence, the geometric scaling technique can be triggered. Efficiency improvements can be found in Section 3. For finite depth, procedures are more complicated. The NP2 deals with finite-depth situations by way of scaling the infinite results by a factor. Hence, it has the same initial and main part as in NP1, except the non-linear energy transfer rate finally multiplies an observed magnitude scaling factor, which is parameterized with a simple function (Herterich and Hasselmann, 198; The Wamdi Group, 1988) R( x) = 1+ 1 x exp - x x 7 4, (1) x=.75k d m 1-2 k n( k, θ) dkdθ =.75 d, (11) n( k, θ) dkdθ where d is the local depth. This method seems easy and fast. Furthermore, the deep-water dispersion relation is still adopted and the geometric scaling also can be enabled. Different from scaling the final transfer rate in NP2, another option for finite depth is NP3, shown in Fig.1c, which uses a depth-dependent coupling coefficient (Herterich and Hasselmann, 198). Differences between the two methods have been discussed in detail based on a JONSWAP spectrum (van Vledder, 26), but effects of the two options in an operational wave model are still unclear (to be discussed later in Section 3). Because a variety of depths are considered now, more pre-calculations are performed in NP3. In the initial part, only the maxdepth is pre-calculated. However, in the main part, a check on whether a pre-calculated file for current -2
5 JIANG et al.: WRT method in MASNUM-WAM depth exists is performed first. If the file exists, the record will be loaded from the database; if not, a new calculation will be performed and the new record added into the database for checking next time. Because of the spectral space transform and the uniform spectral grid adopted, all the parameters needed in the main computation are only changed with depth (the current effects are not considered in this paper). Once a recorded depth is found, the data can be loaded directly without any further checks. After the first traversal of all depths, only a loading of an appropriate record needs to be performed in the remaining loops of NP3, but the pre-calculated files are as numerous as the topographical grid points. NP4 offers a less-files solution by adopting a grid-searching technique. The grid searching technique in the original WRT method finds a nearest pre-calculated depth according to the local one, and a scale factor, based on the local and nearest depths, is generated from Eqs.1 and 11 yielding R( xlocal ) s =. (12) R x ( ) nearest Thus, in the initial part of NP4, more pre-setting depths are pre-calculated. However the more depths are pre-set, the more time is taken in finding the target but the computation time for calculating scale factor cannot be avoided. Initial part Initial part Set depth(1) = 9999 Set depth(1) = 9999 Set configurations Allocate memory space for variables Initial values for spectral grids Initial group speed for depth(1) Set configurations Allocate memory space for variables Initial values for spectral grids Initial group speed for depth(1) Pre-calculation of all loci, interpolation coefficients, coupling coefficients and Jacobian terms for depth(1) Pre-calculation of all loci, interpolation coefficients, coupling coefficients and Jacobian terms for depth(1) Pre-calculated data stored in database in memory Loops for every time step and topographical grid point Pre-calculated data stored in database in memory Loops for every time step and topographical grid point Main part Main part Read all needed data for main computations from memory Read all needed data for main computations from memory Main computaitons based on locus Main computaitons based on locus Snl(h)=R(x) Snl(inf) a Non-linear source term and diagonal term b Non-linear source term and diagonal term Fig.1 Flow diagrams of the four numerical procedures To be continued
6 CHIN. J. OCEANOL. LIMNOL., Fig.1 Continued Initial part Initial part Set depth(1) = 9999 Set configurations Allocate memory space for variables Initial values for spectral grids Initial group speed for depth(1) Set each pre-setting depth(i)=xxxx i=1, 2, 3,, n depth Set configurations Allocate memory space for variables Initial values for spectral grids Initial group speed for depth(1) Pre-calculation of all loci, interpolation coefficients, coupling coefficients and Jacobian terms for depth(1) Pre-calculation of all loci, interpolation coefficients, coupling coefficients and Jacobian terms for each pre-setting depth(i) Pre-calculated data stored in database in memory or on disks Loops for every time step and topographical grid point Pre-calculated data stored in database in memory or on disks Loops for every time step and topographical grid point Main part Main part Adding new calculated depth into the database Calculating all loci, interpolation coefficients, coupling coefficients and Jacobian terms for local depth N Check any pre-calculated results available for local depth in database stored Y Read needed data from the database Search for a NEAREST (depth) record from the database and read all needed data from it Scale factor based on relative distance from NEAREST depth N Check if any precalculated results for current depth exists in database stored Y Read needed data from the database Main computaitons based on locus Main computaitons based on locus Snl(h)=s Snl(hi); hi is the nearest pre-set depth c Non-linear source term and diagonal term d Non-linear source term and diagonal term 3 NUMERICAL EXPERIMENT Some numerical experiments are now described that were performed based on the numerical procedures presented above. All the experiments were run in the C sector of the Rubik s Cube supercomputer, which is supported by the Shanghai Supercomputer Center (SSC), China. The CPUs used in this supercomputer are AMD Barcelona 1.9 GHz, each of which contains 4 processors. 3.1 Idealized experiments Wave growth in deep water For deep-water scenarios, NP1 is adopted. The ST1 experiments are duration-limited, which consider a steady wind blowing over an unlimited ocean for time t with a reasonable initial condition. Unlimited ocean means the advection terms in the energy balance equation are switched off and only one single grid point is relevant. The time-evolving spectra are all recorded for further study. The ST2 experiment is fetch-limited meaning that it considers a wave propagating along a 1-D space and driven by a uniform and steady wind. The time evolution of the total spectral energy along the 1-D fetch is also recorded. The source terms used in the ST1 and ST2 experiments are the wind input, the wave breaking dissipation, and the nonlinear energy transfer. The
7 JIANG et al.: WRT method in MASNUM-WAM Dimensionless spectral energy JONSWAP Donelan85 a 5 m/s 1 m/s 2 m/s 3 m/s 4 m/s Dimensionless peak frequency JONSWAP Donelan85 5 m/s 1 m/s 2 m/s 3 m/s 4 m/s b Dimensionless time Dimensionless time Fig.2 Dimensionless spectral energy and peak frequency growth with dimensionless time simulation time lasts for 3 hours to ensure the wind waves are fully developed in each experiment. The spectral space is set to 4 grids in the k-direction and 36 in the θ-direction. The wave number starts at.7 1 rad/m, corresponding to about.4 Hz in deep water, and increases by the geometric factor of The experiments marked with gs are added to check the efficiency improvement of the geometric scaling technique (the spectral space can be always geometric spaced but the geometric scaling technique is an option to be enabled). There are five wind speeds (5, 1, 2, 3, 4 m/s) employed to drive the model to examine the consistency; the dimensionless parameters are used for further comparison. The dimensionless spectrum energy is ε = σ g/u 1, in which σ 2 is the wave height variance, i.e., the zero-order moment of the spectrum m and U 1 is the wind speed at a 1-m height above the sea surface. Similarly, the dimensionless peak frequency is v=f p U 1 /g, dimensionless wind duration ς=t(g/u 1 ), and 2 dimensionless fetch χ = x(g/ U1) with dimensional f p in Hz, t in seconds, and x in meters. The simulation results are going compared with observed wave growth curves. The JONSWAP fetchlimited wind growth relation (Hasselmann et al., 198) is 7 ε = χ = -.33 ν 3.5 χ. (13) Another denoted Donelan85 (Donelan et al., 1985) is ε = χ = ν 1.85χ (14) Both Eqs.13 and 14 are fetch-limited, which need to be transferred to duration-limited. According to Hwang and Wang (24), during the propagation of the wave energy, functions of time and space are related by ς χ 2πν dς = dχ, (15) R where R is equal to.4 for wind sea. The dimensionless duration and fetch coupled through 2πAf χ ς = χ R ( α f χ + 1) ( α f χ + 1) R χ = ς 2πA f χ α f χ + 1, or ( α f χ + ), (16) 1/ 1 and the parameters in Eq.16 have fixed values A fχ =11.86, α fχ = Figure 2 shows that, with the modified WRT method, the single-grid-point MASNUM-WAM can well simulate the growth of wind-driven sea, because under different steady winds (colored curves) the dimensionless growth produced are similar to the observed results (black curves). The fetch-limited experiments are performed in a 1-D spacing model. The waves are assumed to propagate along a 6-degrees (about 667 km) space and computational grid points are set every.2 degrees. Figure 3 shows the fully developed dimensionless spectral energy along the 1-D fetch compared with measured data (Hasselmann et al., 198; Donelan et al., 1985). Colored curves in Figs.2 and 3 are produced when geometric scaling is enabled; Table 1 shows the time taken if geometric scaling is enabled or not. This indicates that geometric scaling can save about 4% 5% of the computation time in these tests.
8 CHIN. J. OCEANOL. LIMNOL., Dimensionless spectral energy m/s 1 m/s 2 m/s 3 m/s 4 m/s Dimensionless fetch Fig.3 Dimensionless wave growth with dimensionless fetch depth Dimensionless spectral energy JONSWAP Donelan85 2 m 4 m 8 m 16 m 32 m 64 m 128 m 256 m 512 m 124 m ST3-NP3 ST3-NP2gs ST3-NP Dimensionless time Fig.4 Dimensionless wave growth in constant depth Wave growth in constant depth For constant finite depth, NP2 and NP3 are adopted with geometric scaling available in NP2. The performance of the shallow-water schemes in NP2 and NP3 are examined below in a real model. Three groups of tests are made. The ST3_NP2 and ST3_NP3 adopt NP2 and NP3 separately, and ST3_ NP2gs employed the geometric scaling technique. Here, a single-grid-point model like ST1 is employed again, with all the same settings. Tests are performed for different depths, from the very shallow depth of 2 m to 1 24 m, which is acceptable as deep water. Figure 4 shows the modeling results, in which we can figure that, effects of the two shallow-water schemes in NP2 and NP3 are indeed different in a real model. However, the differences become small as depth increases to 1 24 m; there is almost no discrepancy. Nevertheless, the continuously growing curves (solid) of ST3_NP3 seem more reasonable. ST3_NP2 and ST3_NP3 take almost the same simulation time, but geometric scaling saves about half of the time without changing the finally effects (dashed and dotted lines in Fig.4) as indicated in Table 2. Without changes in depth in those tests, the Table 1 Time consumed in ST1 and ST2 experiments (unit: min) Wind speed (m/s) ST1 ST1gs ST2 ST2gs Table 2 Time consumed in ST3 (unit: min) Water depth (m) ST3_NP2 ST3_NP2gs ST3_NP m m m main computation only loads data from memory and no more loci grids are created; hence, a similar time is taken in ST3_NP2 and ST3_NP Modeling of real cases Deep water A real scenario of wind waves in deep-water (ST4) is performed with the parallel version of MASNUM- WAM. The modeling area is located in the South China Sea with a buoy placed at E, N; see Fig.5 (the black triangle). The local water depth is more than 2 1 m and local wind and wave are observed simultaneously. We set the modeling area from 99 E to 125 E, 2 N to 25 N with a resolution of.125 degrees, the spectral grid to 35 grid points in the k-direction, beginning with.7 1 and geometrically increasing by 1.21, and to 12 grid points in the θ-direction covering a full circle. The driving wind is from the NCEP reanalysis product, which conforms with wind observations of the buoy (Fig.6). The modeling time is from July 5th to 26th, 211. During this period, a continuous SSW wind was blowing and the water depth was more than 1 m in the SSW direction of the buoy, which can be classified as deep water.
9 JIANG et al.: WRT method in MASNUM-WAM With the DIA method also included, the experiment is denoted ST4-DIA; ST4-NP1 denotes the MPI version of MASNUM-WAM adopting NP1 with geometric scaling enabled. In both computations, we 25 m 75 N E 125 Fig.5 ST4 modeling area for the South China Sea simulation set the same parameters except the non-linear energy transfer source term, and we employed 128 CPUs to run the models. Denoted ST4-DIAs and ST4-NP1s, two additional computations with 4 grids in the k-direction and 36 in the θ-direction as in the idealized experiments are also conducted to compare computation times. From Fig.7, ST4-NP1 performs well and proves the viability of such implementations in the parallel version of the wave model. ST4-NP1 performs even better than the DIA method in the mean wave period according to the observation data. However, the computation time for ST4-NP1 is much larger by about 54 times that of the DIA model (Table 3). When the calculated quantity expands about 3 times (ST4NP1s), the computation time is up to more than 22 days. However, because of the spectral space transform and geometric scaling technique, it is still much smaller than the computation times of times recorded for WaveWatchIII (Tolman, 29). 14 Observation Model 12 Wind speed (m/s) Wind direction (deg) Date Fig.6 Comparison of wind speed and direction of the NCEP reanalysis product with buoy observations
10 CHIN. J. OCEANOL. LIMNOL., 4 3 Observation NP1 DIA H s (m) Tmean (s) Dir (deg) Date Fig.7 Comparison of simulation results from ST4 with observations Table 3 Experiments settings and time consuming in ST4&5 (unit: min) Tests Spectral grids Finite depth Pre-calculated depths Geometric scaling Grid search Computation time ST4-DIA ST4-NP Yes ST4-DIAs ST4-NP1s Yes ST5-DIA ST5-NP Yes ST5-NP No Yes 4 8 ST5-NP * No No *: with water grid points in ST5, the first traversal in the main computation means pre-calculating all the depths at once. The finite-depth real case is denoted as ST5, and the simulation is performed in an area in the Bohai Sea of China. Connected with the Yellow Sea through 4 N the Bohai Strait, the Bohai Sea is a closed sea area with a mean depth of about 2 m. A buoy is situated in the Bohai Strait with a depth of 33 m ( E, N), and the wave field of this area is influenced mostly by winds coming from the East. The ST5 modeling area is set as E, N with a resolution of 1/12 degrees. The modeling time is from November 5th to 15th, 212. The same spectral grid as in ST4 is adopted. The driven wind field is a hindcasting product of the WRF model and a comparison with observed data shows E Fig.8 ST5 modeling area for the Bohai Sea simulation m
11 JIANG et al.: WRT method in MASNUM-WAM Wind speed (m/s) Observation Model Wind direction (deg) Date Fig.9 Comparison of the WRF numerical hindcasting wind with observational data the viability of such driven wind (Fig.9). For each computation, 48 CPUs were employed. With a finite depth, the situation is a more complicated than deep water. There are four tests taken in ST5. ST5-DIA employs the DIA method which with finite depth adopts the same scale function Eq.1 as in NP2. ST5-NP2 and ST5-NP3 adopt NP2 and NP3 separately, and geometric scaling is enabled in ST5-NP2. For ST5-NP4, the pre-calculated depths are set to 2, 4, 8, 16, 32, 64, and 128 meters, as the max-depth of the modeling area is 88 m. Figure 1 shows the simulation results of the four tests. The wave direction data marks spinning with the time periodically in Fig.1 looks very strange. We checked the observational data and contacted the manufacturer of the buoy, and finally confirmed that there was something amiss with the wave direction sensor of the buoy. Despite the irregularity in the mean direction of the observed wave, the other wave parameters observed were still valid. The simulation results of the tests with DIA, NP3 and NP4 were in good agreement, although ST5-NP2 deviated from the other three, especially in the mean wave period simulation. From Table 3, the test with DIA is clearly the most efficient. The computation time for the exact method in these tests were 1 to 5 times more than the DIA method, which is still smaller than the computation times reported for WaveWatchIII (Tolman, 29) Summary Six groups of 6 experiments were performed. The new numerical procedures were adopted according to the water depth for each of the experiments. For infinite depth, NP1 with the transformed wave number spectral space works very well, and about 4% 5% of the computation time can be saved with the geometric scaling technique. For modeling of real cases with infinite depth, NP1 with geometric scaling only takes about 5 times more time than the DIA method, which is much less than the runtimes for WaveWatchIII. For finite depth, NP2 provides the best efficiency in the exact procedures but gives different results compare with NP3 in the constant finite-depth test and the DIA, NP3 and NP4 in the real-case test. In the real case test, NP2 with geometric scaling takes about 17 times more time than the DIA method, and the well-designed NP3 takes about 245 times more time without geometric scaling. NP2 and NP3, adopting different shallow-water schemes separately, do not produce results in agreement with the idealized experiments (Fig.4). A reasonable explanation for the deviation is that, as for the theoretical analysis (van Vledder, 26), the
12 CHIN. J. OCEANOL. LIMNOL., H s (m) Observation DIA NP2 NP4 NP3 Tmean (s) Dir (deg) Date Fig.1 Comparison of simulating results of ST5 with observational data scaled scheme, as for NP2, distributes more energy into higher frequencies; otherwise, the scheme using depth-dependent coupling coefficient, as in NP3 and NP4, distributes energy more uniformly. The total energy is conservative for non-linear transfer, but a growth limiter (Hersbach and Janssen, 1999) should be imposed later to ensure numerical stability but also to flatten the positive lobe of the spectrum in the equilibrium range. The original conservation is then broken, energy is lost and a smaller significant wave height is generated. The deviation of the energy structure in the spectra could also lead to a prediction of larger wave period, as is seen in Figs.7 and 1. Comparing NP4 with the original grid-searching technique, the new NP3 takes less time and yields the same simulation results. Unlike the original WRT method adopting fixed angular-frequency spectral space, the wave number spectral space is often changed for the local environment. Based on the spectral space transform, all spectral space-dependent loci grids, coefficient, and Jacobian terms are fixed for a certain depth because the spectral space is fixed. Checks for such a record are always passed and all of the pre-calculated information can be accepted as long as a correct depth is found. Hence, it is feasible to pre-calculate all depths once and then load a particular one directly, even though there are as many pre-stored files as topographical grids. Otherwise the grid-searching technique needs to search all of the pre-calculated files downward and upward for every grid point which augments computation demands. 4 CONCLUSION AND DISCUSSION In this paper, we described the transformation of the spectral space in the WRT method and developed four numerical procedures to implement it in the MASNUM wave model for both serial and parallel versions. The transformation from the angularfrequency spectral space to the wave number space formed the foundation, on which geometric-scaling and pre-computation techniques were more efficiently executed in the numerical procedures. Experiments were performed to check the viability and efficiency of the implementation, from the results of which we concluded that: 1) The implementation of the WRT method in MASNUM-WAM for both serial and parallel versions is successful. MASNUM-WAM makes gains in speed useful in numerical experimentation for the purposes of research. 2) The developed numerical procedures which are based on spectral space transformations perform with good efficacy. The geometric scaling technique saves up to 4% 5% of the computation times for infinite
13 JIANG et al.: WRT method in MASNUM-WAM depth conditions, and the numerical procedure for infinite depth only takes about 5 times more time than the DIA method in a real modeling scenario. Precalculating all water depth enables gains over the original grid-searching technique for finite depth conditions, and NP3 takes about 245 times more time than the DIA method. Compared with the computational expense incurred by WaveWatchIII, MASNUM-WAM combined with the WRT method performed more efficiently. Furthermore, other improvements and optimizations for the WRT method have still to be made, such as storing all pre-computed parameters in memory instead of on disk particularly if computer memory is sufficiently large and unused. The spectrum solving in MASNUM-WAM is also influenced by other source terms and parameters. Some calibration and correction can produce more accurate results. The operational use of the WRT method is not recommended because it is still much more computationally intensive than the DIA method when similar results are produced. Nevertheless, exact computation for non-linear energy transfer yields more reasonable spectral shapes, which is necessary for numerical research and helpful for improving and perfecting models. 5 ACKNOWLEDGMENT We thank GAO Shanhong, Ocean University of China, for supplying the WRF hindcasting product to drive the models in our experiments. All of the experiments were performed on the Rubik s Cube supercomputer of Shanghai Supercomputer Center. We thank the supercomputer engineers for their considerate support. Finally, the program codes used in this paper are all based on programs published by Gerbrant van Vledder, ALKYON Hydraulic Consultancy & Research. Some modification was performed under the terms of the GNU General Public License. References Ardhuin F, Herbers T H C, O Reilly W C. 21. A hybrid Eulerian-Lagrangian model for spectral wave evolution with application to bottom friction on the continental shelf. J. Phys. Oceanogr., 31(6): Barnett T P On the generation, dissipation, and prediction of ocean wind waves. J. Geophys. Res., 73(2): Booij N, Haagsma I J G, Holthuijsen L H, Kieftenburg A T M M, Ris R C, van der Westhuysen A J, Zijlema M. 24. SWAN User Manual Cycle III Version Delft University of Technology. Booij N, Ris R C, Holthuijsen L H A third-generation wave model for coastal regions: 1. Model description and validation. J. Geophys. Res., 14(C4): Donelan M A, Hamilton J, Hui W H Directional spectra of wind-generated waves. Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences, 315(1534): Ewing J A A numerical wave prediction method for the North Atlantic Ocean. Dtsch. Hydrogr. Z., 24(6): Hasselmann D E, Dunckel M, Ewing J A Directional wave spectra observed during JONSWAP J. Phys. Oceanogr., 1(8): Hasselmann K On the non-linear energy transfer in a gravity-wave spectrum: part 1. General theory. J. Fluid Mech., 12(4): Hasselmann K. 1963a. On the non-linear energy transfer in a gravity-wave spectrum: part 2. Conservation theorems; wave-particle analogy; irreversibility. J. Fluid Mech., 15(2): Hasselmann K. 1963b. On the non-linear energy transfer in a gravity-wave spectrum: part 3. Evaluation of energy flux and swell-sea interaction for a Neumann spectrum. J. Fluid Mech., 15(3): Hasselmann S, Hasselmann K, Allender J H, Barnett T P Computations and parameterizations of the nonlinear energy transfer in a gravity-wave spectrum. Part II: parameterizations of the nonlinear energy transfer for application in wave models. J. Phys. Oceanogr., 15(11): Hasselmann S, Hasselmann K A symmetrical method of computing the non-linear transfer in a gravity-wave spectrum. Hamburger Geophys, Einzelschriften. 52p. Hasselmann S, Hasselmann K. 1985a. Computations and parameterizations of the nonlinear energy transfer in a gravity-wave spectrum. Part I: a new method for efficient computations of the exact nonlinear transfer integral. J. Phys. Oceanogr., 15(11): Hasselmann S, Hasselmann K. 1985c. The wave model EXACT-NL. In: The SWAMP Group ed. Ocean Wave Modelling. Plenum Press, New York. p Hersbach H, Janssen P A E M Improvement of the short-fetch behavior in the Wave Ocean Model (WAM). Journal of Atmospheric and Oceanic Technology, 16(7): Herterich K, Hasselmann K A similarity relation for the nonlinear energy transfer in a finite-depth gravity-wave spectrum. J. Fluid Mech., 97(1): Hwang P A, Wang D W. 24. Field measurements of durationlimited growth of wind-generated ocean surface waves at young stage of development. J. Phys. Oceanogr., 34(1): Masuda A Nonlinear energy transfer between wind waves. J. Phys. Oceanogr., 1(12): Monbaliu J, Hargreaves J C, Carretero J-C, Gerritsen H,
14 CHIN. J. OCEANOL. LIMNOL., Flather R Wave modelling in the PROMISE project. Coast. Eng., 37(3-4): Resio D T, Perrie W A numerical study of nonlinear energy fluxes due to wave-wave interactions. Part 1: methodology and basic results. J. Fluid Mech., 223: The Wamdi Group The WAM model-a third generation ocean wave prediction model. J. Phys. Oceanogr., 18(12): Tolman H L A third-generation model for wind waves on slowly varying, unsteady, and inhomogeneous depths and currents. J. Phys. Oceanogr., 21(16): Tolman H L. 22. User Manual and System Documentation of WAVEWATCH III Version Technical Report 222, NOAA/NWS/NCEP/MMAB. 133p. Tolman H L. 29. User Manual and System Documentation of WAVEWATCH III TM Version Technical Report 276, NOAA/NWS/NCEP/MMAB. 194p. Tracy B A, Resio D T Theory and Calculation of the Nonlinear Energy Transfer Between Sea Waves in Deep Water. WIS Technical Report 11. US Army Engineer Waterways Experiment Station, Vicksburg, Mississippi. USA. 47p. van Vledder G Ph, Herbers T H C, Jensen R E, Resio D T, Tracy B A. 2. Modelling of non-linear quadruplet wave-wave interactions in operational wave models. In: Proceedings of the 27th International Conference on Coastal Engineering. Sydney. p van Vledder G Ph. 26. The WRT method for the computation of non-linear four-wave interactions in discrete spectral wave models. Coast. Eng., 53(2-3): Wang G S, Qiao F L, Yang Y Z. 27. Study on parallel algorithm for MPI-based LAGFD-WAM numerical wave model. Advances in Marine Science, 25(4): (in Chinese with English abstract) Webb D J Non-linear transfers between sea waves. Deep-Sea Res., 25(3): Yang Y Z, Qiao F L, Zhao W, Teng Y, Yuan Y L. 25. MASNUM ocean wave numerical model in spherical coordinates and its application. Acta Oceanologica Sinica, 27(2): 1-7. (in Chinese with English abstract) Yuan Y L, Hua F, Pan Z D, Sun L T. 1992b. The LAGFD- WAM wave model II. Reginal Characteristic inlaid scheme and its application. Acta Oceanologica Sinica, 14(6): (in Chinese) Yuan Y L, Pan Z D, Hua F, Sun L T. 1992a. The LAGFD- WAM wave model I. Basic physical model. Acta Oceanologica Sinica, 14(5): 1-7. (in Chinese) Yuan Y L, Tung C C, Huang N E Statistical characteristics of breaking waves. In: Phillips O M, Hasselmann K eds. Wave Dynamics and Radio Probing of the Ocean Surface. Springer, New York. p
A Multiscale Nested Modeling Framework to Simulate the Interaction of Surface Gravity Waves with Nonlinear Internal Gravity Waves
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. A Multiscale Nested Modeling Framework to Simulate the Interaction of Surface Gravity Waves with Nonlinear Internal Gravity
More informationQuantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Quantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement Lian Shen Department of Mechanical Engineering
More informationGlobal Stokes Drift and Climate Wave Modeling
Global Stokes Drift and Climate Wave Modeling Adrean Webb University of Colorado, Boulder Department of Applied Mathematics February 20, 2012 In Collaboration with: Research funded by: Baylor Fox-Kemper,
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
1 DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. A Direct Simulation-Based Study of Radiance in a Dynamic Ocean LONG-TERM GOALS Dick K.P. Yue Center for Ocean Engineering
More informationNumerical wave modeling Fundamentals, Model setup and Validation
Chapter 2 Numerical wave modeling Fundamentals, Model setup and Validation 2.1. Introduction This chapter describes the background of numerical wave modeling, the state-ofthe-art wave model MIKE 21 SW,
More informationDevelopment of Real Time Wave Simulation Technique
Development of Real Time Wave Simulation Technique Ben T. Nohara, Member Masami Matsuura, Member Summary The research of ocean s has been important, not changing since the Age of Great Voyages, because
More informationThree-dimensional simulation of floating wave power device Xixi Pan 1, a, Shiming Wang 1, b, Yongcheng Liang 1, c
International Power, Electronics and Materials Engineering Conference (IPEMEC 2015) Three-dimensional simulation of floating wave power device Xixi Pan 1, a, Shiming Wang 1, b, Yongcheng Liang 1, c 1 Department
More informationQuantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Quantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement Lian Shen Department of Mechanical Engineering
More informationCoastal Ocean Modeling & Dynamics
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Coastal Ocean Modeling & Dynamics Roger M. Samelson College of Oceanic and Atmospheric Sciences Oregon State University
More informationDevelopment of the Compliant Mooring Line Model for FLOW-3D
Flow Science Report 08-15 Development of the Compliant Mooring Line Model for FLOW-3D Gengsheng Wei Flow Science, Inc. October 2015 1. Introduction Mooring systems are common in offshore structures, ship
More informationA Diagonal Split-cell Model for the High-order Symplectic FDTD Scheme
PIERS ONLINE, VOL. 2, NO. 6, 2006 715 A Diagonal Split-cell Model for the High-order Symplectic FDTD Scheme Wei Sha, Xianliang Wu, and Mingsheng Chen Key Laboratory of Intelligent Computing & Signal Processing
More information2. Hindcasting of exceptional events having severely damaged maritime structures and for which field records are either incomplete or unavailable.
1/63 TOMAWAC User Manual links from User Manuals, TOMAWAC 1.1 Generals TOMAWAC is a scientific software which models the changes, both in the time and in the spatial domain, of the power spectrum of wind-driven
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
An explicit and conservative remapping strategy for semi-lagrangian advection Sebastian Reich Universität Potsdam, Potsdam, Germany January 17, 2007 Abstract A conservative semi-lagrangian advection scheme
More informationSWAN USER MANUAL. SWAN Cycle III version 41.20
SWAN USER MANUAL SWAN Cycle III version 41.20 SWAN USER MANUAL by : The SWAN team mail address : Delft University of Technology Faculty of Civil Engineering and Geosciences Environmental Fluid Mechanics
More informationSwanOne User Manual V1.3. April 2018
SwanOne User Manual V1.3 April 2018 Version 1.3 (April 2018) 1 Table of Contents 1. Introduction... 3 2. Installation... 3 3. Starting with SwanOne... 4 3.1. File New... 4 4. Input data... 5 4.1. Input
More informationImplicit schemes for wave models
Implicit schemes for wave models Mathieu Dutour Sikirić Rudjer Bo sković Institute, Croatia and Universität Rostock April 17, 2013 I. Wave models Stochastic wave modelling Oceanic models are using grids
More informationQuantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurements
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Quantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurements Dick K.P. Yue Center for Ocean Engineering
More informationModule 1 Lecture Notes 2. Optimization Problem and Model Formulation
Optimization Methods: Introduction and Basic concepts 1 Module 1 Lecture Notes 2 Optimization Problem and Model Formulation Introduction In the previous lecture we studied the evolution of optimization
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
A Direct Simulation-Based Study of Radiance in a Dynamic Ocean Dick K.P. Yue Center for Ocean Engineering Massachusetts Institute of Technology Room 5-321, 77 Massachusetts Ave, Cambridge, MA 02139 phone:
More informationTravel-Time Sensitivity Kernels in Long-Range Propagation
Travel-Time Sensitivity Kernels in Long-Range Propagation Emmanuel Skarsoulis Foundation for Research and Technology Hellas Institute of Applied and Computational Mathematics P.O. Box 1385, GR-71110 Heraklion,
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationThe improvement of MOM4 by adding wave-induced mixing
The improvement of MOM4 by adding wave-induced mixing Qi Shu, Fangli Qiao, Zhenya Song, Changshui Xia and Yongzheng Yang Laboratory of MArine Science and NUmerical Modeling (MASNUM) The First Institute
More informationComputer Project 3. AA Computational Fluid Dyanmics University of Washington. Mishaal Aleem March 17, 2015
Computer Project 3 AA 543 - Computational Fluid Dyanmics University of Washington Mishaal Aleem March 17, 2015 Contents Introduction........................................... 1 3.1 Grid Generator.......................................
More informationGeometric Acoustics in High-Speed Boundary Layers
Accepted for presentation at the 9th International Symposium on Shock Waves. Madison, WI. July -9,. Paper #8 Geometric Acoustics in High-Speed Boundary Layers N. J. Parziale, J. E. Shepherd, and H. G.
More informationWave-ice interactions in wave models like WW3
Wave-ice interactions in wave models like WW3 Will Perrie 1, Hui Shen 1,2, Mike Meylan 3, Bash Toulany 1 1Bedford Institute of Oceanography, Dartmouth NS 2Inst. of Oceanology, China Academy of Sciences,
More information2D numerical simulation of ocean waves
2D numerical simulation of ocean waves Qingjie. Du,*, Y.C. Dennis. Leung Department of Mechanical Engineering, The University of Hong Kong, Hong Kong, China * Corresponding author. Tel: +852 51743593,
More informationCloud Cavitating Flow around an Axisymmetric Projectile in the shallow water
Cloud Cavitating Flow around an Axisymmetric Projectile in the shallow water 1,2 Chang Xu; 1,2 Yiwei Wang*; 1,2 Jian Huang; 1,2 Chenguang Huang 1 Key Laboratory for Mechanics in Fluid Solid Coupling Systems,
More informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More informationAbstract. 1 Introduction
Diffraction of wave spectra by two breakwaters at an angle T. Zographou & C.D. Memos Department of Civil Engineering, National Technical University of Athens, 157 SOZografos, Greece EM ail: memos@hydro.civil.ntua.gr
More informationAutomated Segmentation Using a Fast Implementation of the Chan-Vese Models
Automated Segmentation Using a Fast Implementation of the Chan-Vese Models Huan Xu, and Xiao-Feng Wang,,3 Intelligent Computation Lab, Hefei Institute of Intelligent Machines, Chinese Academy of Science,
More informationLevel-set and ALE Based Topology Optimization Using Nonlinear Programming
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level-set and ALE Based Topology Optimization Using Nonlinear Programming Shintaro Yamasaki
More informationLagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion
Lagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion James F. Price Woods Hole Oceanographic Institution Woods Hole, MA, 02543 July 31, 2006 Summary: This essay
More informationSimulation of Fuel Sloshing Comparative Study
3. LS-DYNA Anwenderforum, Bamberg 2004 Netfreie Verfahren Simulation of Fuel Sloshing Comparative Study Matej Vesenjak 1, Heiner Müllerschön 2, Alexander Hummel 3, Zoran Ren 1 1 University of Maribor,
More informationAdvanced Numerical Methods for Numerical Weather Prediction
Advanced Numerical Methods for Numerical Weather Prediction Francis X. Giraldo Naval Research Laboratory Monterey, CA 93943-5502 phone: (831) 656-4882 fax: (831) 656-4769 e-mail: giraldo@nrlmry.navy.mil
More informationLagrangian methods and Smoothed Particle Hydrodynamics (SPH) Computation in Astrophysics Seminar (Spring 2006) L. J. Dursi
Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Eulerian Grid Methods The methods covered so far in this course use an Eulerian grid: Prescribed coordinates In `lab frame' Fluid elements flow
More informationNumerical Wave Tank Modeling of Hydrodynamics of Permeable Barriers
ICHE 2014, Hamburg - Lehfeldt & Kopmann (eds) - 2014 Bundesanstalt für Wasserbau ISBN 978-3-939230-32-8 Numerical Wave Tank Modeling of Hydrodynamics of Permeable Barriers K. Rajendra & R. Balaji Indian
More informationComputational scheme and effect of complex multi-islands terrain in numerical calculation of waves in Archipelago
Indian Journal of Geo-Marine Sciences Vol. 42(7), November 203, pp. 825-832 Computational scheme and effect of complex multi-islands terrain in numerical calculation of waves in Archipelago Ke-feng Mao,2,Xue-feng
More informationAn explicit feature control approach in structural topology optimization
th World Congress on Structural and Multidisciplinary Optimisation 07 th -2 th, June 205, Sydney Australia An explicit feature control approach in structural topology optimization Weisheng Zhang, Xu Guo
More informationCOMPARISON OF FULL-SCALE MEASUREMENTS WITH CALCULATED MOTION CHARACTERISTICS OF A WEST OF AFRICA FPSO
Proceedings of OMAE3 ND International Conference on Offshore Mechanics and Arctic Engineering June 8 3, 3, Cancun, Mexico OMAE3-378 COMPARISON OF FULL-SCALE MEASUREMENTS WITH CALCULATED MOTION CHARACTERISTICS
More informationAdvances in Simulation for Marine And Offshore Applications. Milovan Perić
Advances in Simulation for Marine And Offshore Applications Milovan Perić Introduction Extensions and enhancements in STAR-CCM+ for marine and offshore applications: Creation of irregular long-crested
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
A Direct Simulation-Based Study of Radiance in a Dynamic Ocean Lian Shen Department of Civil Engineering Johns Hopkins University Baltimore, MD 21218 phone: (410) 516-5033 fax: (410) 516-7473 email: LianShen@jhu.edu
More informationCOMPUTATIONALLY EFFICIENT RAY TRACING ALGORITHM FOR SIMULATION OF TRANSDUCER FIELDS IN ANISOTROPIC MATERIALS
Proceedings of the National Seminar & Exhibition on Non-Destructive Evaluation NDE 2011, December 8-10, 2011 COMPUTATIONALLY EFFICIENT RAY TRACING ALGORITHM FOR SIMULATION OF TRANSDUCER FIELDS IN ANISOTROPIC
More informationRiver inundation modelling for risk analysis
River inundation modelling for risk analysis L. H. C. Chua, F. Merting & K. P. Holz Institute for Bauinformatik, Brandenburg Technical University, Germany Abstract This paper presents the results of an
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 8: 58 63 (2007) Published online 22 May 2007 in Wiley InterScience (www.interscience.wiley.com).151 An explicit and conservative remapping strategy for semi-lagrangian
More informationHigh-Resolution Ocean Wave Estimation
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. High-Resolution Ocean Wave Estimation David Walker SRI International 2100 Commonwealth Boulevard, Third Floor Ann Arbor,
More informationResearch on time optimal trajectory planning of 7-DOF manipulator based on genetic algorithm
Acta Technica 61, No. 4A/2016, 189 200 c 2017 Institute of Thermomechanics CAS, v.v.i. Research on time optimal trajectory planning of 7-DOF manipulator based on genetic algorithm Jianrong Bu 1, Junyan
More informationContinued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data
Continued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data Dick K.P. Yue Center for Ocean Engineering Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge,
More information17. SEISMIC ANALYSIS MODELING TO SATISFY BUILDING CODES
17. SEISMIC ANALYSIS MODELING TO SATISFY BUILDING CODES The Current Building Codes Use the Terminology: Principal Direction without a Unique Definition 17.1 INTRODUCTION { XE "Building Codes" }Currently
More informationGeometric Considerations for Distribution of Sensors in Ad-hoc Sensor Networks
Geometric Considerations for Distribution of Sensors in Ad-hoc Sensor Networks Ted Brown, Deniz Sarioz, Amotz Bar-Noy, Tom LaPorta, Dinesh Verma, Matthew Johnson, Hosam Rowaihy November 20, 2006 1 Introduction
More informationA new 8-node quadrilateral spline finite element
Journal of Computational and Applied Mathematics 195 (2006) 54 65 www.elsevier.com/locate/cam A new 8-node quadrilateral spline finite element Chong-Jun Li, Ren-Hong Wang Institute of Mathematical Sciences,
More informationWAVE PATTERNS, WAVE INDUCED FORCES AND MOMENTS FOR A GRAVITY BASED STRUCTURE PREDICTED USING CFD
Proceedings of the ASME 2011 30th International Conference on Ocean, Offshore and Arctic Engineering OMAE2011 June 19-24, 2011, Rotterdam, The Netherlands OMAE2011-49593 WAVE PATTERNS, WAVE INDUCED FORCES
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More information4.1 QUANTIZATION NOISE
DIGITAL SIGNAL PROCESSING UNIT IV FINITE WORD LENGTH EFFECTS Contents : 4.1 Quantization Noise 4.2 Fixed Point and Floating Point Number Representation 4.3 Truncation and Rounding 4.4 Quantization Noise
More informationCFD-RANS APPLICATIONS IN COMPLEX TERRAIN ANALYSIS. NUMERICAL VS EXPERIMENTAL RESULTS A CASE STUDY: COZZOVALLEFONDI WIND FARM IN SICILY
CFD-RANS APPLICATIONS IN COMPLEX TERRAIN ANALYSIS. NUMERICAL VS EXPERIMENTAL RESULTS A CASE STUDY: COZZOVALLEFONDI WIND FARM IN SICILY J.Maza (*),G.Nicoletti(**), (*) Pisa University, Aerospace Engineering
More informationNumerical Analysis of Shock Tube Problem by using TVD and ACM Schemes
Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics
More informationPossibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics
Possibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics Masanori Hashiguchi 1 1 Keisoku Engineering System Co., Ltd. 1-9-5 Uchikanda, Chiyoda-ku,
More informationGuidelines for proper use of Plate elements
Guidelines for proper use of Plate elements In structural analysis using finite element method, the analysis model is created by dividing the entire structure into finite elements. This procedure is known
More informationApplication and Research of Structure Topology Optimization of Scraper Conveyer with MSC.Nastran.Optishape
COMPUTATIONAL METHODS IN ENGINEERING AND SCIENCE EPMESC X, Aug. 21-23, 2006, Sanya, Hainan, China 2006 Tsinghua University Press & Springer Application and Research of Structure Topology Optimization of
More informationINTERACTION BETWEEN TURBULENT DYNAMICAL PROCESSES AND STATISTICS IN DEFORMED AIR-LIQUID INTERFACES, VIA DNS
INTERACTION BETWEEN TURBULENT DYNAMICAL PROCESSES AND STATISTICS IN DEFORMED AIR-LIQUID INTERFACES, VIA DNS Yoshinobu Yamamoto Department of Nuclear Engineering, Kyoto University Yoshida Sakyo Kyoto, 66-85,
More informationInclusion of Aleatory and Epistemic Uncertainty in Design Optimization
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Inclusion of Aleatory and Epistemic Uncertainty in Design Optimization Sirisha Rangavajhala
More informationPATCH TEST OF HEXAHEDRAL ELEMENT
Annual Report of ADVENTURE Project ADV-99- (999) PATCH TEST OF HEXAHEDRAL ELEMENT Yoshikazu ISHIHARA * and Hirohisa NOGUCHI * * Mitsubishi Research Institute, Inc. e-mail: y-ishi@mri.co.jp * Department
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationThe Spherical Harmonics Discrete Ordinate Method for Atmospheric Radiative Transfer
The Spherical Harmonics Discrete Ordinate Method for Atmospheric Radiative Transfer K. Franklin Evans Program in Atmospheric and Oceanic Sciences University of Colorado, Boulder Computational Methods in
More informationQuantitative flow visualization using the hydraulic analogy
Experiments in Fluids 27 (2000) 165 169 Springer-Verlag 2000 Quantitative flow visualization using the hydraulic analogy S. L. Rani, M. S. Wooldridge 165 Abstract The current work describes the development
More informationTable of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Contents.
Table of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Preface page 2 Part I Fluid Dynamics and Waves 7 1 Elements of fluid dynamics 9 1.1
More informationTurbulent Structure underneath Air-Sea Wavy Interface: Large-Eddy Simulation Mostafa Bakhoday Paskyabi
Turbulent Structure underneath Air-Sea Wavy Interface: Large-Eddy Simulation Mostafa Bakhoday Paskyabi Geophysical Institute, University of Bergen (Mostafa.Bakhoday@uib.no) Introduction Waves Planetary
More informationFull waveform inversion of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC)
of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC) Summary (FWI) is a promising technology for next generation model building. However, it faces with many
More informationSPH: Towards the simulation of wave-body interactions in extreme seas
SPH: Towards the simulation of wave-body interactions in extreme seas Guillaume Oger, Mathieu Doring, Bertrand Alessandrini, and Pierre Ferrant Fluid Mechanics Laboratory (CNRS UMR6598) Ecole Centrale
More informationInteraction of Fluid Simulation Based on PhysX Physics Engine. Huibai Wang, Jianfei Wan, Fengquan Zhang
4th International Conference on Sensors, Measurement and Intelligent Materials (ICSMIM 2015) Interaction of Fluid Simulation Based on PhysX Physics Engine Huibai Wang, Jianfei Wan, Fengquan Zhang College
More informationWETTING PROPERTIES OF STRUCTURED INTERFACES COMPOSED OF SURFACE-ATTACHED SPHERICAL NANOPARTICLES
November 20, 2018 WETTING PROPERTIES OF STRUCTURED INTERFACES COMPOSED OF SURFACE-ATTACHED SPHERICAL NANOPARTICLES Bishal Bhattarai and Nikolai V. Priezjev Department of Mechanical and Materials Engineering
More informationOptimal Segmentation and Understanding of Motion Capture Data
Optimal Segmentation and Understanding of Motion Capture Data Xiang Huang, M.A.Sc Candidate Department of Electrical and Computer Engineering McMaster University Supervisor: Dr. Xiaolin Wu 7 Apr, 2005
More informationQUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD
QUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD Keisuke Yoshida, Tadaharu Ishikawa Dr. Eng., Tokyo Institute
More informationNumerical Simulation Study on Aerodynamic Characteristics of the High Speed Train under Crosswind
2017 2nd International Conference on Industrial Aerodynamics (ICIA 2017) ISBN: 978-1-60595-481-3 Numerical Simulation Study on Aerodynamic Characteristics of the High Speed Train under Crosswind Fan Zhao,
More informationA mass-conservative version of the semi- Lagrangian semi-implicit HIRLAM using Lagrangian vertical coordinates
A mass-conservative version of the semi- Lagrangian semi-implicit HIRLAM using Lagrangian vertical coordinates Peter Hjort Lauritzen Atmospheric Modeling & Predictability Section National Center for Atmospheric
More informationSimultaneous Perturbation Stochastic Approximation Algorithm Combined with Neural Network and Fuzzy Simulation
.--- Simultaneous Perturbation Stochastic Approximation Algorithm Combined with Neural Networ and Fuzzy Simulation Abstract - - - - Keywords: Many optimization problems contain fuzzy information. Possibility
More informationMODELING OF THREE-DIMENSIONAL PROPAGATION ON A COASTAL WEDGE WITH A SEDIMENT SUPPORTING SHEAR
MODELING OF THREE-DIMENSIONAL PROPAGATION ON A COASTAL WEDGE WITH A SEDIMENT SUPPORTING SHEAR Piotr Borejko Vienna University of Technology, Karlsplatz 13/E26/3, 14 Vienna, Austria Fax: + 43 1 588 1 21
More informationCGT 581 G Fluids. Overview. Some terms. Some terms
CGT 581 G Fluids Bedřich Beneš, Ph.D. Purdue University Department of Computer Graphics Technology Overview Some terms Incompressible Navier-Stokes Boundary conditions Lagrange vs. Euler Eulerian approaches
More informationExistence and Dynamics of Bounded Traveling Wave Solutions to Getmanou Equation
Commun. Theor. Phys. 70 (2018) 672 676 Vol. 70, No. 6, December 1, 2018 Existence and Dynamics of Bounded Traveling Wave Solutions to Getmanou Equation Zhen-Shu Wen ( 温振庶 ) School of Mathematical Sciences,
More informationITTC Recommended Procedures and Guidelines. ITTC Quality System Manual. Recommended Procedures and Guidelines
Offshore Structures with Mooring Lines 7.5 Page 1 of 6 ITTC Quality System Manual Recommended Procedure Offshore Structures with Mooring Lines 7.5 Process Control 7.5- Testing and Extrapolation Methods
More informationApplication of a FEA Model for Conformability Calculation of Tip Seal in Compressor
Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 2008 Application of a FEA Model for Conformability Calculation of Tip Seal in Compressor
More informationTruss structural configuration optimization using the linear extended interior penalty function method
ANZIAM J. 46 (E) pp.c1311 C1326, 2006 C1311 Truss structural configuration optimization using the linear extended interior penalty function method Wahyu Kuntjoro Jamaluddin Mahmud (Received 25 October
More informationALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving Objects
Tenth International Conference on Computational Fluid Dynamics (ICCFD10), Barcelona,Spain, July 9-13, 2018 ICCFD10-047 ALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 19, 2013 http://acousticalsociety.org/ ICA 2013 Montreal Montreal, Canada 2-7 June 2013 Underwater Acoustics Session 2aUW: Wave Propagation in a Random Medium
More informationMEASURING SURFACE CURRENTS USING IR CAMERAS. Background. Optical Current Meter 06/10/2010. J.Paul Rinehimer ESS522
J.Paul Rinehimer ESS5 Optical Current Meter 6/1/1 MEASURING SURFACE CURRENTS USING IR CAMERAS Background Most in-situ current measurement techniques are based on sending acoustic pulses and measuring the
More informationSIGNIFICANT WAVE HEIGHT RETRIEVAL FROM SYNTHETIC RADAR IMAGES
Proceedings of the 11 th International Conference on Hydrodynamics (ICHD 2014) October 19 24, 2014, Singapore SIGNIFICANT WAVE HEIGHT RETRIEVAL FROM SYNTHETIC RADAR IMAGES A. P. WIJAYA & E. VAN GROESEN
More informationAurélien Thinat Stéphane Cordier 1, François Cany
SimHydro 2012:New trends in simulation - Hydroinformatics and 3D modeling, 12-14 September 2012, Nice Aurélien Thinat, Stéphane Cordier, François Cany Application of OpenFOAM to the study of wave loads
More informationModeling the Acoustic Scattering from Axially Symmetric Fluid, Elastic, and Poroelastic Objects due to Nonsymmetric Forcing Using COMSOL Multiphysics
Modeling the Acoustic Scattering from Axially Symmetric Fluid, Elastic, and Poroelastic Objects due to Nonsymmetric Forcing Using COMSOL Multiphysics Anthony L. Bonomo *1 and Marcia J. Isakson 1 1 Applied
More informationAdaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited
Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Summary We present a new method for performing full-waveform inversion that appears
More informationMinimal Test Cost Feature Selection with Positive Region Constraint
Minimal Test Cost Feature Selection with Positive Region Constraint Jiabin Liu 1,2,FanMin 2,, Shujiao Liao 2, and William Zhu 2 1 Department of Computer Science, Sichuan University for Nationalities, Kangding
More informationNUMERICAL ANALYSIS OF WIND EFFECT ON HIGH-DENSITY BUILDING AERAS
NUMERICAL ANALYSIS OF WIND EFFECT ON HIGH-DENSITY BUILDING AERAS Bin ZHAO, Ying LI, Xianting LI and Qisen YAN Department of Thermal Engineering, Tsinghua University Beijing, 100084, P.R. China ABSTRACT
More informationCELLULAR AUTOMATA IN MATHEMATICAL MODELING JOSH KANTOR. 1. History
CELLULAR AUTOMATA IN MATHEMATICAL MODELING JOSH KANTOR 1. History Cellular automata were initially conceived of in 1948 by John von Neumann who was searching for ways of modeling evolution. He was trying
More informationAvailable online at ScienceDirect. Energy Procedia 69 (2015 )
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 69 (2015 ) 1885 1894 International Conference on Concentrating Solar Power and Chemical Energy Systems, SolarPACES 2014 Heliostat
More informationLATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS
14 th European Conference on Mixing Warszawa, 10-13 September 2012 LATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS Felix Muggli a, Laurent Chatagny a, Jonas Lätt b a Sulzer Markets & Technology
More informationInversion of Array Induction Logs and Its Application
2007 Petroleum Science Vol.4 No.3 Inversion of Array Induction Logs and Its Application Gao Jie 1, Zhao Aibin 2, Peng Fei 1 and Li Hongqi 1 (1. School of Resources and Information, China University of
More informationThis is an author-deposited version published in: Eprints ID: 4362
This is an author-deposited version published in: http://oatao.univ-toulouse.fr/ Eprints ID: 4362 To cite this document: CHIKHAOUI Oussama, GRESSIER Jérémie, GRONDIN Gilles. Assessment of the Spectral
More informationStable simulations of illumination patterns caused by focusing of sunlight by water waves
Stable simulations of illumination patterns caused by focusing of sunlight by water waves Sjoerd de Ridder ABSTRACT Illumination patterns of underwater sunlight have fascinated various researchers in the
More informationINTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012
INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012 Copyright 2010 All rights reserved Integrated Publishing services Research article ISSN 0976 4399 Efficiency and performances
More informationFull Plane STWAVE: SMS Graphical Interface
Full Plane STWAVE: SMS Graphical Interface by Jane McKee Smith and Alan Zundel PURPOSE: The purpose of this Coastal and Hydraulics Engineering Technical Note (CHETN) is to describe the graphical interface
More informationLearning the Three Factors of a Non-overlapping Multi-camera Network Topology
Learning the Three Factors of a Non-overlapping Multi-camera Network Topology Xiaotang Chen, Kaiqi Huang, and Tieniu Tan National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy
More informationElement Quality Metrics for Higher-Order Bernstein Bézier Elements
Element Quality Metrics for Higher-Order Bernstein Bézier Elements Luke Engvall and John A. Evans Abstract In this note, we review the interpolation theory for curvilinear finite elements originally derived
More information