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.

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.

A spatiotemporal evolution model of a short-circuit arc to a secondary arc based on the improved charge simulation method

The initial shape of the secondary arc considerably influences its subsequent shape.To establish the model for the arcing time of the secondary arc and modify the single-phase reclosing sequence,theoretical and experimental analysis of the evolution process of the short-circuit arc to the secondary arc is critical.In this study,an improved charge simulation method was used to develop the internal-space electric-field model of the short-circuit arc.The intensity of the electric field was used as an independent variable to describe the initial shape of the secondary arc.A secondary arc evolution model was developed based on this model.Moreover,the accuracy of the model was evaluated by comparison with physical experimental results.When the secondary arc current increased,the arcing time and dispersion increased.There is an overall trend of increasing arc length with increasing arcing time.Nevertheless,there is a reduction in arc length during arc ignition due to short circuits between the arc columns.Furthermore,the arcing time decreased in the range of 0°-90°as the angle between the wind direction and the x-axis increased.This work investigated the method by which short-circuit arcs evolve into secondary arcs.The results can be used to develop the secondary arc evolution model and to provide both a technical and theoretical basis for secondary arc suppression.

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.

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.

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.

Physics-based seismic analysis of ancient wood structure:fault-to-structure simulation

Based on the domain reduction method,this study employs an SEM-FEM hybrid workflow which integrates the advantages of the spectral element method(SEM)for flexible and highly efficient simulation of seismic wave propagation in a three-dimensional(3D)regional-scale geophysics model and the finite element method(FEM)for fine simulation of structural response including soil-structure interaction,and performs a physics-based simulation from initial fault rupture on an ancient wood structure.After verification of the hybrid workflow,a large-scale model of an ancient wood structure in the Beijing area,The Tower of Buddhist Incense,is established and its responses under the 1665 Tongxian earthquake and the 1730 Yiheyuan earthquake are simulated.The results from the simulated ground motion and seismic response of the wood structure under the two earthquakes demonstrate that this hybrid workflow can be employed to efficiently provide insight into the relationships between geophysical parameters and the structural response,and is of great significance toward accurate input for seismic simulation of structures under specific site and fault conditions.

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.

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.

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 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.

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.

Class of modified parallel combined methods of real-time numerical simulation for a stiff system

A class of modified parallel combined methods of real-time numerical simulation are presented for a stiff dynamic system. By combining the parallelism across the system with the parallelism across the method, and relaxing the dependence of stage value computation on sampling time of input function, a class of modified real-time parallel combined methods are constructed. Stiff and nonstiff subsystems are solved in parallel on a parallel computer by a parallel Rosen-brock method and a parallel RK method, respectively. Their order conditions and convergences are discussed. The numerical simulation experiments show that this class of modified algorithms can get high speed and efficiency.

Harmonic balance method with alternating frequency/time domain technique for nonlinear dynamical system with fractional exponential

Comparisons of the common methods for obtaining the periodic responses show that the harmonic balance method with alternating frequency/time （HB-AFT） do- main technique has some advantages in dealing with nonlinear problems of fractional exponential models. By the HB-AFT method, a rigid rotor supported by ball bearings with nonlinearity of Hertz contact and ball passage vibrations is considered. With the aid of the Floquet theory, the movement characteristics of interval stability are deeply studied. Besides, a simple strategy to determine the monodromy matrix is proposed for the stability analysis.

A time domain induced polarization relaxation time spectrum inversion method based on a damping factor and residual correction

Relaxation time spectra （RTS） derived from time domain induced polarization data （TDIP） are helpful to assess oil reservoir pore structures. However, due to the sensitivity to the signal-to-noise ratio （SNR）, the inversion accuracy of the traditional singular value decomposition （SVD） inversion method reduces with a decrease of SNR. In order to enhance the inversion accuracy and improve robustness of the inversion method to the SNR, an improved inversion method, based on damping factor and spectrum component residual correction, is proposed in this study. The numerical inversion results show that the oscillation of the RTS derived from the SVD method increased with a decrease of SNR, which makes it impossible to get accurate inversion components. However, the SNR has little influence on inversion components of the improved method, and the RTS has high inversion accuracy and robustness. Moreover, RTS derived from core sample data is basically in accord with the pore-size distribution curve, and the RTS derived from the actual induced polarization logging data is smooth and continuous, which indicates that the improved method is practicable.