Fully nonlinear modeling of radiated waves generated by floating flared structures

The nonlinear radiated waves generated by a structure in forced motion, are simulated numerically based on the potential theory. A fully nonlinear numerical model is developed by using a higher-order boundary element method （HOBEM）. In this model, the instantaneous body position and the transient free surface are updated at each time step. A Lagrangian technique is employed as the time marching scheme on the free surface. The mesh regridding and interpolation methods are adopted to deal with the possible numerical instability. Several auxiliary functions are proposed to calculate the wave loads indirectly, instead of directly predicting the temporal derivative of the velocity potential. Numerical experiments are carried out to simulate the heave motions of a submerged sphere in infinite water depth, the heave and pitch motions of a truncated flared cylinder in finite depth. The results are verified against the published numerical results to ensure the effectiveness of the proposed model. Moreover, a series of higher harmonic waves and force components are obtained by the Fourier transformation to investigate the nonlinear effect of oscillation frequency. The difference among fully nonlinear, body-nonlinear and linear results is analyzed. It is found that the nonlinearity due to free surface and body surface has significant influences on the numerical results of the radiated waves and forces.

Fully Nonlinear Simulations of Wave Resonance by An Array of Cylinders in Vertical Motions

The finite element method （FEM） is employed to analyze the resonant oscillations of the liquid confined within multiple or an array of floating bodies with fully nonlinear boundary conditions on the free surface and the body surface in two dimensions. The velocity potentials at each time step are obtained through the FEM with 8-node quadratic shape functions. The finite element linear system is solved by the conjugate gradient （CG） method with a symmetric successive overelaxlation （SSOR） preconditioner. The waves at the open boundary are absorbed by the combination of the damping zone method and the Sommerfeld-Orlanski equation. Numerical examples are given by an array of floating wedge- shaped cylinders and rectangular cylinders. Results are provided for heave motions including wave elevations, profiles and hydrodynamic forces. Comparisons are made in several cases with the results obtained from the second order solution in the time domain. It is found that the wave amplitude in the middle region of the array is larger than those in other places, and the hydrodynamic force on a cylinder increases with the cylinder closing to the middle of the array.

An Improved Nearshore Wave Breaking Model Based on the Fully Nonlinear Boussinesq Equations

This paper aims to propose an improved numerical model for wave breaking in the nearshore region based on the fully nonlinear form of Boussinesq equations. The model uses the κ equation turbulence scheme to determine the eddy viscosity in the Boussinesq equations. To calculate the turbulence production term in the equation, a new formula is derived based on the concept of surface roller. By use of this formula, the turbulence production in the one-equation turbulence scheme is directly related to the difference between the water particle velocity and the wave celerity. The model is verified by Hansen and Svendsen's experimental data (1979) in terms of wave height and setup and setdown. The comparison between the model and experimental results of wave height and setup and setdown shows satisfactory agreement. The modeled turbulence energy decreases as waves attenuate in the surf zone. The modeled production term peaks at the breaking point and decreases as waves propagate shoreward. It is also suggested that both convection and diffusion play their important roles in the transport of turbulence energy immediately after wave breaking. When waves approach to the shoreline, the production and dissipation of turbulence energy are almost balanced. By use of the slot technique for the simulation of the movable shoreline boundary, wave runup in the swash zone is well simulated by the present model.

Fully Nonlinear Simulation for Fluid/Structure Impact：A Review

This paper presents a review of the work on fluid/structure impact based on inviscid and imcompressible liquid and irrotational flow. The focus is on the velocity potential theory together with boundary element method （BEM）. Fully nonlinear boundary conditions are imposed on the unknown free surface and the wetted surface of the moving body. The review includes （1） vertical and oblique water entry of a body at constant or a prescribed varying speed, as well as free fall motion, （2） liquid droplets or column impact as well as wave impact on a body, （3） similarity solution of an expanding body. It covers two dimensional （2D）, axisymmetric and three dimensional （3D） cases. Key techniques used in the numerical simulation are outlined, including mesh generation on the multivalued free surface, the stretched coordinate system for expanding domain, the auxiliary function method for decoupling the mutual dependence of the pressure and the body motion, and treatment for the jet or the thin liquid film developed during impact.

