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 two-dimensional interior ballistics model of solid propellant electrothermal-chem ical launch with discharge rod plasma generator

Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-dimensional interior ballistics model of the solid propellant ETC gun(2D-IB-SPETCG)is presented to describe the process of the ETC launch.Both calculated pressure and projectile muzzle velocity accord well with the experimental results.The feasibility of the 2D-IB-SPETCG model is proved.Depending on the experimental data and initial parameters,detailed distribution of the ballistics parameters can be simulated.With the distribution of pressure and temperature of the gas phase and the propellant,the influence of plasma during the ignition process can be analyzed.Because of the radial flowing plasma,the propellant in the area of the DRPG is ignited within 0.01 ms,while all propellant in the chamber is ignited within 0.09 ms.The radial ignition delay time is much less than the axial delay time.During the ignition process,the radial pressure difference is less than 5 MPa at the place 0.025 m away from the breech.The radial ignition uniformity is proved.The temperature of the gas increases from several thousand K(conventional ignition)to several ten thousand K(plasma ignition).Compare the distribution of the density and temperature of the gas,we know that low density and high temperature gas appears near the exits of the DRPG,while high density and low temperature gas appears at the wall near the breech.The simulation of the 2D-IB-SPETCG model is an effective way to investigate the interior ballistics process of the ETC launch.The 2D-IB-SPETC model can be used for prediction and improvement of experiments.

Device simulation of quasi-two-dimensional perovskite/silicon tandem solar cells towards 30%-efficiency

Perovskite/silicon(Si) tandem solar cells have been recognized as the next-generation photovoltaic technology with efficiency over 30% and low cost. However, the intrinsic instability of traditional three-dimensional(3D) hybrid perovskite seriously hinders the lifetimes of tandem devices. In this work, the quasi-two-dimensional(2D)(BA)_(2)(MA)_(n-1)Pbn I_(3n+1)(n = 1, 2, 3, 4, 5)(where MA denotes methylammonium and BA represents butylammonium), with senior stability and wider bandgap, are first used as an absorber of semitransparent top perovskite solar cells(PSCs) to construct a fourterminal(4T) tandem devices with a bottom Si-heterojunction cell. The device model is established by Silvaco Atlas based on experimental parameters. Simulation results show that in the optimized tandem device, the top cell(n = 4) obtains a power conversion efficiency(PCE) of 17.39% and the Si bottom cell shows a PCE of 11.44%, thus an overall PCE of 28.83%. Furthermore, by introducing a 90-nm lithium fluoride(LiF) anti-reflection layer to reduce the surface reflection loss, the current density(J_(sc)) of the top cell is enhanced from 15.56 m A/cm^(2) to 17.09 m A/cm^(2), the corresponding PCE reaches 19.05%, and the tandem PCE increases to 30.58%. Simultaneously, in the cases of n = 3, 4, and 5, all the tandem PCEs exceed the limiting theoretical efficiency of Si cells. Therefore, the 4T quasi-2D perovskite/Si devices provide a more cost-effective tandem strategy and long-term stability solutions.

Two-Dimensional Fluid Simulation of Collisional Plasma Sheath over rf Powered Electrode with Cylindrical Hole

The characteristics of collisional radio-frequency (rf) sheath dynamics over an elec-trode with a cylindrical hole is simulated by means of a self-consistent model which consists of two-dimensional time-dependent fluid equations coupled with Poisson equation. In addition, an equivalent-circuit model is coupled to the fluid equations in order to self-consistently determine re-lationship between the instantaneous potential at the rf-biased electrode and the sheath thickness. Two-dimensional profiles of the potential, the ion fluid velocity, and the distributions of the ion and electron densities within the sheath are computed under various discharge conditions, such as the discharge powers and the gas pressures. The results show that the existence of the cylindrical hole on the electrode significantly affects the sheath structure and generates a potential trap in the horizontal direction, which is particularly strong when the sheath thickness is comparable to the depth of the hole. Moreover, it is found that the collisional effects have a significant influence on the sheath characteristics.

Two-Dimensional Simulation of a Dual Frequency Sheath Near an Electrode with a Cylindrical Hole

