Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φn-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated. 相似文献
1 IntroductionThe El Ni’o atmospheric physics oscillation is anabnormal phenomenon involved in the tropical Pacificocean- atmosphere interactions. Studies on the El Ni(ophenomenon are very attractive (Lin et al., 2000; Linet al., 2001, 2002; Wang, 2001; Feng et al., 2001;Feng et al., 2002; Liu et al., 2002; Wei and Chen,2003; Xie et al., 2002; Zhu et al., 2002; Pu et al.,2003; Gu et al., 2004; Yu and Liu, 2004). Lin and Mo(2004), Mo and Lin (2004) and Mo et al. (2004) alsostudied a s… 相似文献
A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed. 相似文献
A new method to determine semi-analytical solutions of one-dimensional contaminant transport problem with nonlinear sorption is described. This method is based on operator splitting approach where the convective transport is solved exactly and the diffusive transport by finite volume method. The exact solutions for all sorption isotherms of Freundlich and Langmuir type are presented for the case of piecewise constant initial profile and zero diffusion. Very precise numerical results for transport with small diffusion can be obtained even for larger time steps (e.g., when the Courant-Friedrichs-Lewy (CFL) condition failed). 相似文献
This paper presents results recently obtained for generating site-specific ground motions needed for design of critical facilities. The general approach followed in developing these ground motions using either deterministic or probabilistic criteria is specification of motions for rock outcrop or very firm soil conditions followed by adjustments for site-specific conditions. Central issues in this process include development of appropriate attenuation relations and their uncertainties, differences in expected motions between Western and Eastern North America, and incorporation of site-specific adjustments that maintain the same hazard level as the control motions, while incorporating uncertainties in local dynamic material properties. For tectonically active regions, such as the Western United States (WUS), sufficient strong motion data exist to constrain empirical attenuation relations for M up to about 7 and for distances greater than about 10–15 km. Motions for larger magnitudes and closer distances are largely driven by extrapolations of empirical relations and uncertainties need to be substantially increased for these cases.
For the Eastern United States (CEUS), due to the paucity of strong motion data for cratonic regions worldwide, estimation of strong ground motions for engineering design is based entirely on calibrated models. The models are usually calibrated and validated in the WUS where sufficient strong motion data are available and then recalibrated for applications to the CEUS. Recalibration generally entails revising parameters based on available CEUS ground motion data as well as indirect inferences through intensity observations. Known differences in model parameters such as crustal structure between WUS and CEUS are generally accommodated as well. These procedures are examined and discussed. 相似文献