A Fully Nonlinear HOBEM with the Domain Decomposition Method for Simulation of Wave Propagation and Diffraction

A higher-order boundary element method(HOBEM) for simulating the fully nonlinear regular wave propagation and diffraction around a fixed vertical circular cylinder is investigated. The domain decomposition method with continuity conditions enforced on the interfaces between the adjacent sub-domains is implemented for reducing the computational cost. By adjusting the algorithm of iterative procedure on the interfaces, four types of coupling strategies are established, that is, Dirchlet/Dirchlet-Neumman/Neumman(D/D-N/N), Dirchlet-Neumman(D-N),Neumman-Dirchlet(N-D) and Mixed Dirchlet-Neumman/Neumman-Dirchlet(Mixed D-N/N-D). Numerical simulations indicate that the domain decomposition methods can provide accurate results compared with that of the single domain method. According to the comparisons of computational efficiency, the D/D-N/N coupling strategy is recommended for the wave propagation problem. As for the wave-body interaction problem, the Mixed D-N/N-D coupling strategy can obtain the highest computational efficiency.

Numerical Investigations on Transient Behaviours of Two 3-D Freely Floating Structures by Using a Fully Nonlinear Method

Two floating structures in close proximity are very commonly seen in offshore engineering. They are often subjected to steep waves and, therefore, the transient effects on their hydrodynamic features are of great concem. This paper uses the quasi arbitrary Lagrangian Eulerian finite element method （QALE-FEM）, based on the fully nonlinear potential theory （FNPT）, to numerically investigate the interaction between two 3-D floating structures, which undergo motions with 6 degrees of freedom （DOFs）, and are subjected to waves with different incident angles. The transient behaviours of floating structures, the effect of the accompanied structures, and the nonlinearity on the motion of and the wave loads on the structures are the main focuses of the study. The investigation reveals an important transient effects causing considerably larger structure motion than that in steady state. The results also indicate that the accompanied structure in close proximity enhances the interaction between different motion modes and results in stronger nonlinearity causing 2hal-order component to be of similar significance to the fundamental one.

Numerical Comparison for Focused Wave Propagation Between the Fully Nonlinear Potential Flow and the Viscous Fluid Flow Models

Numerical simulations on focused wave propagation are carried out by using three types of numerical models,including the linear potential flow,the nonlinear potential flow and the viscous fluid flow models.The wave-wave interaction of the focused wave group with different frequency bands and input wave amplitudes is examined,by which the influence of free surface nonlinearity and fluid viscosity on the related phenomenon of focused wave is investigated.The significant influence of free surface nonlinearity on the characteristics of focused wave can be observed,including the increased focused wave crest,delayed focused time and downstream shift of focused position with the increase of input amplitude.It can plot the evident difference between the results of the nonlinear potential flow and linear potential flow models.However,only a little discrepancy between the nonlinear potential flow and viscous fluid flow models can be observed,implying the insignificant effect of fluid viscosity on focused wave behavior.Therefore,the nonlinear potential flow model is recommended for simulating the non-breaking focused wave problem in this study.

A STUDY ON GRADIENT BLOW UP FOR VISCOSITY SOLUTIONS OF FULLY NONLINEAR,UNIFORMLY ELLIPTIC EQUATIONS

We investigate sharp conditions for boundary and interior gradient estimates of continuous viscosity solutions to fully nonlinear, uniformly elliptic equations under Dirichlet boundary conditions. When these conditions are violated, there can be blow up of the gradient in the interior or on the boundary of the domain. In particular we de- rive sharp results on local and global Lipschitz continuity of continuous viscosity solutions under more general growth conditions than before. Lipschitz regularity near the boundary allows us to predict when the Dirichlet condition is satisfied in a classical and not just in a viscosity sense, where detachment can occur. Another consequence is this： if interior gra- dient blow up occurs, Perron-type solutions can in general become discontinuous, so that the Dirichlet problem can become unsolvable in the class of continuous viscosity solutions.

DEFORMING METRICS WITH POSITIVE CURVATURE BY A FULLY NONLINEAR FLOW

By studying a fully nonlinear flow deforming conformal metrics on a compact and connected manifold, we prove the long time existence and the exponential convergence of the solutions of the flow for any initial metric g0 with the Schouten tensor Ag0 ∈ Γk.

