Wave propagation across fluid-solid interfaces with LBM-LSM coupling schemes

Seismic wave propagation in fluid-solid coupled media is currently a popular topic. However, traditional wave equation-based simulation methods have to consider complex boundary conditions at the fluid-solid interface. To address this challenge, we propose a novel numerical scheme that integrates the lattice Boltzmann method(LBM) and lattice spring model(LSM). In this scheme, LBM simulates viscoacoustic wave propagation in the fluid area and LSM simulates elastic wave propagation in the solid area. We also introduce three different LBM-LSM coupling strategies, a standard bounce back scheme, a specular reflection scheme, and a hybrid scheme, to describe wave propagation across fluid-solid boundaries. To demonstrate the accuracy of these LBM-LSM coupling schemes, we simulate wave propagation in a two-layer model containing a fluid-solid interface. We place excitation sources in the fluid layer and the solid layer respectively, to observe the wave phenomena when seismic waves propagate to interface from different sides. The simulated results by LBM-LSM are compared with the reference wavefields obtained by the finite difference method(FDM) and the analytical solution(ANA).Our LBM-LSM coupling scheme was verified effective, as the relative errors between the LBM-LSM solutions and reference solutions were within an acceptable range, sometimes around 1.00%. The coupled LBM-LSM scheme is further used to model seismic wavefields across a more realistic rugged seabed,which reveals the potential applications of the coupled LBM-LSM scheme in marine seismic imaging techniques, such as reverse-time migration and full-waveform inversion. The method also has potential applications in simulating wave propagation in complex two-and multi-phase media.

Application of the generalized quasi-complementary energy principle to the fluid-solid coupling problem

The fluid-solid coupling theory, an interdisciplinary science between hydrodynamics and solid mechanics, is an important tool for response analysis and direct design of structures in naval architecture and ocean engineering. By applying the corresponding relations between generalized forces and generalized displacements, convolutions were performed between the basic equations of elasto-dynamics in the primary space and corresponding virtual quantities. The results were integrated and then added algebraically. In light of the fact that body forces and surface forces are both follower forces, the generalized quasi-complementary energy principle with two kinds of variables for an initial value problem is established in non-conservative systems. Using the generalized quasi-complementary energy principle to deal with the fluid-solid coupling problem and to analyze the dynamic response of structures, a method for using two kinds of variables simultaneously for calculation of force and displacement was derived.

Fluid-solid coupling numerical simulation of charge process in variable-mass thermodynamic system

Abstract： A joint solution model of variabk：-mass flow in two-phase region and fluid-solid coupling heat transfer, concerned about the charge process of variable-mass thermodynamic system, is built up and calculated by the finite element method （FEM）. The results are basically consistent with relative experimental data. The calculated average heat transfer coefficient reaches 1.7~105 W/（m2. K）. When the equal percentage valve is used, the system needs the minimum requirements of valve control, but brings the highest construction cost. With the： decrease of initial steam pressure, the heat transfer intensity also weakens but the steam flow increases. With the initial water filling coefficient increasing or the temperature of steam supply decreasing, the amount of accumulative steam flow increases with the growth of steam pressure. When the pressure of steam supply drops, the steam flow gradient increases during the maximum opening period of control valve, and causes the maximum steam flow to increase.

Fully fluid-solid coupling dynamic model for seismic response of underground structures in saturated soils

The seismic response characteristics of underground structures in saturated soils are investigated.A fully fluid-solid coupling dynamic model is developed and implemented into ABAQUS with a user-defined element to simulate the dynamic behavior of saturated soils.The accuracy of the model is validated using a classic example in literature.The performance of the model is verified by its application on simulating the seismic response characteristics of a subway station built in saturated soils.The merits of the model are demonstrated by comparing the difference of the seismic response of an underground structure in saturated soils between using the fully coupling model and a single-phase medium model.The study finds that the fully coupling model developed herein can simulate the dynamic response characteristics of the underground structures in saturated soils with high accuracy.The seismic response of the underground structure tends to be underestimated by using the single-phase medium model compared with using the fully coupling model,which provides a weaker confining action to the underground structure.

