A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface i...A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface immersed boundary(IB)method,which is attractive for simulating moving-boundary flows with large deformations.The adaptive mesh refinement technique is employed to reduce the computational cost while retain the desired resolution.The dynamic response of the parachute is solved with the finite element approach.The canopy and cables of the parachute system are modeled with the hyperelastic material.A tether force is introduced to impose rigidity constraints for the parachute system.The accuracy and reliability of the present framework is validated by simulating inflation of a constrained square plate.Application of the present framework on several canonical cases further demonstrates its versatility for simulation of parachute inflation.展开更多
We perform a comprehensive study of the electron-doped t-t′-J model on cylinders with density matrix renormalization group(DMRG).We conduct a systematic study on the finite-size and boundary condition effects on t-t...We perform a comprehensive study of the electron-doped t-t′-J model on cylinders with density matrix renormalization group(DMRG).We conduct a systematic study on the finite-size and boundary condition effects on t-t′-J model on cylinders.Periodic and anti-periodic boundary conditions are implemented along the circumference direction,with the system’s width extending up to as large as 8 lattice units.We study doping levels of 1/6,1/8,and 1/12,which represent the most interesting region in the phase diagram of electron-doped cuprates.We find that for width-4 and width-6 systems,the ground state for fixed doping switches between anti-ferromagnetic Neel state and stripe state under different boundary conditions and system widths,indicating the presence of large finite size effect in the t-t′-J model.We also have a careful analysis of the d-wave pairing correlations which also change quantitatively with boundary conditions and widths of the system.However,the pairing correlations are enhanced when the system becomes wider for all dopings,suggesting the existence of possible long-range superconducting order in the thermodynamic limit.The width-8 results are found to be dependent on the starting state in the DMRG calculation for the kept states we can reach.For the width-8 system,only Neel(stripe)state can be stabilized in DMRG calculation for 1/12(1/6)doping,while both stripe and Neel states are stable in the DMRG sweep for 1/8 doping,regardless of the boundary conditions.These results indicate that 1/8 doping is likely to lie on the boundary of a phase transition between the Neel phase with lower doping and the stripe phase with higher doping,consistent with the previous study.The sensitivity of the ground state on boundary conditions and size observed for narrow systems is similar to that found in the t′-Hubbard model,where the t′term introduces frustration and makes the stripe state fragile.The study of different boundary conditions provides a useful tool to check the finite size effect in the future DMRG calculations.展开更多
The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analy...The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analytical one. And then examples are given to demonstrate the solution of several problems of the irregular object scattering. The method can not only produce clear physical pictures, but can efficiently handle many complicated scattering problems.展开更多
The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-eleme...The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.展开更多
By establishing equivalent fixed point theorem, the boundary value problems of p Laplace equations with finite time delay are studied. It’s the first time that the functional differential equation is discussed w...By establishing equivalent fixed point theorem, the boundary value problems of p Laplace equations with finite time delay are studied. It’s the first time that the functional differential equation is discussed with p Laplacian. The topological degree and fixed point theorem on cone are used to prove the existence of solution and positive solution. The conditions are all easy to check.展开更多
In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered gri...In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition （ABC） to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.展开更多
An implicit finite difference method is developed for a one-dimensional frac- tional percolation equation （FPE） with the Dirichlet and fractional boundary conditions. The stability and convergence are discussed for ...An implicit finite difference method is developed for a one-dimensional frac- tional percolation equation （FPE） with the Dirichlet and fractional boundary conditions. The stability and convergence are discussed for two special cases, i.e., a continued seep- age flow with a monotone percolation coefficient and a seepage flow with the fractional Neumann boundary condition. The accuracy and efficiency of the method are checked with two numerical examples.展开更多
In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference me...In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.展开更多
By using cone theory and the MSnch fixed theorem combined with a monotone iterative technique, we investigate the existence of positive solutions for systems of second- order nonlinear singular differential equations ...By using cone theory and the MSnch fixed theorem combined with a monotone iterative technique, we investigate the existence of positive solutions for systems of second- order nonlinear singular differential equations with integral boundary conditions on infinite interval and establish the existence theorem of positive solutions and iterative sequence for approximating the positive solutions. The results in this paper improve some known results.展开更多
An increment-dimensional scaled boundary finite element method (ID-SBFEM) is developed to solve the transient temperature field.To improve the accuracy of SBFEM,the effect of high frequency factor on dynamic stiffness...An increment-dimensional scaled boundary finite element method (ID-SBFEM) is developed to solve the transient temperature field.To improve the accuracy of SBFEM,the effect of high frequency factor on dynamic stiffness is considered,and the first-order continued fraction technique is used.After the derivation,the SBFE equations are obtained,and the dimensions of thermal conduction,the thermal capacity matrix and the vector of the right side term in the equations are doubled.An example is presented to illustrate the feasibility and good accuracy of the proposed method.展开更多
The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique combining the advantage of the finite element method （FEM） and the boundary element method （BEM） with its unique properties....The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique combining the advantage of the finite element method （FEM） and the boundary element method （BEM） with its unique properties. In this paper, the SBFEM is used for computing wave passing submerged breakwaters, and the reflection coeffcient and transmission coefficient are given for the case of wave passing by a rectangular submerged breakwater, a rigid submerged barrier breakwater and a trapezium submerged breakwater in a constant water depth. The results are compared with the analytical solution and experimental results. Good agreement is obtained. Through comparison with the results using the dual boundary element method （DBEM）, it is found that the SBFEM can obtain higher accuracy with fewer elements. Many submerged breakwaters with different dimensions are computed by the SBFEM, and the changing character of the reflection coeffcient and the transmission coefficient are given in the current study.展开更多
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 ...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.展开更多
Time-resolved particle image velocimetry（TRPIV） experiments are performed to investigate the coherent structure＇s performance of riblets in a turbulent boundary layer（TBL） at a friction Reynolds number of 185. To...Time-resolved particle image velocimetry（TRPIV） experiments are performed to investigate the coherent structure＇s performance of riblets in a turbulent boundary layer（TBL） at a friction Reynolds number of 185. To visualize the energetic large-scale coherent structures（CSs） over a smooth surface and riblets, the proper orthogonal decomposition（POD） and finite-time Lyapunov exponent（FTLE） are used to identify the CSs in the TBL. Spatial-temporal correlation is implemented to obtain the characters and transport properties of typical CSs in the FTLE fields. The results demonstrate that the generic flow structures, such as hairpin-like vortices, are also observed in the boundary layer flow over the riblets, consistent with its smooth counterpart. Low-order POD modes are more sensitive to the riblets in comparison with the high-order ones,and the wall-normal movement of the most energy-containing structures are suppressed over riblets. The spatial correlation analysis of the FTLE fields indicates that the evolution process of the hairpin vortex over riblets are inhibited. An apparent decrease of the convection velocity over riblets is noted, which is believed to reduce the ejection/sweep motions associated with high shear stress from the viscous sublayer. These reductions exhibit inhibition of momentum transfer among the structures near the wall in the TBL flows.展开更多
In this paper, the problem of a crack perpendicular to and terminating at an interface in bimaterial structure with finite boundaries is investigated. The dislocation simulation method and boundary collocation approac...In this paper, the problem of a crack perpendicular to and terminating at an interface in bimaterial structure with finite boundaries is investigated. The dislocation simulation method and boundary collocation approach are used to derive and solve the basic equations. Two kinds of loading form are considered when the crack lies in a softer or a stiffer material, one is an ideal loading and the other one fits to the practical experiment loading. Complete solutions of the stress field including the T stress are obtained as well as the stress intensity factors. Influences of T stress on the stress field ahead of the crack tip are studied. Finite boundary effects on the stress intensity factors are emphasized. Comparisons with the problem presented by Chen et al. (Int. J. Solids and Structure, 2003, 40, 2731–2755) are discussed also.展开更多
Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,compara...Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,comparatively fewer works have been devoted to study the effect of adsorption boundary layer(ABL)in throats based on the digital rock method.By considering an ABL,we investigate its effects on fluid flow.We build digital rock model based on computed tomography technology.Then,microscopic pore structures are extracted with watershed segmentation and pore geometries are meshed through Delaunay triangulation approach.Finally,using the meshed digital simulation model and finite element method,we investigate the effects of viscosity and thickness of ABL on microscale flow.Our results demonstrate that viscosity and thickness of ABL are major factors that significantly hinder fluid flow in throats.展开更多
The numerical solutions for uncertain viscoelastic problems have important theo- retical and practical significance. The paper develops a new approach by combining the scaled boundary finite element method （SBFEM） a...The numerical solutions for uncertain viscoelastic problems have important theo- retical and practical significance. The paper develops a new approach by combining the scaled boundary finite element method （SBFEM） and fuzzy arithmetic. For the viscoelastic problems with zero uncertainty, the SBFEM and the temporally piecewise adaptive algorithm is employed in the space domain and the time domain, respectively, in order to provide an accurate semi- analytical boundary-based approach and to ensure the accuracy of discretization in the time domain with different sizes of time step at the same time. The fuzzy arithmetic is used to address the uncertainty analysis of viscoelastic material parameters, and the transformation method is used for computation with the advantages of effectively avoiding overestimation and reducing the computational costs. Numerical examples are provided to test the performance of the proposed method. By comparing with the analytical solutions and the Monte Carlo method, satisfactory results are achieved.展开更多
Axisymmetric fundamental solutions that are applied in the consolidation calculations of a finite clay layer with impeded boundaries were derived. Laplace and Hankel integral transforms were utilized with respect to t...Axisymmetric fundamental solutions that are applied in the consolidation calculations of a finite clay layer with impeded boundaries were derived. Laplace and Hankel integral transforms were utilized with respect to time and radial coordinates, respectively in the analysis. The derivation of fundamental solutions considers two boundary value problems involving unit point loading and ring loading in the vertical. The solutions are extended to circular distributed and strip distributed normal load. The computation and analysis of settlements, vertical total stress and excess pore pressure in the consolidation layer subject to circular loading are presented.展开更多
A general and efficient parallel approach is proposed for the first time to parallelize the hybrid finiteelement-boundary-integral-multi-level fast multipole algorithm （FE-BI-MLFMA）. Among many algorithms of FE-BI-M...A general and efficient parallel approach is proposed for the first time to parallelize the hybrid finiteelement-boundary-integral-multi-level fast multipole algorithm （FE-BI-MLFMA）. Among many algorithms of FE-BI-MLFMA, the decomposition algorithm （DA） is chosen as a basis for the parallelization of FE-BI-MLFMA because of its distinct numerical characteristics suitable for parallelization. On the basis of the DA, the parallelization of FE-BI-MLFMA is carried out by employing the parallelized multi-frontal method for the matrix from the finiteelement method and the parallelized MLFMA for the matrix from the boundary integral method respectively. The programming and numerical experiments of the proposed parallel approach are carried out in the high perfor- mance computing platform CEMS-Liuhui. Numerical experiments demonstrate that FE-BI-MLFMA is efficiently parallelized and its computational capacity is greatly improved without losing accuracy, efficiency, and generality.展开更多
The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special fe...The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors （SIFs） directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.展开更多
The coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and ex...The coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.A kind of characteristic finite difference schemes is put forward,from which optimal order estimates in l~2 norm are derived for the error in the approximate solutions.The research is important both theoretically and practically for the model analysis in the field,the model numerical method and software development.展开更多
基金supported by the Open Project of Key Laboratory of Aerospace EDLA,CASC(No.EDL19092208)。
文摘A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface immersed boundary(IB)method,which is attractive for simulating moving-boundary flows with large deformations.The adaptive mesh refinement technique is employed to reduce the computational cost while retain the desired resolution.The dynamic response of the parachute is solved with the finite element approach.The canopy and cables of the parachute system are modeled with the hyperelastic material.A tether force is introduced to impose rigidity constraints for the parachute system.The accuracy and reliability of the present framework is validated by simulating inflation of a constrained square plate.Application of the present framework on several canonical cases further demonstrates its versatility for simulation of parachute inflation.
基金supported by the National Key Research and Development Program of China(Grant No.2022YFA1405400)the Innovation Program for Quantum Science and Technology(Grant No.2021ZD0301902),the National Natural Science Foundation of China(Grant No.12274290)the sponsor-ship from Yangyang Development Fund.
文摘We perform a comprehensive study of the electron-doped t-t′-J model on cylinders with density matrix renormalization group(DMRG).We conduct a systematic study on the finite-size and boundary condition effects on t-t′-J model on cylinders.Periodic and anti-periodic boundary conditions are implemented along the circumference direction,with the system’s width extending up to as large as 8 lattice units.We study doping levels of 1/6,1/8,and 1/12,which represent the most interesting region in the phase diagram of electron-doped cuprates.We find that for width-4 and width-6 systems,the ground state for fixed doping switches between anti-ferromagnetic Neel state and stripe state under different boundary conditions and system widths,indicating the presence of large finite size effect in the t-t′-J model.We also have a careful analysis of the d-wave pairing correlations which also change quantitatively with boundary conditions and widths of the system.However,the pairing correlations are enhanced when the system becomes wider for all dopings,suggesting the existence of possible long-range superconducting order in the thermodynamic limit.The width-8 results are found to be dependent on the starting state in the DMRG calculation for the kept states we can reach.For the width-8 system,only Neel(stripe)state can be stabilized in DMRG calculation for 1/12(1/6)doping,while both stripe and Neel states are stable in the DMRG sweep for 1/8 doping,regardless of the boundary conditions.These results indicate that 1/8 doping is likely to lie on the boundary of a phase transition between the Neel phase with lower doping and the stripe phase with higher doping,consistent with the previous study.The sensitivity of the ground state on boundary conditions and size observed for narrow systems is similar to that found in the t′-Hubbard model,where the t′term introduces frustration and makes the stripe state fragile.The study of different boundary conditions provides a useful tool to check the finite size effect in the future DMRG calculations.
文摘The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analytical one. And then examples are given to demonstrate the solution of several problems of the irregular object scattering. The method can not only produce clear physical pictures, but can efficiently handle many complicated scattering problems.
基金sponsored by the National Natural Science Foundation of China Research(Grant No.41274138)the Science Foundation of China University of Petroleum(Beijing)(No.KYJJ2012-05-02)
文摘The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.
文摘By establishing equivalent fixed point theorem, the boundary value problems of p Laplace equations with finite time delay are studied. It’s the first time that the functional differential equation is discussed with p Laplacian. The topological degree and fixed point theorem on cone are used to prove the existence of solution and positive solution. The conditions are all easy to check.
基金supported by the National Natural Science Foundation of China(NSFC)(Grant No. 41074100)the Program for New Century Excellent Talents in University of Ministry of Education of China(Grant No. NCET-10-0812)
文摘In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition （ABC） to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.
基金supported by the National Natural Science Foundation of China(Nos.11171193 and11371229)the Natural Science Foundation of Shandong Province(No.ZR2014AM033)the Science and Technology Development Project of Shandong Province(No.2012GGB01198)
文摘An implicit finite difference method is developed for a one-dimensional frac- tional percolation equation （FPE） with the Dirichlet and fractional boundary conditions. The stability and convergence are discussed for two special cases, i.e., a continued seep- age flow with a monotone percolation coefficient and a seepage flow with the fractional Neumann boundary condition. The accuracy and efficiency of the method are checked with two numerical examples.
基金heprojectissupportedbyNNSFofChina (No .1 9972 0 39) .
文摘In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.
基金SuppoSed by the NSF of Anhui Provincial Education Depaxtment(KJ2012A265,KJ2012B187)
文摘By using cone theory and the MSnch fixed theorem combined with a monotone iterative technique, we investigate the existence of positive solutions for systems of second- order nonlinear singular differential equations with integral boundary conditions on infinite interval and establish the existence theorem of positive solutions and iterative sequence for approximating the positive solutions. The results in this paper improve some known results.
基金supported by the Innovation Training Project for Students in NUAA(No.2016C-X0010-129)the Key Laboratory of Aircraft Environment Control and Life Support(NUAA),Ministry of Industry and Information Technology
文摘An increment-dimensional scaled boundary finite element method (ID-SBFEM) is developed to solve the transient temperature field.To improve the accuracy of SBFEM,the effect of high frequency factor on dynamic stiffness is considered,and the first-order continued fraction technique is used.After the derivation,the SBFE equations are obtained,and the dimensions of thermal conduction,the thermal capacity matrix and the vector of the right side term in the equations are doubled.An example is presented to illustrate the feasibility and good accuracy of the proposed method.
基金This research wasfinanciallysupported bythe National Natural Science Foundation of China(Grant No.50639030)a Programfor Changjiang ScholarsInnovative Research Teamin Dalian University of Technology(Grant No.IRTO420)
文摘The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique combining the advantage of the finite element method （FEM） and the boundary element method （BEM） with its unique properties. In this paper, the SBFEM is used for computing wave passing submerged breakwaters, and the reflection coeffcient and transmission coefficient are given for the case of wave passing by a rectangular submerged breakwater, a rigid submerged barrier breakwater and a trapezium submerged breakwater in a constant water depth. The results are compared with the analytical solution and experimental results. Good agreement is obtained. Through comparison with the results using the dual boundary element method （DBEM）, it is found that the SBFEM can obtain higher accuracy with fewer elements. Many submerged breakwaters with different dimensions are computed by the SBFEM, and the changing character of the reflection coeffcient and the transmission coefficient are given in the current study.
基金the State Key Program of National Natural Science of China under Grant No.51138001Science Fund for Creative Research Groups of the National Natural Science Foundation of China under Grant No.51121005Open Research Fund Program of State key Laboratory of Hydro science and Engineering under Grant No.shlhse-2010-C-03
文摘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.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11332006,11732010,11572221,and 11502066)the Natural Science Foundation of Tianjin City(Grant No.18JCQNJC5100)
文摘Time-resolved particle image velocimetry（TRPIV） experiments are performed to investigate the coherent structure＇s performance of riblets in a turbulent boundary layer（TBL） at a friction Reynolds number of 185. To visualize the energetic large-scale coherent structures（CSs） over a smooth surface and riblets, the proper orthogonal decomposition（POD） and finite-time Lyapunov exponent（FTLE） are used to identify the CSs in the TBL. Spatial-temporal correlation is implemented to obtain the characters and transport properties of typical CSs in the FTLE fields. The results demonstrate that the generic flow structures, such as hairpin-like vortices, are also observed in the boundary layer flow over the riblets, consistent with its smooth counterpart. Low-order POD modes are more sensitive to the riblets in comparison with the high-order ones,and the wall-normal movement of the most energy-containing structures are suppressed over riblets. The spatial correlation analysis of the FTLE fields indicates that the evolution process of the hairpin vortex over riblets are inhibited. An apparent decrease of the convection velocity over riblets is noted, which is believed to reduce the ejection/sweep motions associated with high shear stress from the viscous sublayer. These reductions exhibit inhibition of momentum transfer among the structures near the wall in the TBL flows.
基金The project supported by the National Natural Science Foundation of China(10202023 and 10272103)the Key Project of CAS(KJCX2-SW-L2).
文摘In this paper, the problem of a crack perpendicular to and terminating at an interface in bimaterial structure with finite boundaries is investigated. The dislocation simulation method and boundary collocation approach are used to derive and solve the basic equations. Two kinds of loading form are considered when the crack lies in a softer or a stiffer material, one is an ideal loading and the other one fits to the practical experiment loading. Complete solutions of the stress field including the T stress are obtained as well as the stress intensity factors. Influences of T stress on the stress field ahead of the crack tip are studied. Finite boundary effects on the stress intensity factors are emphasized. Comparisons with the problem presented by Chen et al. (Int. J. Solids and Structure, 2003, 40, 2731–2755) are discussed also.
基金National Natural Science Foundation of China(No.51674280,51774308,51704033,51722406,51950410591)Shandong Provincial Natural Science Foundation(ZR2019JQ21,JQ201808)+3 种基金the Fundamental Research Funds for the Central Universities(No.20CX02113A)National Science and Technology Major Project(2016ZX05014-000407)Program for Changjiang Scholars and Innovative Research Team in University(IRT_16R69)PetroChina Innovation Foundation(No.2018D-5007-0210)。
文摘Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,comparatively fewer works have been devoted to study the effect of adsorption boundary layer(ABL)in throats based on the digital rock method.By considering an ABL,we investigate its effects on fluid flow.We build digital rock model based on computed tomography technology.Then,microscopic pore structures are extracted with watershed segmentation and pore geometries are meshed through Delaunay triangulation approach.Finally,using the meshed digital simulation model and finite element method,we investigate the effects of viscosity and thickness of ABL on microscale flow.Our results demonstrate that viscosity and thickness of ABL are major factors that significantly hinder fluid flow in throats.
文摘The numerical solutions for uncertain viscoelastic problems have important theo- retical and practical significance. The paper develops a new approach by combining the scaled boundary finite element method （SBFEM） and fuzzy arithmetic. For the viscoelastic problems with zero uncertainty, the SBFEM and the temporally piecewise adaptive algorithm is employed in the space domain and the time domain, respectively, in order to provide an accurate semi- analytical boundary-based approach and to ensure the accuracy of discretization in the time domain with different sizes of time step at the same time. The fuzzy arithmetic is used to address the uncertainty analysis of viscoelastic material parameters, and the transformation method is used for computation with the advantages of effectively avoiding overestimation and reducing the computational costs. Numerical examples are provided to test the performance of the proposed method. By comparing with the analytical solutions and the Monte Carlo method, satisfactory results are achieved.
文摘Axisymmetric fundamental solutions that are applied in the consolidation calculations of a finite clay layer with impeded boundaries were derived. Laplace and Hankel integral transforms were utilized with respect to time and radial coordinates, respectively in the analysis. The derivation of fundamental solutions considers two boundary value problems involving unit point loading and ring loading in the vertical. The solutions are extended to circular distributed and strip distributed normal load. The computation and analysis of settlements, vertical total stress and excess pore pressure in the consolidation layer subject to circular loading are presented.
文摘A general and efficient parallel approach is proposed for the first time to parallelize the hybrid finiteelement-boundary-integral-multi-level fast multipole algorithm （FE-BI-MLFMA）. Among many algorithms of FE-BI-MLFMA, the decomposition algorithm （DA） is chosen as a basis for the parallelization of FE-BI-MLFMA because of its distinct numerical characteristics suitable for parallelization. On the basis of the DA, the parallelization of FE-BI-MLFMA is carried out by employing the parallelized multi-frontal method for the matrix from the finiteelement method and the parallelized MLFMA for the matrix from the boundary integral method respectively. The programming and numerical experiments of the proposed parallel approach are carried out in the high perfor- mance computing platform CEMS-Liuhui. Numerical experiments demonstrate that FE-BI-MLFMA is efficiently parallelized and its computational capacity is greatly improved without losing accuracy, efficiency, and generality.
基金The project supported by the National Natural Science Foundation of China (50579081)the Australian Research Council (DP0452681)The English text was polished by Keren Wang
文摘The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors （SIFs） directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.
基金the Major State Basic Research Program of China(No.G19990328)the National Tackling Key Problem Program(No.20050200069)+1 种基金the National Natural Science Foundation of China(Nos.10771124,10372052)the Ph.D.Programs Foundation of Ministry of Education of China(No.20030422047)
文摘The coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.A kind of characteristic finite difference schemes is put forward,from which optimal order estimates in l~2 norm are derived for the error in the approximate solutions.The research is important both theoretically and practically for the model analysis in the field,the model numerical method and software development.