The characteristics of a collisional dual frequency （DF） sheath near an electrode with a cylindrical hole are studied by utilizing a two-dimensional model which includes time-dependent fluid equations coupled with the Poisson equation and an equivalent-circuit model, The effects of the gas pressure on the two-dimensional profiles of the potential, electric field, ion fluid velocity in a DF sheath are investigated. The simulation results show that the cylindrical hole on the electrode has a significant influence on the DF sheath structure, i.e., the sheath profile tends to wrap around the contour of the hole feature. Moreover, it is shown that the structure of the DF sheath is different from that of a single frequency （SF） sheath because the profile of the DF sheath is modulated by the combination of the high and low frequency sources. In addition the characteristics of the DF sheath are obviously affected by the collisional effects in the DF sheath.

Numerical simulation for the two-dimensional nonlinear shallow water waves

This study deals with the general numerical model to simulate the two-dimensional tidal flow, flooding wave (long wave) and shallow water waves (short wave). The foundational model is based on nonlinear Boussinesq equations. Numerical method for modelling the short waves is investigated in detail. The forces, such as Coriolis forces, wind stress, atmosphere and bottom friction, are considered. A two-dimensional implicit difference scheme of Boussinesq equations is proposed. The low-reflection outflow open boundary is suggested. By means of this model,both velocity fields of circulation current in a channel with step expansion and the wave diffraction behind a semi-infinite breakwater are computed, and the results are satisfactory.

Effect of Parallel-Plate Geometry on Mode Transition Behavior in Argon Microplasmas: Two-Dimensional Simulation

A two-dimensional self-consistent fluid model is employed to investigate radio-frequency process parameters on the plasma properties in Ar microdischarges. The neutral gas density and temperature balance equations are taken into account. We mainly investigate the effect of the electrode gap on the spatial distribution of the electron density and electron temperature profiles, due to a mode transition from the regime（secondary electrons emission is responsible for the significant ionization） to the regime（sheath oscillations and bulk electrons are responsible for sustaining discharge） induced by a sudden decrease of electron density and electron temperature.The pressure, radio-frequency sources frequency and voltage effects on the electron density are also elaborately investigated.

Forward Dynamics Analysis of the 6-PUS Mechanism Based on Platform-Legs Composite Simulation

The dynamics analysis plays an important role for the control, simulation and optimization of the parallel manipulators. Normally, the Stewart type manipulators have a platform and several legs. The inverse dynamics can be solved efficiently by taking the advantage of such structural characteristics. However, for the forward dynamics analysis, this structural decomposition still faces challenges from both modeling and computation. In this paper, an efficient approach is proposed for the forward dynamics of the 6-PUS manipulator based on the platform-legs composite simulation. By composite method, the dynamics modeling of the parallel manipulator is separated into the forward dynamics of the platform and the kineto-statics of the legs. The global simulation model can be constructed by connecting the predefined platform model and leg models according to the manipulator＇s topology. Thus, the global simulation can be decomposed into the independent calculations of purely algebraic equations and ordinary differential equations （ODEs）, the computational cost can be reduced and the stability of the simulation can be improved. For the purpose of solving the manipulator＇s forward dynamics accurately, the algebraic-loop problem is discussed and a closed form algorithm is proposed. A numerical example of the 6-PUS manipulator is given to demonstrate the effectiveness of the proposed approach. The example results show that the modeling efficiency can be improved and the simulation stability can be ensured for decomposing the system equations into purely algebraic equations and ODEs.

An agent-based simulation with NetLogo platform to evaluate forward osmosis process(PRO Mode)

Forward osmosis(FO), as an emerging technology, is influenced by different factors such as operating conditions,module characteristics, and membrane properties. The general aim of this study was to develop a suitable(flexible,comprehensive, and convenient to use) computational tool which is able to simulate osmosis through an asymmetric membrane oriented in pressure retarded osmosis(PRO) mode in a wide variety of scenarios. For this purpose, an agent-based model was created in NetLogo platform, which is an easy-to-use application environment with graphical visualization abilities and well suited for modeling a complex system evolving over time. The simulation results were validated with empirical data obtained from literature and a great agreement was observed. The effect of various parameters on process performance was investigated in terms of temperature,cross-flow velocity, length of the module, pure water permeability coefficient, and structural parameter of the membrane. Results demonstrated that the increase in all parameters, except structural parameter of the membrane and the length of module led to the increase of average water flux. Moreover, nine different draw solutes were selected in order to assess the influence of net bulk osmotic pressure difference between the draw solution(DS) and feed solution(FS)(known as the driving force of FO process) on water flux. Based on the findings of this paper, the performance of FO process(PRO mode) can be efficiently evaluated using the NetL ogo platform.

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.

