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.

Stabilized Multi-domain Simulation Algorithms and Their Application in Simulation Platform for Forging Manipulator

Most researches focused on the analytical stabilized algorithm for the modular simulation of single domain, e.g., pure mechanical systems. Only little work has been performed on the problem of multi-domain simulation stability influenced by algebraic loops. In this paper, the algebraic loop problem is studied by a composite simulation method to reveal the internal relationship between simulation stability and system topologies and simulation unit models. A stability criterion of multi-domain composite simulation is established, and two algebraic loop compensation algorithms are proposed using numerical iteration and approximate function in multi-domain simulation. The numerical stabilized algorithm is the Newton method for the solution of the set of nonlinear equations, and it is used here in simulation of the system composed of mechanical system and hydraulic system. The approximate stabilized algorithm is the construction of response surface for inputs and outputs of unknown unit model, and it is utilized here in simulation of the system composed of forging system, mechanical and hydraulic system. The effectiveness of the algorithms is verified by a case study of multi-domain simulation for forging system composed of thermoplastic deformation of workpieces, mechanical system and hydraulic system of a manipulator. The system dynamics simulation results show that curves of motion and force are continuous and convergent. This paper presents two algorithms, which are applied to virtual reality simulation of forging process in a simulation platform for a manipulator, and play a key role in simulation efficiency and stability.

Design of Multi-domain Collaborative Simulation System for Dual Clutch Transmission Based on CAN Bus

A multi-domain collaborative simulation(MDCS) system for dual clutch transmission(DCT) was presented based on controller area network(CAN) bus.An interface card of CAN bus was designed,in which MDCS subsystems were linked as the nodes according to the interface mode of MDCS.A DCT simulation model was established based on Matlab/Simdriveline,whose running process was accurately controlled by the designed control system.The playback system of vehicle state(VPS) was proposed whose input was the road-test data,with a real vehicle test environment for the development of transmission control unit(TCU) being provided.A DCT kinematic system model was set up,and the running status of DCT parts could be displayed in real time.The functions of MDCS were verified based on the extra-urban driving cycle(EUDC) and the vehicle road-test data respectively.The results show the functions of MDCS are accomplished,and the unified supporting platform for the development of TCU is achieved by MDCS.

Solving Strategy for Repetitious Simulation of Multi-domain Physical System Based on Modelica

Repetitious simulation after modifying parameters of multi-domain physical system based on Modelica often appears in model experiment and optimization design. At present, the solvers based on Modelica need calculate all the coupled blocks during every simulation run after updating parameters. Based on discussing scale decomposition methods of simulation model, subdivision solving strategy and minimum solving strategy are put forward to improve the efficiency of repetitious simulation, by which the numerical solution of the simulation model can be achieved by only calculating the solving sequence influenced by altered parameters. A simplified model of aircraft is used to demonstrate the efficiency of the strategies presented.

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.

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

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.

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.

Phase field simulation of the 180°domain-switching process in PbTiO_3 single crystal under an antiparallel electric field

The process of 180°domain switching in PbTiO_3 single crystal under an antiparallel electric field was investigated by the three-dimensional phase field simulation,especially the effect of electric field on the type and duration of domain switching.It is found that the polarization reversal of domains takes place under an antiparallel electric field in PbTiO_3 single crystal.The results of the phase field simulation indicate that there is only 90°domain switching under a weak electric field.With the rise of the electric field,180°domain switching appears.If the electric field is strengthened further,90°domain switching disappears and the duration of domain switching is shortened.

Physics-informed neural network-based petroleum reservoir simulation with sparse data using domain decomposition

Recent advances in deep learning have expanded new possibilities for fluid flow simulation in petroleum reservoirs.However,the predominant approach in existing research is to train neural networks using high-fidelity numerical simulation data.This presents a significant challenge because the sole source of authentic wellbore production data for training is sparse.In response to this challenge,this work introduces a novel architecture called physics-informed neural network based on domain decomposition(PINN-DD),aiming to effectively utilize the sparse production data of wells for reservoir simulation with large-scale systems.To harness the capabilities of physics-informed neural networks(PINNs)in handling small-scale spatial-temporal domain while addressing the challenges of large-scale systems with sparse labeled data,the computational domain is divided into two distinct sub-domains:the well-containing and the well-free sub-domain.Moreover,the two sub-domains and the interface are rigorously constrained by the governing equations,data matching,and boundary conditions.The accuracy of the proposed method is evaluated on two problems,and its performance is compared against state-of-the-art PINNs through numerical analysis as a benchmark.The results demonstrate the superiority of PINN-DD in handling large-scale reservoir simulation with limited data and show its potential to outperform conventional PINNs in such scenarios.

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.

Numerical Simulation of a Domain-Tessellation Pattern on a Spherical Surface Using a Phase Field Model

A numerical simulation scheme is proposed to analyze domain tessellation and pattern formation on a spherical surface using the phase-field method. A multi-phase-field model is adopted to represent domain growth, and the finite-difference method (FDM) is used for numerical integration. The lattice points for the FDM are distributed regularly on a spherical surface so that a mostly regular triangular domain division is realized. First, a conventional diffusion process is simulated using this lattice to confirm its validity. The multi-phase-field equation is then applied, and pattern formation processes under various initial conditions are simulated. Unlike pattern formation on a flat plane, where the regular hexagonal domains are always stable, certain different patterns are generated. Specifically, characteristic stable patterns are obtained when the number of domains, n, is 6, 8, or 12;for instance, a regular pentagonal domain division pattern is generated for n = 12, which corresponds to a regular dodecahedron.