FLUID-SOLID COUPLING MATHEMATICAL MODEL OF CONTAMINANT TRANSPORT IN UNSATURATED ZONE AND ITS ASYMPTOTICAL SOLUTION

The process of contaminant transport is a problem of multicomponent and multiphase flow in unsaturated zone. Under the presupposition that gas existence affects water transport, a coupled mathematical model of contaminant transport in unsaturated zone has been established based on fluid_solid interaction mechanics theory. The asymptotical solutions to the nonlinear coupling mathematical model were accomplished by the perturbation and integral transformation method. The distribution law of pore pressure, pore water velocity and contaminant concentration in unsaturated zone has been presented under the conditions of with coupling and without coupling gas phase. An example problem was used to provide a quantitative verification and validation of the model. The asymptotical solution was compared with Faust model solution. The comparison results show reasonable agreement between asymptotical solution and Faust solution, and the gas effect and media deformation has a large impact on the contaminant transport. The theoretical basis is provided for forecasting contaminant transport and the determination of the relationship among pressure_saturation_permeability in laboratory.

Algorithm Selection Method Based on Coupling Strength for Partitioned Analysis of Structure-Piezoelectric-Circuit Coupling

In this study, we propose an algorithm selection method based on coupling strength for the partitioned analysis ofstructure-piezoelectric-circuit coupling, which includes two types of coupling or inverse and direct piezoelectriccoupling and direct piezoelectric and circuit coupling. In the proposed method, implicit and explicit formulationsare used for strong and weak coupling, respectively. Three feasible partitioned algorithms are generated, namely(1) a strongly coupled algorithm that uses a fully implicit formulation for both types of coupling, (2) a weaklycoupled algorithm that uses a fully explicit formulation for both types of coupling, and (3) a partially stronglycoupled and partially weakly coupled algorithm that uses an implicit formulation and an explicit formulation forthe two types of coupling, respectively.Numerical examples using a piezoelectric energy harvester,which is a typicalstructure-piezoelectric-circuit coupling problem, demonstrate that the proposed method selects the most costeffectivealgorithm.

Fluid-solid coupling model for studying wellbore instability in drilling of gas hydrate bearing sediments

As the oil or gas exploration and development activities in deep and ultra- deep waters become more and more, encountering gas hydrate bearing sediments （HBS） is almost inevitable. The variation in temperature and pressure can destabilize gas hydrate in nearby formation around the borehole, which may reduce the strength of the formation and result in wellbore instability. A non-isothermal, transient, two-phase, and fluid-solid coupling mathematical model is proposed to simulate the complex stability performance of a wellbore drilled in HBS. In the model, the phase transition of hydrate dissociation, the heat exchange between drilling fluid and formation, the change of mechanical and petrophysical properties, the gas-water two-phase seepage, and its interaction with rock deformation are considered. A finite element simulator is developed, and the impact of drilling mud on wellbore instability in HBS is simulated. Results indicate that the re- duction in pressure and the increase in temperature of the drilling fluid can accelerate hydrate decomposition and lead to mechanical properties getting worse tremendously. The cohesion decreases by 25% when the hydrate totally dissociates in HBS. This easily causes the wellbore instability accordingly. In the first two hours after the formation is drilled, the regions of hydrate dissociation and wellbore instability extend quickly. Then, with the soaking time of drilling fluid increasing, the regions enlarge little. Choosing the low temperature drilling fluid and increasing the drilling mud pressure appropriately can benefit the wellbore stability of HBS. The established model turns out to be an efficient tool in numerical studies of the hydrate dissociation behavior and wellbore stability of HBS.

Numerical simulation on fault water-inrush based on fluid-solid coupling theory