Numerical simulation and response analysis of microspherical focused logging in inclined micro-fractured formation

Tight and unconventional reservoirs have become the focus with the progress of petroleum exploration and development.Micro-fractures in these reservoirs can effectively improve reservoir permeability,and well-developed micro-fractures can serve to directly improve productivity.Compared with the centered electrical well logging method,the Micro Spherical Focused Logging(MSFL)is more suitable for microfracture identification due to its high resolution and near borehole wall measuring method.In this study,an anisotropic model is used to depict micro-fractured formation.First,a forward model with microfractured formation,borehole,logging instrument and surrounding rock is established.Subsequently,MSFL responses under different micro-fracture porosity,resistivity,dip angle and borehole radius,are calculated based on the finite element method(FEM).Finally,the MSFL responses under different microfracture parameters are analyzed with the response laws clarified.

3D seismic forward modeling from the multiphysical inversion at the Ketzin CO_(2) storage site

From June 2008 to August 2013,approximately 67 kt of CO_(2) was injected into a deep saline formation at the Ketzin pilot CO_(2) storage site.During injection,3D seismic surveys have been performed to monitor the migration of sequestered CO_(2).Seismic monitoring results are limited by the acquisition and signal-to-noise ratio of the acquired data.The multiphysical reservoir simulation provides information regarding the CO_(2) fluid behavior,and the approximated model should be calibrated with the monitoring results.In this work,property models are delivered from the multiphysical model during 3D repeated seismic surveys.The simulated seismic data based on the models are compared with the real data,and the results validate the effectiveness of the multiphysical inversion method.Time-lapse analysis shows the trend of CO_(2) migration during and after injection.

Different human-simulated characteristics of the driver models in the forward simulation of hybrid electric vehicles

Based on the control theories of PID, fuzzy logic and expert PID, the driver models are built and applied in the forward simulation for hybrid electric vehicles （HEV）. The impact to the vehicle speed tracking and the fuel economy is compared among the different driver models. The different human-simulated characteristics of the driver models are emphatically analyzed. The analysis results indicate that the driver models based on PID, simple fuzzy logic and expert PID are corresponding to the handling characteristics of different drives. The driver models of different human-simulated characteristics bring the handling divergence of drivers with different driving level and habit to the HEV forward simulation, and that is significant to the all-around verification and validation of the control strategy for HEV. System simulation results of different driver models validate the impact of driver models to the dynamic and fuel economy performance of HEV.

Fast forward modeling of gravity anomalies for two-dimensional bodies of arbitrary shape and density distribution

A fast and high precision spatial domain algorithm is presented for forward modeling of two-dimensional(2D)body gravity anomalies of arbitrary shape and density distribution.The new algorithm takes advantage of the convolution properties of the expression for 2D gravity anomalies,uses a rectangular cell as a grid subdivision unit,and then 2D bodies with irregular cross-sections are approximated by a combination of 2D bodies with a rectangular cross section.The closed-form expression is used to calculate the gravitational anomalies of the combination of 2D bodies with a rectangular cross section.To improve computing effi ciency,the new algorithm uses a fast algorithm for the implementation of the Toeplitz matrix and vector multiplication.The synthetic 2D models with rectangular and circular cross-sections and constant and variable densities are designed to evaluate the computational accuracy and speed of the new algorithm.The experiment results show that the computation costs less than 6 s for a grid subdivision with 10000×10000 elements.Compared to the traditional forward modeling methods,the proposed method significantly improved computational effi ciency while guaranteeing computational accuracy.

Forward simulation of array laterolog resistivity in anisotropic formation and its application

