A new method to predict the ultimate strength of fiber reinforced composites under arbitrary load condition is introduced. The micromechanics strength theory is used to perform the final failure prediction of composit...A new method to predict the ultimate strength of fiber reinforced composites under arbitrary load condition is introduced. The micromechanics strength theory is used to perform the final failure prediction of composite laminates. The theory is based on unit cell analytic model which can provide the ply composite material properties by only using the constituent fiber and matrix properties and the laminate geometric parameters without knowing any experimental information of the laminates. To show that this method is suitable for predicting the strength of composite laminates, the micromechanics strength theory is ranked by comparing it with all the micro-level and the best two macro-level theories chosen from the World Wide Failure Exercise. The results show that this method can be used for predicting strength of any composite laminates and provide a direct reference for composite optimum design.展开更多
The aim of the paper is to discover the general creep mechanisms for the short fiber reinforcement matrix composites (MMCs) under uniaxial stress states and to build a relationship between the macroscopic steady creep...The aim of the paper is to discover the general creep mechanisms for the short fiber reinforcement matrix composites (MMCs) under uniaxial stress states and to build a relationship between the macroscopic steady creep behavior and the material micro geometric parameters. The unit cell models were used to calculate the macroscopic creep behavior with different micro geometric parameters of fibers on different loading directions. The influence of the geometric parameters of the fibers and loading directions on the macroscopic creep behavior had been obtained, and described quantitatively. The matrix/fiber interface had been considered by a third layer, matrix/fiber interlayer, in the unit cells with different creep properties and thickness. Based on the numerical results of the unit cell models, a statistic model had been presented for the plane randomly-distributed-fiber MMCs. The fiber breakage had been taken into account in the statistic model for it starts experimentally early in the creep life. With the distribution of the geometric parameters of the fibers, the results of the statistic model agree well with the experiments. With the statistic model, the influence of the geometric parameters and the breakage of the fibers as well as the properties and thickness of, the interlayer on the macroscopic steady creep rate have been discussed.展开更多
As one of Chinese minority languages,Tibetan speech recognition technology was not researched upon as extensively as Chinese and English were until recently.This,along with the relatively small Tibetan corpus,has resu...As one of Chinese minority languages,Tibetan speech recognition technology was not researched upon as extensively as Chinese and English were until recently.This,along with the relatively small Tibetan corpus,has resulted in an unsatisfying performance of Tibetan speech recognition based on an end-to-end model.This paper aims to achieve an accurate Tibetan speech recognition using a small amount of Tibetan training data.We demonstrate effective methods of Tibetan end-to-end speech recognition via cross-language transfer learning from three aspects:modeling unit selection,transfer learning method,and source language selection.Experimental results show that the Chinese-Tibetan multi-language learning method using multilanguage character set as the modeling unit yields the best performance on Tibetan Character Error Rate(CER)at 27.3%,which is reduced by 26.1%compared to the language-specific model.And our method also achieves the 2.2%higher accuracy using less amount of data compared with the method using Tibetan multi-dialect transfer learning under the same model structure and data set.展开更多
The damage effects of water sorption on the mechanical properties of the hydroxyapatite particle reinforced Bis-GMA/TEGDMA copolymer (HA/Bis-GMA/TEGDMA) h-ave been predicted using 3D finite cell models. The plastici...The damage effects of water sorption on the mechanical properties of the hydroxyapatite particle reinforced Bis-GMA/TEGDMA copolymer (HA/Bis-GMA/TEGDMA) h-ave been predicted using 3D finite cell models. The plasticizer effect on the polymer matrix was considered as a variation of its Young's modulus. Three different cell models were used to determine the influence of varying particle contents, interphase strength and moisture concentration on the debonding damage. The stress distribution pattern has been examined and the stress transfer mode clarified. The Young's modulus and fracture strength of the Bis-GMA/TEGDMA composite were also predicted using the model with and without consideration of the damage. ine Iormer results with consideration of the debonding damage are in good agreement with existing literature experimental data. The shielding effect of our proposed model and an alternative approach were discussed. The FCC cell model has also been extended to predict the critical load for the damaged and the undamaged composite subject to the 3-point flexural test.展开更多
The electric energy injection from a pulsed power supply to a planar type of dielectric barrier discharge(DBD) reactor at atmospheric pressure was studied. Relations of the energy injection with barrier materials, b...The electric energy injection from a pulsed power supply to a planar type of dielectric barrier discharge(DBD) reactor at atmospheric pressure was studied. Relations of the energy injection with barrier materials, barrier thickness, peak voltage, gap distance, electrode area,and operation temperature were experimentally investigated. The energy injection is a function of relative permittivity, barrier thickness, peak voltage, gap distance, and electrode area. The influence of operation temperature on energy injection is slight in the range of 27-300℃ but becomes obvious in the range of 300-500℃. A model was established using which the energy injection can be easily predicted.展开更多
The growth habit of the basic magnesium oxysulfate whisker was investigated based on the theoreticalmodelof anion coordination polyhedron growth units.It is found that typicalbasic magnesium oxysulfate whisker growth ...The growth habit of the basic magnesium oxysulfate whisker was investigated based on the theoreticalmodelof anion coordination polyhedron growth units.It is found that typicalbasic magnesium oxysulfate whisker growth is consistent with anion tetrahedralcoordination incorporation rules.The growth units of basic magnesium oxysulfate whiskers are [Mg-(OH)_4]^(2-) and HSO_4^-.[Mg-(OH)_4]^(2-) is the favorable growth unit and whisker growth is in the direction of the [Mg-(OH)_4]^(2-) combination.A plurality of [Mg-(OH)_4]^(2-) s combine and become a larger dimensionalgrowth unit in a one-dimensionaldirection.Then HSO_4^- and larger dimensionalgrowth units connect as basic magnesium sulfate whiskers,according to the structuralcharacteristics of the basic magnesium sulfate whisker,which can guide the synthesis of magnesium hydroxide whisker.展开更多
According to the existing method including testing the frequency and establishing the relationship between moisture content and frequency, a corresponding instrument was designed. In order to further improve the accur...According to the existing method including testing the frequency and establishing the relationship between moisture content and frequency, a corresponding instrument was designed. In order to further improve the accuracy and rapidity of the system, a new approach to describe the relationship between the measurement error and the temperature was proposed. The error band could be obtained and divided into several parts(based on the range of temperature) to indicate the error value that should compensate the grain moisture content for the changes in temperature. By calculating the error band at the maximum and the minimum operating temperatures, as well as by determining the error compensation value from the error band based on the measurement moisture content, the final effective result was derived.展开更多
To design the amorphous hydrogen storage alloy efficiently, the maximum hydrogen capacities for Zr - Ni amorphous alloy were calculated. Based on the Rhomb Unit Structure Model(RUSM) for amorphous alloy and the experi...To design the amorphous hydrogen storage alloy efficiently, the maximum hydrogen capacities for Zr - Ni amorphous alloy were calculated. Based on the Rhomb Unit Structure Model(RUSM) for amorphous alloy and the experimental result that hydrogen atoms exist in 3Zr1Ni and 4Zr tetrahedron interstices in Zr-Ni amolphous alloy, the numbers of 3Zr-1Ni and 4Zr tetrahedron interstices in a RUSM were calculated which correspond to the hydrogen capacity. The two extremum Zr distribution states were calculated, such as highly heterogeneous Zr distribution and homogeneous Zr distribution. The calculated curves of hydrogen capacity with different Zr contents at two states indicate that the hydrogen capacity increases with increasing Zr content and reaches its maximum when Zr is 75%. The theoretical maximum hydrogen capacity for Zr - Ni amorphous alloy is 2.0 (H/M). Meanwhile, the hydrogen capacity of heterogenous Zr distribution alloy is higher than that of homogenous one at the same Zr content. The experimental results prove the calculated results reasonable, and accordingly, the experimental results that the distribution of Zr atom in amorphous alloy occur heterogeneous after a few hydrogen absorption desorption cycles can be explained.展开更多
In this paper the influence of a knot on the structure of a polymethylene (PM) strand in the tensile process is investigated by using the steered molecular dynamics (SMD) method. The gradual increasing of end-to-e...In this paper the influence of a knot on the structure of a polymethylene (PM) strand in the tensile process is investigated by using the steered molecular dynamics (SMD) method. The gradual increasing of end-to-end distance, R, results in a tighter knot and a more stretched contour. That the break in a knotted rope almost invariably occurs at a point just outside the 'entrance' to the knot, which has been shown in a good many experiments, is further theoretically verified in this paper through the calculation of some structural and thermodynamic parameters. Moreover, it is found that the analyses on bond length, torsion angle and strain energy can facilitate to the study of the localization and the size of a knot in the tensile process. The symmetries of torsion angles, bond lengths and bond angles in the knot result in the whole symmetry of the knot in microstructure, thereby adapting itself to the strain applied. Additionally, the statistical property of the force-dependent average knot size illuminates in detail the change in size of a knot with force f, and therefore the minimum size of the knot in the restriction of the potentials considered in this work for a PM chain is deduced. At the same time, the difference in response to uniaxial strain, between a knotted PM strand and an unknotted one is also investigated. The force-extension profile is easily obtained from the simulation. As expected, for a given f, the knotted chain has an R significantly smaller than that of an unknotted polymer. However, the scaled difference becomes less pronounced for larger values of N, and the results for longer chains approach those of the unknotted chains.展开更多
From the perspective of growth units, the growth mechanism of Mg2(OH)2CO3.3H2O whisker is investigated in this paper. Results show that the growth morphology of Mg2(OH)2CO33H2O whisker is consistent with the model...From the perspective of growth units, the growth mechanism of Mg2(OH)2CO3.3H2O whisker is investigated in this paper. Results show that the growth morphology of Mg2(OH)2CO33H2O whisker is consistent with the model of anion coordination polyhedron growth units. The growth solution Raman shift of Mg2(OH)2CO;3H2O was monitored using Raman spectroscopy. The growth units are [Mg-(OH)4]2- and H2COv The growth process of Mg2(OH)2COf3H2O whisker is as follows: growth unit [Mg-(OH)4]2- first incorporates into the larger dimension [Mg-(OH)4]2-, then the [Mg-(OH)4]2-n combines with H2CO3 into a linear skeleton Mg2(OH)2CO3 in the same line. Mg2(OH)2CO3 combines with H2O by hydrogen bonds and ultimately transforms into Mg2(OH)2COf3H2O whisker. Magnesium carbonate whiskers have a layered structure, each of which is made of magnesium, carbon, oxygen, with H2O in between each layer. When skeletons are superimposed within the same plane as a parallelepiped one, they grow into solid cuboid-shaped whiskers. When the parallelepiped skeletons planes combine with each other through the cascading links, they grow into hollow cylindrical whiskers.展开更多
This work is concerned with the estimation from the outside of effective yield strength for the stage II consolidated material package of axisymmetric solid particles. Once an appropriate simple representative axisymm...This work is concerned with the estimation from the outside of effective yield strength for the stage II consolidated material package of axisymmetric solid particles. Once an appropriate simple representative axisymmetric unit cell is chosen, the kinematical approach of the yield design homogenization method is used in order to obtain external estimates which has been found depending on the loading history (isostatic and closed die compactions) as well as on the relative density of the material powder. For comparison purpose, finite element simulations that describe the behavior of spherical elastic plastic particles uniformly distributed inside the material powder are carried out.展开更多
For steady-state heat conduction,a new variational functional for a unit cell of composites with periodic microstructures is constructed by considering the quasi-periodicity of the temperature field and in the periodi...For steady-state heat conduction,a new variational functional for a unit cell of composites with periodic microstructures is constructed by considering the quasi-periodicity of the temperature field and in the periodicity of the heat flux fields. Then by combining with the eigenfunction expansion of complex potential which satisfies the fiber-matrix interface conditions, an eigenfunction expansion-variational method (EEVM)based on a unit cell is developed. The effective transverse thermal conductivities of doubly-periodic fiber reinforced composites are calculated, and the first-order approximation formula for the square and hexagonal arrays is presented, which is convenient for engineering application. The numerical results show a good convergency of the presented method,even though the fiber volume fraction is relatively high. Comparisons with the existing analytical and experimental results are made to demonstrate the accuracy and validity of the first-order approximation formula for the hexagonal array.展开更多
With the unstructured grid, the Finite Volume Coastal Ocean Model(FVCOM) is converted from its original FORTRAN code to a Compute Unified Device Architecture(CUDA) C code, and optimized on the Graphic Processor U...With the unstructured grid, the Finite Volume Coastal Ocean Model(FVCOM) is converted from its original FORTRAN code to a Compute Unified Device Architecture(CUDA) C code, and optimized on the Graphic Processor Unit(GPU). The proposed GPU-FVCOM is tested against analytical solutions for two standard cases in a rectangular basin, a tide induced flow and a wind induced circulation. It is then applied to the Ningbo's coastal water area to simulate the tidal motion and analyze the flow field and the vertical tide velocity structure. The simulation results agree with the measured data quite well. The accelerated performance of the proposed 3-D model reaches 30 times of that of a single thread program, and the GPU-FVCOM implemented on a Tesla k20 device is faster than on a workstation with 20 CPU cores, which shows that the GPU-FVCOM is efficient for solving large scale sea area and high resolution engineering problems.展开更多
基金Funded in part by China Postdoctoral Science Foundation(No.2013M541574)Hi-tech Research and Development Program of China(No.2012AA040209)
文摘A new method to predict the ultimate strength of fiber reinforced composites under arbitrary load condition is introduced. The micromechanics strength theory is used to perform the final failure prediction of composite laminates. The theory is based on unit cell analytic model which can provide the ply composite material properties by only using the constituent fiber and matrix properties and the laminate geometric parameters without knowing any experimental information of the laminates. To show that this method is suitable for predicting the strength of composite laminates, the micromechanics strength theory is ranked by comparing it with all the micro-level and the best two macro-level theories chosen from the World Wide Failure Exercise. The results show that this method can be used for predicting strength of any composite laminates and provide a direct reference for composite optimum design.
文摘The aim of the paper is to discover the general creep mechanisms for the short fiber reinforcement matrix composites (MMCs) under uniaxial stress states and to build a relationship between the macroscopic steady creep behavior and the material micro geometric parameters. The unit cell models were used to calculate the macroscopic creep behavior with different micro geometric parameters of fibers on different loading directions. The influence of the geometric parameters of the fibers and loading directions on the macroscopic creep behavior had been obtained, and described quantitatively. The matrix/fiber interface had been considered by a third layer, matrix/fiber interlayer, in the unit cells with different creep properties and thickness. Based on the numerical results of the unit cell models, a statistic model had been presented for the plane randomly-distributed-fiber MMCs. The fiber breakage had been taken into account in the statistic model for it starts experimentally early in the creep life. With the distribution of the geometric parameters of the fibers, the results of the statistic model agree well with the experiments. With the statistic model, the influence of the geometric parameters and the breakage of the fibers as well as the properties and thickness of, the interlayer on the macroscopic steady creep rate have been discussed.
基金This work was supported by three projects.Zhao Y received the Grant with Nos.61976236 and 2020MDJC06Bi X J received the Grant with No.20&ZD279.
文摘As one of Chinese minority languages,Tibetan speech recognition technology was not researched upon as extensively as Chinese and English were until recently.This,along with the relatively small Tibetan corpus,has resulted in an unsatisfying performance of Tibetan speech recognition based on an end-to-end model.This paper aims to achieve an accurate Tibetan speech recognition using a small amount of Tibetan training data.We demonstrate effective methods of Tibetan end-to-end speech recognition via cross-language transfer learning from three aspects:modeling unit selection,transfer learning method,and source language selection.Experimental results show that the Chinese-Tibetan multi-language learning method using multilanguage character set as the modeling unit yields the best performance on Tibetan Character Error Rate(CER)at 27.3%,which is reduced by 26.1%compared to the language-specific model.And our method also achieves the 2.2%higher accuracy using less amount of data compared with the method using Tibetan multi-dialect transfer learning under the same model structure and data set.
基金Project supported by the Research Committee of the Hong Kong Polytechnic University (No.G-U042).
文摘The damage effects of water sorption on the mechanical properties of the hydroxyapatite particle reinforced Bis-GMA/TEGDMA copolymer (HA/Bis-GMA/TEGDMA) h-ave been predicted using 3D finite cell models. The plasticizer effect on the polymer matrix was considered as a variation of its Young's modulus. Three different cell models were used to determine the influence of varying particle contents, interphase strength and moisture concentration on the debonding damage. The stress distribution pattern has been examined and the stress transfer mode clarified. The Young's modulus and fracture strength of the Bis-GMA/TEGDMA composite were also predicted using the model with and without consideration of the damage. ine Iormer results with consideration of the debonding damage are in good agreement with existing literature experimental data. The shielding effect of our proposed model and an alternative approach were discussed. The FCC cell model has also been extended to predict the critical load for the damaged and the undamaged composite subject to the 3-point flexural test.
基金supported by National Natural Science Foundation of China(No.11575159)National Natural Science Foundation of China(No.51206146)+1 种基金Zhejiang Provincial Natural Science Foundation of China(No.LY13B070004)Program for Zhejiang Leading Team of S&T Innovation(No.2013TD07)
文摘The electric energy injection from a pulsed power supply to a planar type of dielectric barrier discharge(DBD) reactor at atmospheric pressure was studied. Relations of the energy injection with barrier materials, barrier thickness, peak voltage, gap distance, electrode area,and operation temperature were experimentally investigated. The energy injection is a function of relative permittivity, barrier thickness, peak voltage, gap distance, and electrode area. The influence of operation temperature on energy injection is slight in the range of 27-300℃ but becomes obvious in the range of 300-500℃. A model was established using which the energy injection can be easily predicted.
基金Funded by the National Natural Science Foundation of China(No.51272207)
文摘The growth habit of the basic magnesium oxysulfate whisker was investigated based on the theoreticalmodelof anion coordination polyhedron growth units.It is found that typicalbasic magnesium oxysulfate whisker growth is consistent with anion tetrahedralcoordination incorporation rules.The growth units of basic magnesium oxysulfate whiskers are [Mg-(OH)_4]^(2-) and HSO_4^-.[Mg-(OH)_4]^(2-) is the favorable growth unit and whisker growth is in the direction of the [Mg-(OH)_4]^(2-) combination.A plurality of [Mg-(OH)_4]^(2-) s combine and become a larger dimensionalgrowth unit in a one-dimensionaldirection.Then HSO_4^- and larger dimensionalgrowth units connect as basic magnesium sulfate whiskers,according to the structuralcharacteristics of the basic magnesium sulfate whisker,which can guide the synthesis of magnesium hydroxide whisker.
基金Supported by the National Natural Science Foundation of China(51275145)
文摘According to the existing method including testing the frequency and establishing the relationship between moisture content and frequency, a corresponding instrument was designed. In order to further improve the accuracy and rapidity of the system, a new approach to describe the relationship between the measurement error and the temperature was proposed. The error band could be obtained and divided into several parts(based on the range of temperature) to indicate the error value that should compensate the grain moisture content for the changes in temperature. By calculating the error band at the maximum and the minimum operating temperatures, as well as by determining the error compensation value from the error band based on the measurement moisture content, the final effective result was derived.
基金Supported by Foundation of National Nature Science(59872024)
文摘To design the amorphous hydrogen storage alloy efficiently, the maximum hydrogen capacities for Zr - Ni amorphous alloy were calculated. Based on the Rhomb Unit Structure Model(RUSM) for amorphous alloy and the experimental result that hydrogen atoms exist in 3Zr1Ni and 4Zr tetrahedron interstices in Zr-Ni amolphous alloy, the numbers of 3Zr-1Ni and 4Zr tetrahedron interstices in a RUSM were calculated which correspond to the hydrogen capacity. The two extremum Zr distribution states were calculated, such as highly heterogeneous Zr distribution and homogeneous Zr distribution. The calculated curves of hydrogen capacity with different Zr contents at two states indicate that the hydrogen capacity increases with increasing Zr content and reaches its maximum when Zr is 75%. The theoretical maximum hydrogen capacity for Zr - Ni amorphous alloy is 2.0 (H/M). Meanwhile, the hydrogen capacity of heterogenous Zr distribution alloy is higher than that of homogenous one at the same Zr content. The experimental results prove the calculated results reasonable, and accordingly, the experimental results that the distribution of Zr atom in amorphous alloy occur heterogeneous after a few hydrogen absorption desorption cycles can be explained.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 20274040,20574052 and 20774066)the Program for New Century Excellent Talents in University,China (Grant No NCET-05-0538)the Natural Science Foundation of Zhejiang Province,China (Grant No R404047)
文摘In this paper the influence of a knot on the structure of a polymethylene (PM) strand in the tensile process is investigated by using the steered molecular dynamics (SMD) method. The gradual increasing of end-to-end distance, R, results in a tighter knot and a more stretched contour. That the break in a knotted rope almost invariably occurs at a point just outside the 'entrance' to the knot, which has been shown in a good many experiments, is further theoretically verified in this paper through the calculation of some structural and thermodynamic parameters. Moreover, it is found that the analyses on bond length, torsion angle and strain energy can facilitate to the study of the localization and the size of a knot in the tensile process. The symmetries of torsion angles, bond lengths and bond angles in the knot result in the whole symmetry of the knot in microstructure, thereby adapting itself to the strain applied. Additionally, the statistical property of the force-dependent average knot size illuminates in detail the change in size of a knot with force f, and therefore the minimum size of the knot in the restriction of the potentials considered in this work for a PM chain is deduced. At the same time, the difference in response to uniaxial strain, between a knotted PM strand and an unknotted one is also investigated. The force-extension profile is easily obtained from the simulation. As expected, for a given f, the knotted chain has an R significantly smaller than that of an unknotted polymer. However, the scaled difference becomes less pronounced for larger values of N, and the results for longer chains approach those of the unknotted chains.
基金Funded by the National Natural Science Foundation of China(No.51272207)
文摘From the perspective of growth units, the growth mechanism of Mg2(OH)2CO3.3H2O whisker is investigated in this paper. Results show that the growth morphology of Mg2(OH)2CO33H2O whisker is consistent with the model of anion coordination polyhedron growth units. The growth solution Raman shift of Mg2(OH)2CO;3H2O was monitored using Raman spectroscopy. The growth units are [Mg-(OH)4]2- and H2COv The growth process of Mg2(OH)2COf3H2O whisker is as follows: growth unit [Mg-(OH)4]2- first incorporates into the larger dimension [Mg-(OH)4]2-, then the [Mg-(OH)4]2-n combines with H2CO3 into a linear skeleton Mg2(OH)2CO3 in the same line. Mg2(OH)2CO3 combines with H2O by hydrogen bonds and ultimately transforms into Mg2(OH)2COf3H2O whisker. Magnesium carbonate whiskers have a layered structure, each of which is made of magnesium, carbon, oxygen, with H2O in between each layer. When skeletons are superimposed within the same plane as a parallelepiped one, they grow into solid cuboid-shaped whiskers. When the parallelepiped skeletons planes combine with each other through the cascading links, they grow into hollow cylindrical whiskers.
文摘This work is concerned with the estimation from the outside of effective yield strength for the stage II consolidated material package of axisymmetric solid particles. Once an appropriate simple representative axisymmetric unit cell is chosen, the kinematical approach of the yield design homogenization method is used in order to obtain external estimates which has been found depending on the loading history (isostatic and closed die compactions) as well as on the relative density of the material powder. For comparison purpose, finite element simulations that describe the behavior of spherical elastic plastic particles uniformly distributed inside the material powder are carried out.
基金National Natural Science Foundation of China(90716002)The Open Fund of LNM
文摘For steady-state heat conduction,a new variational functional for a unit cell of composites with periodic microstructures is constructed by considering the quasi-periodicity of the temperature field and in the periodicity of the heat flux fields. Then by combining with the eigenfunction expansion of complex potential which satisfies the fiber-matrix interface conditions, an eigenfunction expansion-variational method (EEVM)based on a unit cell is developed. The effective transverse thermal conductivities of doubly-periodic fiber reinforced composites are calculated, and the first-order approximation formula for the square and hexagonal arrays is presented, which is convenient for engineering application. The numerical results show a good convergency of the presented method,even though the fiber volume fraction is relatively high. Comparisons with the existing analytical and experimental results are made to demonstrate the accuracy and validity of the first-order approximation formula for the hexagonal array.
基金Project supported by the National Natural Science Foundation of China(Grant No.51279028,51479175)the Public Science and Technology Research Funds Projects of Ocean(Grant No.201405025)
文摘With the unstructured grid, the Finite Volume Coastal Ocean Model(FVCOM) is converted from its original FORTRAN code to a Compute Unified Device Architecture(CUDA) C code, and optimized on the Graphic Processor Unit(GPU). The proposed GPU-FVCOM is tested against analytical solutions for two standard cases in a rectangular basin, a tide induced flow and a wind induced circulation. It is then applied to the Ningbo's coastal water area to simulate the tidal motion and analyze the flow field and the vertical tide velocity structure. The simulation results agree with the measured data quite well. The accelerated performance of the proposed 3-D model reaches 30 times of that of a single thread program, and the GPU-FVCOM implemented on a Tesla k20 device is faster than on a workstation with 20 CPU cores, which shows that the GPU-FVCOM is efficient for solving large scale sea area and high resolution engineering problems.