A hybrid contact approach for modeling soil-structure interaction using the material point method

The grid-based multi-velocity field technique has become increasingly popular for simulating the Material Point Method(MPM)in contact problems.However,this traditional technique has some shortcomings,such as(1)early contact and contact penetration can occur when the contact conditions are unsuitable,and(2)the method is not available for contact problems involving rigid-nonrigid materials,which can cause numerical instability.This study presents a new hybrid contact approach for the MPM to address these limitations to simulate the soil and structure interactions.The approach combines the advantages of point-point and point-segment contacts to implement contact detection,satisfying the impenetrability condition and smoothing the corner contact problem.The proposed approach is first validated through a disk test on an inclined slope.Then,several typical cases,such as granular collapse,bearing capacity,and deformation of a flexible retaining wall,are simulated to demonstrate the robustness of the proposed approach compared with FEM or analytical solutions.Finally,the proposed method is used to simulate the impact of sand flow on a deformable structure.The results show that the proposed contact approach can well describe the phenomenon of soil-structure interaction problems.

Parallelization strategies for resolved simulations of fluid-structure-particle interactions

Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.

Modeling pipe-soil interaction under vertical downward relative offset using B-spline material point method

To analyze the pipeline response under permanent ground deformation,the evolution of resistance acting on the pipe during the vertical downward offset is an essential ingredient.However,the efficient simulation of pipe penetration into soil is challenging for the conventional finite element(FE)method due to the large deformation of the surrounding soils.In this study,the B-spline material point method(MPM)is employed to investigate the pipe-soil interaction during the downward movement of rigid pipes buried in medium and dense sand.To describe the density-and stress-dependent behaviors of sand,the J2-deformation type model with state-dependent dilatancy is adopted.The effectiveness of the model is demonstrated by element tests and biaxial compression tests.Afterwards,the pipe penetration process is simulated,and the numerical outcomes are compared with the physical model tests.The effects of pipe size and burial depth are investigated with an emphasis on the mobilization of the soil resistance and the failure mechanisms.The simulation results indicate that the bearing capacity formulas given in the guidelines can provide essentially reasonable estimates for the ultimate force acting on buried pipes,and the recommended value of yield displacement may be underestimated to a certain extent.

Experimental Verification of the Theory of Non-Force (Information) Interaction

The work is devoted to the demonstration of the possibility of applying the formulas of information handling obtained in the theory of non-force interaction for the natural language processing. These formulas were obtained in computer experiments in modelling the movement and interaction of material objects by changing the amount of information that triggers this movement. The hypothesis, objective and tasks of the experimental research were defined. The methods and software tools were developed to conduct the experiments. To compare different results of the simulation of the processes in a human brain during speech production, there was a range of methods proposed to calculate the estimate of sequence of fragments of natural language texts including the methods based on linear approximation. The experiments confirmed that the formulas of information handling obtained in the theory of non-force interaction reflect the processes of language formation. It is shown that the offered approach can successfully be used to create systems of reactive artificial intelligence machines. Experimental and, presented in this work, practical results constitute that the non-force (informational) interaction formulae are generally valid.

Modeling One Dimensional Two-Cell Model with Tumor Interaction Using Krylov Subspace Methods

A brain tumor occurs when abnormal cells grow, sometimes very rapidly, into an abnormal mass of tissue. The tumor can infect normal tissue, so there is an interaction between healthy and infected cell. The aim of this paper is to propose some efficient and accurate numerical methods for the computational solution of one-dimensional continuous basic models for the growth and control of brain tumors. After computing the analytical solution, we construct approximations of the solution to the problem using a standard second order finite difference method for space discretization and the Crank-Nicolson method for time discretization. Then, we investigate the convergence behavior of Conjugate gradient and generalized minimum residual as Krylov subspace methods to solve the tridiagonal toeplitz matrix system derived.

An extended multiple-support response spectrum method incorporating fluid-structure interaction for seismic analysis of deep-water bridges