Fully Nonlinear Shallow Water Waves Simulation Using Green-Naghdi Theory

Green-Naghdi （G-N） theory is a fully nonlinear theory for water waves. Some researchers call it a fully nonlinear Boussinesq model. Different degrees of complexity of G-N theory are distinguished by ＂levels＂ where the higher the level, the more complicated and presumably more accurate the theory is. In the research presented here a comparison was made between two different levels of G-N theory, specifically level II and level III G-N restricted theories. A linear analytical solution for level III G-N restricted theory was given. Waves on a planar beach and shoaling waves were both simulated with these two G-N theories. It was shown for the first time that level III G-N restricted theory can also be used to predict fluid velocity in shallow water. A level III G-N restricted theory is recommended instead of a level II G-N restricted theory when simulating fullv nonlinear shallow water waves.

POSITIVE RADIAL SOLUTIONS OF FULLY NONLINEAR ELLIPTIC EQUATIONS IN R^n

By the Schauder-Tychonoff fixed-point theorem, we inyestigate the existence and asymptotic behavior of positive radial solutions of fully nonlinear elliptic equations in Re. We give some sufficient conditions to guarantee the existence of bounded and unbounded radial solutions and consider the nonexistence of positive solution in Rn.

Existence and Uniqueness to a Fully Nonlinear Version of the Loewner–Nirenberg Problem Dedicated to Celebrate the Sixtieth Anniversary of USTC

We consider the problem of finding on a given Euclidean domainof dimension n≥3 a complete conformally flat metric whose Schouten curvature A satisfies some equations of the form f(λ(−A))=1.This generalizes a problem considered by Loewner and Nirenberg for the scalar curvature.We prove the existence and uniqueness of such metric when the boundary δΩ is a smooth bounded hypersurface(of codimension one).When δΩ contains a compact smooth submanifold ∑ of higher codimension with δΩ\∑ being compact,we also give a‘sharp’condition for the divergence to infinity of the conformal factor near ∑ in terms of the codimension.

FULLY NONLINEAR WAVE COMPUTATIONS FOR ARBITRARY FLOATING BODIES USING THE DELTA METHOD

Fully nonlinear water wave problems are solved using Eulerian-Lagrangian timestepping methods in conjunction with a desingularized approach to solve the mixed boundary valueproblem that arises at each time step. In the desingularized approach, the singularities generatingthe flow field are outside the fluid domain. This allows the singularity distribution to be replacedby isolated Rankine sources with the corresponding reduction in computational complexity andcomputer time. A moving boundary technique is applied to e-liminate the reflection waves fromlimited computational boundaries. Examples of the use of the method in three-dimensions are givenfor the exciting forces acting on a modified wigley hull and Series 60 are presented. The numericalresults show good agreements with those of experiments.

A fully nonlinear and weakly dispersive water wave model for simulating the propagation,interaction,and transformation of solitary waves

This paper presents the development of a theoretical model of fully nonlinear and weakly dispersive(FNWD)waves and numerical techniques for simulating the propagation,interaction,and transformation of solitary waves.Using the standard expansion method and without the limit of small nonlinear parameter defined as the ratio of the wave height versus water depth,a set of model equations describing the FNWD waves in a domain of moderately varying bottom topography are formulated.Exact solitary wave solutions satisfying the FNWD equations are also derived.Numerically,a time-accurate and stabilized finite-element code to solve the governing equations is developed for wave simulations.The solitary wave solutions of FNWD,weakly nonlinear and weakly dispersive(WNWD),and Laplace equations based models in terms of wave profile and phase speed are compared to examine their related features and differences.Investigations on the overtaking collision of two unidirectional solitary waves of different amplitudes,i.e.,ax and a2 where a1>a2,are carried out using both the FNWD and WNWD water wave models.Selected cases by running the FNWD and WNWD models are performed to identify the critical values of a1/a2 for forming a flattened merging wave peak,which is the condition used to determine if the stronger wave is to pass through the weaker one or both waves are to remain separated during the encountering process.It is interesting to note the critical values of a1/a2 obtained from the FNWD and WNWD models are found to be different and greater than the value of 3 proposed by Wu through the theoretical analysis of the Korteweg-de Vries(KdV)equations.Finally,the phenomena of wave splitting and nonlinear focusing of a solitary wave propagating over a three-dimensional semicircular shoal are simulated.The results obtained from both the FNWD and WNWD models showing the fission process of separating a main solitary wave into multiple waves of decreasing amplitudes are presented,compared,and discussed.

