A 1D time-domain method for in-plane wave motions in a layered half-space

A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the layered half-space is discretized on the basis of the propagation characteristic of elastic wave according to the Snell law. Then, the finite element method with lumped mass and the central difference method are incorporated to establish 2D wave motion equations, which can be transformed into 1D equations by discretization principle and explicit finite element method. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained, and the wave motions in layered half-space are finally determined based on the characteristic of traveling wave. Both the theoretical analysis and the numerical results demonstrate that the proposed method has high accuracy and good stability.

Second-Order Wave Diffraction Around 3-D Bodies by A Time-Domain Method

A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.

Time-Domain Analysis of the Relative Motions Between Side-by-Side FLNG and LNGC Under Oblique Waves

Strong hydrodynamic interactions during the side-by-side offloading operation between floating liquefied natural gas(FLNG) and liquefied natural gas carrier(LNGC) can induce high risks of collision. The weather vane effect of a single-point mooring system normally results in the satisfactory hydrodynamic performance of the side-by-side configuration in head seas. Nevertheless, the changes in wave directions in real sea conditions can significantly influence the relative motions. This article studies the relative motions of the side-by-side system by using the theoretical analysis method and the numerical calculation method. Based on the three-dimensional potential theory modified by artificial damping-lid method, the frequency-domain hydrodynamic coefficients can be improved to calculate the retardation functions for the multi-body problem. An in-house code is then developed to perform the time-domain simulation of two vessels, through which the relative motions are subsequently obtained. A range of oblique waves are chosen for the extensive calculation of relative motions between the two vessels, which are further analyzed in terms of the phase shift of motion responses induced by specific resonant wave patterns. Investigation results show that wave directions have a significant influence on the relative sway, roll, and yaw motions. Under the circumstance that the absolute phase shift between the roll motions of two vessels approaches 180°, stronger relative motions are induced when LNGC is on the weather side.Moreover, the gap water resonances at high frequencies tend to cause the dangerous opposed oscillation of two vessels in the sway and yaw modes, whereas FLNG reduces the gap water resonances and relative motions when located on the weather side.

Time-domain Wave Propagation in Dispersive Media

The equation of time-domain wave propagation in dispersive media and the explicit beam propagation method are presented in this paper.This method is demonstrated by the short optical pulses in a directional coupler with second order dispersive effect and shows to be in full agreement with former references.This method is simple,easy and practical.

Time-Domain Nonlinear Wave-Current Interaction with A Steep Wave Riser Considering Internal Flow Effect

The nonlinear dynamic response induced by the wave-current interaction on a deepwater steep wave riser(SWR)is numerically investigated based on a three-dimensional(3 D)time-domain finite element method(FEM).The governing equation considering internal flow is established in the global coordinate system.The whole SWR consists of three segments:the decline segment,buoyancy segment and hang-off segment,in which the buoyancy segment is wrapped by several buoyancy modules in the middle section,leading to the arch bend and sag bend.A Newmark-β iterative scheme is adopted for the accurate analysis to solve the governing equation and update the dynamic response at each time step.The proposed method is verified through the published results for the dynamic response of steel catenary riser(SCR)and static configuration of steel lazy wave riser(SLWR).Simulations are executed to study the influence of wave height,current velocity/direction,internal flow density/velocity and top-end pressure on the tension,configuration and bending moment of the SWR.The results indicate that the influence of the current on the configuration and mechanical behavior of the SWR is greater than that of the wave,especially in the middle section.With increasing current velocity,the suspending height of the middle section drops,meanwhile,its bending moment decreases accordingly,but the tension increases significantly.For a fixed external load,the increasing internal flow density induces the amplification of the tension at the hang-off segment and the mitigation at the decline segment,while the opposite trend occurs at the bending moment.

A BOUNDARY ELEMENT METHOD FOR STEADY SHIP WAVE-MAKING POTENTIAL PROBLEMS