About 75% water-inrush accidents in China are caused by geological structure such as faults, therefore, it is necessary to investigate the water-inrush mechanism of faults to provide references for the mining activity above confined water. In this paper, based on the fluid-solid coupling theory, we built the stress-seepage coupling model for rock, then we combined with an example of water-inrush caused by fault, studied the water-inrush mechanism by using the numerical software COMSOL Mutiphysics, analyzed the change rule of shear stress, vertical stress, plastic area and water pressure for stope with a fault, and estimated the water-inrush risk at the different distances between working faces and the fault. The numerical simula- tion results indicate that： （1） the water-inrush risk will grow as the decrease of the distance between working face and the fault; （2） the failure mode of the rock in floor with fault is shear failure; （3） the rock between water-containing fault and working face failure is the reason for water-inrush.

Simulation Analysis of Torsion Beam Hydroforming Based on the Fluid-Solid Coupling Method

Hydroformed parts are widely used in industrial automotive parts because of their higher stiffness and fatigue strength and reduced weight relative to their corresponding cast and welded parts.This paper reports a hydraulicforming experimental platform for rectangular tube fittings that was constructed to conduct an experiment on the hydraulic forming of rectangular tube fittings.A finite element model was established on the basis of the fluid–solid coupling method and simulation analysis.The correctness of the simulation analysis and the feasibility of the fluid–solid coupling method for hydraulic forming simulation analysis were verified by comparing the experimental results with the simulation results.On the basis of the simulation analysis of the hydraulic process of the torsion beam using the fluid–solid coupling method,a sliding mold suitable for the hydroforming of torsion beams was designed for its structural characteristics.The effects of fluid characteristics,shaping pressure,axial feed rate,and friction coefficient on the wall thicknesses of torsions beams during formation were investigated.Fluid movement speed was related to tube deformation.Shaping pressure had a significant effect on rounded corners and straight edges.The axial feed speed was increased,and the uneven distribution of wall thicknesses was effectively improved.Although the friction coefficient had a nonsignificant effect on the wall thickness of the ladder-shaped region,it had a significant influence on a large deformation of wall thickness in the V-shaped area.In this paper,a method of fluid-solid coupling simulation analysis and sliding die is proposed to study the high pressure forming law in torsion beam.

Discrete modulational instability and bright localized spin wave modes in easy-axis weak ferromagnetic spin chains involving the next-nearest-neighbor coupling

We report a theoretical work on the properties of modulational instability and bright type nonlinear localized modes in one-dimensional easy-axis weak ferromagnetic spin lattices involving next-nearest-neighbor couplings.With a linear stability analysis, we calculate the growth rates of the modulational instability, and plot the instability regions.When the strength of the next-nearest-neighbor coupling is large enough, two new asymmetric modulational instability regions appear near the boundary of the first Brillouin zone.Furthermore, analytical forms of the bright nonlinear localized modes are constructed by means of a quasi-discreteness approach.The influence of the next-nearest-neighbor coupling on the Brillouin zone center mode and boundary mode are discussed.In particular, we discover a reversal phenomenon of the propagation direction of the Brillouin zone boundary mode.

Control schemes for synchronizing two subnetworks with weak couplings

In this paper, we propose a well-designed network model with a parameter and study full and partial synchronization of the network model based on the stability analysis. The network model is composed of a star-coupled subnetwork and a globally coupled subnetwork. By analyzing the special coupling configuration, three control schemes are obtained for synchronizing the network model. Further analysis indicates that even if the inner couplings in each subnetwork are very weak, two of the control schemes are still valid. In particular, if the outer coupling weight parameter 0 is larger than （n2 - 2n）/4, or the subnetwork size n is larger than 02, the two subnetworks with weak inner couplings can achieve synchronization. In addition, the synchronizability is independent of the network size in case of 0 〈 0 〈 n/（n ＋ 1 ）. Finally, we carry out some numerical simulations to confirm the validity of the obtained control schemes. It is worth noting that the main idea of this paper also applies to any network consisting of a dense subnetwork and a sparse network.