The effects of ground motion spatial variability(GMSV)or fluid-structure interaction(FSI)on the seismic responses of deep-water bridges have been extensively examined.However,there are few studies on the seismic performance of bridges considering GMSV and FSI effects simultaneously.In this study,the original multiple-support response spectrum(MSRS)method is extended to consider FSI effect for seismic analysis of deep-water bridges.The solution of hydrodynamic pressure on a pier is obtained using the radiation wave theory,and the FSI-MSRS formulation is derived according to the random vibration theory.The influence of FSI effect on the related coefficients is analyzed.A five-span steel-concrete continuous beam bridge is adopted to conduct the numerical simulations.Different load conditions are designed to investigate the variation of the bridge responses when considering the GMSV and FSI effects.The results indicate that the incoherence effect and wave passage effect decrease the bridge responses with a maximum percentage of 86%,while the FSI effect increases the responses with a maximum percentage of 26%.The GMSV and FSI effects should be included in the seismic design of deep-water bridges.

Design of N-11-Azaartemisinins Potentially Active against Plasmodium falciparum by Combined Molecular Electrostatic Potential, Ligand-Receptor Interaction and Models Built with Supervised Machine Learning Methods

N-11-azaartemisinins potentially active against Plasmodium falciparum are designed by combining molecular electrostatic potential (MEP), ligand-receptor interaction, and models built with supervised machine learning methods (PCA, HCA, KNN, SIMCA, and SDA). The optimization of molecular structures was performed using the B3LYP/6-31G* approach. MEP maps and ligand-receptor interactions were used to investigate key structural features required for biological activities and likely interactions between N-11-azaartemisinins and heme, respectively. The supervised machine learning methods allowed the separation of the investigated compounds into two classes: cha and cla, with the properties εLUMO+1 (one level above lowest unoccupied molecular orbital energy), d(C6-C5) (distance between C6 and C5 atoms in ligands), and TSA (total surface area) responsible for the classification. The insights extracted from the investigation developed and the chemical intuition enabled the design of sixteen new N-11-azaartemisinins (prediction set), moreover, models built with supervised machine learning methods were applied to this prediction set. The result of this application showed twelve new promising N-11-azaartemisinins for synthesis and biological evaluation.

Strength weakening and its micromechanism in water–rock interaction,a short review in laboratory tests

Water–rock interaction(WRI)is a topic of interest in geology and geotechnical engineering.Many geological hazards and engineering safety problems are severe under the WRI.This study focuses on the water weakening of rock strength and its infuencing factors(water content,immersion time,and wetting–drying cycles).The strength of the rock mass decreases to varying degrees with water content,immersion time,and wetting–drying cycles depending on the rock mass type and mineral composition.The corresponding acoustic emission count and intensity and infrared radiation intensity also weaken accordingly.WRI enhances the plasticity of rock mass and reduces its brittleness.Various microscopic methods for studying the pore characterization and weakening mechanism of the WRI were compared and analyzed.Various methods should be adopted to study the pore evolution of WRI comprehensively.Microscopic methods are used to study the weakening mechanism of WRI.In future work,the mechanical parameters of rocks weakened under long-term water immersion(over years)should be considered,and more attention should be paid to how the laboratory scale is applied to the engineering scale.

Land-Air Interaction over Arid/Semi-arid Areas in China and Its Impact on the East Asian Summer Monsoon. Part I:Calibration of the Land Surface Model (BATS)Using Multicriteria Methods

To improve the land surface simulation in the arid and semi-arid areas of northern China, the observational data from two field experiments in Dunhuang and Tongyu are used to optimize the parameters in the land surface model, BATS, through calibration with the multicriteria method. Sensitivity analysis to the parameters in Dunhuang and Tongyu indicates that different parameters need to be calibrated in two sites with different environmental and climate regimes. Comparison of observed sensible heat flux, latent heat flux, and ground surface temperature with the simulated ones shows the simulations with the optimized parameters have been substantially improved. Especially, the holistic simulations with the calibration of the parameter values are much closer to the observations in the arid region （Dunhuang）, and the energy partition with the calibrated parameters can also be simulated well in the semi-arid region （Tongyu）. Whole results demonstrate that the parameter calibration of the land surface model is important when the model is to be used to investigate the land-air interaction.

Seismic wave input method for three-dimensional soil-structure dynamic interaction analysis based on the substructure of artificial boundaries