A boundary element method for three-dimensional steady ship wave-making potential problems is established with the Rankine source function as its fundamental solution. In the treatment of the linearized free surface condition, one-sided, upstream finite difference operator (FDO) is used to suppress the upstream waves, and the equation of the disturbance velocity is established so that the first order FDO can be used in place of the second order FDO. Compared with the method with the second order FDO, the current method gives better precision and stability. Numerical examples are presented for verification.

Making Waves

Cities begin to hire chiefs to protect rivers and lakes from pollution IN addition to his work as chief of the Communist Party of China(CPC)Committee of Yixing,a city in Jiangsu Province in east China,Wang Zhongsu was given an extra task.In 2011,he was appointed chief

Propagation characteristics of oblique incidence terahertz wave through non-uniform plasma

The propagation characteristics of oblique incidence terahertz(THz) waves through non-uniform plasma are investigated by the shift-operator finite-difference time-domain(SO-FDTD) method combined with the phase matching condition.The electron density distribution of the non-uniform plasma is assumed to be in a Gaussian profile. Validation of the present method is performed by comparing the results with those obtained by an analytical method for a homogeneous plasma slab.Then the effects of parameters of THz wave and plasma layer on the propagation properties are analyzed. It is found that the transmission coefficients greatly depend on the incident angle as well as on the thickness of the plasma, while the polarization of the incident wave has little influence on the propagation process in the range of frequency considered in this paper. The results confirm that the THz wave can pass through the plasma sheath effectively under certain conditions,which makes it a potential candidate to overcome the ionization blackout problem.

Numerical investigation of influence of wave directionality on the water resonance at a narrow gap between two rectangular barges

A three-dimensional time-domain potential flow model with second-order nonlinearity was applied to simulate the wave resonance in a gap between two side-by-side rectangular barges. In the model, the velocity potential was decomposed into the incident potential and unknown scattered potential which was obtained by solving the boundary integral equation. The fourth-order predict-correct method was applied to enforce the free surface conditions in the time integration. The influence of the wave direction on the first and second-order gap surface elevations was investigated. The results reveal that the incident wave angle does not affect the resonant wave frequency and the maximum surface elevation at resonance always occurs at the middle location along the gap. However, the corresponding maximum wave surface elevation at resonance varies with the incident wave angle. The location of the maximum wave elevation shifts either upstream or downstream along the gap, depending on the relative magnitude of incident wave frequency to the resonant frequency.

Electromagnetic Waves from TNT Explosions

Experimental observations indicate that electromagnetic (EM) radiation is emitted after the detonation of high explosives (HE) charges. The movement of ionized atoms, particles and electrons seems to be the underlying cause. Expansion of the detonation products (DP) drives a strong (~1 kb) shock in surrounding air. This forms an intense thermal wave (T ~11,000 K) with duration of ~20 microseconds. Such temperatures create significant ionization of the air. According to Ohm’s Law, movement of ionized patches generates current;and according to the Biot-Savart Law, such currents induce electric and magnetic fields. We investigate these effects through numerical simulations of TNT explosions. A high-order Godunov scheme is used to integrate the one-dimensional conservation laws of gasdynamics. An extremely fine grid (10 microns) was needed to get converged temperature and conductivity profiles. The gasdynamic solution provided a source current, which was fed into a time-domain Green’s function code to predict three-dimensional electromagnetic waves emanating from the TNT explosion. This analysis clearly demonstrates one mechanism—the Boronin current—as the source of EM emissions from TNT explosions, but other mechanisms are also possible.

Time-history responses on the surface by regularly distributed enormous embedded cavities:Incident SH-waves

The time-history responses of the surface were obtained for a linear elastic half-plane including regularly distributed enormous embedded circular cavities subjected to propagating obliquely incident plane SH-waves. An advanced numerical approach named half-plane time-domain boundary element method（BEM）, which only located the meshes around the cavities, was used to create the model. By establishing the modified boundary integral equation（BIE）independently for each cavity and forming the matrices, the final coupled equation was solved step-by-step in the timedomain to obtain the boundary values. The responses were developed for a half-plane with 512 cavities. The amplification patterns were also obtained to illustrate the frequencydomain responses for some cases. According to the results,the presence of enormous cavities affects the scattering and diffraction of the waves arrived to the surface. The introduced method can be recommended for geotechnical/mechanical engineers to model structures in the fields of earthquake engineering and composite materials.