Splitting of One-Phonon Coupling States in the Weakly Coupling Model

The splitting of one-phonon coupling states in the weakly coupling model for odd-A nuclei is studied. Matrix elements are given explicitly, and an analytical expression for the energy levels is given in a perturbation approach. An application is given to ^(93)Nb and ^(51)V. A comparison between the weakly coupling model and the shell model treatment for ^(51)V is discussed.

Numerical simulation of aluminum holding furnace with fluid-solid coupled heat transfer

To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mathematics models of aluminum holding furnace in the premixed combustion processing were established based on mass conservation,moment conservation,momentum conservation,energy conservation and chemistry species conservation.Computational results agree well with the test data of the typical condition.The maximum combustion temperature is 1 850 K.The average temperature of the molten aluminum is 1 158 K,and the maximum temperature difference is about 240 K.The average temperature increases 0.3 ℃ while the temperature of combustion air increases 1 ℃.The optimal excess air ratio is 1.25-1.30.

Numerical investigation on permeability evolution behavior of rock by an improved flow-coupling algorithm in particle flow code

Permeability is a vital property of rock mass, which is highly affected by tectonic stress and human engineering activities. A comprehensive monitoring of pore pressure and flow rate distributions inside the rock mass is very important to elucidate the permeability evolution mechanisms, which is difficult to realize in laboratory, but easy to be achieved in numerical simulations. Therefore, the particle flow code （PFC）, a discrete element method, is used to simulate permeability behaviors of rock materials in this study. Owe to the limitation of the existed solid-fluid coupling algorithm in PFC, an improved flow-coupling algorithm is presented to better reflect the preferential flow in rock fractures. The comparative analysis is conducted between original and improved algorithm when simulating rock permeability evolution during triaxial compression, showing that the improved algorithm can better describe the experimental phenomenon. Furthermore, the evolution of pore pressure and flow rate distribution during the flow process are analyzed by using the improved algorithm. It is concluded that during the steady flow process in the fractured specimen, the pore pressure and flow rate both prefer transmitting through the fractures rather than rock matrix. Based on the results, fractures are divided into the following three types： I） fractures link to both the inlet and outlet, II） fractures only link to the inlet, and III） fractures only link to the outlet. The type I fracture is always the preferential propagating path for both the pore pressure and flow rate. For type II fractures, the pore pressure increases and then becomes steady. However, the flow rate increases first and begins to decrease after the flow reaches the stop end of the fracture and finally vanishes. There is no obvious pore pressure or flow rate concentration within type III fractures.

Fluid solid coupling model based on endochronic damage for roller compacted concrete dam

According to the characteristics of thin-layer rolling and pouting construction technology and the complicated mechanical behavior of the roller compacted concrete dam （RCCD） construction interface, a constitutive model of endochronic damage was established based on the endochronic theory and damage mechanics. The proposed model abandons the traditional concept of elastic-plastic yield surface and can better reflect the real behavior of rolled control concrete. Basic equations were proposed for the fluid-solid coupling analysis, and the relationships among the corresponding key physical parameters were also put forward. One three-dimensional finite element method （FEM） program was obtained by studying the FEM type of the seepage-stress coupling intersection of the RCCD. The method was applied to an actual project, and the results show that the fluid-solid interaction influences dam deformation and dam abutment stability, which is in accordance with practice. Therefore, this model provides a new method for revealing the mechanical behavior of RCCD under the coupling field.

Flow Field and Temperature Field of Water-Cooling-Type Magnetic Coupling

