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.

Combined hybrid energy storage system and transmission grid model for peak shaving based on time series operation simulation

This study proposes a combined hybrid energy storage system(HESS) and transmission grid(TG) model, and a corresponding time series operation simulation(TSOS) model is established to relieve the peak-shaving pressure of power systems under the integration of renewable energy. First, a linear model for the optimal operation of the HESS is established, which considers the different power-efficiency characteristics of the pumped storage system, electrochemical storage system, and a new type of liquid compressed air energy storage. Second, a TSOS simulation model for peak shaving is built to maximize the power entering the grid from the wind farms and HESS. Based on the proposed model, this study considers the transmission capacity of a TG. By adding the power-flow constraints of the TG, a TSOS-based HESS and TG combination model for peak shaving is established. Finally, the improved IEEE-39 and IEEE-118 bus systems were considered as examples to verify the effectiveness and feasibility of the proposed model.

Parareal in Time Simulation of Morphological Transformation in Cubic Alloys with Spatially Dependent Composition

In this paper,a reducedmorphological transformation model with spatially dependent composition and elastic modulus is considered.The parareal in time algorithm introduced by Lions et al.is developed for longer-time simulation.The fine solver is based on a second-order scheme in reciprocal space,and the coarse solver is based on a multi-model backward Euler scheme,which is fast and less expensive.Numerical simulations concerning the composition with a randomnoise and a discontinuous curve are performed.Some microstructure characteristics at very low temperature are obtained by a variable temperature technique.

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.

Real Time Math Simulation of Contact Interaction during Spacecraft Docking and Berthing

Contact reactions of guide surfaces of assembly interfaces lead to the decreasing of theirs lateral and angular misalignments. The focus of this paper is the development of algorithms for computation of guide surfaces contact forces with acceptable engineering accuracy for real time simulation of assembly operations. Therefore, each complex guide surface is described as a set of contacting elements. Each contacting element for one＇s part can be represented by a finite set of geometric primitives which geometry is described by low order algebraic equations. So contact conditions and geometric parameters for all pairs of primitives are determined by analytical expressions. Math models are developed for two classes of contact interaction. The first class includes all cases when each contacting surface has several degrees of freedom of motion. Therefore, contact reactions introduced into differential equations of motion are calculated by using contacting elements penetrations, stiffness and damping parameters. The second class corresponds to all cases when one of contacting surfaces has insignificant inertia and only one degree of freedom of relative displacement counteracted by a spring. Here contact reactions are calculated from spring tension with any practical accuracy. This is very useful in some practical applications. Presented algorithms provide real time simulation together with some approaches for reduction of redundant comnutations.

Efficiency analysis of numerical integrations for finite element substructure in real-time hybrid simulation

Finite element（FE） is a powerful tool and has been applied by investigators to real-time hybrid simulations（RTHSs）. This study focuses on the computational efficiency, including the computational time and accuracy, of numerical integrations in solving FE numerical substructure in RTHSs. First, sparse matrix storage schemes are adopted to decrease the computational time of FE numerical substructure. In this way, the task execution time（TET） decreases such that the scale of the numerical substructure model increases. Subsequently, several commonly used explicit numerical integration algorithms, including the central difference method（CDM）, the Newmark explicit method, the Chang method and the Gui-λ method, are comprehensively compared to evaluate their computational time in solving FE numerical substructure. CDM is better than the other explicit integration algorithms when the damping matrix is diagonal, while the Gui-λ（λ = 4） method is advantageous when the damping matrix is non-diagonal. Finally, the effect of time delay on the computational accuracy of RTHSs is investigated by simulating structure-foundation systems. Simulation results show that the influences of time delay on the displacement response become obvious with the mass ratio increasing, and delay compensation methods may reduce the relative error of the displacement peak value to less than 5% even under the large time-step and large time delay.

Subgrid-scale contributions to Lagrangian time correlations in isotropic turbulence

