A Time-Domain Numerical Simulation for Free Motion Responses of Two Ships Advancing in Head Waves

The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems within a time domain framework,the free water surface needs to simultaneously satisfy both the kinematic and dynamic boundary conditions of the free water surface.This provides conditions for adding artificial damping layers.Using the Runge−Kutta method to solve equations related to time.An upwind differential scheme is used in the present method to deal with the convection terms on the free surface to prevent waves upstream.Through the comparison with the available experimental data and other numerical methods,the present method is proved to have good mesh convergence,and satisfactory results can be obtained.The constant panel method is applied to calculate the hydrodynamic interaction responses of two parallel ships advancing in head waves.Numerical simulations are conducted on the effects of forward speed,different longitudinal and lateral distances on the motion response of two modified Wigley ships in head waves.Then further investigations are conducted on the effects of different ship types on the motion response.

Simulation of pavement roughness based on time domain model

In order to describe pavement roughness more intuitively and effectively, a method of pavement roughness simulation, i.e., the stochastic sinusoidal wave, is introduced. The method is based on the primary idea that pavement roughness is denoted as the sum of numerous sines or cosines with stochastic phases, and uses the discrete spectrum to approach the target stochastic process. It is a discrete numerical method used to simulate pavement roughness. According to a given pavement power spectral density （PSD） coefficient, under the condition that the character of displacement frequency based on the time domain model is in accordance with the given pavement surface spectrum, the pavement roughness is optimized to stochastic equivalent vibrations by computer simulation, and the curves that describe pavement roughness under each grade are obtained. The results show that the stochastic sinusoidal wave is suitable for simulation of measured pavement surface spectra based on the time domain model. The method of the stochastic sinusoidal wave is important to the research on vehicle ride comfort due to its rigorous mathematical derivation, extensive application range and intuitive simulation curve. Finally, a roughness index defined as the nominal roughness index （NRI） is introduced, and it has correlation with the PSD coefficient.

Time Domain Simulation of a One Line Failure for a DP-assisted Mooring System

This paper focuses on the research of a semi-submersible platform equipped with a DP-assisted mooring system. Based on the working principles of the DP-assisted mooring system and the model of the platform motion, a time domain simulation program is applied to analyze the impact, in the case of one line failure, on the platform motion, power consumption of the thrusters and the tension of the mooring lines. The results show that, under the 10-year wind dominant, a one line failure will have little impact on the tension of the mooring lines. When the failure line is windward, the power consumption will increase greatly with a weakened position of accuracy. However when the failure line is leeward, the power consumption will be reduced with a partly strengthened oosition of accuracy.

Time-domain simulation for water wave radiation by floating structures (Part A)

Direct time-domain simulation of floating structures has advantages:it can calculate wave pressure fields and forces directly; and it is useful for coupled analysis of floating structures with a mooring system. A time-domain boundary integral equation method is presented to simulate three-dimensional water wave radiation problems. A stable form of the integration free-surface boundary condition (IFBC) is used to update velocity potentials on the free surface. A multi-transmitting formula (MTF) method with an artificial speed is introduced to the artificial radiation boundary (ARB). The method was applied to simulate a semi-spherical liquefied natural gas (LNG) carrier and a semi-submersible undergoing specified harmonic motion. Numerical parameters such as the form of the ARB, and the time and space discretization related to this method are discussed. It was found that a good agreement can be obtained when artificial speed is between 0.6 and 1.6 times the phase velocity of water waves in the MTF method. A simulation can be done for a long period of time by this method without problems of instability, and the method is also accurate and computationally efficient.

Time Domain Simulation of Transient Responses of Very Large Floating Structures Under Unsteady External Loads

A time domain finite element method （FEM） for the analysis of transient elastic response of a very large floating structure （VLFS） subjected to arbitrary time-dependent external loads is presented. This method is developed directly in time domain and the hydrodynamic problem is formulated based on linear, inviscid and slightly compressible fluid theory and the structural response is analyzed on the thin plate assumption. The time domain finite element procedure herein is validated by comparing numerical results with available experimental data. Finally, the transient elastic response of a pontoon-type VLFS under the landing of an airplane is computed by the proposed time domain FEM. The time histories of the applied force and the position and velocity of an airplane during landing are modeled with data from a Boeing 747-400 jumbo jet.

Finite difference time domain method forward simulation of complex geoelectricity ground penetrating radar model

The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.

An Approximation to Energy Dissipation in Time Domain Simulation of Sloshing Waves Based on Linear Potential Theory

This paper proposes a new approximation to energy dissipation in time domain simulation of sloshing waves by use of a linear potential theory. The boundary value problem is solved by the NURBS （non-uniform rational B-spline） higher-order panel method, in which a time-domain Green function is employed. The energy dissipation is modeled by changing the boundary condition on solid boundaries. Model experiments are carried out in a partially filled rectangular tank with forced horizontal motion. Sloshing-induced internal pressures and horizontal force obtained numerically and experimentally are compared with each other. It is observed that the present energy dissipation approximation can help produce a fair agreement between experimental forces and those of numerical simulations.

