Fluid structure interaction simulation method and apparatus, and computer-readable storage medium

Patent metadata

Field	Value
Publication number	US-10114911-B2
Application number	US-201113104495-A
Country	US
Kind code	B2
Filing date	May 10, 2011
Priority date	May 24, 2010
Publication date	Oct 30, 2018
Grant date	Oct 30, 2018

Abstract

A fluid structure interaction simulation method includes a graph information forming process to form graph information of nodes obtained by discretizing a computing region for each of a fluid and a structure that are represented by meshes, and a main time development loop process to simulate a physical phenomenon. The loop process includes arranging IMEs (Interaction Mediating Elements) that move with a displacement of the structure, on a boundary of the structure, defining, within the IME, correcting functions of a pressure and a velocity of the fluid that interact with the pressure and the velocity of the fluid and the displacement of the structure, and executing a simulation based on the correcting functions, in a state in which the meshes of the fluid are mismatched to the meshes of the structure.

First claim

What is claimed is:

1. A fluid structure interaction simulation method to perform a process for simulating mechanical interaction between a fluid and a structure in a biomechanical system, executed by a processor included in a fluid structure interaction simulation apparatus, the process comprising:
a graph information forming process, executed by the processor, to form graph information of nodes obtained by discretising a computing region for each of a fluid and a structure in the biomechanical system that are represented by meshes; and
a main time development loop process, executed by the processor, to simulate the mechanical interaction between the fluid and the structure,
wherein the main time development loop process includes:
arranging interaction mediating elements that move with a displacement of the structure, on a boundary of the structure and the fluid;
defining, within the interaction mediating element, a correcting function of a pressure and a correcting function of a velocity of the fluid that interact with the pressure and the velocity of the fluid and the displacement of the structure, to mediate interaction of the displacement of the structure and the pressure and the velocity of the fluid;
executing a simulation based on the correcting functions, in a state in which the meshes of the fluid are mismatched to the meshes of the structure,

wherein the correcting functions are Lagrangian multipliers having, as a constraint condition, a relation satisfied by an incompressibility condition of the fluid and a relation satisfied by a viscous operator of a Newtonian fluid, at a boundary part between the fluid and the structure, based on Green's theorem, thereby compensating for the mismatched meshes of the fluid and the structure,
wherein the main time development loop process includes:
a first process to create both of a matrix and an equivalent nodal force of each element of the structure by using a time discretization technique and to create both of a matrix and an equivalent nodal force of each element of the fluid by using the time discretization technique;
a second process to create a matrix representing interaction of the velocity, the pressure, and the correcting functions for each interaction mediating element;
a third process to create a matrix equation with respect to all degrees of freedom by restructuring a contribution from the matrix of each element, based on the matrix of each element and the matrix representing the interaction;
a fourth process to solve and evaluate a variation of dependent variables, including the pressure, the velocity, and the displacement, of the fluid and the structure every time a Newton-Raphson loop is executed in the main time development loop process; and
a fifth process to display the variation of the dependent variables and the evaluation thereof.

2. The method as claimed in claim 1 , wherein the main time development loop process includes:
searching for a background fluid element having an overlap with the interaction mediating element, based on a change in an interaction pair formed by the interaction mediating element and a degree of freedom of the fluid as the interaction mediating element moves;
determining a fluid node that is a member of the fluid element as belonging to the interaction pair; and
executing the simulation based on the interaction pair.

3. The method as claimed in claim 1 , wherein the second process includes:
computing an interaction matrix of the correcting function of the pressure, and the velocity of the fluid;
computing an interaction matrix of the correcting function of the velocity, and the velocity of the fluid;
computing an interaction matrix of the correcting function of the pressure, and a displacement velocity of the structure; and
computing an interaction matrix of the correcting function of the velocity, and the displacement velocity of the structure.

4. The method as claimed in claim 1 ,
wherein the correcting functions have a non-zero value only within the interaction mediating element; and
wherein the main time development loop process applies a condition to each of a region having a value, a function value and a derived function on the structure, and a function value and a derived function on an element boundary, so that the correcting functions represent a discontinuity of the pressure and a velocity gradient while mediating the interaction of the fluid and the structure.

5. An apparatus comprising:
a storage configured to store a program; and
a processor, coupled to the storage, and configured to perform a fluid structure interaction simulation process for simulating mechanical interaction between a fluid and a structure in a biomechanical system by executing the program,
wherein the fluid structure interaction simulation process includes:
forming graph information of nodes obtained by discretising a computing region for each of a fluid and a structure that are represented by meshes; and
simulating the mechanical interaction between the fluid and the structure, wherein the simulating includes:
arranging interaction mediating elements that move with a displacement of the structure, on a boundary of the structure and the fluid;
defining, within the interaction mediating element, a correcting function of a pressure and a correcting function of a velocity of the fluid that interact with the pressure and the velocity of the fluid and the displacement of the structure to mediate interaction of the displacement of the structure and the pressure and the velocity of the fluid;
executing a simulation based on the correcting functions, in a state in which the meshes of the fluid are mismatched to the meshes of the structure,

wherein the correcting functions are Lagrangian multipliers having, as a constraint condition, a relation satisfied by an incompressibility condition of the fluid and a relation satisfied by a viscous operator of a Newtonian fluid, at a boundary part between the fluid and the structure, based on Green's theorem, thereby compensating for the mismatched meshes of the fluid and the structure,
wherein the simulating further includes:
first creating both of a matrix and an equivalent nodal force of each element of the structure by using a time discretization technique and to create both of a matrix and an equivalent nodal force of each element of the fluid by using the time discretization technique;
second creating a matrix representing interaction of the velocity, the pressure, and the correcting functions for each interaction mediating element;
creating a matrix equation with respect to all degrees of freedom by restructuring a contribution from the matrix of each element, based on the matrix of each element and the matrix representing the interaction;
solving and evaluating a variation of dependent variables, including the pressure, the velocity, and the displacement, of the fluid and the structure every time a Newton-Raphson loop is executed in the simulating; and
displaying the variation of the dependent variables and the evaluation thereof.

6. The apparatus as claimed in claim 5 , wherein the simulating further includes:
searching for a background fluid element having an overlap with the interaction mediating element, based on a change in an interaction pair formed by the interaction mediating element and a degree of freedom of the fluid as the interaction mediating element moves;
determining a fluid node that is a member of the fluid element as belonging to the interaction pair; and
executing the simulation based on the interaction pair.

7. The apparatus as claimed in claim 5 , wherein the second creating includes:
computing an interaction matrix of the correcting function of the pressure, and the velocity of the fluid;

Assignees

• Yamazaki Takashi
• Kadooka Yoshimasa
• Hisada Toshiaki
• Sugiura Seiryo
• Washio Takumi
• Okada Junichi
• Fujitsu Ltd
• Univ Tokyo

Inventors

• Yamazaki Takashi
• Kadooka Yoshimasa
• Hisada Toshiaki
• Sugiura Seiryo
• Washio Takumi
• Okada Junichi

Classifications

• G06F2111/10

  Numerical modelling · CPC title

• G06F30/20

  Design optimisation, verification or simulation (optimisation, verification or simulation of circuit designs G06F30/30) · CPC title

• G06F30/23 Primary

  using finite element methods [FEM] or finite difference methods [FDM] · CPC title

• G06F17/5018

  Physics · mapped topic

• G06F17/5009 Primary

  Physics · mapped topic