A fully nonlinear approach for efficient ship-wave simulation

This paper presents an efficient time-domain method for simulating nonlinear ship waves.The proposed method,implemented in an earth-fixed coordinate system,integrates a compact boundary element domain within a high-order spectral layer,enabling accurate modeling of both near-field and far-field ship waves.An overset mesh method and an attention mechanism are employed to track the moving ship.The effectiveness of the method is validated through simulations of Wigley and Series 60 ships sailing at various speeds.The numerical results,including the nonlinear wave run-up at the ship bow,surface pressure distribution on the hull,and the ship resistance,agree well with experimental data and published numerical results,confirming that the method is capable of accurately simulating the nonlinear ship waves.

A priori estimates for classical solutions of fully nonlinear elliptic equations

For the fully nonlinear uniformly elliptic equation F(D2u) = 0, it is well known that the viscosity solutions are C2,α if the nonlinear operator F is convex (or concave). In this paper, we study the classical solutions for the fully nonlinear elliptic equation where the nonlinear operator F is locally C1,β a.e. for any 0 < β < 1. We will prove that the classical solutions u are C2,α. Moreover, the C2,α norm of u depends on n,F and the continuous modulus of D2u.

Fully Nonlinear Parabolic Equations and the Dini Condition

Interior regularity results for viscosity solutions of fully nonlinear uniformly parabolic equations under the Dini condition, which improve and generalize a result due to Kovats, are obtained by the use of the approximation lemma.

MIXED INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS FOR ONE-DIMENSIONAL FULLY NONLINEAR SECOND ORDER ELLIPTIC AND PARABOLIC EQUATIONS

This paper is concerned with developing accurate and efficient numerical methods for one-dimensional fully nonlinear second order elliptic and parabolic partial differential equations （PDEs）. In the paper we present a general framework for constructing high order interior penalty discontinuous Galerkin （IP-DG） methods for approximating viscosity solutions of these fully nonlinear PDEs. In order to capture discontinuities of the second order derivative uxx of the solution u, three independent functions p1,p2 and p3 are introduced to represent numerical derivatives using various one-sided limits. The proposed DG frame- work, which is based on a nonstandard mixed formulation of the underlying PDE, embeds a nonlinear problem into a mostly linear system of equations where the nonlinearity has been modified to include multiple values of the second order derivative uxz. The proposed framework extends a companion finite difference framework developed by the authors in [9] and allows for the approximation of fully nonlinear PDEs using high order polynomials and non-uniform meshes. In addition to the nonstandard mixed formulation setting, another main idea is to replace the fully nonlinear differential operator by a numerical operator which is consistent with the differential operator and satisfies certain monotonicity （called g-monotonicity） properties. To ensure such a g-monotonicity, the crux of the construction is to introduce the numerical moment, which plays a critical role in the proposed DG frame- work. The g-monotonicity gives the DG methods the ability to select the mathematically ＂correct＂ solution （i.e., the viscosity solution） among all possible solutions. Moreover, the g-monotonicity allows for the possible development of more efficient nonlinear solvers as the special nonlinearity of the algebraic systems can be explored to decouple the equations. This paper also presents and analyzes numerical results for several numerical test problems which are used to guage the accuracy and efficiency of the proposed DG methods.

The Neumann Problem for a Class of Fully Nonlinear Elliptic Partial Differential Equations

In this paper,we establish global C2 estimates to the Neumann problem for a class of fully nonlinear elliptic equations.As an application,we prove the existence and uniqueness of k-admissible solutions to the Neumann problems.

Global Solution for Fully Nonlinear Parabolic Equations in One-Dimensional Space

We discuss the existence of global classical solution for the uniformly parabolic equation ■ut=a(x,t,u,ux,uxx)+b（x,t,u,ux）,（x,t）∈（-1,1）×(0,T], u（±1,t）=0,u（x,0）=■（x）, where a is strongly nonlinear with respect to uxxand ■ is not necessarily small.We also deal with nonuniform case.