At present, the water-cooling simulation of the water-cooled magnetic coupler is based on the water-cooled motor and the hydraulic coupler, which cannot accurately characterize the temperature distribution of the rotating watercooled coupling of the coupler. Focusing on rotating water cooling radiating, the present paper proposes simulating the water cooling temperature field as well as the flow field through the method of combining fluid-solid coupled heat transfer and MRF(Multiphase Reference Frame). In addition, taking an 800 kW magnetic coupling as an example, the paper optimizes the shape, number, cooling water inlet speed? and so on? of the cooling channel. Considering factors such as the complete machine’s temperature, and drag torque, it is proved that the cooling e ect is best when there are 36 involute curved channels and when the inlet speed is 3 m/s. Further, through experiments, the actual temperature values at six di erent positions when 50 kW and 70 kW thermal losses di er are measured. The measured values agree with the simulation results, proving the correctness of the proposed method. Further, data have been collected during the entire experimental procedure? and the variation in the coupling’s temperature is analyzed in depth, with the objective of laying a foundation for the estimation of the inner temperature rise as well as for the optimization of the structural design.

Nonlinear coupling analysis of coal seam floor during mining based on FLAC^3D

Based on the hydro-geological conditions of 1028 mining face in Suntuan Coal Mine, mining seepage strain mechanism of seam floor was simulated by a nonlinear coupling method, which applied fluid-solid coupling analysis module of FLAC^3D. The results indicate that the permeability coefficient of adjoining rock changes a lot due to mining. The maximum value reaches 1 379.9 times to the original value, where it is at immediate roof of the mined-out area. According to the analysis on the seepage field, mining does not destroy water resistance of the floor aquiclude. The mining fissure does not conduct lime-stone aquifer, and it is less likely to form damage. The plastic zone does not exactly correspond to the seepage area, and the scope of the altered seepage area is much larger than the plastic zone.

Importance of Three-Dimensional Piezoelectric Coupling Modeling in Quantitative Analysis of Piezoelectric Actuators

This paper demonstrates the importance of three-dimensional(3-D)piezoelectric coupling in the electromechanical behavior of piezoelectric devices using three-dimensional finite element analyses based on weak and strong coupling models for a thin cantilevered piezoelectric bimorph actuator.It is found that there is a significant difference between the strong and weak coupling solutions given by coupling direct and inverse piezoelectric effects(i.e.,piezoelectric coupling effect).In addition,there is significant longitudinal bending caused by the constraint of the inverse piezoelectric effect in the width direction at the fixed end(i.e.,3-D effect).Hence,modeling of these effects or 3-D piezoelectric coupling modeling is an electromechanical basis for the piezoelectric devices,which contributes to the accurate prediction of their behavior.

Fully coupled fluid-solid productivity numerical simulation of multistage fractured horizontal well in tight oil reservoirs

A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite element method was used to determine the numerical solution and the accuracy of the model was verified. On this basis, the model was used to simulate productivity of multistage fractured horizontal wells in tight oil reservoirs. The results show that during the production of tight oil wells, the reservoir region close to artificial fractures deteriorated in physical properties significantly, e.g. the aperture and conductivity of artificial fractures dropped by 52.12% and 89.02% respectively. The simulations of 3000-day production of a horizontal well in tight oil reservoir showed that the predicted productivity by the uncoupled model had an error of 38.30% from that by the fully-coupled model. Apparently, ignoring the influence of fluid-solid interaction effect led to serious deviations of the productivity prediction results. The productivity of horizontal well in tight oil reservoir was most sensitive to the start-up pressure gradient, and second most sensitive to the opening of artificial fractures. Enhancing the initial conductivity of artificial fractures was helpful to improve the productivity of tight oil wells. The influence of conductivity, spacing, number and length of artificial fractures should be considered comprehensively in fracturing design. Increasing the number of artificial fractures unilaterally could not achieve the expected increase in production.

Spin-star environment assisted entanglement generation in weakly coupled bipartite systems

We study the entanglement evolution in a weakly coupled bipartite system with a large energy level difference under the influence of spin-star environments. The subsystems can be coupled to a pure state or a thermal equilibrium state spin-star environment. Our results show that, in the case of the coupling strength being less than the energy level difference of the subsystems （weakly coupled）, the spin-star environment can always be used to assist the entanglement generation of the bipartite system.