The application of large-eddy simulation （LES） to particle-laden turbulence raises such a fundamental question as whether the LES with a subgrid scale （SGS） model can correctly predict Lagrangian time correlations （LTCs）. Most of the currently existing SGS models are constructed based on the energy budget equations. Therefore, they are able to correctly predict energy spectra, but they may not ensure the correct prediction on the LTCs. Previous researches investigated the effect of the SGS modeling on the Eulerian time correlations. This paper is devoted to study the LTCs in LES. A direct numerical simulation （DNS） and the LES with a spectral eddy viscosity model are performed for isotropic turbulence and the LTCs are calculated using the passive vector method. Both a priori and a posteriori tests are carried out. It is observed that the subgrid;scale contributions to the LTCs cannot be simply ignored and the LES overpredicts the LTCs than the DNS. It is concluded from the straining hypothesis that an accurate prediction of enstrophy spectra is most critical to the prediction of the LTCs.

Probabilistic prediction of expected ground condition and construction time and costs in road tunnels

Ground condition and construction （excavation and support） time and costs are the key factors in decision-making during planning and design phases of a tunnel project. An innovative methodology for probabilistic estimation of ground condition and construction time and costs is proposed, which is an integration of the ground prediction approach based on Markov process, and the time and cost variance analysis based on Monte-Carlo （MC） simulation. The former provides the probabilistic description of ground classification along tunnel alignment according to the geological information revealed from geological profile and boreholes. The latter provides the probabilistic description of the expected construction time and costs for each operation according to the survey feedbacks from experts. Then an engineering application to Hamro tunnel is presented to demonstrate how the ground condition and the construction time and costs are estimated in a probabilistic way. In most items, in order to estimate the data needed for this methodology, a number of questionnaires are distributed among the tunneling experts and finally the mean values of the respondents are applied. These facilitate both the owners and the contractors to be aware of the risk that they should carry before construction, and are useful for both tendering and bidding.

Digitalization of mine operations:Scenarios to benefit in real-time truck dispatching

One of the key factors in a profitable open-pit mine is the efficiency of the waste disposal system. Using GPS-technology, the truck-dispatching decisions can be made in real-time but the chosen strategy has a crucial role. Therefore, finding the optimal dispatching strategy for truck-shovel operations is extremely important. Dispatching strategies have been reported in the literature, but the comparison of these strategies is still missing. This paper illustrates the differences between the strategies by conducting a stochastic simulation study based on the data gathered from an actual mine. The findings underline the importance of the global vision in dispatching decisions.

Interval Sampling Method for Steady-State Simulation Output Analysis

A new algorithm namely the interval sampling method, applicable to the analysisof steady-state simulation output is proposed. This algorithm uses the time series analysisto carry out conrrelation analysis of the steady-state simulation output so as to obtain theobservation data which are actually uncorrelated in nature. On the basis of theseuncorrelated data gathered, some satisfactory deductions cam be made on the data under re search. A comparison between batch means method and the interval sampling method hasbeen performed by taking the M/M/l queuing system as an example. The results attestedthat the interval sampling method is mere accurate than the batch means method.

Kinetic effects of nanosecond discharge on ignition delay time

The effects of nanosecond discharge on ignition characteristics of a stoichiometric methane–air mixture without inert diluent gas were studied by numerical simulation at 0.1 MPa and an initial temperature of 1300 K. A modified non-equilibrium plasma kinetic model was developed to simulate the temporal evolution of particles produced during nanosecond discharge and its afterglow. As important roles in ignition, path fluxes of O and H radicals were analyzed in detail. Different strength of E/N and different discharge duration were applied to the discharge process in this study. And the results presented that a deposited energy of 1–30 m J·cm^(-3) could dramatically reduce the ignition delay time. Furthermore, temperature and radicals analysis was conducted to investigate the effect of non-equilibrium plasma on production of intermediate radicals. Finally, sensitivity analysis was employed to have further understanding on ignition chemistries of the mixture under nanosecond discharge.