Domain ontology driven composable simulation framework

The simulation composability is one of important development directions in simulation domain. The key issue is the semantic composability. Currently there is no feasible approach to realize the effective semantic representation and composition of simulation components. Based on domain knowledge and Web Ontology Language (OWL), this paper proposes a composable simulation framework, which includes conceptual model semantics, model components semantics, model framework semantics, and simulation scenario semantics. Additionally, all the semantics are utilized in the model components development process, the simulation system development process, and the simulation system execution process respectively. The consistency checking among those semantics is also proposed. The detailed mapping processes between different semantic models can help to build the domain ontology driven composable simulation system.

Some of Advances on the Error Analysis in the Frequency Domain for a Digital Simulation Model

This paper surveys a number of recent advances in the error analysis in the frequency domain for a digital simulation model. It is emphasized to discuss the errors in characteristic roots and transfer funcnon of the digital simulation model, the frequency domain errors of the data transfers between thesimulation submodels, and some compensation methods for the errors. Some of the questions to be answered are also presented.

Parallel Simulation of 3D Wave Propagation by Domain Decomposition

In order to perform large scale numerical simulation of wave propagation in 3D heterogeneous multiscale viscoelastic media, Finite Difference technique and its parallel implementation based on domain decomposition is used. A couple of typical statements of borehole geophysics are dealt with—sonic log and cross well measurements. Both of them are essentially multiscales, which claims to take into account heterogeneities of very different sizes in order to provide reliable results of simulations. Locally refined spatial grids help us to avoid the use of redundantly tiny grid cells in a target area, but cause some troubles with uniform load of Processor Units involved in computations. We present results of scalability tests together with results of numerical simulations for both statements performed for some realistic models.

MD Simulations of the P53 oncoprotein structure: the effect of the Arg273→His mutation on the DNA binding domain

A comparative molecular dynamics (MD) simulation study was performed on the p53 oncoprotein to investigate the effect of the Arg273His (R273H) mutation on the p53→DNA Binding Domain (DBD). The two p53 dimer structures of the wild-type and mutant Arg273His (R273H) were simulated with the same thermodynamic and environmental parameters. The obtained results demonstrate that the induced Arg273His mutation has a considerable effect on the p53→DNA close contact interaction and changes the picture of hydrogen formation. The Arg273His mutation, in some cases, destroys the existing native hydrogen bond, but, in other cases, forms a strong p53→DNA hydrogen bond, which is not proper for the native protein. The MD simulation results illustrate some molecular mechanism of the conformational changes of the Arg273His key amino acid residue in the p53→DNA binding domain, which might be important for the understanding of the physiological functioning of the p53 protein and the origin of cancer.

Three-dimensional numerical simulation of geothermal fi eld in space-wavenumber domain

Large-scale,fine,and efficient numerical simulation of a geothermal field plays an important role in geothermal energy development.Confronted with the problem of large computation and high storage requirements for complex underground models in a three-dimensional(3-D)numerical simulation of a geothermal fi eld,a mixed space-wavenumber domain 3-D numerical simulation algorithm is proposed in this paper.According to the superposition principle of temperature field,the geothermal field is decomposed into background and abnormal temperature fi elds for calculation.The uniform layered model is used to solve the background field.When the abnormal field is solved,the horizontal two-dimensional(2-D)Fourier transform is used to transform the 3-D diff erential equation satisfi ed by an abnormal field into a series of one-dimensional ordinary differential equations with diff erent wavenumbers,which greatly reduces the calculation and storage.The unit division of an ordinary diff erential equation is fl exible,and the calculation amount is small.The algorithm fully takes advantage of the effi ciency of the Fourier transform and the quickness of the catch-up method to solve linear equations with a fixed bandwidth,which effectively improves the computational efficiency.Compared with the COMSOL Multiphysics professional simulation finite element software,the time consumption and memory requirements of the algorithm proposed in this paper are reduced by multiple orders of magnitude in terms of ensuring accuracy and the same mesh division.The more the number of calculated nodes is,the more obvious is the advantage.We design models to study the thermal conductivity,heat fl ux boundary,regional tectonic morphology,and topographic relief of the geothermal fi eld distribution.A 3-D geophysical model is developed based on topographic elevation data,geothermal geology,and geophysical exploration data in the Qiabuqia area of Gonghe Basin,Qinghai Province,China.Numerical simulation of the geothermal fi eld in this area is realized,which shows that the algorithm is suitable for precise and effi cient simulation of an arbitrary complex terrain and geological conditions.

Simulation on Domain Rotation Path and Magnetostriction of Terfenol-D Alloy

Based on a simplified domain rotation model, the rotation path of internal domains and cor- responding magnetostriction of <112>, <111> oriented single Tedenol-D crystals under com- pressive prestresses have been simulated. Comparisons with results of experiment and other calculation have been made. Results of simulation showed that the <111> oriented single crystal has better low-field magnetostriction properties than the <112> oriented one. Under a compressive prestress of 10 MPa, up to 2300×10-6 saturation magnetostriction of <111> oriented crystal has been obtained at 800 Oe, while for <112> oriented. 1600×10-6 maximum magnetostriction has been reached at 1000 Oe.

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.