The method of inputting the seismic wave determines the accuracy of the simulation of soil-structure dynamic interaction. The wave method is a commonly used approach for seismic wave input, which converts the incident wave into equivalent loads on the cutoff boundaries. The wave method has high precision, but the implementation is complicated, especially for three-dimensional models. By deducing another form of equivalent input seismic loads in the fi nite element model, a new seismic wave input method is proposed. In the new method, by imposing the displacements of the free wave fi eld on the nodes of the substructure composed of elements that contain artifi cial boundaries, the equivalent input seismic loads are obtained through dynamic analysis of the substructure. Subsequently, the equivalent input seismic loads are imposed on the artifi cial boundary nodes to complete the seismic wave input and perform seismic analysis of the soil-structure dynamic interaction model. Compared with the wave method, the new method is simplifi ed by avoiding the complex processes of calculating the equivalent input seismic loads. The validity of the new method is verifi ed by the dynamic analysis numerical examples of the homogeneous and layered half space under vertical and oblique incident seismic waves.

NUMERICAL SOLUTION OF FLUID-STRUCTURE INTERACTION IN LIQUID-FILLED PIPES BY METHOD OF CHARACTERISTICS

Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenomenon of FSI due to friction coupling and Poisson coupling being taken into account, is utilized to describe the FSI of fluid-filled pipe system. Terse compatibility equations are educed by the method of characteristics （MOC） to describe the fluid-filled pipe system. To shorten computing time needed to get the solutions under the condition of keeping accuracy requirement, two steps are adopted, firstly the time step Δt and divided number of the straight pipe are optimized, sec-ondly the mesh spacing Δz close to boundary is subdivided in several submeshes automatically ac-cording to the speed gradient of fluid. The mathematical model and arithmetic are validated by com-parisons between simulation solutions of two straight pipe systems and experiment known from lit-erature.

Numerical simulation of unsteady interaction of centrifugal impeller with its diffuser using Lagrangian discrete vortex method

In this study, an advanced Lagrangian vortex- boundary element method is applied to simulate the unsteady impeller-diffuser interactions in a diffuser pump not only for design but also for off-design considerations. In velocity calculations based on the Biot-Savart law we do not have to grid large portions of the flow field and the calculation points are concentrated in the regions where vorticity is present. Lagrangian representation of the evolving vorticity field is well suited to moving boundaries. An integral pressure equation shows that the pressure distribution can be estimated directly from the instantaneous velocity and vorticity field. The numerical results are compared with the experimental data and the comparisons show that the method used in this study can provide us insight into the complicated unsteady impeller-diffuser interaction phenomena in a diffuser pump.

A SEMI-ANALYTICAL AND SEMI-NUMERICAL METHOD FOR SOLVING 2-D SOUND-STRUCTURE INTERACTION PROBLEMS

Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harmonic excitation.Within any integration segment, as long as its length is small enough,along the circumferential curvilinear coordinate,the non- homogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper.For the exterior fluid domain,the multi-circular virtual source simulation technique is adopted.The source density distributed on each virtual circular curve may be ex- panded as the Fourier's series.Combining with the inverse fast Fourier transformation,a higher accuracy and efficiency method for solving 2-D exterior Helmholtz's problems is presented in this paper.In the aspect of solution to the coupling equations,the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper.According to the superposition princi- ple and compatibility conditions at the interface between the elastic ring and fluid,the algebraic equation of system can be directly constructed by using the least square approximation.Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented.Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.

Dynamic interaction numerical models in the time domain based on the high performance scaled boundary finite element method

Consideration of structure-foundation-soil dynamic interaction is a basic requirement in the evaluation of the seismic safety of nuclear power facilities. An efficient and accurate dynamic interaction numerical model in the time domain has become an important topic of current research. In this study, the scaled boundary finite element method （SBFEM） is improved for use as an effective numerical approach with good application prospects. This method has several advantages, including dimensionality reduction, accuracy of the radial analytical solution, and unlike other boundary element methods, it does not require a fundamental solution. This study focuses on establishing a high performance scaled boundary finite element interaction analysis model in the time domain based on the acceleration unit-impulse response matrix, in which several new solution techniques, such as a dimensionless method to solve the interaction force, are applied to improve the numerical stability of the actual soil parameters and reduce the amount of calculation. Finally, the feasibility of the time domain methods are illustrated by the response of the nuclear power structure and the accuracy of the algorithms are dynamically verified by comparison with the refinement of a large-scale viscoelastic soil model.