Comparison of delay compensation methods for real-time hybrid simulation using frequency-domain evaluation index

The delay compensation method plays an essential role in maintaining the stability and achieving accurate real-time hybrid simulation results. The effectiveness of various compensation methods in different test scenarios, however, needs to be quantitatively evaluated. In this study, four compensation methods （i.e., the polynomial extrapolation, the linear acceleration extrapolation, the inverse compensation and the adaptive inverse compensation） are selected and compared experimentally using a frequency evaluation index （FEI） method. The effectiveness of the FEI method is first verified through comparison with the discrete transfer fimction approach for compensation methods assuming constant delay. Incomparable advantage is further demonstrated for the FEI method when applied to adaptive compensation methods, where the discrete transfer function approach is difficult to implement. Both numerical simulation and laboratory tests with predefined displacements are conducted using sinusoidal signals and random signals as inputs. Findings from numerical simulation and experimental results demonstrate that the FEI method is an efficient and effective approach to compare the performance of different compensation methods, especially for those requiring adaptation of compensation parameters.

NONLINEAR DRIFT MOTION SIMULATION OF MULTI-CHAIN MOORING BUOY IN TIME DOMAIN

This paper deals with the nonlinear effect of the drift motion of multi-chain mooring buoys. The buoy's motion in time domain is determined for the case that the wave and mooring force are nonlinear. The Kotorayama's method of hydrodynamic effect in a single mooring chain is expanded to multi-mooring chains. The time history of drift motion of the mooring buoy in regular waves and wave groups is calculated. The relation between the drift motion and the wave height or difference frequency is discussed. It can be shown that the effect of the hydrodynamic force acting on the mooring chain has remarkable influence on the total drift motion of the mooring buoy. In a wave group, its amplitude is mainly controled by the wave height and has little relation with difference frequency.

Time-Domain Higher-Order Boundary Element Method for Simulating High Forward-Speed Ship Motions in Waves

The hydrodynamic performance of a high forward-speed ship in obliquely propagating waves is numerically examined to assess both free motions and wave field in comparison with a low forward-speed ship.This numerical model is based on the time-domain potential flow theory and higher-order boundary element method,where an analytical expression is completely expanded to determine the base-unsteady coupling flow imposed on the moving condition of the ship.The ship in the numerical model may possess different advancing speeds,i.e.stationary,low speed,and high speed.The role of the water depth,wave height,wave period,and incident wave angle is analyzed by means of the accurate numerical model.It is found that the resonant motions of the high forward-speed ship are triggered by comparison with the stationary one.More specifically,a higher forward speed generates a V-shaped wave region with a larger elevation,which induces stronger resonant motions corresponding to larger wave periods.The shoaling effect is adverse to the motion of the low-speed ship,but is beneficial to the resonant motion of the high-speed ship.When waves obliquely propagate toward the ship,the V-shaped wave region would be broken due to the coupling effect between roll and pitch motions.It is also demonstrated that the maximum heave motion occurs in beam seas for stationary cases but occurs in head waves for high speeds.However,the variation of the pitch motion with period is hardly affected by wave incident angles.

Inversion of Seabed Geotechnical Properties in the Arctic Chukchi Deep Sea Basin Based on Time Domain Adaptive Search Matching Algorithm

The chirp sub-bottom profiler,for its high resolution,easy accessibility and cost-effectiveness,has been widely used in acoustic detection.In this paper,the acoustic impedance and grain size compositions were obtained based on the chirp sub-bottom profiler data collected in the Chukchi Plateau area during the 11th Arctic Expedition of China.The time-domain adaptive search matching algorithm was used and validated on our established theoretical model.The misfit between the inversion result and the theoretical model is less than 0.067%.The grain size was calculated according to the empirical relationship between the acoustic impedance and the grain size of the sediment.The average acoustic impedance of sub-seafloor strata is 2.5026×10^(6) kg(s m^(2))^(-1)and the average grain size(θvalue)of the seafloor surface sediment is 7.1498,indicating the predominant occurrence of very fine silt sediment in the study area.Comparison of the inversion results and the laboratory measurements of nearby borehole samples shows that they are in general agreement.

Multi-Time Scale Operation and Simulation Strategy of the Park Based on Model Predictive Control

Due to the impact of source-load prediction power errors and uncertainties,the actual operation of the park will have a wide range of fluctuations compared with the expected state,resulting in its inability to achieve the expected economy.This paper constructs an operating simulation model of the park power grid operation considering demand response and proposes a multi-time scale operating simulation method that combines day-ahead optimization and model predictive control(MPC).In the day-ahead stage,an operating simulation plan that comprehensively considers the user’s side comfort and operating costs is proposed with a long-term time scale of 15 min.In order to cope with power fluctuations of photovoltaic,wind turbine and conventional load,MPC is used to track and roll correct the day-ahead operating simulation plan in the intra-day stage to meet the actual operating operation status of the park.Finally,the validity and economy of the operating simulation strategy are verified through the analysis of arithmetic examples.

