Method for the numerical simulation of a physical phenomenon with a preferential direction
First Claim
1. A computation method for simulating behavior of a flow interacting with an object, comprising the steps of:
- developing a simulated numerical representation in N dimensions, N≧
3, composed of a plurality of approximated values at a multitude of points in at least a part of space where a flow interacts with an object, the approximated values being of a physical parameter u of the flow to which is associated a velocity field {right arrow over (a)} which determines a preferential direction by executing a numerical scheme wherein at least one spatial pth derivative Dp, p≧
1, of the parameter u is approximated at the points of the part of space, through the steps ofi) constructing a discrete N-dimensional grid constructed by N families of coordinate lines and using the constructed grid for the part of space,ii) computing, in at least one point P of the grid, called the point of computation, an approximated value of DpA of Dp with an error ε
n, by using values uS of the parameter in a collection of grid points, called the stencil S, and computational functions, evaluated with the values uS, the computational functions depending on the numerical framework in which Dp is expressed, andiii) choosing the computational functions for the approximated value DpA in such a way that the approximated value DpA is optimized for the preferential direction,wherein the stencil S contains at least one point situated outside all the coordinate lines passing through the point of computation P, and the stencil S contains at least a first point and a second point, the first point being defined by N first coordinate lines of the N families of lines, the second point being defined by N second coordinate lines of the N families of lines, and for at least one family Nƒ
of the coordinate lines, the first coordinate line belonging to the family Nƒ
is different from and not adjacent to the second coordinate line belonging to the same family Nƒ
;
outputting the simulated numerical representation that simulates, for the part of space, the behavior of the flow interacting with the object; and
determining the behavior of the flow interacting with the object by using the output simulated numerical representation.
0 Assignments
0 Petitions
Accused Products
Abstract
The invention provides a method of simulating behavior of a flow interacting with an object. The invention improves the accuracy of the approximation of the spatial pth derivative Dp, p≧1, and therefore reduces the cost of the simulation. It includes in the approximation values of the parameter at points which do not lay on grid lines passing through the point of computation P, and by using these values to optimize the approximation. The use of these additional points depends on a preferential direction, such as determined by the advection direction of the flow. The points used in the approximation extend beyond a unit cube on the computational grid. The simulation is in three or more dimensions. The numerical simulation produces output which can be used in the design or optimization of the object which interacts with the flow. In certain cases, it is the output by itself which is important.
-
Citations
52 Claims
-
1. A computation method for simulating behavior of a flow interacting with an object, comprising the steps of:
-
developing a simulated numerical representation in N dimensions, N≧
3, composed of a plurality of approximated values at a multitude of points in at least a part of space where a flow interacts with an object, the approximated values being of a physical parameter u of the flow to which is associated a velocity field {right arrow over (a)} which determines a preferential direction by executing a numerical scheme wherein at least one spatial pth derivative Dp, p≧
1, of the parameter u is approximated at the points of the part of space, through the steps ofi) constructing a discrete N-dimensional grid constructed by N families of coordinate lines and using the constructed grid for the part of space, ii) computing, in at least one point P of the grid, called the point of computation, an approximated value of DpA of Dp with an error ε
n, by using values uS of the parameter in a collection of grid points, called the stencil S, and computational functions, evaluated with the values uS, the computational functions depending on the numerical framework in which Dp is expressed, andiii) choosing the computational functions for the approximated value DpA in such a way that the approximated value DpA is optimized for the preferential direction, wherein the stencil S contains at least one point situated outside all the coordinate lines passing through the point of computation P, and the stencil S contains at least a first point and a second point, the first point being defined by N first coordinate lines of the N families of lines, the second point being defined by N second coordinate lines of the N families of lines, and for at least one family Nƒ
of the coordinate lines, the first coordinate line belonging to the family Nƒ
is different from and not adjacent to the second coordinate line belonging to the same family Nƒ
;outputting the simulated numerical representation that simulates, for the part of space, the behavior of the flow interacting with the object; and determining the behavior of the flow interacting with the object by using the output simulated numerical representation. - View Dependent Claims (2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52)
-
2. The method of claim 1, wherein for varying preferential directions, the approximation DpA depends continuously on the values uS.
-
3. The method of claim 1, wherein the computing step of DpA comprises the steps of:
-
providing a local basis B({right arrow over (e)}1,{right arrow over (e)}2,{right arrow over (e)}3, . . . ) of curvilinear coordinates which has the unit vector {right arrow over (e)}1 along the preferential direction, and choosing the computational functions so that a contribution to the error ε
n of at least one pure or one mixed derivative as expressed in the local basis B is minimized, while using as a formulation of the values uS of the parameter at each of the points of the stencil S, a truncated Taylor series expansion with respect to the point of computation P with an error, called the truncation error ε
S.
-
-
4. The method of claim 3, wherein the computational functions are individual coefficients CS and the approximated value DpA is a linear combination of values uS of the parameter.
-
5. The method of claim 3, wherein the computing step of DpA comprises the steps of:
-
using an integral formulation for the computation of the derivative Dp with computational functions which are fluxes through a control volume, computing in at least one volume an approximated value DpA of Dp with an error ε
n, where the approximated value is a function of the flux formulation employed.
-
-
6. The method of claim 3, wherein the computing step of DpA comprises the steps of:
-
subdividing the part of space into elements, containing nodes at which the approximation of the physical parameter u is stored, defining basis computational functions φ
, called interpolation functions, on the elements, where the interpolation functions are used to approximate the physical value u on the element,computing the integral of the derivative Dp on the element, using a test computational function ψ
, called weighting function,computing in at least one element an approximated value Dpφ
,ψ
of Dp with an error ε
n where the approximated value is a function of the interpolation functions φ and
the weighting function ψ
.
-
-
7. The method of claim 3, wherein the computing step of DpA comprises the steps of:
-
subdividing the part of space into elements, containing nodes at which the approximation of the physical parameter u is stored, defining basis computational functions φ
, called interpolation functions, on the elements, where the interpolation functions are used to approximate the physical value u on the element,computing the integral Iel of the derivative Dp on volumes, and distributing parts α
iIel to nodes i, where α
i represent computational functions called distribution coefficients.
-
-
8. The method of claim 3, wherein the computing step of DpA comprises the steps of:
-
subdividing the part of space into elements, containing nodes at which the approximation of the physical parameter u is stored, defining basis computational functions φ
, called interpolation functions, on the elements, where the interpolation functions are used to approximate the physical value u on the element,computing fluxes ƒ
at the surfaces of volumes, and distributing parts α
iƒ
to nodes i, where α
i represent computational functions called distribution coefficients.
-
-
9. The method of claim 1, wherein the computing step of DpA comprises the steps of:
-
using a representation of the numerical solutions in Fourier components, and choosing in the approximated value DpA the computational functions in such a way that the Fourier components are optimized for certain directions which are related to the velocity {right arrow over (a)}, while using the values uS of the parameter at each of the points of the stencil S in the Fourier representation.
-
-
10. The method of claim 9, wherein the computational functions are individual coefficients CS and the approximated value DpA is a linear combination of values uS of the parameter.
-
11. The method of claim 9, wherein the computing step of DpA comprises the steps of:
-
using an integral formulation for the computation of the derivative Dp with computational functions which are fluxes through a control volume, computing in at least one volume an approximated value DpA of Dp with an error ε
n, where the approximated value is a function of the flux formulation employed.
-
-
12. The method of claim 9, wherein the computing step of DpA comprises the steps of:
-
subdividing the part of space into elements, containing nodes at which the approximation of the physical parameter u is stored, defining basis computational functions φ
, called interpolation functions, on the elements, where the interpolation functions are used to approximate the physical value u on the element,computing the integral of the derivative Dp on the element, using a test computational function ψ
, called weighting function,computing in at least one element an approximated value Dpφ
,ψ
of Dp with an error ε
n, where the approximated value is a function of the interpolation functions φ and
the weighting function ψ
.
-
-
13. The method of claim 9, wherein the computing step of DpA comprises the steps of:
-
subdividing the part of space into elements, containing nodes at which the approximation of the physical parameter u is stored, defining basis computational functions φ
, called interpolation functions, on the elements, where the interpolation functions are used to approximate the physical value u on the element,computing the integral Iel of the derivative Dp on volumes, and distributing parts α
iIel to nodes i, where α
i represent computational functions called distribution coefficients.
-
-
14. The method of claim 9, wherein the computing step of DpA comprises the steps of:
-
subdividing the part of space into elements, containing nodes at which the approximation of the physical parameter u is stored, defining basis computational functions φ
, called interpolation functions, on the elements, where the interpolation functions are used to approximate the physical value u on the element,computing fluxes ƒ
at the surfaces of volumes, and distributing parts α
iƒ
to nodes i, where α
i represent computational functions called distribution coefficients.
-
-
15. The method of claim 1, wherein the N-dimensional grid is expressed in a coordinate system which is chosen from the group consisting of:
- rectangular coordinates, spherical coordinates, cylindrical coordinates, parabolic cylindrical coordinates, paraboloidal coordinates, elliptic cylindrical coordinates, prolate spheroidal coordinates, oblate spheroidal coordinates, bipolar coordinates, toroidal coordinates, conical coordinates, confocal ellipsoidal coordinates and confocal paraboloidal coordinates.
-
16. The method of claim 1, wherein the N-dimensional grid is chosen from the group consisting of:
- a grid with non-uniform mesh spacing, a grid which is moving, a grid which is deforming, a grid which is rotating, and a grid which is staggered, and any combination thereof.
-
17. The method of claim 1, wherein the spatial pth derivative Dp is a pure derivative
-
p u ∂ e i p .
-
-
18. The method of claim 1, wherein p=1.
-
19. The method of claim 1, wherein p=2.
-
20. The method of claim 1, wherein the spatial pth derivative Dp is a mixed derivative
-
p u ∂ e 1 p 1 ∂ e 2 p 2 … with p1+p2+ . . . =p.
-
-
21. The method of claim 3, computing an approximated value D1A of D1, where the approximated value is denoted by
-
… ( u y ) i , j , k , … ( u z ) i , j , k , … ⋮ ) , ( 27 ) wherein wherein a−
m≠
0, an≠
0, m and n are given integers, m+n>
0, and m+n>
1 if m*n=0, the terms Tx, Ty, Tz, . . . represent the degrees of freedom which are used in the optimization of the approximated value D1A, and where indices (i, j, k, . . . ) define the point of computation P on the N-dimensional grid, and Δ
x, Δ
y, Δ
z, . . . denote the mesh spacings of the N-dimensional grid in each coordinate direction.
-
-
22. The method of claim 21, computing an approximated value of D1A of
-
u ∂ e 1 by a discretization of order M, in which the terms with M2+M3+ . . . =M+1 are eliminated in the approximated value D1A.
-
-
23. The method of claim 21, computing an approximated value D1A of
-
u ∂ e 1 by a discretization of order M, in which the terms with M1+M2+M3+ . . . =M+1 and M1<
M+1 are eliminated in the approximated value D1A in the case that {right arrow over (e)}1 is along the x-axis or along diagonals.
-
-
24. The method of claim 1, computing an approximated value D1A of D1, where the approximated value is denoted by
-
… ( u y ) i , j , k , … ( u z ) i , j , k , … ⋮ ) , ( 29 ) wherein wherein a−
m≠
0, an≠
0, m and n are given integers, m+n>
0, and m+n>
1 if m*n=0, and where indices (i, j, k, . . . ) define the point of computation P on the N-dimensional grid, and Δ
x, Δ
y, Δ
z, . . . denote the mesh spacings of the N-dimensional grid in each coordinate direction.
-
-
25. The method of claim 21, wherein
-
… = 1 Δ x { 1 2 ( u i + 1 , j , k , … - u i - 1 , j , k , … ) + Tx } , ( u y ) i , j , k , … = 1 Δ y { 1 2 u i + 1 , j + 1 , k , … - u i + 1 , j , k , … + u i - 1 , j , k , … - u i - 1 , j - 1 , k , … ) + Ty } , ( u z ) i , j , k , … = 1 Δ z { 1 2 ( u i + 1 , j + 1 , k + 1 , … - u i + 1 , j + 1 , k , … + u i - 1 , j - 1 , k , … - u i - 1 , j - 1 , k - 1 , … ) + Tz } , ⋮ ( 31 )
-
-
26. The method of claim 21, wherein
-
… = 1 Δ x { 1 12 ( u i - 2 , j , k , … - 8 u i - 1 , j , k , … + 8 u i + 1 , j , k , … - u i + 2 , j , k , … ) + Tx } , ( u y ) i , j , k , … = 1 Δ y { 1 12 ( u i - 2 , j - 2 , k , … - u i - 2 , j , k , … - 8 u i - 1 , j - 1 , k , … + 8 u i - 1 , j , k , … + 8 u i + 1 , j + 1 , k , … - 8 u i + 1 , j , k , … - u i + 2 , j + 2 , k , … + u i + 2 , j , k , … ) + Ty } , ( u z ) i , j , k , … = 1 Δ z { 1 12 ( u i - 2 , j - 2 , k - 2 , … - u i - 2 , j - 2 , k , … - 8 u i - 1 , j - 1 , k - 1 , … + 8 u i - 1 , j - 1 , k , … + 8 u i + 1 , j + 1 , k + 1 , … - 8 u i + 1 , j + 1 , k , … - u i + 2 , j + 2 , k + 2 , … + u i + 2 , j + 2 , k , … ) + Tz } , ⋮ ( 32 )
-
-
27. The method of claim 21, wherein
-
… = 1 Δ x ( 3 2 u i , j , k , … - 2 u i - 1 , j , k , … + 1 2 u i - 2 , j , k , … + Tx ) , ( u y ) i , j , k , … = 1 Δ y ( - 2 ( u i - 1 , j - 1 , k , … - u i - 1 , j , k , … ) + 1 2 ( u i - 2 , j - 2 , k , … - u i - 2 , j , k , … ) + Ty ) , ( u z ) i , j , k , … = 1 Δ z ( - 2 ( u i - 1 , j - 1 , k - 1 , … - u i - 1 , j - 1 , k , … ) + 1 2 ( u i - 2 , j - 2 , k - 2 , … - u i - 2 , j - 2 , k , … ) + Tz ) , ⋮ ( 33 )
-
-
28. A method according to claim 1, in three dimensions, obtaining the approximation DpA with order M of the derivative
-
p u ∂ e 1 p 1 ∂ e 2 p 2 ∂ e 3 p 3 with p1+p2+p3=p on a grid of given extent from the output of a computer program.
-
-
29. A method according to claim 28, wherein in the approximation DpA by a discretization of order M terms are eliminated, which are the terms
-
M + 1 u ∂ e 2 M 2 ∂ e 3 M 3 … with M2+M3+ . . . =M+1, or the terms with M1+M2+M3+ . . . =M+1 and M1<
M+1 in the case that {right arrow over (e)}1 is along the x-axis or along diagonals.
-
-
30. A method according to claim 28, wherein the approximation of DpA is a first order accurate discretization.
-
31. The method of claim 3, computing in three dimensions an approximation of the second derivative
-
2 u ∂ e 1 2 , where the approximation D2A is expressed in the terms uxx, uyy, uzz, uxy, uyz and uzx which are given by wherein the terms Txx, Txy, Txz, Tyy, Tyz and Tzz represent the degrees of freedom which are used in the optimization of the approximated value D2A and where indices (i, j, k) define the point of computation P on the three-dimensional grid, and Δ
x, Δ
y, Δ
z denote the mesh spacings of the three-dimensional grid in each coordinate direction.
-
-
32. The method according to claim 1, wherein
DpA=Σ - nLnDp,nA where Ln are constants and each Dp,nA is a function of values uS of the parameter in a collection of grid points, called stencil Sn, with individual computation functions, which depend on the numerical framework in which Dp,nA is expressed, and wherein in the approximation Dp,nA, the computation functions are chosen in such a way that the approximated value Dp,nA is optimized for the preferential direction.
-
33. The method of claim 1, wherein
DpA=Σ - nLnDp,nA where Ln are limiting functions of the values uS of the stencil S, and at least one Dp,nA is a function of values uS of the parameter in a collection of grid points, called stencil Sn, with individual computation functions, which depend on the numerical framework in which Dp,nA is expressed, and wherein in the approximation Dp,nA, the computation functions are chosen in such a way that said approximated value Dp,nA is optimized for the preferential direction.
-
34. The method of claim 1, wherein the stencil S is chosen from the group consisting of:
- upwind discretization stencils, centered discretization stencils, and discretization stencils which are a combination of at least one upwind discretization stencil and at least one centered discretization stencil.
-
35. The method of claim 1, wherein the numerical discretization is a non-linear discretization.
-
36. The method of claim 1, wherein the numerical scheme is chosen from the group consisting of:
- the Lax-Wendroff scheme, the Lax-Friedrich scheme, the MacCormack scheme, the leap-frog scheme, the Crank-Nicholson scheme, the Stone-Brian scheme, the box scheme, Henn'"'"'s scheme, the QUICK scheme, the κ
scheme, the Flux Corrected Transport (FTC) scheme, the family of ENO schemes, schemes in the class of the Piecewise Parabolic Method (PPM), multi-level schemes, and schemes obtained with the fractional step method and variants thereof.
- the Lax-Wendroff scheme, the Lax-Friedrich scheme, the MacCormack scheme, the leap-frog scheme, the Crank-Nicholson scheme, the Stone-Brian scheme, the box scheme, Henn'"'"'s scheme, the QUICK scheme, the κ
-
37. The method of claim 1, wherein the numerical scheme includes the discretization of a plurality of equations.
-
38. The method of claim 1 for the numerical simulation of physical phenomena which are modeled by the Navier-Stokes equations with equation(s) of state.
-
39. The method of claim 1 for the numerical simulation of physical phenomena which are modeled by the Euler equations with equation(s) of state.
-
40. The method of claim 1 for the numerical simulation of physical phenomena which are modeled by the magneto-hydrodynamic equations with equation(s) of state.
-
41. The method of claim 1 in combination with at least one model for the physical phenomena chosen from the group consisting of:
- equations to model turbulence, equations to model chemical reactions, equations to model electromagnetic phenomena, equations to model multiphase flow and equations to model multiphysics phenomena.
-
42. The method of claim 1 for the numerical simulation of physical phenomena in combination with at least one acceleration technique chosen from the group consisting of:
- local time-stepping, multi-grid, GMRES and preconditioning.
-
43. The method of claim 1 for the simulation of a physical phenomenon which includes a material flow.
-
44. The method of claim 1 for the simulation of a physical phenomenon which includes a material object interacting with material flow.
-
45. The method of claim 1, wherein the object is a vehicle.
-
46. The method of claim 1, wherein the object includes a rotating blade.
-
47. The method of claim 1, wherein the object is at least a part of the atmosphere of the earth.
-
48. The method of claim 1, wherein the flow includes oil.
-
49. The method of claim 1, wherein the flow includes combustion.
-
50. A data processing system programmed to implement a method according to claim 1.
-
51. A computer program that can be loaded in a data processing system so as to implement a method according to claim 1.
-
52. A digital storage computer-readable medium containing a stored computer program that is configured to be loadable in a data processing system so as to control the data processing system to implement a method according to claim 1.
-
2. The method of claim 1, wherein for varying preferential directions, the approximation DpA depends continuously on the values uS.
-
Specification
- Resources
-
Current AssigneeRobert Struijs
-
Original AssigneeRobert Struijs
-
InventorsStruijs, Robert
-
Primary Examiner(s)FREJD, RUSSELL WARREN
-
Application NumberUS10/779,685Publication NumberTime in Patent Office1,231 DaysField of Search703/2, 703/9, 345/473US Class Current703/2CPC Class CodesG06F 2111/10 Numerical modellingG06F 30/23 using finite element method...