The Boundary Element Method for Pile-Soil Interaction

In this paper, several mathmatical models for the pile- soil interaction are outlined. The Boundary Element Method is one of the very effective methods for the reasonable models of elasticity and elastoplasticity. The major of this paper is concerned with the Boundary Element Method for the pile-soil interaction, including general methods and calculating formulation of static and dynamic analysis of the pile and pile groups. Some results of analysis are also given.

OpenIFEM:A High Performance Modular Open-Source Software of the Immersed Finite Element Method for Fluid-Structure Interactions

We present a high performance modularly-built open-source software-OpenIFEM.OpenIFEM is a C++implementation of the modified immersed finite element method(mIFEM)to solve fluid-structure interaction(FSI)problems.This software is modularly built to perform multiple tasks including fluid dynamics(incompressible and slightly compressible fluid models),linear and nonlinear solid mechanics,and fully coupled fluid-structure interactions.Most of open-source software packages are restricted to certain discretization methods;some are under-tested,under-documented,and lack modularity as well as extensibility.OpenIFEM is designed and built to include a set of generic classes for users to adapt so that any fluid and solid solvers can be coupled through the FSI algorithm.In addition,the package utilizes well-developed and tested libraries.It also comes with standard test cases that serve as software and algorithm validation.The software can be built on cross-platform,i.e.,Linux,Windows,and Mac OS,using CMake.Efficient parallelization is also implemented for high-performance computing for large-sized problems.OpenIFEM is documented using Doxygen and publicly available to download on GitHub.It is expected to benefit the future development of FSI algorithms and be applied to a variety of FSI applications.

Galaxy Interactions in Filaments and Sheets:Insights from EAGLE Simulations

We study the color and star formation rates of paired galaxies in filaments and sheets using the EAGLE simulations.We find that the major pairs with pair separation<50 kpc are bluer and more star-forming in filamentary environments compared to those hosted in sheet-like environments.This trend reverses beyond a pair separation of~50 kpc.The interacting pairs with larger separations(>50 kpc)in filaments are on average redder and low-star-forming compared to those embedded in sheets.The galaxies in filaments and sheets may have different stellar mass and cold gas mass distributions.Using a KS test,we find that for paired galaxies with pair separation<50 kpc,there are no significant differences in these properties in sheets and filaments.The filaments transport gas toward the cluster of galaxies.Some earlier studies find preferential alignment of galaxy pairs with the filament axis.Such alignment of galaxy pairs may lead to different gas accretion efficiency in galaxies residing in filaments and sheets.We propose that the enhancement of star formation rate at smaller pair separation in filaments is caused by the alignment of galaxy pairs.A recent study with SDSS data reports the same findings.The confirmation of these results by the EAGLE simulations suggests that the hydrodynamical simulations are powerful theoretical tools for studying galaxy formation and evolution in the cosmic web.

ALE Fractional Step Finite Element Method for Fluid-Structure Nonlinear Interaction Problem

A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as the free surface motion, the arbitrary Lagrangian-Eulerian formulation is employed as the basis of the finite element spatial discretization. For numerical integration in time, the fraction,step method is used. This method is useful because one can use the same linear interpolation function for both velocity and pressure. The method is applied to the nonlinear interaction of a structure and a tuned liquid damper. All computations are performed with a personal computer.

Numerical method of studying nonlinear interactions between long waves and multiple short waves

Although the nonlinear interactions between a single short gravity wave and a long wave can be solved analytically, the solution is less tractable in more general cases involving multiple short waves. In this work we present a numerical method of studying nonlinear interactions between a long wave and multiple short harmonic waves in infinitely deep water. Specifically, this method is applied to the calculation of the temporal and spatial evolutions of the surface elevations in which a given long wave interacts with several short harmonic waves. Another important application of our method is to quantitatively analyse the nonlinear interactions between an arbitrary short wave train and another short wave train. From simulation results, we obtain that the mechanism for the nonlinear interactions between one short wave train and another short wave train （expressed as wave train 2） leads to the energy focusing of the other short wave train （expressed as wave train 3）. This mechanism occurs on wave components with a narrow frequency bandwidth, whose frequencies are near that of wave train 3.