Air Gap Prediction for Floating Bodies Using a 3D Numerical Wave Tank Approach

Computations for air gap response of a semisubmersible platform based on a 3D numerical wave tank approach are presented.The developed method is in time domain and can consider nonlinearities associated with incident wave and hydrostatic forces exactly in determining the body response, but the interaction hydrodynamics of radiation and diffraction are based on simplified linearization assumptions. The incident wave can be defined by any suitable wave theory and here defined by a fully nonlinear numerical wave model. After verifying the present computations results in its degenerated linearized version against the usual linear 3D Green function–based frequency-domain results for air gap predictions, systematic comparative studies are undertaken between linear and the approximate nonlinear solutions. It is found that nonlinear computations can yield considerably conservative predictions as compared to fully linear calculations, amounting to a difference of up to 30%–40% in the minimum air gap in steep ambient incident waves at high and moderate frequencies.

Study on Wave Added Resistance of a Deep-V Hybrid Monohull Based on Panel Method

In this paper, a panel method based on three dimensional potential flow theory is used to study the problem of wave added resistance. The time-domain motion response of Wigely III ship in head waves is calculated by AQWA, and then the wave added resistance of ship is obtained by near-field pressure integration method. By comparing the calculated results with the experimental data in literature, it is shown that the variation trend and peak value are in good agreement, and the accuracy and efficiency meet the research requirements. Based on the above mentioned method, the wave added resistance of a deep-V hybrid monohull in head waves is studied. The motions and wave added resistances of the deep-V hybrid monohull and the deep-V original ship advancing in head waves with various forward speed and wave frequencies are calculated and analyzed. The results show that the longitudinal motion response of the deep-V hybrid monohull is effectively suppressed and the wave added resistance is obviously reduced, the new type of ship has good engineering application prospects. The present method provides an approach of satisfactory accuracy and efficiency to predict wave added resistance of ships voyaging in waves.

Application of A Fully Nonlinear Higher-Order Element Method for Modelling Three-Dimensional Wave Entry of A Cone

Fully nonlinear water entry of a cone into waves with gravity effect has been analyzed based on a three-dimensional(3D)higher-order boundary method(HOBEM).The total velocity potential at the initial time is divided into the incident and scattering components.In the subsequent time steps,the solution of the velocity potential is defined as a whole through instantaneous boundary conditions.Based on the image theory,a modified Green function is applied to establish the integral equations so that only half of the calculation domain is considered and the seabed can be excluded.The free surface elevation is tracked along a given azimuth plane in the polar coordinate system,while the horizontal motion of the water particle is updated by using a segment-spring analogy method,which redistributes nodes and maintains mesh connectivity according to linear stiffness.An auxiliary function is applied to solve the pressure distribution,instead of directly calculating time derivative of the velocity potential.The high accuracy of the present numerical method is achieved through a detailed convergence study and comparison with results in the literature.Simulations are emphatically performed to examine the effects of gravity,wave nonlinearity,entry location,and oblique entry.

The transmission enhancement of a THz pulse through an Ag/Ag_2 O layer detected by terahertz time-domain spectroscopy

This paper reports a new way to detect the enhanced transmission of a THz electromagnetic wave through an Ag/Ag2O layer by THz-TDS （time-domain spectroscopy）. As the THz beam illuminates the sub-wavelength Ag particles gained by Ag2O thermal decomposition, the evanescent wave is generated. The evanescent wave is coupled by a 500μm-GaAs substrate, which attaches behind the Ag/Ag2O layer, and then it transmits to the far field to be detected. The experimental results indicate that the transmitting amplitude is enhanced, as well as the frequent shifting and spectra broadening.

Assessment of right ventricular afterload by pressure waveform analysis in acute pulmonary hypertension

