After a brief review of studies on artificial boundaries in dynamic soil-structure interaction, a three-dimensional viscous-spring artificial boundary （VSAB） in the time domain is developed in this paper. First, the...After a brief review of studies on artificial boundaries in dynamic soil-structure interaction, a three-dimensional viscous-spring artificial boundary （VSAB） in the time domain is developed in this paper. First, the 3D VSAB equations in the normal and tangential directions are derived based on the elastic wave motion theory. Secondly, a numerical simulation technique of wave motion equations along with the VSAB condition in the time domain is studied. Finally, numerical examples of some classical elastic wave motion problems are presented and the results are compared with the associated theoretical solutions, demonstrating that high precision and adequate stability can be achieved by using the proposed 3D VSAB. The proposed 3D VSAB can be conveniently incorporated in the general finite element program, which is commonly used to study dynamic soil-structure interaction problems.展开更多
A method is presented to calculate the resistance of a high-speed displacement ship taking the effect of sinkage and trim and viscosity of fluid into account.A free surface flow field is evaluated by solving Reynolds ...A method is presented to calculate the resistance of a high-speed displacement ship taking the effect of sinkage and trim and viscosity of fluid into account.A free surface flow field is evaluated by solving Reynolds averaged Navier-Stokes(RANS) equations with volume of fluid(VoF) method.The sinkage and trim are computed by equating the vertical force and pitching moment to the hydrostatic restoring force and moment.The software Fluent,Maxsurf and MATLAB are used to implement this method.With dynamic mesh being used,the position of a ship is updated by the motion of "ship plus boundary layer" grid zone.The hull factors are introduced for fast calculating the running attitude of a ship.The method has been applied to the ship model INSEAN2340 for different Froude numbers and is found to be efficient for evaluating the flow field,resistance,sinkage and trim.展开更多
The hydrodynamic analysis of multi-floating bodies is important and widely used in marine engineering. In this study, we systematically simulated the wave diffraction problem of a fixed vertical four-cylinder structur...The hydrodynamic analysis of multi-floating bodies is important and widely used in marine engineering. In this study, we systematically simulated the wave diffraction problem of a fixed vertical four-cylinder structure in regular waves in the time domain in a viscous numerical wave tank. The hydrodynamic interaction of waves with a bottom-mounted structure consisting of four vertical cylinders arranged at the corners of a square shows a complicated interference phenomenon. In this paper, we illustrate and analyze the run-up around the structure and the corresponding wave forces. To investigate the viscous effect on the near-trapping phenomenon, we pay particular attention to investigating the waves near-trapped inside the four-cylinder structure,and make a comparative study of the viscous-and inviscid-flow solutions with the experimental measurements. The results show that the maximum wave elevation occurs on the inner side of the leeside cylinder, and that the wave elevations on the outer side of the cylinders are lower than those on the inner side. We can conclude that viscosity has an obvious damping effect on wave elevations inside the structure. The cylinders show a tendency to drift apart from each other when the near-trapping phenomenon occurs.展开更多
The three-dimensional numerical manifold method(3DNMM) method is further enriched to simulate wave propagation across homogeneous/jointed rock masses. For the purpose of minimizing negative effects from artificial bou...The three-dimensional numerical manifold method(3DNMM) method is further enriched to simulate wave propagation across homogeneous/jointed rock masses. For the purpose of minimizing negative effects from artificial boundaries, a viscous nonreflecting boundary, which can effectively absorb the energy of a wave, is firstly adopted to enrich 3DNMM. Then, to simulate the elastic recovery property of an infinite problem domain, a viscoelastic boundary, which is developed from the viscous nonreflecting boundary, is further adopted to enrich 3DNMM. Finally, to eliminate the noise caused by scattered waves, a force input method which can input the incident wave correctly is incorporated into 3DNMM. Five typical numerical tests on P/S-wave propagation across jointed/homogeneous rock masses are conducted to validate the enriched 3DNMM. Numerical results indicate that wave propagation problems within homogeneous and jointed rock masses can be correctly and reliably modeled with the enriched 3DNMM.展开更多
In this paper,it is proved that the weak solution to the Cauchy problem for the scalar viscous conservation law,with nonlinear viscosity,different far field states and periodic perturbations,not only exists globally i...In this paper,it is proved that the weak solution to the Cauchy problem for the scalar viscous conservation law,with nonlinear viscosity,different far field states and periodic perturbations,not only exists globally in time,but also converges towards the viscous shock wave of the corresponding Riemann problem as time goes to infinity.Furthermore,the decay rate is shown.The proof is given by a technical energy method.展开更多
In this paper, the in-house multifunction solver naoe-FOAM-SJTU is applied to study the resistance and wave-making performance of a high-speed catamaran sailing at different velocity in calm water. The volume of fluid...In this paper, the in-house multifunction solver naoe-FOAM-SJTU is applied to study the resistance and wave-making performance of a high-speed catamaran sailing at different velocity in calm water. The volume of fluid（VOF） method is used to capture the free interface and the finite volume method（FVM） is adopted as the discretization scheme. The hull model is fixed with initial trim and sinkage. The numerical results of the presented paper agree very well with the measurement data of model test. Wave making and vortex field are well simulated to analyze the hydrodynamic performance of a catamaran.展开更多
This paper is concerned with an ideal polytropic model of non-viscous and heatconductive gas in a one-dimensional half space. We focus our attention on the outflow problem when the flow velocity on the boundary is neg...This paper is concerned with an ideal polytropic model of non-viscous and heatconductive gas in a one-dimensional half space. We focus our attention on the outflow problem when the flow velocity on the boundary is negative and we prove the stability of the viscous shock wave and its superposition with the boundary layer under some smallness conditions.Our waves occur in the subsonic area. The intrinsic properties of our system are more challenging in mathematical analysis, however, in the subsonic area, the lack of a boundary condition on the density provides us with a special manner for defining the shift for the viscous shock wave, and helps us to construct the asymptotic profiles successfully. New weighted energy estimates are introduced and the perturbations on the boundary are handled by some subtle estimates.展开更多
To evaluate the trapped wave energy and energy loss, the problem of wave scattering by twin fixed vertical surface- piercing plates over a stepped bottom is numerically simulated using the open source package OpenFOAM...To evaluate the trapped wave energy and energy loss, the problem of wave scattering by twin fixed vertical surface- piercing plates over a stepped bottom is numerically simulated using the open source package OpenFOAM and the associated toolbox waves2Foam. The volume of fluid (VOF) method was employed to capture the free surface in the time domain. The validation of the present numerical model was performed by comparing with both the analytical and experimental results. The effects of the spacing between two plates and the configuration of stepped bottom on the hydrodynamic characteristics, such as reflection and transmission coefficients, viscous dissipation ratio, and relative wave height between the plates (termed as trapped wave energy), were examined. Moreover, the nonlinear effects of the incident wave height on the hydrodynamic characteristics were addressed as well. The results show that the step configuration can be tuned for efficient-performance of wave damping, and the optimum configurations of the step length B, the step height h1 and the spacing b, separately equaling λ/4, 3h/4, and 0.05h (λ and h are the wavelength and the water depth, respectively), are recommended for the trapping of wave energy.展开更多
A design of offshore floating structure is mainly based on the extreme response analysis due to the forces experienced. The extreme response can induce the negative air gap response and potential impact to the deck bo...A design of offshore floating structure is mainly based on the extreme response analysis due to the forces experienced. The extreme response can induce the negative air gap response and potential impact to the deck bottom of floating structure. It is important to predict the slamming load in order to check the strength of local structures which withstand the wave slamming. In recent years, studies of the effects of wind load on air gap response and slamming load are ignored. When the platform suffers the extreme wave, the wind is also harsh.Moreover, the wind load can affect the motion response of the platform. The wind load cannot be simulated easily by model test in towing tank whereas it can be simulated accurately in wind tunnel test. Though the model test results are not accurate enough for air gap and slamming load evaluation due to the loss of wind effect, they can be used as a good basis for tuning the radiation damping and viscous drag in numerical simulation. This paper aims at presenting the sensitivity analysis results of wave slamming load with respect to the wind load for the design of semi-submersible platform. As an example of semi-submersible drilling platform design, the wind tunnel test has been carried out, and the sea-keeping model test is also performed in towing tank, while the wind load effect is ignored. According to the model test results, a numerical model is tuned and validated by ANSYS AQWA. Sensitivity analysis studies of the relative velocity between water particle and platform surface and the wave slamming load with respect to the wind load are performed in time domain by the tuned numerical model.Five simulation cases about the presented platform are simulated based on the results of wind tunnel tests and sea-keeping tests. The sensitivity analysis results are valuable for the floating platform design.展开更多
This work attempts to extend the fundamental theory for classic gas dynamics to viscous compressible flow,of which aeroacoustics will naturally be a special branch.As a continuation of Part I.Unbounded fluid(Mao et al...This work attempts to extend the fundamental theory for classic gas dynamics to viscous compressible flow,of which aeroacoustics will naturally be a special branch.As a continuation of Part I.Unbounded fluid(Mao et al.,2022),this paper studies the source of longitudinal field at solid boundary,caused by the on-wall kinematic and viscous dynamic coupling of longitudinal and transverse processes.We find that at this situation the easiest choice for the two independent thermodynamic variables is the dimensionless pressure P and temperature T.The two-level structure of boundary dynamics of longitudinal field is obtained by applying the continuity equation and its normal derivative to the surface.We show that the boundary dilatation flux represents faithfully the boundary production of vortex sound and entropy sound,and the mutual generation mechanism of the longitudinal and transverse fields on the boundary does not occur symmetrically"at the samc level,but appears along a zigzag route.At the first level,it is the pressure gradient that generates vorticity unidirectionally;while at the second level,it is the vorticity that generates dilatation unidirectionally.展开更多
Either potential flow or viscous flow based model may be flawed for numerical wave simulations.The two-way coupling of potential and viscous flow models with the domain decomposition utilizing respective strengths has...Either potential flow or viscous flow based model may be flawed for numerical wave simulations.The two-way coupling of potential and viscous flow models with the domain decomposition utilizing respective strengths has been a trending research topic.In contrast to existing literatures in which closed source potential models were used,the widely used open source OceanWave3D,OpenFOAM-v2012 are used in the present research.An innovative overlapping two-way coupling strategy is developed utilizing the ghost points in OceanWave3D.To guarantee computational stability,a relaxation zone used both for outlet damping and data transfer is built over the overlapping region in OceanWave3D.The free surface elevation in the relaxation zone is directly probed in OpenFOAM while the velocity potential is indirectly built upon its temporal variation which is calculated by the free surface boundary condition using the probed velocity.Strong coupling is achieved based on the fourth-order Runge-Kutta(RK)algorithm.Both two-and three-dimensional tests including linear,nonlinear,irregular,and multi-directional irregular waves,are conducted.The effectiveness of the coupling procedure in bidirectional data transfer is fully demonstrated,and the model is validated to be accurate and efficient,thus providing a competitive alternative for ocean wave simulations.展开更多
基金National Natural Science Foundation of ChinaUnder Grant No.50478014Special Funds for Major State Basic Research Project Under Grant No.2002CB412706Research Funds from National Civil Defense Oficce of Chinafor the Tenth Five-year Plan。
文摘After a brief review of studies on artificial boundaries in dynamic soil-structure interaction, a three-dimensional viscous-spring artificial boundary （VSAB） in the time domain is developed in this paper. First, the 3D VSAB equations in the normal and tangential directions are derived based on the elastic wave motion theory. Secondly, a numerical simulation technique of wave motion equations along with the VSAB condition in the time domain is studied. Finally, numerical examples of some classical elastic wave motion problems are presented and the results are compared with the associated theoretical solutions, demonstrating that high precision and adequate stability can be achieved by using the proposed 3D VSAB. The proposed 3D VSAB can be conveniently incorporated in the general finite element program, which is commonly used to study dynamic soil-structure interaction problems.
基金the National Natural Science Foundation of China (No.50879090)the Advanced Research Program of GAD of the P.L.A (No.7131005)
文摘A method is presented to calculate the resistance of a high-speed displacement ship taking the effect of sinkage and trim and viscosity of fluid into account.A free surface flow field is evaluated by solving Reynolds averaged Navier-Stokes(RANS) equations with volume of fluid(VoF) method.The sinkage and trim are computed by equating the vertical force and pitching moment to the hydrostatic restoring force and moment.The software Fluent,Maxsurf and MATLAB are used to implement this method.With dynamic mesh being used,the position of a ship is updated by the motion of "ship plus boundary layer" grid zone.The hull factors are introduced for fast calculating the running attitude of a ship.The method has been applied to the ship model INSEAN2340 for different Froude numbers and is found to be efficient for evaluating the flow field,resistance,sinkage and trim.
基金the National Natural Science Foundation of China (11502059)the Natural Scientific Research Innovation Foundation in Harbin Institute of Technology (HIT.NSRIF.201726)+1 种基金the Weihai ScienceTechnology Development Planning (2014DXGJ10).
文摘The hydrodynamic analysis of multi-floating bodies is important and widely used in marine engineering. In this study, we systematically simulated the wave diffraction problem of a fixed vertical four-cylinder structure in regular waves in the time domain in a viscous numerical wave tank. The hydrodynamic interaction of waves with a bottom-mounted structure consisting of four vertical cylinders arranged at the corners of a square shows a complicated interference phenomenon. In this paper, we illustrate and analyze the run-up around the structure and the corresponding wave forces. To investigate the viscous effect on the near-trapping phenomenon, we pay particular attention to investigating the waves near-trapped inside the four-cylinder structure,and make a comparative study of the viscous-and inviscid-flow solutions with the experimental measurements. The results show that the maximum wave elevation occurs on the inner side of the leeside cylinder, and that the wave elevations on the outer side of the cylinders are lower than those on the inner side. We can conclude that viscosity has an obvious damping effect on wave elevations inside the structure. The cylinders show a tendency to drift apart from each other when the near-trapping phenomenon occurs.
基金supported by the Youth Innovation Promotion Association CAS(Grant No. 2020327)the National Natural Science Foundation of China(Grant Nos. 12202024, 52130905, 12272393, and 12072357)。
文摘The three-dimensional numerical manifold method(3DNMM) method is further enriched to simulate wave propagation across homogeneous/jointed rock masses. For the purpose of minimizing negative effects from artificial boundaries, a viscous nonreflecting boundary, which can effectively absorb the energy of a wave, is firstly adopted to enrich 3DNMM. Then, to simulate the elastic recovery property of an infinite problem domain, a viscoelastic boundary, which is developed from the viscous nonreflecting boundary, is further adopted to enrich 3DNMM. Finally, to eliminate the noise caused by scattered waves, a force input method which can input the incident wave correctly is incorporated into 3DNMM. Five typical numerical tests on P/S-wave propagation across jointed/homogeneous rock masses are conducted to validate the enriched 3DNMM. Numerical results indicate that wave propagation problems within homogeneous and jointed rock masses can be correctly and reliably modeled with the enriched 3DNMM.
文摘In this paper,it is proved that the weak solution to the Cauchy problem for the scalar viscous conservation law,with nonlinear viscosity,different far field states and periodic perturbations,not only exists globally in time,but also converges towards the viscous shock wave of the corresponding Riemann problem as time goes to infinity.Furthermore,the decay rate is shown.The proof is given by a technical energy method.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11072154,51379125)the National Key Basic Research Development Program of China(973 Program,Grant No.2013CB036103)
文摘In this paper, the in-house multifunction solver naoe-FOAM-SJTU is applied to study the resistance and wave-making performance of a high-speed catamaran sailing at different velocity in calm water. The volume of fluid（VOF） method is used to capture the free interface and the finite volume method（FVM） is adopted as the discretization scheme. The hull model is fixed with initial trim and sinkage. The numerical results of the presented paper agree very well with the measurement data of model test. Wave making and vortex field are well simulated to analyze the hydrodynamic performance of a catamaran.
基金the Natural Science Foundation of China(11871388)。
文摘This paper is concerned with an ideal polytropic model of non-viscous and heatconductive gas in a one-dimensional half space. We focus our attention on the outflow problem when the flow velocity on the boundary is negative and we prove the stability of the viscous shock wave and its superposition with the boundary layer under some smallness conditions.Our waves occur in the subsonic area. The intrinsic properties of our system are more challenging in mathematical analysis, however, in the subsonic area, the lack of a boundary condition on the density provides us with a special manner for defining the shift for the viscous shock wave, and helps us to construct the asymptotic profiles successfully. New weighted energy estimates are introduced and the perturbations on the boundary are handled by some subtle estimates.
基金partially supported by the National Natural Science Foundation of China(Grant Nos.11702244,51679014,51809209)the Joint Fund of Zhoushan City and Zhejiang University(Grant No.2017C82223)the Open Fund of Hunan Provincial Key Laboratory of Key Technology on Hydropower Development(Grant No.PKLHD201707)
文摘To evaluate the trapped wave energy and energy loss, the problem of wave scattering by twin fixed vertical surface- piercing plates over a stepped bottom is numerically simulated using the open source package OpenFOAM and the associated toolbox waves2Foam. The volume of fluid (VOF) method was employed to capture the free surface in the time domain. The validation of the present numerical model was performed by comparing with both the analytical and experimental results. The effects of the spacing between two plates and the configuration of stepped bottom on the hydrodynamic characteristics, such as reflection and transmission coefficients, viscous dissipation ratio, and relative wave height between the plates (termed as trapped wave energy), were examined. Moreover, the nonlinear effects of the incident wave height on the hydrodynamic characteristics were addressed as well. The results show that the step configuration can be tuned for efficient-performance of wave damping, and the optimum configurations of the step length B, the step height h1 and the spacing b, separately equaling λ/4, 3h/4, and 0.05h (λ and h are the wavelength and the water depth, respectively), are recommended for the trapping of wave energy.
基金the National Natural Science Foundation of China(No.11072149)
文摘A design of offshore floating structure is mainly based on the extreme response analysis due to the forces experienced. The extreme response can induce the negative air gap response and potential impact to the deck bottom of floating structure. It is important to predict the slamming load in order to check the strength of local structures which withstand the wave slamming. In recent years, studies of the effects of wind load on air gap response and slamming load are ignored. When the platform suffers the extreme wave, the wind is also harsh.Moreover, the wind load can affect the motion response of the platform. The wind load cannot be simulated easily by model test in towing tank whereas it can be simulated accurately in wind tunnel test. Though the model test results are not accurate enough for air gap and slamming load evaluation due to the loss of wind effect, they can be used as a good basis for tuning the radiation damping and viscous drag in numerical simulation. This paper aims at presenting the sensitivity analysis results of wave slamming load with respect to the wind load for the design of semi-submersible platform. As an example of semi-submersible drilling platform design, the wind tunnel test has been carried out, and the sea-keeping model test is also performed in towing tank, while the wind load effect is ignored. According to the model test results, a numerical model is tuned and validated by ANSYS AQWA. Sensitivity analysis studies of the relative velocity between water particle and platform surface and the wave slamming load with respect to the wind load are performed in time domain by the tuned numerical model.Five simulation cases about the presented platform are simulated based on the results of wind tunnel tests and sea-keeping tests. The sensitivity analysis results are valuable for the floating platform design.
基金supported by the National Natural Science Foundation of China(Grant Nos.12102365,91752202,11472016,11621202,and 12272371).
文摘This work attempts to extend the fundamental theory for classic gas dynamics to viscous compressible flow,of which aeroacoustics will naturally be a special branch.As a continuation of Part I.Unbounded fluid(Mao et al.,2022),this paper studies the source of longitudinal field at solid boundary,caused by the on-wall kinematic and viscous dynamic coupling of longitudinal and transverse processes.We find that at this situation the easiest choice for the two independent thermodynamic variables is the dimensionless pressure P and temperature T.The two-level structure of boundary dynamics of longitudinal field is obtained by applying the continuity equation and its normal derivative to the surface.We show that the boundary dilatation flux represents faithfully the boundary production of vortex sound and entropy sound,and the mutual generation mechanism of the longitudinal and transverse fields on the boundary does not occur symmetrically"at the samc level,but appears along a zigzag route.At the first level,it is the pressure gradient that generates vorticity unidirectionally;while at the second level,it is the vorticity that generates dilatation unidirectionally.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.52101324,52131102,51879159 and 52131102)the National Key Research and Development Program of China(Grant No.2019YFB1704200).
文摘Either potential flow or viscous flow based model may be flawed for numerical wave simulations.The two-way coupling of potential and viscous flow models with the domain decomposition utilizing respective strengths has been a trending research topic.In contrast to existing literatures in which closed source potential models were used,the widely used open source OceanWave3D,OpenFOAM-v2012 are used in the present research.An innovative overlapping two-way coupling strategy is developed utilizing the ghost points in OceanWave3D.To guarantee computational stability,a relaxation zone used both for outlet damping and data transfer is built over the overlapping region in OceanWave3D.The free surface elevation in the relaxation zone is directly probed in OpenFOAM while the velocity potential is indirectly built upon its temporal variation which is calculated by the free surface boundary condition using the probed velocity.Strong coupling is achieved based on the fourth-order Runge-Kutta(RK)algorithm.Both two-and three-dimensional tests including linear,nonlinear,irregular,and multi-directional irregular waves,are conducted.The effectiveness of the coupling procedure in bidirectional data transfer is fully demonstrated,and the model is validated to be accurate and efficient,thus providing a competitive alternative for ocean wave simulations.