Hydrodynamic Characteristics of Three-Bucket Jacket Foundation for Offshore Wind Turbines During the Lowering Process

The three-bucket jacket foundation is a new type of foundation for offshore wind turbine that has the advantages of fast construction speed and suitability for deep water. The study of the hoisting and launching process is of great significance to ensure construction safety in actual projects. In this paper, a new launching technology is proposed that is based on the foundation of the three-bucket jacket for offshore wind turbine. A complete time domain simulation of the launching process of three-bucket jacket foundation is carried out by a theoretical analysis combined with hydrodynamic software Moses. At the same time, the effects of different initial air storage and sea conditions on the motion response of the structure and the hoisting cable tension are studied. The results show that the motion response of the structure is the highest when it is lowered to 1.5 times the bucket height. The natural period of each degree of freedom of the structure increases with the increase of the lowering depth. The structural motion response and the hoisting cable tension vary greatly in the early phases of Stages Ⅰ and Ⅲ, smaller in Stage Ⅱ, and gradually stabilize in the middle and late phases of Stage Ⅲ.

Implementation of Distributed Control System for Rice Mill Using C#

The paper presents the distributed control system for rice mill using C#language.The real-time manufacturing system can be implemented by utilizing the signal from the real time control units that is more operative than other old-fashioned control systems in the extent of modern industrial days.The software-based Distributed Control System(DCS)is novel fashionable than any other conventional control systems in the state-ofthe-art manufacturing developments.This research study emphasizes on the implementation of the DCS-based rice mill using visual C#.net.The Industrial Ethernet(IE)is realized between the top level controller for the operator and the controlled station for the remote devices.The model of client-server approach is more appropriate for the automation and manufacturing research purposes.In this study,the computer graphical simulation of the complete control development is depicted in real-time status quo by visual C#language under Visual Studio 2008 software.The parallel ports in the computers of remote terminal level and the master terminal level controllers have been interconnected with port interface coding by visual C#program.

Optimal Error Estimates of the Semi-Discrete Local Discontinuous Galerkin Method and Exponential Time Differencing Schemes for the Thin Film Epitaxy Problem without Slope Selection

In this paper,we prove the optimal error estimates in L2 norm of the semidiscrete local discontinuous Galerkin(LDG)method for the thin film epitaxy problem without slope selection.To relax the severe time step restriction of explicit time marching methods,we employ a class of exponential time differencing(ETD)schemes for time integration,which is based on a linear convex splitting principle.Numerical experiments of the accuracy and long time simulations are given to show the efficiency and capability of the proposed numerical schemes.

SIMULATION OF TANK SLOSHING BASED ON OPENFOAM AND COUPLING WITH SHIP MOTIONS IN TIME DOMAIN

Tank sloshing in ship cargo is excited by ship motions, which induces impact load on tank wall and then affects the ship motion. Wave forces acting on ship hull and the retardation function are solved by using three-dimensional frequency domain theory and an impulse response function method based on the potential flow theory, and global ship motion is examined coupling with nonlinear tank sloshing which is simulated by viscous flow theory. Based on the open source Computational Fluid Dynamics （CFD） development platform Open Field Operation and Manipulation （OpenFOAM）, numerical calculation of ship motion coupled with tank sloshing is achieved and the corresponding numerical simulation and validation are carried out. With this method, the interactions of wave, ship body and tank sloshing are completely taken into consideration. This method has quite high efficiency for it takes advantage of potential flow theory for outer flow field and viscous flow theory for inside tank sloshing respectively. The numerical and experimental results of the ship motion agree well with each other.

Dynamic Response Analysis of the Equivalent Water Depth Truncated Point of the Catenary Mooring Line