Time-Domain Simulation for Coupled Motions of Three Barges Moored Side-by-Side in Floatover Operation

Simulating the coupled motions of multiple bodies in the time domain is a complex problem because of the strong hydrodynamic interactions and coupled effect of various mechanical connectors. In this study, we investigate the hydrodynamic responses of three barges moored side-by-side in a floatover operation in the frequency and time domains. In the frequency domain, the damping lid method is adopted to improve the overestimated hydrodynamic coefficients calculated from conventional potential flow theory. A time-domain computing program based on potential flow theory and impulse theory is compiled for analyses that consider multibody hydrodynamic interactions and mechanical effects from lines and fenders. Correspondingly, an experiment is carried out for comparison with the numerical results. All statistics, time series, and power density spectra from decay and irregular wave tests are in a fairly good agreement.

Modeling and Simulations in Symmetrical Supercapacitors Using Time Domain Mathematical Expressions

This study presents the deduction of time domain mathematical equations to simulate the curve of the charging process of a symmetrical electrochemical supercapacitor with activated carbon electrodes fed by a source of constant electric potential in time ε and the curve of the discharge process through two fixed resistors. The first resistor RCo is a control that aims to prevent sudden variations in the intensity of the electric current i1(t) present at the terminals of the electrochemical supercapacitor at the beginning of the charging process. The second resistor is the internal resistance RA of the ammeter used in the calculation of the intensity of the electric current i1(t) over time in the charging and discharging processes. The mathematical equations generated were based on a 2R(C + kUC(t)) electrical circuit model and allowed to simulate the effects of the potential-dependent capacitance (kUC(t)) on the charge and discharge curves and hence on the calculated values of the fixed capacitance C, the equivalent series resistance (ESR), the equivalent parallel resistance (EPR) and the electrical potential dependent capacitance index k.

Nonlinear Time-Domain Theory for the Simulation of Moored Floating Body Motion

Nonlinear wave loads can induce low-frequency and high-frequency resonance motions of a moored platform in deep water. For the analysis of the nonlinear response of an offshore platform under the action of irregular waves, the most widely used method in practice is the Cummins method, in which the second-order exciting forces in the time domain are computed by a two-term Volterra series model based on incident waves, first-order body motion response, and quadratic transfer functions(QTFs). QTFs are bichromatic waves acting on a body and are computed in the frequency domain in advance. For moving bodies, QTFs are related to the first-order body response, which is to be determined in the simulation process of body motion response but is unknown in the computation procedure of QTFs. In solving this problem, Teng and Cong(2017) proposed a method to divide the QTFs into different components,which are unrelated to the body response. With the application of the new QTF components, a modified Cummins method can be developed for the simulation of the nonlinear response of a moored floating platform. This paper presents a review of the theory.

Optical simulation of in-plane-switching blue phase liquid crystal display using the finite-difference time-domain method

The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell＇s equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.

A Full Dynamic Voltage Stability Research Based on Time-domain Simulation

The voltage stability is substantially a dynamic stability, but the primary method which is more mature and engineering practical to analyze the stability of voltage is still static analysis. The time-domain simulation is an important measure in research of complex power grid. With the development of full dynamic simulation technology, the research of dynamic voltage stability by using full dynamic simulation program which is based on time-domain simulation can be carried out. This paper uses full dynamic simulation program in dynamic voltage stability research, lays special stress on research in how generator over-excitation limiter functioned and influence in dynamic voltage stability research, and raise 2 methods and steps to figure out dynamic stable voltage in both over-excitation counted and not counted. The simulation results of examples indicate the correctness and effectiveness of these methods, and also fully verify the function and influence of generator over-excitation limiter in full dynamic voltage stability research.

Temperature and Pressure Effects on Terahertz Time-Domain Spectroscopy of Crystalline Methedrine by Molecular Dynamics Simulation

In this study, the terahertz time-domain spectroscopy (THz-TDS) of crystalline methedrine, which is one of the illegal drugs, is performed using molecular dynamics simulation by the Fourier transform of time derivative auto-correlation functions of the dipole moment. In order to accurately detect the drugs from samples, it is necessary to build a complete database for terahertz spectra under different external conditions from theoretical calculation, which are hardly obtained from the experiments directly. Our results show remarkable consistency with the available experimental data in the frequency range of 10 - 100 cm-1 indicating that the presented method has significant capability to simulate terahertz spectra at various conditions. We investigated the effects of temperature and pressure on THz-TDS by simulating the system at temperature range between 78.4 K and 400 K at pressures up to 100 atm. Results show the spectral features of THz-TDS both in intensity and profile are highly sensitive to the variation of temperature and with a lower magnitude to the variation of pressure. The vanishing, rebuilding and shifting of spectral peaks are due to the complex mechanisms such as the anharmonicity, shifting in the vibration energy levels, formation and destruction of hydrogen-binding and the deformation of the potential energy surface during the environment changing. This improved our understanding for complicated THz-TDS of crystalline methedrine and would be useful for assignment of the practical measurements.

Viscoacoustic prestack reverse time migration based onthe optimal time-space domain high-order finite-difference method

Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.