One of the critical issues in numerical simulation of fluid-structure interaction problems is inaccuracy of the solutions,especially for flows past a stationary thin elastic structure where large deformations occur.Hi...One of the critical issues in numerical simulation of fluid-structure interaction problems is inaccuracy of the solutions,especially for flows past a stationary thin elastic structure where large deformations occur.High resolution is required to capture the flow characteristics near the fluid-structure interface to enhance accuracy of the solutions within proximity of the thin deformable body.Hence,in this work,an algorithm is developed to simulate fluid-structure interactions of moving deformable structures with very thin thicknesses.In this algorithm,adaptive mesh refinement(AMR)is integrated with immersed boundary finite element method(IBFEM)with two-stage pressure-velocity corrections.Despite successive interpolation of the flow field by IBM,the governing equations were solved using a fixed structured mesh,which significantly reduces the computational time associated with mesh reconstruction.The cut-cell IBM is used to predict the body forces while FEM is used to predict deformation of the thin elastic structure in order to integrate the motions of the fluid and solid at the interface.AMR is used to discretize the governing equations and obtain solutions that efficiently capture the thin boundary layer at the fluid-solid interface.The AMR-IBFEM algorithm is first verified by comparing the drag coefficient,lift coefficient,and Strouhal number for a benchmark case(laminar flow past a circular cylinder at Re=100)and the results showed good agreement with those of other researchers.The algorithm is then used to simulate 2-D laminar flows past stationary and moving thin structures positioned perpendicular to the freestream direction.The results also showed good agreement with those obtained from the arbitrary Lagrangian-Eulerian(ALE)algorithm for elastic thin boundaries.It is concluded that the AMR-IBFEM algorithm is capable of predicting the characteristics of laminar flow past an elastic structure with acceptable accuracy(error of-0.02%)with only-1%of the computational time for simulations with full mesh refinement.展开更多
This paper is concerned with the weak Galerkin mixed finite element method(WG-MFEM)for the second-order hyperbolic acoustic wave equation in velocity-pressure formulation.In this formulation,the original second-order ...This paper is concerned with the weak Galerkin mixed finite element method(WG-MFEM)for the second-order hyperbolic acoustic wave equation in velocity-pressure formulation.In this formulation,the original second-order differential equation in time and space is reduced to first-order differential equations by introducing the velocity and pressure variables.We employ the usual discontinuous piecewise-polynomials of degree k0 for the pressure and k+1 for the velocity.Furthermore,the normal component of the pressure on the interface of elements is enhanced by polynomials of degree k+1.The time derivative is approximated by the backward Euler difference.We show the stability of the semi-discrete and fullydiscrete schemes,and obtain the suboptimal order error estimates for the velocity and pressure variables.Numerical experiment confirms our theoretical analysis.展开更多
文摘One of the critical issues in numerical simulation of fluid-structure interaction problems is inaccuracy of the solutions,especially for flows past a stationary thin elastic structure where large deformations occur.High resolution is required to capture the flow characteristics near the fluid-structure interface to enhance accuracy of the solutions within proximity of the thin deformable body.Hence,in this work,an algorithm is developed to simulate fluid-structure interactions of moving deformable structures with very thin thicknesses.In this algorithm,adaptive mesh refinement(AMR)is integrated with immersed boundary finite element method(IBFEM)with two-stage pressure-velocity corrections.Despite successive interpolation of the flow field by IBM,the governing equations were solved using a fixed structured mesh,which significantly reduces the computational time associated with mesh reconstruction.The cut-cell IBM is used to predict the body forces while FEM is used to predict deformation of the thin elastic structure in order to integrate the motions of the fluid and solid at the interface.AMR is used to discretize the governing equations and obtain solutions that efficiently capture the thin boundary layer at the fluid-solid interface.The AMR-IBFEM algorithm is first verified by comparing the drag coefficient,lift coefficient,and Strouhal number for a benchmark case(laminar flow past a circular cylinder at Re=100)and the results showed good agreement with those of other researchers.The algorithm is then used to simulate 2-D laminar flows past stationary and moving thin structures positioned perpendicular to the freestream direction.The results also showed good agreement with those obtained from the arbitrary Lagrangian-Eulerian(ALE)algorithm for elastic thin boundaries.It is concluded that the AMR-IBFEM algorithm is capable of predicting the characteristics of laminar flow past an elastic structure with acceptable accuracy(error of-0.02%)with only-1%of the computational time for simulations with full mesh refinement.
基金supported by the National Natural Science Foundation of China(No.11971337)the Key Fund Project of Sichuan Provincial Department of Education(No.18ZA0276).
文摘This paper is concerned with the weak Galerkin mixed finite element method(WG-MFEM)for the second-order hyperbolic acoustic wave equation in velocity-pressure formulation.In this formulation,the original second-order differential equation in time and space is reduced to first-order differential equations by introducing the velocity and pressure variables.We employ the usual discontinuous piecewise-polynomials of degree k0 for the pressure and k+1 for the velocity.Furthermore,the normal component of the pressure on the interface of elements is enhanced by polynomials of degree k+1.The time derivative is approximated by the backward Euler difference.We show the stability of the semi-discrete and fullydiscrete schemes,and obtain the suboptimal order error estimates for the velocity and pressure variables.Numerical experiment confirms our theoretical analysis.