The real-time computer-controlled actuators are used to connect the truncated parts of moorings and risers in the active hybrid model testing system. This must be able to work in model-scale real time, based on feedback input from the floater motions. Thus, mooring line dynamics and damping effects are artificially simulated in real time, based on a computer-based model of the problem. In consideration of the nonlinear characteristics of the sea platform catenary mooring line, the equations of the mooring line motion are formulated by using the lumped-mass method and the dynamic response of several points on the mooring line is investigated by the time and frequency domain analysis method. The dynamic response of the representative point on the mooring line is analyzed under the condition of two different corresponding upper endpoint movements namely sine wave excitation and random wave excitation. The corresponding laws of the dynamic response between the equivalent water depth truncated points at different locations and the upper endpoint are obtained, which can provide technical support for further study of the active hybrid model test.

Nonlinear Numerical Analysis of Vortex-Induced Vibration of A Three-Dimensional Deepwater Steep Wave Riser with Large Deformation Features

The cross-flow(CF)vortex-induced vibration(VIV)of a deepwater steep wave riser(SWR)subjected to uniform or shear flow loads is investigated numerically.The model is based on a three-dimensional(3D)nonlinear elastic rod theory coupled with a wake oscillator model.In this numerical simulation,the nonlinear motion equations of the riser with large deformation features are established in a global coordinate system to avoid the transformation between global and local coordinate systems,and are discretized with the time-domain finite element method(FEM).A wakeoscillator model is employed to study the vortex shedding,and the lift force generated by the wake flow is described in a van der Pol equation.A Newmark-βiterative scheme is used to solve their coupling equation for the VIV response of the SWR.The developed model is validated against the existing experimental results for the VIV response of the top-tension riser(TTR).Then,the numerical simulations are executed to determine VIV characteristics of the SWR.The effects of both flow velocity and the spanwise length of the flow field on the drag coefficient in the inline(IL)direction and the lift coefficient in the CF direction are investigated systematically.The results indicate that compared with TTR,the low frequency and multi-modal vibration are the main components of the SWR due to the large deformation and flexible characteristics.For shear flow,the multi-frequency resonance dominates the VIV response of the SWR,especially at the hang-off segment.

Hydration Process and Crack Tendency of Concrete Based on Resistivity and Restrained Shrinkage Crack

Hydration process, crack potential and setting time of concrete grade C30, C40 and C50 were monitored by using a non-contact electrical resistivity apparatus, a novel plastic ring mould and penetration resistance methods, respectively. The results show the highest resistivity of C30 at the early stage until a point when C50 accelerated and overtook the others. It has been experimentally confirmed that the crossing point of C30 and C50 corresponds to the final setting time of C50. From resistivity derivative curve, four different stages were observed upon which the hydration process is classified; these are dissolution, induction, acceleration and deceleration periods. Consequently, restrained shrinkage crack and setting time results demonstrated that C50 set and cracked the earliest. The cracking time of all the samples occurred within a reasonable experimental period thus the novel plastic ring is a convenient method for predicting concrete＇s crack potential. The highest inflection time（t_i） obtained from resistivity curve and the final setting time（t_f） were used with crack time（t_c） in coming up with mathematical models for the prediction of concrete＇s cracking age for the range of concrete grade considered. Finally, an ANSYS numerical simulation supports the experimental findings in terms of the earliest crack age of C50 and the crack location.

Wind-induced vibration of single-layer reticulated shell structures

Aiming at the dynamic response of reticulated shell structures under wind load,systematic parameter analyses on wind-induced responses of Kiewitt6-6 type single-layer spherical reticulated shell structures and three-way grid single-layer cylindrical reticulated shell structures were performed with the random simulation method in time domain,including geometric parameters,structural parameters and aerodynamic parameters.Moreover,a wind-induced vibration coefficient was obtained,which can be a reference to the wind-resistance design of reticulated shell structures.The results indicate that the geometric parameters are the most important factor influencing wind-induced responses of the reticulated shell structures;the wind-induced vibration coeffi-cient is 3.0-3.2 for the spherical reticulated shell structures and that is 2.8-3.0 for the cylindrical reticula-ted shell structures,which shows that the wind-induced vibration coefficients of these two kinds of space frames are well-proportioned.