This study presents an effective numerical model for the dynamic response of poroelastic seabed under wave action with enhanced performance. The spatial discretization is based on the Element-Free Galerkin （EFG） met...This study presents an effective numerical model for the dynamic response of poroelastic seabed under wave action with enhanced performance. The spatial discretization is based on the Element-Free Galerkin （EFG） method and the time integration based on the GN11 scheme. A stability strategy that adopts a smaller number of nodes for the pore water pressure compared with those for the displacements of the soil skeleton is suggested to resolve the similar difficulty as encountered in the finite element method for a problem with mixed formulation when the pore water is incompressible and the soil skeleton impervious. The accuracy of the numerical model is verified through applying it to a typical case with critical permeability. Good agreement between computational and analytical solutions is obtained.展开更多
The multiscale hybrid-mixed(MHM)method is applied to the numerical approximation of two-dimensional matrix fluid flow in porous media with fractures.The two-dimensional fluid flow in the reservoir and the one-dimensio...The multiscale hybrid-mixed(MHM)method is applied to the numerical approximation of two-dimensional matrix fluid flow in porous media with fractures.The two-dimensional fluid flow in the reservoir and the one-dimensional flow in the discrete fractures are approximated using mixed finite elements.The coupling of the two-dimensional matrix flow with the one-dimensional fracture flow is enforced using the pressure of the one-dimensional flow as a Lagrange multiplier to express the conservation of fluid transfer between the fracture flow and the divergence of the one-dimensional fracture flux.A zero-dimensional pressure(point element)is used to express conservation of mass where fractures intersect.The issuing simulation is then reduced using the MHM method leading to accurate results with a very reduced number of global equations.A general system was developed where fracture geometries and conductivities are specified in an input file and meshes are generated using the public domain mesh generator GMsh.Several test cases illustrate the effectiveness of the proposed approach by comparing the multiscale results with direct simulations.展开更多
基金supported by the National Natural Science Foundation of China (Grant No. 10772099)
文摘This study presents an effective numerical model for the dynamic response of poroelastic seabed under wave action with enhanced performance. The spatial discretization is based on the Element-Free Galerkin （EFG） method and the time integration based on the GN11 scheme. A stability strategy that adopts a smaller number of nodes for the pore water pressure compared with those for the displacements of the soil skeleton is suggested to resolve the similar difficulty as encountered in the finite element method for a problem with mixed formulation when the pore water is incompressible and the soil skeleton impervious. The accuracy of the numerical model is verified through applying it to a typical case with critical permeability. Good agreement between computational and analytical solutions is obtained.
文摘The multiscale hybrid-mixed(MHM)method is applied to the numerical approximation of two-dimensional matrix fluid flow in porous media with fractures.The two-dimensional fluid flow in the reservoir and the one-dimensional flow in the discrete fractures are approximated using mixed finite elements.The coupling of the two-dimensional matrix flow with the one-dimensional fracture flow is enforced using the pressure of the one-dimensional flow as a Lagrange multiplier to express the conservation of fluid transfer between the fracture flow and the divergence of the one-dimensional fracture flux.A zero-dimensional pressure(point element)is used to express conservation of mass where fractures intersect.The issuing simulation is then reduced using the MHM method leading to accurate results with a very reduced number of global equations.A general system was developed where fracture geometries and conductivities are specified in an input file and meshes are generated using the public domain mesh generator GMsh.Several test cases illustrate the effectiveness of the proposed approach by comparing the multiscale results with direct simulations.