A scaling-down experiment system of array laterolog resistivity was developed, and a corresponding formation model was built by 3 D finite element numerical method to study the effect of different factors on the logging response quantitatively. The error between the experimental and numerical results was less than 5%, validating the reliability of the numerical simulation method. The single factor analysis of the formation relative dip, resistivity anisotropy and drilling fluid invasion was carried out by numerical simulation method, and the results show that:(1) The increase of relative dip can lead to the increase of formation resistivity, but the increasing value is relatively small, and the values of five array resistivity curves will reverse when the relative dip angle reaches a certain degree.(2) The increase of anisotropy coefficient λ can also cause the formation resistivity to rise, and the resistivity will increase by about 10% when λ increases from 1.0 to 1.5 in vertical wells.(3) Drilling fluid invasion has a more significant effect on the logging response than the former two factors. The order of the five curves will change due to drilling fluid invasion in anisotropic formation and the change rule is contrary to resistivity anisotropy. Taking the logging data of the Yingxi oilfield in the Qaidam Basin as an example, an anisotropic formation model considering drilling fluid invasion was built, and the numerical simulation results from the above methods were basically consistent with the logging data, which verified the accuracy of the method again. The results of this study lay a theoretical foundation for multiple-parameter inversion in anisotropic formation under complex well conditions.

Dominant pulse simulation of near fault ground motions

In this study, a new mathematical model is developed composed of two parts, including harmonic and polynomial expressions for simulating the dominant velocity pulse of near fault ground motions. Based on a proposed velocity function, the corresponding expressions for the ground acceleration and displacement time histories are also derived. The proposed model is then fitted using some selected pulse-like near fault ground motions in the Next Generation Attenuation （NGA） project library. The new model is not only simple in form but also simulates the long-period portion of actual velocity near fault records with a high level of precision. It is shown that the proposed model-based elastic response spectra are compatible with the near fault records in the neighborhood of the prevailing frequency of the pulse. The results indicate that the proposed model adequately simulates the components of the time histories. Finally, the energy of the proposed pulse was compared with the energy of the actual record to confirm the compatibility.

Numerical Simulation of Viscous Flow Around a Rolling Cylinder with Ship-Like Section

Two-dimensional unsteady incompressible viscous flow around a rolling cylinder with ship-like section is numerically simulated by employing the computational scheme previously developed by the authors, in which the continuity and momentum equations are satisfied simultaneously at each time step for oscillating flow. The numerical results show that the motion of vortices around a rolling ship hull is cyclical. It is found that the location of the vortices is very similar to the existing experimental result. Using these simulation results, we can calculate the roll damping of ships including viscous effects.

Diagnostic problems in two-dimensional shear wave elastography of the liver

Two-dimensional shear wave elastography(2D-SWE)is used in the clinical setting for observation of the liver.Unfortunately,a wide spectrum of artifactual images are frequently encountered in 2D-SWE,the precise mechanisms of which remain incompletely understood.This review was designed to present many of the artifactual images seen in 2D-SWE of the liver and to analyze them by computer simulation models that support clinical observations.Our computer simulations yielded the following suggestions:(1)When performing 2D-SWE in patients with chronic hepatic disease,especially liver cirrhosis,it is recommended to measure shear wave values through the least irregular hepatic surface;(2)The most useful 2D-SWE in patients with focal lesion will detect lesions that are poorly visible on B-mode ultrasound and will differentiate true tumors from pseudo-tumors(e.g.,irregular fatty change);and(3)Measurement of shear wave values in the area posterior to a focal lesion must be avoided.

Using TOUGH2 numerical simulation to analyse the geothermal formation in Guide basin,China

The Guide sedimentary basin is located in the northeastern part of Qinghai-Xizang Plateau,which is rich in geothermal resources.However,exploitation of the geothermal resources has so far been limited,because of limited understanding of the resources quantity and storage gained from scientific researches.In this study,using a typical cross section across the basin and taking into account its geothermal and geological conditions,a new waterheat coupled model was built and associated modelling was done by the software TOUGH2.During modelling process,the accuracy and applicability of the model was confirmed through the calibration of relevant parameters for modelling the heat and water transport and the formation of geothermal reservoir across the basin,with particular focus on the Neogene geothermal field.Results show that the groundwater that flows from the basin margins to the center is heated by the Neogene and Paleogene sedimentary rocks with high geothermal gradients.Since the east-west extending fault F1 is conductive,it acts as preferential flow paths which on one hand provide additional and rapid flows to the thermal reservoir;and on the other hand,cool down the thermal water to a certain extent due to the infiltration of shallower water sources in the vicinity of the fault.Furthermore,the estimated geothermal resources quantity is close to that of previous studies.In comparison with the Paleogene rock formations,the Neogene geothermal reservoir shows a better nature in terms of water content,aquifer permeability and resources exploitability,although the resource quantity of the Paleogene reservoir is considerable.