We considermulti-physics computationswhere theNavier-Stokes equations of compressible fluid flow on some parts of the computational domain are coupled to the equations of elasticity on other parts of the computational...We considermulti-physics computationswhere theNavier-Stokes equations of compressible fluid flow on some parts of the computational domain are coupled to the equations of elasticity on other parts of the computational domain.The different subdomains are separated by well-defined interfaces.We consider time accurate computations resolving all time scales.For such computations,explicit time stepping is very efficient.We address the issue of discrete interface conditions between the two domains of different physics that do not lead to instability,or to a significant reduction of the stable time step size.Finding such interface conditions is non-trivial.We discretize the problem with high order centered difference approximations with summation by parts boundary closure.We derive L2 stable interface conditions for the linearized one dimensional discretized problem.Furthermore,we generalize the interface conditions to the full non-linear equations and numerically demonstrate their stable and accurate performance on a simple model problem.The energy stable interface conditions derived here through symmetrization of the equations contain the interface conditions derived through normal mode analysis by Banks and Sj¨ogreen in[8]as a special case.展开更多
基金supported by Lawrence Livermore National Laboratory under the auspices of the U.S.Department of Energy through contract number DE-AC52-07NA27344.
文摘We considermulti-physics computationswhere theNavier-Stokes equations of compressible fluid flow on some parts of the computational domain are coupled to the equations of elasticity on other parts of the computational domain.The different subdomains are separated by well-defined interfaces.We consider time accurate computations resolving all time scales.For such computations,explicit time stepping is very efficient.We address the issue of discrete interface conditions between the two domains of different physics that do not lead to instability,or to a significant reduction of the stable time step size.Finding such interface conditions is non-trivial.We discretize the problem with high order centered difference approximations with summation by parts boundary closure.We derive L2 stable interface conditions for the linearized one dimensional discretized problem.Furthermore,we generalize the interface conditions to the full non-linear equations and numerically demonstrate their stable and accurate performance on a simple model problem.The energy stable interface conditions derived here through symmetrization of the equations contain the interface conditions derived through normal mode analysis by Banks and Sj¨ogreen in[8]as a special case.