AIM: To characterize hydraulic right ventricle (RV) af- terload by pulmonary arterial pressure waveform analy- sis in an acute pulmonary hypertension (PH) model. METHODS: Pulmonary artery (PA) flow and pressure were recorded in six anesthetized sheep. Acute iso- baric PH was induced by phenylephrine (active) and PA mechanical constriction (passive). We estimated the amplitude of the forward and reflected pressure waves according to the inflection point. In most cases the in- flection pressure was smooth, thus the inflection point was defined as the time at which the first derivative ofpulmonary arterial pressure reached its first minimum. We calculated the input and characteristic (Z C , time- domain Li method) impedances, the capacitance index (stroke volume/pulse pressure), the augmentation index (AI) (reflected pressure/pulse pressure), the frac- tional pulse pressure (pulse pressure/mean pressure) and the wasted energy generated by the RV due to wave reflection during ejection (E W ). RESULTS: Pulse pressure, fractional pulse pressure, AI and Z C increased and capacitance index decreased during passive PH with respect to control (P < 0.05). In contrast, Z C and the capacitance index did not change and E W and the AI decreased during active PH. Pulse pressure correlated with E W and Z C and the AI was cor- related with E W (r > 0.6, P < 0.05). CONCLUSION: PA pressure waveform analysis al- lows the quantification of the dynamic RV afterload. Prospective clinical studies will be necessary to validate this time-domain approach to evaluate the dynamic RV afterload in chronic PH.

DYNAMIC INTERACTION OF PLANE WAVES WITH A UNILATERALLY FRICTIONALLY CONSTRAINED INCLUSION-TIME DOMAIN BOUNDARY ELEMENT ANALYSIS

A 2D time domain boundary element method(BEM)is developed to solve the transient scattering of plane waves by a unilaterally frictionally constrained inclusion.Coulomb friction is assumed along the contact interface.The incident wave is assumed strong enough so that localized slip and separation take place along the interface.The present problem is in effect a nonlinear boundary value problem since the mixed boundary conditions involve unknown intervals (slip,separation and stick regions).In order to determine the unknown intervals,an iterative technique is developed.As an example,we consider the scattering of a circular cylinder embedded in an infinite solid.

Nonlinear Analysis of Bidirectional Vortex-Induced Vibration of A Deepwater Steep Wave Riser Subjected to Oblique Currents

An improved three-dimensional(3D)time-domain couple model is established in this paper to simulate the bidirectional vortex-induced vibration(VIV)of a deepwater steep wave riser(SWR)subjected to oblique currents.In this model,the nonlinear motion equations of the riser are established in the global coordinate system based on the slender rod theory with the finite element method.Van der Pol equations are used to describe the lift forces induced by the x-and y-direction current components,respectively.The coupled equations at each time step are solved by a Newmark-βiterative scheme for the SWR VIV.The present model is verified by comparison with the published experimental results for a top-tension riser.Then,a series of simulations are executed to determine the influences of the oblique angle/velocity of the current,different top-end positions and the length of the buoyancy segment on the VIV displacement,oscillating frequency as well as hydrodynamic coefficients of the SWR.The results demonstrate that there exists a coupled resonant VIV corresponding to x-direction and y-direction,respectively.However,the effective frequency is almost identical between the vibrations at the hang-off segment along x and y directions.The addition of the buoyancy modules in the middle of the SWR has a beneficial impact on the lift force of three segments and simultaneously limits the VIV response,especially at the decline segment and the hang-off segments.Additionally,the incident current direction significantly affects the motion trajectory of the SWR which mainly includes the fusiform and rectangle shapes.

激波/湍流边界层干扰中的自适应控制技术

从激波/湍流边界层干扰机理以及流动控制的迫切需求入手,从自适应涡流发生器、自适应鼓包、自适应微射流以及自适应次流循环四个方面对激波/湍流边界层干扰中的自适应控制技术研究进展进行了总结。分析认为,结合AI技术发展自适应流动控制技术,加速控制方式智能化,可作为新一代高超声速飞行器宽速域飞行的重要技术手段。具体来说,就是通过调节外加激励对高超声速飞行器不同区域实现局部流动加/减速、气动热防护、气动控制等功能,根据流场参数建立控制反馈回路,自适应调整局部流场结构,以满足工程实际需求。