Method for simulating forces applied to a wing in a fluid flow

Abstract

A method implemented by a computer for the real-time simulation of a force and/or of a moment applied to a wing in a fluid flow, wherein the body is previously modelled by at least two scalar and vector point sources of a velocity field, the method including at least one of the steps: assessing at least one element selected among a force and a moment, through a linear function dependent on the main velocity vector, V, of the flow; at least one flow rate, referred to as the scalar mass flow rate, As; and at least one flow rate, referred to as the vector mass flow rate, AR, associated with the sources, on at least one point corresponding to one of the sources.

Claims

1. A method implemented by a computer for real-time simulation of at least one element chosen from a force and a moment, applied to a body defined by a volume in a fluid flow having a main velocity vector representing the velocity of said flow at infinity, in which said body is modelled beforehand by at least two point sources of a velocity field, including at least one scalar source, generating a radial velocity field and associated with a scalar mass flow rate λ.sub.S, and at least one vector source—, generating a solenoidal velocity field, associated with a vector mass flow rate λ.sub.R, said sources being arranged at barycentre inside said volume, said method comprising the steps of: (a) simulating a superposition of the contributions of the fluid flows at least of one main velocity flow, of one said scalar source associated with a scalar mass flow rate λ.sub.S and of one said vector source associated with a vector mass flow rate λ.sub.R; and (b) evaluating at least one element chosen from a force and a moment, by a linear function dependent on said main velocity vector, of at least one said scalar mass flow rate λ.sub.S and of at least one said vector mass flow rate λ.sub.R, at least one point corresponding to one said source, said body being a wing.

2. The method according to claim 1, in which said simulation is a far-field simulation.

3. The method according to claim 2, in which said body is modelled by fewer than one hundred of said sources.

4. The method according to claim 1, in which two of said scalar sources are arranged so as to form a dipole.

5. The method according to claim 4, in which the scalar mass flow rates λ.sub.S associated with the two said sources forming said dipole are strictly different and have opposite signs.

6. The method according to claim 1, in which at least two contrarotating vector sources are arranged inside a volume defining said body, associated vector mass flow rates λ.sub.R of which have opposite directions.

7. The method according to claim 6, in which said vector mass flow rates are oriented parallel to said direction of a main velocity vector, the two said contrarotating vector sources forming a direction normal to the direction of said main velocity vector.

8. The method according to claim 7, in which the sources are arranged along a line corresponding to a simple wing.

9. The method according to claim 1, in which one element representative of the environment of the body is modelled beforehand by at least one element chosen from a scalar external source associated with an external scalar mass flow rate λ.sub.S′ and from one vector external source associated with an external vector mass flow rate λ.sub.R′, and in which: the contributions of the fluid flows of each of said sources are also superposed during said step (a); and in said step (b) said analytic linear function is evaluated also as a function of at least one said external mass flow rate λ.sub.S′ and at least one external vector mass flow rate λ.sub.R′.

10. The method according to claim 9, in which: said body is close to the ground, said ground coinciding locally with a plane, said sources being arranged beforehand in a first half-space a limit of which is formed by said plane, the first half-space being above the ground; and at least one mirror scalar external source associated with a mirror external scalar mass flow rate and at least one said mirror vector external source associated with a mirror external vector mass flow rate are arranged in the other half-space a limit of which is formed by said plane, the other half-space being below the ground plane, and in which: the contributions of the fluid flows of said mirror scalar external source and mirror vector external source are also superposed during said step (a); and in said step (b) said analytic linear function is evaluated also as a function of at least one said mirror external mass flow rate and of at least one said mirror external vector mass flow rate.

11. A computer program product stored on a non-transitory recording medium suitable for implementing a method according to claim 1.

12. A non-transitory recording medium on which the program according to claim 11 is stored.

13. A flight simulator comprising: at least one calculator configured in order to receive at least first data acquired by a user, a memory containing second data representative of the fluid flow around said body, and in that said calculator is configured in order to implement a method according to claim 1.

14. A method for the construction of a body model in a fluid flow having at least one main velocity vector representing the velocity at infinity, said body being associated with at least two aerodynamic parameters chosen from a drag force coefficient, a lift coefficient, a roll coefficient, a pitch coefficient and a yaw coefficient, each coefficient being determined in one or more known spatial coefficients of said body, in said fluid flow, said method comprising at least the steps consisting of: (a) arranging in the simulated space at least two point sources of fluid, including at least one scalar source of fluid, associated with an unknown scalar mass flow rate λ.sub.S, and at least one vector source of fluid, associated with an unknown vector mass flow rate λ.sub.R, said sources being arranged inside said volume defined by said body; (b) determining the value of each said scalar mass flow rate λ.sub.S and each said vector mass flow rate λ.sub.R as a function of a part of the boundary conditions imposed by said at least one main velocity vector and as a linear function of said at least one main velocity vector, said body being a wing.

15. The method according to claim 14 in which in said step (b), the value of each said scalar mass flow rate λ.sub.S and each said vector mass flow rate λ.sub.R is determined by the distributed point source method.

Description

DESCRIPTION OF THE DRAWINGS

(1) The invention will be better understood, and other advantages, details and characteristics thereof will become apparent during the following explanatory description, given by way of example with reference to the attached drawings in which:

(2) FIG. 1 shows the arrangement in space of sources S according to the invention:

(3) FIG. 2 shows a method for the construction of a model according to an embodiment of the invention;

(4) FIG. 3 shows the construction of a model according to the invention adapted to the simulation of an interaction between a wing and an external event:

(5) FIG. 4 shows a wing A in a configuration able to be affected by a ground effect;

(6) FIG. 5 shows the construction of a model according to the invention adapted to a simulation of the ground effect:

(7) FIG. 6 shows a simulation of the lift coefficient according to the invention:

(8) FIG. 7 shows a model construction according to the invention making it possible to simulate the Magnus effect and a shape drag force;

(9) FIG. 8 shows a cross section of an infinite cylinder, corresponding to a body or to a wing:

(10) FIG. 9 shows a simulation of a fluid flow around a wing A, the shape of which is an infinite cylinder simulated by a model according to the invention in which two scalar sources are arranged as a dipole;

(11) FIG. 10 shows drag force coefficients and lift coefficients evaluated comparatively according to a method of the prior art and according to the invention;

(12) FIG. 11 shows a wing A and an arrangement of sources S according to the invention,

(13) FIG. 12 is a photograph of a simulator SIM according to the invention.

DETAILED DESCRIPTION

(14) The invention comprises a step of calculation resolution by the distributed point source method or DPSM, described in patent applications WO2015177364 and WO2007071735. As a preliminary, this method requires knowledge of the equations governing the development of physical variables in the different simulated environments and their particular solution in the case of a point source (Green's function). The principle thereof can be compared to methods of the integral boundary type, singularity methods, or other boundary element methods (BEM). This method only requires meshing of the surfaces or of the physical interfaces between the objects of the simulation. The DPSM method is based on a spatial distribution of point sources, arranged on either side of active surfaces of the objects of the simulation. This semi-analytical calculation method relies on a superposition of emission (or “luminous”) points called sources, the weights of which, called associated values, are determined so as to satisfy the set of conditions at the boundaries of a problem. The principle consists of substituting for the objects present in the system, layers or surfaces comprising point sources situated on either side of an interface. The arrangement of the sources is associated with a mesh of test points situated on the interfaces. These sources are intended to substitute for the physical variables (fields, potential, pressure etc.) present in the simulation object and are calculated so as to verify the boundary conditions on the test points distributed over interfaces.

(15) The DPSM method uses the equality C=M.sub.i.Math.λ, in which C is a vector comprising the boundary conditions of the domain of a simulation. M.sub.i is an intercoupling matrix between sources S and test points, the test points being arranged in a mesh at the interfaces, and the sources S arranged on either side of the test points representing interfaces in prior art embodiments. λ is a vector of the values associated with the sources, unknown at the time of construction of a model. In all of the embodiments of the invention, the elements of λ are mass flow rates. The matrix M.sub.i can be inverted so as to determine λ. The value of a physical variable simulated at a point P of the space, or of the spatial domain of simulation, is calculated or evaluated by using the product M.sub.p.Math.λ (M.sub.p in this case representing an observation matrix).

(16) The inventors discovered that it is possible to use a part of this method so as to evaluate or simulate forces or moments applied to a wing: by calculating macroscopic variables. It is possible to use jointly, integral formulations such as the rotational theorem (Stokes' theorem) and/or the divergence theorem (Green-Ostrogradski theorem). In the embodiments of the invention, these integral equations are used in order to calculate forces or moments; by expressing forces or moments under far-field hypotheses. These hypotheses make it possible not to produce a mesh of the objects of the simulation, but to use representative points, as one or more barycentres. In this way, it is possible to calculate forces by using fewer than 100 sources S and preferentially fewer than 10 sources S, while a mesh produced in the prior art requires at least more than one thousand sources S, making the calculation more complex and extending the calculation processing time.

(17) FIG. 1 shows the arrangement in space of sources S according to the invention. In particular, FIG. 1 corresponds to the case in which a body C corresponds to a wing A. Generally, the different embodiments described can be implemented equally well for a body C or a wing A. In embodiments of the invention, a model is constructed of a wing A (or aircraft) in a fluid flow FLU. The fluid flow FLU is described at least by a main velocity vector V, shown by a bold arrow in FIG. 1. In embodiments of the invention, several vectors V can be considered: V can be different for each of the sources.

(18) The body C or the wing A can be described by a volume. This volume is for example defined by the interface between the wing A and the fluid flow FLU.

(19) By “scalar source S” is meant a point source S resulting in a radial velocity field. Each scalar source is associated with a scalar mass flow rate λ.sub.S. In embodiments of the invention, the dimension of the scalar mass flow rate λ.sub.S associated with a source S is homogenous at a mass flow rate, expressed in kg.Math.s.sup.−1.

(20) By “vector source S” is meant a point source S resulting in a solenoidal zero divergence velocity field. Each vector source is associated with a vector mass flow rate λ.sub.R. In embodiments of the invention, the dimension of a vector mass flow rate λ.sub.R is homogenous at a mass flow rate, expressed in kg.Math.s.sup.−1.

(21) By “far-field simulation” or “simulation produced in far-field” is meant a simulation proceeding to a phenomenological decomposition of the drag and/or the lift of a body C or of a wing A to be evaluated, as opposed to a “near-field simulation” corresponding to an evaluation of a drag and/or of a lift by integration of forces of pressure and/or of friction at the interface defining the volume of said body C or of said wing A.

(22) In FIG. 1, a set of four point sources S (S1, S2, S3 and S4) is arranged inside the volume defined by the wing A. Generally, at least two sources S can be arranged inside the volume defined by the wing A. In FIG. 1, two sources S are arranged at the centre of the wing A: the source S1 which is a scalar source S associated with an unknown scalar mass flow rate λ.sub.S and the source S2, which is a vector source, associated with an unknown vector mass flow rate λ.sub.R.

(23) Generally, at least two sources S are arranged inside the volume, including a source S called scalar source and a source S called vector source. Two other sources are shown in FIG. 1; S3 and S4 are sources arranged at two ends of the wing A. S3 and S4 are contrarotating vector sources S, i.e. the vector mass flow rates λ.sub.R that are respectively associated therewith have opposite directions. In the invention, the vector mass flow rates λ.sub.R associated with the sources S3 and S4 can be oriented parallel to the direction of the main velocity vector V, the two sources S forming a direction normal to the direction of the main velocity vector V. The advantages of a modelling shown in FIG. 1 are described hereinafter.

(24) FIG. 2 shows a method for the construction of a model according to an embodiment of the invention. Generally, during the construction of a model according to the invention, the body C or the wing A is associated with at least two aerodynamic parameters. These parameters can be chosen from a drag force coefficient, a lift coefficient, a roll coefficient, a pitch coefficient and a yaw coefficient. These coefficients correspond to parameters of the wing A determined in the state of the art for known spatial and aerodynamic configurations of the wing A, for example determined experimentally beforehand, in a known fluid flow FLU. These parameters can correspond to parameters of the nomographs used by the nomograph pointers in the embodiments of the prior art.

(25) The sources S (S1, S2, S3, and S4) shown in FIG. 2 correspond to the four sources shown in FIG. 1: the modelling comprises two sources S1 and S2 arranged at the centre (S1 being a scalar source S, associated with the scalar mass flow rate λ.sub.S and S2 being a vector source S, associated with the vector mass flow rate λ.sub.Ry. The modelling also comprises two other vector sources S, called boundary sources. S3 and S4. In the embodiments of the invention, the sources S3 and S4 are contrarotating and are respectively associated with the vector mass flow rates λ.sub.MTx and λ.sub.MBx. The source S3 (called starboard rotational boundary source) and the source S4 (called port rotational boundary source) are arranged on either side of the sources S1 and S2, at a distance b from the sources S1 and S2 (2.b corresponding to the span of the wing A). In the case shown in FIG. 2, the directions of the vector sources S3 and S4 are normal to the main velocity vector V shown. This condition is for example not verified in a modelling in which a wind component is a crosswind.

(26) Generally, for a wing A with surface area S, having a drag force coefficient C.sub.x and range coefficient C.sub.z, it is possible to derive the following lift force F.sub.Z and the drag force F.sub.X:
F.sub.z=½ρ.Math.V.sup.2.Math.S.Math.C.sub.z and  (1)
F.sub.X=½ρ.Math.V.sup.2.Math.S.Math.C.sub.X.  (2)
where ρ is the mass density of the fluid FLU in Kg.Math.m.sup.−3. The vector mass flow rate λ.sub.R can be written in the form λ.sub.R=ρ.Math.(ω.Math.dS).Math.dl, ω being the vorticity in s.sup.−1 and dl a vector length element.

(27) By analogy with electromagnetic forces, the lift force F.sub.Z and the drag force F.sub.X can be expressed as follows:
F.sub.Z=V∧λ.sub.R  (3)
F.sub.X=V.Math.λ.sub.S  (4)

(28) The equation (3), corresponding to the source S2, can be written in matrix form:

(29) $\begin{array}{cc}\stackrel{.fwdarw.}{F}=\left(\begin{array}{c}{f}_x\\ f_y\\ f_z\end{array}\right)=\left(\begin{array}{c}{v}_y.Math.{\lambda }_{\mathrm{Rz}}-v_z.Math.{\lambda }_{\mathrm{Ry}}\\ v_z.Math.{\lambda }_{\mathrm{Rx}}-v_x.Math.{\lambda }_{\mathrm{Rz}}\\ v_x.Math.{\lambda }_{\mathrm{Ry}}-v_y.Math.{\lambda }_{\mathrm{Rx}}\end{array}\right)=\left(\begin{array}{ccc}0& -v_z& +v_y\\ +v_z& 0& -v_x\\ -v_y& +v_x& 0\end{array}\right).Math.\left(\begin{array}{c}{\lambda }_{\mathrm{Rx}}\\ {\lambda }_{\mathrm{Ry}}\\ {\lambda }_{\mathrm{Rz}}\end{array}\right)& \left(5\right)\end{array}$

(30) In the modelling shown in FIG. 2, the main velocity vector V is orientated along the axis x shown, such that V.sub.y=V.sub.z=0. In the example shown, non-limitatively, the vector mass flow rate associated with S2 is oriented along y, the components along x and z being zero. In this case:

(31) $\begin{array}{cc}\left(\begin{array}{c}{f}_x\\ f_y\\ f_z\end{array}\right)=\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& -v_x\\ 0& +v_x& 0\end{array}\right).Math.\left(\begin{array}{c}0\\ {\lambda }_{\mathrm{Ry}}\\ 0\end{array}\right)=\left(\begin{array}{c}0\\ 0\\ v_x{\lambda }_{\mathrm{Ry}}\end{array}\right)& \left(6\right)\end{array}$

(32) The lift force in the equation (6) is reduced to f.sub.z=v.sub.x.Math.λ.sub.Ry.

(33) The equation (4), corresponding to the source S1, can be written in matrix form:

(34) $\begin{array}{cc}\stackrel{.fwdarw.}{F}=\left(\begin{array}{c}{f}_x\\ f_y\\ f_z\end{array}\right)=\left(\begin{array}{c}{v}_x\\ v_y\\ v_z\end{array}\right).Math.{\lambda }_S=\left(\begin{array}{c}{v}_x.Math.{\lambda }_S\\ v_y.Math.{\lambda }_S\\ v_z.Math.{\lambda }_S\end{array}\right)=\left(\begin{array}{ccc}{v}_x& 0& 0\\ 0& v_y& 0\\ 0& 0& v_z\end{array}\right).Math.\left(\begin{array}{c}{\lambda }_S\\ {\lambda }_S\\ {\lambda }_S\end{array}\right)& \left(7\right)\end{array}$

(35) The drag force can be reduced in this modelling example to f.sub.x=v.sub.x.Math.λ.sub.S: the force F has a single x (drag force) component oriented along x if V has only one component along x.

(36) The preceding equations show that, considering the sources S1 and S2, it is possible to model forces parallel and/or normal to the fluid flow FLU having main velocity vector V by using respectively the properties of the scalar and rotational sources S. The resulting force applied at the point corresponding to S1 and S2, associated with {λ.sub.Rr, λ.sub.Ry, λ.sub.Rx, λ.sub.S} is:

(37) $\begin{array}{cc}\stackrel{.fwdarw.}{F}=\left(\begin{array}{ccc}0& -v_z& +v_y\\ +v_z& 0& -v_x\\ -v_y& +v_x& 0\end{array}\right).Math.\left(\begin{array}{c}{\lambda }_{\mathrm{Rx}}\\ {\lambda }_{\mathrm{Ry}}\\ {\lambda }_{\mathrm{Rz}}\end{array}\right)+\left(\begin{array}{ccc}{v}_x& 0& 0\\ 0& v_y& 0\\ 0& 0& v_z\end{array}\right).Math.\left(\begin{array}{c}{\lambda }_S\\ {\lambda }_S\\ {\lambda }_S\end{array}\right)& \left(8\right)\end{array}$$\begin{array}{cc}i.e.:& \end{array}$$\begin{array}{cc}\stackrel{.fwdarw.}{F}=\left(\begin{array}{ccc}{\lambda }_S& {\lambda }_{\mathrm{Rz}}& -{\lambda }_{\mathrm{Ry}}\\ -{\lambda }_{\mathrm{Rz}}& {\lambda }_S& {\lambda }_{\mathrm{Rx}}\\ {\lambda }_{\mathrm{Ry}}& -{\lambda }_{\mathrm{Rx}}& {\lambda }_S\end{array}\right).Math.\left(\begin{array}{c}{v}_x\\ v_y\\ v_z\end{array}\right)& \left(9\right)\end{array}$

(38) It is also possible to derive the analytical expression of the velocities created by a scalar source S (isotropic. Green's function) associated with a scalar mass flow rate λ.sub.S. For example a scalar source S is considered, with coordinate r.sub.0, radiating at a coordinate point r. The relative coordinate R thereof can be defined as follows:

(39) $\begin{array}{cc}\stackrel{.fwdarw.}{R}=\stackrel{.fwdarw.}{r}-{\stackrel{.fwdarw.}{r}}_0=\left(\begin{array}{c}x-x_0\\ y-y_0\\ z-z_0\end{array}\right)=\left(\begin{array}{c}{r}_x\\ r_y\\ r_z\end{array}\right)& \left(14\right)\end{array}$$\begin{array}{cc}\mathrm{with}:& \end{array}$$\begin{array}{cc}R={\left(r_x^2+r_y^2+r_z^2\right)}^{\frac{1}{2}}& \left(15\right)\end{array}$

(40) The velocity radiated at the coordinate point r is then given by:

(41) $\begin{array}{cc}\stackrel{.fwdarw.}{V}=\left(\begin{array}{c}{v}_x\\ v_y\\ v_z\end{array}\right)=\frac{1}{4.Math.\pi .Math.\rho .Math.R^3}.Math.\left(\begin{array}{ccc}+r_x& 0& 0\\ 0& +r_y& 0\\ 0& 0& +r_z\end{array}\right).Math.\left(\begin{array}{c}{\lambda }_S\\ {\lambda }_S\\ {\lambda }_S\end{array}\right)& \left(16\right)\end{array}$

(42) The scalar source S associated with the mass flow rate λ.sub.S radiates isotropically.

(43) Similarly, it is possible to consider a vector source S with coordinate r.sub.0, associated with a vector mass flow rate λ.sub.R, and derive the expression of the velocity radiated at a coordinate point r. The equations (14) and (15) at still verified and:

(44) $\begin{array}{cc}\stackrel{.fwdarw.}{V}=\left(\begin{array}{c}{v}_x\\ v_y\\ v_z\end{array}\right)=\frac{{\stackrel{.fwdarw.}{\lambda }}_R\wedge \stackrel{.fwdarw.}{R}}{4.Math.\pi .Math.R^3}=\left(\begin{array}{c}{\lambda }_{\mathrm{Rx}}\\ {\lambda }_{\mathrm{Ry}}\\ {\lambda }_{\mathrm{Rz}}\end{array}\right)\wedge \left(\begin{array}{c}{r}_x\\ r_y\\ r_z\end{array}\right)& \left(17\right)\end{array}$
and therefore:

(45) $\begin{array}{cc}\stackrel{.fwdarw.}{V}=\left(\begin{array}{c}{v}_x\\ v_y\\ v_z\end{array}\right)=\frac{1}{4.Math.\pi .Math.\rho .Math.R^3}.Math.\left(\begin{array}{c}{\lambda }_{\mathrm{Ry}}.Math.r_z-{\lambda }_{\mathrm{Ry}}.Math.r_y\\ {\lambda }_{\mathrm{Rz}}.Math.r_x-{\lambda }_{\mathrm{Rx}}.Math.r_z\\ {\lambda }_{\mathrm{Rx}}.Math.r_y-{\lambda }_{\mathrm{Ry}}.Math.r_x\end{array}\right)& \left(18\right)\end{array}$

(46) If r.sub.0 is the coordinate of a first vector source S with respect to a second vector source, the first vector source radiates the following velocity on the second vector source:

(47) $\begin{array}{cc}\stackrel{.fwdarw.}{V}=\left(\begin{array}{c}{v}_x\\ v_y\\ v_z\end{array}\right)=\frac{1}{4.Math.\pi .Math.\rho .Math.R^3}.Math.\left(\begin{array}{ccc}0& +r_z& -r_y\\ -r_z& 0& +r_x\\ +r_y& -r_x& 0\end{array}\right).Math.\left(\begin{array}{c}{\lambda }_{\mathrm{Rx}}\\ {\lambda }_{\mathrm{Ry}}\\ {\lambda }_{\mathrm{Rz}}\end{array}\right)& \left(19\right)\end{array}$
and radiates the following force:

(48) 0$\begin{array}{cc}\stackrel{.fwdarw.}{F}=\left(\begin{array}{c}{f}_x\\ f_y\\ f_z\end{array}\right)=\left(\begin{array}{ccc}0& +{\lambda }_{R2z}& -{\lambda }_{R2y}\\ -{\lambda }_{R2z}& 0& +{\lambda }_{R2x}\\ +{\lambda }_{R2y}& -{\lambda }_{R2x}& 0\end{array}\right).Math.\left(\begin{array}{c}{v}_x\\ v_y\\ v_z\end{array}\right)& \left(20\right)\end{array}$

(49) In the case where λ.sub.R only comprises one x component (the case corresponding physically to wingtip vortices, shown by S3 and S4 in FIG. 1), the first source S results in a force, at the point of sources S1 and S2, at the centre of the wing, comprising only a single z component.

(50) In the case of the starboard wingtip vortex, (corresponding to the source S3) λ.sub.Rx is positive, r.sub.y is negative, and V.sub.z is negative.

(51) In the case of the port wingtip vortex, (corresponding to the source S4) λ.sub.RX is negative, r.sub.y is positive, and V.sub.z is negative.

(52) The velocity induced by the vector source S3 on the torque of sources S1 and S2 is therefore:

(53) $\begin{array}{cc}{v}_z=-\frac{1}{4.Math.\pi .Math.\rho .Math.b^2}.Math.{\lambda }_{\mathrm{MTx}}& \left(21\right)\end{array}$

(54) Similarly, the velocity induced by the source S4 on the torque of sources S1 and S2 is:

(55) $\begin{array}{cc}{v}_z=+\frac{1}{4.Math.\pi .Math.\rho .Math.b^2}.Math.{\lambda }_{\mathrm{MBx}}& \left(22\right)\end{array}$

(56) The different components of the forces induced by the sources S3 and S4 on S1 can be deduced therefrom:

(57) $\begin{array}{cc}{\stackrel{.fwdarw.}{F}}_z=-\frac{{\lambda }_S}{4.Math.\pi .Math.\rho .Math.b^2}.Math.\left({\lambda }_{\mathrm{MTx}}-{\lambda }_{\mathrm{MBx}}\right)& \left(23\right)\end{array}$
and on S2:

(58) $\begin{array}{cc}{\stackrel{.fwdarw.}{F}}_x=+\frac{{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho .Math.b^2}.Math.\left({\lambda }_{\mathrm{MTx}}-{\lambda }_{\mathrm{MBx}}\right)& \left(24\right)\end{array}$

(59) The mass flow rates associated with the different sources are proportional to the velocity:
λ.sub.Ry=k.sub.Ry.Math.V  (25)
λ.sub.S=k.sub.S.Math.V  (26)
λ.sub.MTx=k.sub.MTx.Math.V  (27)
λ.sub.MBx=k.sub.MBx.Math.V  (28)

(60) Symmetrical flight of the wing A results in the condition:
k.sub.MBx=−k.sub.MTx  (29)

(61) Preservation of the vorticity results in:
k.sub.MTx=k.sub.Ry=−k.sub.MBx  (30)

(62) The lift can then be written in the two following ways:

(63) $\begin{array}{cc}{\stackrel{.fwdarw.}{F}}_z=k_{\mathrm{Ry}}.Math.V_x^2-\frac{k_S}{4.Math.\pi .Math.\rho .Math.b^2}.Math.\left(2.Math.k_{\mathrm{MTx}}\right).Math.V_x^2=\frac{1}{2}.Math.\rho .Math.V_x^2.Math.S.Math.C_z& \left(31\right)\end{array}$

(64) Similarly, the drag force can be written in the two following ways:

(65) $\begin{array}{cc}{\stackrel{.fwdarw.}{F}}_x=k_S.Math.V_x^2-\frac{k_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho .Math.b^2}.Math.\left(2.Math.k_{\mathrm{MTx}}\right).Math.V_x^2=\frac{1}{2}.Math.\rho .Math.V_x^2.Math.S.Math.C_x& \left(32\right)\end{array}$

(66) In the model shown, the sources S3 and S4, modelling the wingtip vortices, result in a reduction in the lift (induced incidence effect) and an increase in the drag force (induced drag force effect). In the absence of the sources S3 and S4, the terms corresponding to a typical lift and an atypical drag force would be found.

(67) With the following notation:
B=4.Math.π.Math.ρ.Math.b.sup.2 and  (33)
A=½.Math.ρ.Math.S  (34)

(68) The expressions for the lift and drag force can be written respectively as follows:

(69) $\begin{array}{cc}{\stackrel{.fwdarw.}{F}}_z=k_{\mathrm{Ry}}.Math.\left[1-\frac{k_S}{B}.Math.\left(2\right)\right].Math.V_x^2=A.Math.V_x^2.Math.C_z\mathrm{and}& \left(35\right)\\ {\stackrel{.fwdarw.}{F}}_x=1.Math.\left[k_S-\frac{k_{\mathrm{Ry}}^2}{B}.Math.\left(2\right)\right].Math.V_x^2=A.Math.V_x^2.Math.C_x.& \left(36\right)\end{array}$

(70) The lift and drag force coefficients can then be expressed as follows:

(71) $\begin{array}{cc}{C}_z=\frac{1}{A}.Math.k_{\mathrm{Ry}}.Math.\left[1-2.Math.\frac{k_S}{B}\right]\mathrm{and}& \left(37\right)\\ C_x=\frac{1}{A}.Math.\left[k_S+2.Math.\frac{k_{\mathrm{Ry}}^2}{B}\right]& \left(38\right)\end{array}$
and therefore the expression of k.sub.Ry:

(72) $\begin{array}{cc}{k}_{\mathrm{Ry}}=\frac{A.Math.C_z}{\left[1-2.Math.k_s/B\right]}& \left(39\right)\end{array}$
The expression for the drag force coefficient becomes:

(73) 0$\begin{array}{cc}{C}_x=\frac{1}{A}.Math.\left[k_S+\frac{2.Math.A^2.Math.{\mathrm{Cz}}^2}{B.Math.{\left(1-2.Math.k_S/B\right)}^2}\right]& \left(40\right)\end{array}$

(74) It is possible to solve the equation (40) numerically so as to identify k.sub.s and subsequently the different values for the other mass flow rates associated with the other sources of the model. Generally, a step of a method according to the invention consists of determining the value of each scalar mass flow rate λ.sub.S and of each vector mass flow rate λ.sub.R as a function of a part, or the totality, of the boundary conditions imposed by at least one main velocity vector V and of at least one aerodynamic parameter, such as the drag force coefficient and the lift coefficient in the example shown in FIG. 2, in a given configuration. Advantageously, the values of each scalar mass flow rate λ.sub.S and of each vector mass flow rate λ.sub.R can be determined within the general framework of a resolution according to the distributed point source method (DPSM).

(75) By simplifying the equation (40) to the first order, it is possible to write:

(76) $\begin{array}{cc}{C}_x=\frac{k_S}{A}+\frac{2.Math.A.Math.{\mathrm{Cz}}^2}{B}=\frac{k_S}{A}+\frac{2.Math.\left(\frac{1}{2}.Math.\rho .Math.S\right).Math.{\mathrm{Cz}}^2}{4.Math.\pi .Math.\rho .Math.b^2}& \left(41\right)\end{array}$
and therefore:

(77) $\begin{array}{cc}{C}_x=\frac{k_S}{A}+\frac{S.Math.{\mathrm{Cz}}^2}{4.Math.\pi .Math.b^2}=\frac{k_S}{A}+\frac{{\mathrm{Cz}}^2}{\pi .Math.\left(4.Math.b^2/S\right)}& \left(42\right)\end{array}$

(78) The aspect ratio All is defined by the ratio of the span to the chord of a wing A, i.e. in the model shown in FIG. 2:
All=(2.Math.b)/(S/2.Math.b)=4.Math.b.sup.2/S  (43)

(79) The equation (42) becomes:

(80) $\begin{array}{cc}{C}_x=\frac{k_S}{A}+\frac{{\mathrm{Cz}}^2}{\pi .Math.\left(4.Math.b^2/S\right)}={\mathrm{Cx}}_0+\frac{{\mathrm{Cz}}^2}{\pi .Math.\left(\mathrm{All}\right)}& \left(44\right)\end{array}$

(81) The equation (44) makes it possible to identify a term similar to the induced drag force of Prandtl's theory, added to C.sub.x0.

(82) Similarly, it is possible to express C.sub.x as a function of C.sub.z:

(83) $\begin{array}{cc}{C}_x=\frac{1}{A}.Math.\left[k_S+2.Math.\frac{k_{\mathrm{Ry}}^2}{B}\right]& \left(45\right)\\ C_z=\frac{1}{A}.Math.k_{\mathrm{Ry}}.Math.\left[1-2.Math.\frac{k_S}{B}\right]& \left(46\right)\\ k_S=A.Math.C_x-\frac{2.Math.k_{\mathrm{Ry}}^2}{B}& \left(47\right)\\ C_z=\frac{1}{A}.Math.k_{\mathrm{Ry}}.Math.\left[1-\frac{2.Math.A.Math.\mathrm{Cx}}{B}+\frac{4.Math.k_{\mathrm{Ry}}^2}{B^2}\right]& \left(48\right)\end{array}$

(84) On the basis of a dimensional analysis, it is possible to show that:

(85) $\begin{array}{cc}4.Math.\frac{k_{\mathrm{Ry}}^2}{B^2}\cong {\left(2.Math.\frac{A}{B}\right)}^2\cong {\left(\frac{1}{\pi .Math.\mathrm{All}}\right)}^2<<\left(2.Math.\frac{A}{B}\right).Math.\mathrm{Cx}& \left(49\right)\end{array}$

(86) It is possible to simplify the equation (48) to the first order as follows:

(87) $\begin{array}{cc}{C}_z=\frac{1}{A}.Math.k_{\mathrm{Ry}}.Math.\left[1-\frac{2.Math.A.Math.\mathrm{Cx}}{B}\right]=\frac{k_{\mathrm{Ry}}}{A}\left[1-\frac{2.Math.\left(\frac{1}{2}.Math.\rho .Math.S\right).Math.\mathrm{Cx}}{4.Math.\pi .Math.\rho .Math.b^2}\right]& \left(50\right)\\ i.e.\text{:}& \\ C_z=\frac{k_{\mathrm{Ry}}}{A}\left[1-\frac{S.Math.\mathrm{Cx}}{4.Math.\pi .Math.b^2}\right]=\frac{k_{\mathrm{Ry}}}{A}\left[1-\frac{\mathrm{Cx}}{\pi .Math.\left(4.Math.b^2/S\right)}\right],& \left(51\right)\\ C_z=\frac{k_{\mathrm{Ry}}}{A}\left[1-\frac{\mathrm{Cx}}{\pi .Math.\left(4.Math.b^2/S\right)}\right]={\mathrm{Cz}}_0\left[1-\frac{\mathrm{Cx}}{\pi .Math.\mathrm{All}}\right]& \left(52\right)\end{array}$

(88) It is possible to identify a corrective term in the equation (42), dependent on C.sub.x, which reduces the lift C.sub.z0.

(89) Another embodiment of the invention is a method, implemented by a computer in real-time, for simulating at least one force F or a moment M, applied to a body C, modelled by a volume described above, and using a model the construction of which has been described above. i.e. a model in which the fluid flow FLU comprises at least one main velocity vector V representing a velocity of the flow at infinity, in which the body C is modelled beforehand by at least two point sources S of a velocity field, including at least one scalar source S, generating a radial velocity field and associated with a scalar mass flow rate λ.sub.S, and at least one vector source S, generating a solenoidal velocity field, associated with a vector mass flow rate λ.sub.R. It is possible to calculate, beforehand, for at least one configuration of the body C, the values of each scalar mass flow rate λ.sub.S and of each vector mass flow rate λ.sub.R as described above. An embodiment of the invention comprises at least one step consisting of: simulating a superposition of the contributions of the fluid flows FLU of at least one main velocity flow V, of a scalar source S associated with a scalar mass flow rate λ.sub.S and of a vector source S associated with a vector mass flow rate λ.sub.R, and/or evaluating at least one element chosen from a force F and a moment M, by a linear function dependent on the main velocity vector V, of at least one scalar mass flow rate λ.sub.S and at least one said vector mass flow rate λ.sub.R, at at least one point corresponding to one said source S. In an embodiment in FIG. 2, this step can for example consist of calculating one of the forces, or the two forces corresponding to the equations (34) and/or (35). This embodiment of the invention makes it possible to solve at least a part of the technical problems raised: the inventors discovered that by using a method according to the invention, it is possible to simulate a force and/or a moment by evaluating a simple linear expression. This calculation can be carried out by a computer in less than 6 ms, and makes it possible to use this model in order to simulate the flight of a wing A by using the calculation power of an ordinary computer at the time of the invention. i.e. one with performance that is not exceptional at the time of the invention.

(90) Advantageously, the simulation described above is a simulation carried out in far field.

(91) Another embodiment of the invention is a computer program product, stored on a recording medium suitable for implementing a method described above. Another embodiment of the invention is a recording medium on which such a program is stored.

(92) Another embodiment of the invention is a flight simulator SIM, comprising a calculator, for example incorporated in a computer, the calculator being configured in order to implement a method according to the invention described above. Such a flight simulator SIM can for example be configured in order to receive first data, acquired by a user of the simulator. A memory can comprise second data representative of the fluid flow FLU, updated in real time. These second data can also comprise values of the different mass flow rates λ.sub.S and λ.sub.R, calculated beforehand by a method for the construction of the model according to the invention. They can also comprise nomographs. In this case, a simulator SIM can choose to evaluate forces applied to a wing A below a height defined beforehand for the wing A in the simulation. Below this height, these forces can be evaluated by nomograph pointers according to the prior art. In a variant, a calculator of a second computer connected to a first computer of the simulator can be responsible for the simulation of a force F or of a moment 9M sequentially with the simulation carried out according to the state of the art by the first computer.

(93) FIG. 3 shows the construction of a model according to the invention adapted to a simulation of the effect of an element of the environment external to the body C or to the wing A. Generally, the inventors discovered that it was possible to model the effect of an external element on a model of an existing body C, by adding to the models described above, at least one scalar external source S′, associated with a mirror scalar mass flow rate λ.sub.S′. Advantageously, it is also possible to add at least one vector external source S′, associated with a mirror vector mass flow rate λ.sub.R′, each of the external sources S′ being arranged outside the volume of the body C or of the wing A. FIG. 3 shows the construction of a model adapted to a simulation of the external effect, comprising sources S of an embodiment shown in FIG. 2, and external sources S′ modelling a premises or a building. More generally, different systems of external sources S′ can be arranged so as to simulate for example a naval vessel (in the case of landing at sea), buildings on land, the ground, or even particular crosswinds. In this embodiment of the invention, the contributions of the fluid flows FLU of each source S′ are also superposed during step (a) of the method, and in step (b) of the method an analytic linear function is evaluated, also as a function of at least one external mass flow rate λ.sub.S′ and/or of at least one external vector mass flow rate λ.sub.R′. An embodiment of such a model is given in the remainder of the description, in particular in modelling the ground effect.

(94) FIG. 4 shows a wing A in a configuration able to be affected by a ground effect. Experimentally, this configuration modifies the lift and the drag force of a wing A, close to the ground, with respect to a configuration of the wing A at a cruising height. In FIG. 4, it is assumed that the ground coincides locally with a plane GND. The wing A is shown above the ground GND, in a first half-space a limit of which is formed by the plane GND. Conventionally, the ground effect on a wing A is simulated, experimentally or numerically, by replacing the ground GND by a mirror wing A′, symmetrical with the wing A with respect to the plane GND.

(95) FIG. 5 shows the construction of a model according to the invention adapted to a simulation of the ground effect. Generally, the inventors discovered that it was possible to model the ground effect on a model of a wing, adding to the models described above at least one mirror scalar external source S′, associated with a mirror scalar mass flow rate λ.sub.S′ and at least one mirror vector external source S′, associated with a mirror vector mass flow rate λ.sub.R′, each of the mirror external sources S′ being arranged in a second half-space defined by the plane GND and different from the first half-space, in which the sources S are arranged. FIG. 5 shows the construction of a model adapted to a simulation of the ground effect, comprising sources S of an embodiment shown in FIG. 2, and mirror external sources S′ symmetrical with the sources S with respect to the plane GND, each of the sources S′ being separated from the symmetrical source S with respect to the plane GND by a distance a=2.h.

(96) The sources S1′, S2′, S3′ and S4′ are respectively symmetrical, with respect to the plane GND, with the sources S1, S2, S3 and S4, and are respectively associated with the mass flow rates λ.sub.S,2, λ.sub.Ry,2, λ.sub.MTx,2 and λ.sub.MBx,2. The conditions of symmetry also define:
λ.sub.MBx,2=−λ.sub.MTx,2=λ.sub.MTx=−λ.sub.MBx=λ.sub.Ry  (53)
λ.sub.Ry,2=−λ.sub.Ry and  (54)
λ.sub.S,2=λ.sub.S.  (55)

(97) It is possible to calculate the influence of the model of the mirror wing A′ comprising the mirror external sources S′ on the forces exerted on the wing A′. In the following calculation, by “centre” is meant the point of sources S1 and S2. The vector R, with norm equal to a, between the sources S1′ and S2′ and S1 and S2 is defined by:

(98) $\begin{array}{cc}\stackrel{.fwdarw.}{R}=\left(\begin{array}{c}{r}_x\\ r_y\\ r_z\end{array}\right)=\left(\begin{array}{c}0\\ 0\\ a\end{array}\right)& \left(56\right)\end{array}$

(99) The contribution of the vector source S2′ can be expressed as follows:

(100) $\begin{array}{cc}{\stackrel{.fwdarw.}{V}}_{R,\mathrm{Centre},2}=\left(\begin{array}{c}{v}_x\\ v_y\\ v_z\end{array}\right)=\frac{1}{4.Math.\pi .Math.\rho .Math.R^3}.Math.\left(\begin{array}{ccc}0& +r_z& -r_y\\ -r_z& 0& +r_x\\ +r_y& -r_x& 0\end{array}\right).Math.\left(\begin{array}{c}{\lambda }_{\mathrm{Rx},2}\\ {\lambda }_{{\mathrm{Ry}}_{,2}}\\ {\lambda }_{{\mathrm{Rz}}_{,2}}\end{array}\right)& \left(57\right)\\ i.e.\text{:}& \\ {\stackrel{.fwdarw.}{V}}_{R,\mathrm{Centre},2}=\frac{1}{4.Math.\pi .Math.\rho .Math.R^3}.Math.\left(\begin{array}{ccc}0& a& 0\\ -a& 0& 0\\ 0& 0& 0\end{array}\right).Math.\left(\begin{array}{c}0\\ {\lambda }_{{\mathrm{Ry}}_{,2}}\\ 0\end{array}\right)=\frac{1}{4.Math.\pi .Math.\rho .Math.a^3}\left(\begin{array}{c}a.Math.{\lambda }_{{\mathrm{Ry}}_{,2}}\\ 0\\ 0\end{array}\right)=\frac{-{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho .Math.a^2}\stackrel{.fwdarw.}{e_x}& \left(58\right)\end{array}$

(101) The contribution of the scalar source S1′ can be expressed as follows:

(102) $\begin{array}{cc}{\stackrel{.fwdarw.}{V}}_{S,\mathrm{Centre},2}=\left(\begin{array}{c}{v}_x\\ v_y\\ v_z\end{array}\right)=\frac{1}{4.Math.\pi .Math.\rho .Math.R^3}.Math.\left(\begin{array}{ccc}{r}_x& 0& 0\\ 0& r_y& 0_x\\ 0& 0& r_z\end{array}\right).Math.\left(\begin{array}{c}{\lambda }_{S,2}\\ {\lambda }_{S_{,2}}\\ {\lambda }_{S_{,2}}\end{array}\right)& \left(59\right)\end{array}$$\begin{array}{cc}i.e.:& \end{array}$$\begin{array}{cc}{\stackrel{.fwdarw.}{V}}_{S,\mathrm{Centre},2}=\frac{1}{4.Math.\pi .Math.\rho .Math.a^3}.Math.\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 0\\ 0& 0& a\end{array}\right).Math.\left(\begin{array}{c}{\lambda }_{S,2}\\ {\lambda }_{S_{,2}}\\ {\lambda }_{S_{,2}}\end{array}\right)=\left(\begin{array}{c}0\\ 0\\ \frac{{\lambda }_{S,2}}{4.Math.\pi .Math.\rho .Math.a^2}\end{array}\right)=\frac{{\lambda }_S}{4.Math.\pi .Math.\rho .Math.a^2}\stackrel{.fwdarw.}{e_z}& \left(60\right)\end{array}$

(103) The mirror external sources S3′ and S4′, corresponding to the image wingtip vortices, make the following contribution:

(104) 0$\begin{array}{cc}{\stackrel{.fwdarw.}{V}}_{\mathrm{MT},\mathrm{Centre}}=\left(\begin{array}{c}0\\ 0\\ \frac{-{\lambda }_{\mathrm{MTx}}}{4.Math.\pi .Math.\rho .Math.b^2}\end{array}\right)=\frac{-{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho .Math.b^2}\stackrel{.fwdarw.}{e_z}& \left(61\right)\end{array}$

(105) The sum of the contributions of the sources S and the mirror external sources S′ at the centre, corresponding to the source S2 can be expressed by:

(106) $\begin{array}{cc}{v}_x=\frac{-{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho .Math.a^2},& \left(62\right)\end{array}$$\begin{array}{cc}{v}_y=0\mathrm{and}& \left(63\right)\end{array}$$\begin{array}{cc}{v}_z=\frac{{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho }\left(\frac{-2}{b^2}+\frac{2.Math.b}{{\sqrt{a^2+b^2}}^3}\right)+\frac{{\lambda }_S}{4.Math.\pi .Math.\rho .Math.a^2}& \left(64\right)\end{array}$

(107) It is possible to calculate the forces linked to the velocities induced on the source S1 on the basis of the equation (7):

(108) $\begin{array}{cc}{F}_x=\frac{-{\lambda }_S.Math.{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho .Math.a^2},& \left(65\right)\end{array}$$\begin{array}{cc}{F}_y=0\mathrm{and}& \left(66\right)\end{array}$$\begin{array}{cc}{F}_z=\frac{{\lambda }_S.Math.{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho }\left(\frac{-2}{b^2}+\frac{2.Math.b}{{\sqrt{a^2+b^2}}^3}\right)+\frac{{\lambda }_S.Math.{\lambda }_S}{4.Math.\pi .Math.\rho .Math.a^2}& \left(67\right)\end{array}$

(109) Similarly, it is possible to calculate the forces linked to the velocities induced on the source S2 on the basis of the equation (5):

(110) $\begin{array}{cc}{F}_x=\frac{{\lambda }_{\mathrm{Ry}}.Math.{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho }\left(\frac{-2}{b^2}+\frac{2.Math.b}{{\sqrt{a^2+b^2}}^3}\right)-\frac{{\lambda }_S.Math.{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho .Math.a^2}& \left(68\right)\end{array}$$\begin{array}{cc}{F}_y=0& \left(69\right)\end{array}$$\begin{array}{cc}{F}_z=\frac{-{\lambda }_{\mathrm{Ry}}.Math.{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho .Math.a^2}& \left(70\right)\end{array}$

(111) A sum of the forces applied at the point of sources S1 and S2 of the wing A can be written in the form:

(112) $\begin{array}{cc}{F}_x=-\frac{{\lambda }_{\mathrm{Ry}}.Math.{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho }\left(\frac{-2}{b^2}+\frac{2.Math.b}{{\sqrt{a^2+b^2}}^3}\right)-\frac{2.Math.{\lambda }_S.Math.{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho .Math.a^2}& \left(71\right)\end{array}$$\begin{array}{cc}{F}_y=0& \left(72\right)\end{array}$$\begin{array}{cc}{F}_x=\frac{{\lambda }_S.Math.{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho }\left(\frac{-2}{b^2}+\frac{2.Math.b}{{\sqrt{a^2+b^2}}^3}\right)+\frac{{\lambda }_S.Math.{\lambda }_S}{4.Math.\pi .Math.\rho .Math.a^2}-\frac{{\lambda }_{\mathrm{Ry}}.Math.{\lambda }_{\mathrm{Ry}}}{4.Math.\pi .Math.\rho .Math.a^2}& \left(73\right)\end{array}$

(113) It is for example possible to choose a condition a=b (the ground effect exists when the height of the aircraft is less than one third of the span, corresponding to h=(2.Math.b)/4 and

(114) ${\lambda }_S=\frac{{\lambda }_{\mathrm{Ry}}}{5}$
hypothetically). By defining

(115) $K^2=.Math.\frac{{\lambda }_{\mathrm{Ry}}^2}{4.Math.\pi .Math.\rho b^2}.Math.,$
it is possible to summarize the different contributions of the forces using the following table:

(116) TABLE-US-00001 S3, S4 S1’ S2’ S3’, S4’ F.sub.x F = K.sup.2 (=1) $F=\frac{K^2}{5}$ (=0.2) $F=\frac{K^2}{5}$ (=0.2) $F=\frac{K^2}{{\left(\sqrt{2}\right)}^3}$ (=0.35) F.sub.y e Zero zero zero F.sub.z 0$F=\frac{K^2}{5}$ (=0.2) $F=\frac{K^2}{25}$ (=0.04) F = K.sup.2 (=1) $F=\frac{K^2}{5.{\left(\sqrt{2}\right)}^3}$ (=0.07)

(117) The majority of the variations in the lift and drag force originate from the two rotational sources situated at the centre of the wing and of the image wing. The variation in the drag force is much greater than the variation in the lift when the aircraft approaches the ground. Close to the ground, one wing A appears experimentally to glide, as the drag force is the variable that reduces the most, essentially as a result of the reduction in the effect of the wingtip vortices of the wing A approaching the wingtip vortices of the mirror wing A′. In a method for the simulation of a force according to the invention, in a configuration of the wing A close to the ground, the contributions of the fluid flows FLU of the sources S′ are also superposed during a step of simulation of the fluid FLU around the body C, and in a step of evaluation of a force and/or a moment, a linear function is evaluated, also as a function of at least one mirror mass flow rate λ.sub.S of at least one mirror vector mass flow rate λ.sub.R′.

(118) FIG. 6 shows a simulation of the lift coefficient according to the invention. The data from a simulator SIM contain the values of different aerodynamic parameters of the wing A, for example the values of Cx and Cz for different angles of incidence of the wing A (in data tables). Based on these data, the values of Cx and Cz can be amended as a function of the height of the aircraft with respect to the ground, and of the results of the evaluation of the aerodynamic parameters simulated according to the invention. The data thus evaluated are weighted, so as to be adjusted to a real behaviour of a wing A.

(119) The diagram in FIG. 6 shows the drag force coefficients C.sub.z of a flight simulator according to the prior art (curves (b)) and the drag force coefficients C.sub.z evaluated according to the invention (curves (a)), using a model comprising mirror external sources, described above. The C.sub.z values are shown as a function of altitude, for altitude or height of the wing A and its angle of incidence.

(120) FIG. 7 shows the construction of a model according to the invention making it possible to simulate the Magnus effect and a shape drag force.

(121) Theoretically and generally, considering a cylinder in rotation, with an angular velocity ω.sub.cy1, of infinite length, the cylinder being considered according to its circular cross section, Poisson equations make it possible to derive the potential of the velocities in polar coordinates:

(122) $\begin{array}{cc}\varphi \left(r,\theta \right)=\frac{\cos \left(\theta \right)}{2}\left(r+\frac{1}{r}\right)-\frac{r}{2\pi }\theta .& \left(74\right)\end{array}$

(123) The streamlines can be described by:

(124) $\begin{array}{cc}\psi \left(r,\theta \right)=\frac{\sin \left(\theta \right)}{2}\left(r-\frac{1}{r}\right)+\frac{r}{2\pi }\ln \left(r\right),& \left(75\right)\end{array}$
and the radial and angular velocities, respectively by:

(125) $\begin{array}{cc}{v}_r=\frac{\partial \varphi }{\partial r}=\frac{\cos \left(\theta \right)}{2}\left(1-\frac{1}{r^2}\right)& \left(76\right)\end{array}$$\begin{array}{cc}{v}_{\theta }=\frac{1}{r}\frac{\partial \varphi }{\partial \theta }=-\frac{\sin \left(\theta \right)}{2}\left(1+\frac{1}{r^2}\right)-\frac{r}{2\pi .Math.r}.& \left(77\right)\end{array}$
In this case the lift force is:
F.sub.p=−ρv.sub.0Γ  (78)

(126) In this theoretical framework, it is possible to express a lift force, but not a drag force. Outside this theoretical framework, as the cylinder is infinite, there is no induced drag force. There is however a shape drag force and a friction drag force.

(127) In an embodiment of the invention, a model is constructed in which two sources S are arranged so as to form a dipole. In embodiments of the invention, this dipole can be oriented, as shown in FIG. 7, in the direction of the main velocity vector V. Two sources S arranged in this manner are shown in FIG. 7 by the sources S5 and S6. These two sources are respectively associated with the scalar mass flow rates λ.sub.S5 and λ.sub.S6, and placed on either side of the point corresponding to the sources S1 and S2, and along an axis of the same direction as the main velocity vector V. The two sources S forming the dipoles can be associated with scalar mass flow rates λ.sub.S that are different and have opposite signs.

(128) In variants of the invention, the scalar mass flow rates λ.sub.S5 and λ.sub.S6 associated with the sources S5 and S6 can be equal or different.

(129) In a variant of the invention, when the mass flow rates λ.sub.S5 and λ.sub.S6 are different, a supplementary condition is necessary in order to make it possible to solve the systems of equations allowing the different scalar and/or vector mass flow rates associated with the sources S to be calculated. A condition corresponding to a zero normal flow velocity at a point of the interface of a virtual wing A or a virtual body C can be added to the modelling.

(130) FIG. 7 shows the construction of a model comprising a dipole of scalar sources S. According to the equations (1) and (2) it is possible to derive:

(131) $\begin{array}{cc}{\lambda }_R{V}_{\infty }=\frac{1}{2}\rho V_{\infty }{C}_{Z}S& \left(79\right)\end{array}$$\begin{array}{cc}{\lambda }_S{V}_{\infty }=\frac{1}{2}\rho V_{\infty }{C}_{X}S& \left(80\right)\end{array}$$\begin{array}{cc}i.e.:& \end{array}$$\begin{array}{cc}\{\begin{array}{c}{V}_{\infty }{C}_Z=\frac{1}{\frac{1}{2}\rho .Math.S}.Math.{\lambda }_R\\ V_{\infty }{C}_X=\frac{1}{\frac{1}{2}\rho .Math.S}.Math.{\lambda }_S\end{array}& \left(81\right),\left(82\right)\end{array}$

(132) A intercoupling matrix can be identified by deriving the following equation:

(133) $\begin{array}{cc}\left(\begin{array}{c}{V}_{\infty }{C}_Z\\ V_{\infty }{C}_X\end{array}\right)=\frac{1}{\frac{1}{2}\rho .Math.S.Math.}\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right).Math.\left(\begin{array}{c}{\lambda }_R\\ {\lambda }_S\end{array}\right)=\left(\begin{array}{cc}\frac{1}{\frac{1}{2}\rho .Math.S.Math.}& 0\\ 0& \frac{1}{\frac{1}{2}\rho .Math.S.Math.}\end{array}\right).Math.\left(\begin{array}{c}{\lambda }_R\\ {\lambda }_S\end{array}\right)& \left(83\right)\end{array}$

(134) FIG. 8 shows a circular cross section of an infinite cylinder, corresponding to a body C or to a wing A. Two scalar sources S5 and S6 are arranged in the same way as shown in FIG. 7. The two sources are separated by a distance d and the cylinder has a radius R. It is possible to derive an equation corresponding to a zero normal velocity condition at the interface of a body C or of a wing A.

(135) $\begin{array}{cc}{V}^n=0=\frac{\left(R-\frac{d}{2}\right).Math.{\lambda }_{S1}}{4.Math.\pi .Math.\rho .Math.{\left(R-\frac{d}{2}\right)}^3}+\frac{\left(R+\frac{d}{2}\right).Math.{\lambda }_{S2}}{4.Math.\pi .Math.\rho .Math.{\left(R+\frac{d}{2}\right)}^3}=\frac{{\lambda }_{S1}}{4.Math.\pi .Math.\rho .Math.{\left(R-\frac{d}{2}\right)}^2}+\frac{{\lambda }_{S2}}{4.Math.\pi .Math.\rho .Math.{\left(R+\frac{d}{2}\right)}^2}& \left(84\right)\end{array}$
A intercoupling matrix can be identified in the following equation:

(136) $\begin{array}{cc}\left(\begin{array}{c}0\\ V_{\infty }{C}_x\\ V_{\infty }{C}_z\end{array}\right)=\left(\begin{array}{ccc}\frac{1}{\frac{1}{2}\rho .Math.S}& 0& 0\\ 0& \frac{1}{\frac{1}{2}\rho .Math.S}& \frac{1}{\frac{1}{2}\rho .Math.S}\\ \frac{1}{4.Math.\pi .Math.\rho .Math.{R^2\left(1-\frac{d}{2R}\right)}^2}& \frac{1}{4.Math.\pi .Math.\rho .Math.{R^2\left(1+\frac{d}{2R}\right)}^2}& 0\end{array}\right).Math.\left(\begin{array}{c}{\lambda }_R\\ {\lambda }_{S1}\\ {\lambda }_{S2}\end{array}\right)=\left(\begin{array}{cccc}\frac{1}{\frac{1}{2}\rho .Math.S}|& \begin{array}{c}1\\ 0\end{array}& \begin{array}{c}0\\ 1\end{array}& \begin{array}{c}0\\ 1\end{array}\\ \frac{1}{4.Math.\pi .Math.\rho .Math.R^2}|& {\left(1-\frac{d}{2R}\right)}^{-2}& {\left(1+\frac{d}{2R}\right)}^{-2}& 0\end{array}\right).Math.\left(\begin{array}{c}{\lambda }_R\\ {\lambda }_{S1}\\ {\lambda }_{S2}\end{array}\right)& \left(85\right)\end{array}$

(137) FIG. 9 shows a simulation of a fluid flow FLU around a wing A, the shape of which is an infinite cylinder simulated by a model according to the invention in which two scalar sources S5 and S6 are arranged as a dipole. The scalar sources S5 and S6 are arranged as described in the embodiments of the invention corresponding to FIGS. 6 and 7. The different grey levels show the value of the fluid velocity FLU (the black corresponding to zero velocity and the white corresponding to maximum velocity). The black lines comprising multiple arrows correspond to streamlines. The single black arrows, in panels B. C and D in FIG. 9 correspond to the resultant of the forces calculated by a model constructed according to the invention and simulated according to the invention.

(138) The panel A in FIG. 9 shows a simulation according to the invention, using a model constructed according to the invention, comprising sources S1, S5 and S6 described above, associated respectively with the scalar λ.sub.S5, λ.sub.S6 and vector λ.sub.Ry mass flow rates. In this case, λ.sub.Ry=0 and λ.sub.S5=λ.sub.S6. The resultant cylinder of the simulation is fixed and without drag force. In this case, the lift and the drag force are zero: this case corresponding to d'Alembert's paradox.

(139) The panel B in FIG. 9 shows a simulation according to the invention, using a model constructed according to the invention, comprising sources S1, S5 and S6 described above, associated respectively with the scalar λ.sub.S5, λ.sub.S6 and vector λ.sub.Ry mass flow rates. In this case, λ.sub.Ry≠0 and λ.sub.S5=λ.sub.S6. The resultant cylinder of the simulation is rotational about its main axis, without drag force. In this case, the lift is non-zero and the drag force is zero: this case corresponds to the Magnus effect.

(140) The panel C in FIG. 9 shows a simulation according to the invention, using a model constructed according to the invention, comprising sources S1, S5 and S6 described above, associated respectively with the scalar λ.sub.S5, λ.sub.S6 and vector λ.sub.Ry mass flow rates. In this case, λ.sub.Ry=0 and λ.sub.S5≠λ.sub.S6. The resultant cylinder of the simulation is fixed and involves drag force. In this case, the lift is zero and the drag force is non-zero: this case corresponds to a shape drag force.

(141) The panel D in FIG. 9 shows a simulation according to the invention, using a model constructed according to the invention, comprising sources S1, S5 and S6 described above, associated respectively with the scalar λ.sub.S5, λ.sub.S6 and vector λ.sub.Ry mass flow rates. In this case. λ.sub.Ry≠0 and λ.sub.S5≠λ.sub.S6. The resultant cylinder of the simulation is rotational about its main axis, and involves a drag force. In this case, the lift is zero and the drag force is non-zero: this case approaches a real case.

(142) FIG. 10 shows drag force coefficients and lift coefficients comparatively evaluated according to a method of the prior art and according to the invention.

(143) The panel A in FIG. 10 shows on the one hand, drag force coefficients C.sub.X evaluated according to a method implemented by a flight simulator, according to a method of the prior art using nomograph pointers. These values of the drag force coefficients C.sub.X are shown by the continuous curves, for example the curve (a). Moreover, the panel A in FIG. 10 shows drag force coefficients C.sub.X evaluated according to the invention, and corresponding to the curves marked with stars, for example the curve (b). The C.sub.X evaluations according to the invention are obtained by weighting the scalar mass flow rates by a factor 0.25 and weighting the vector mass flow rates by a factor 0.1 so as to adjust the evaluations to the real behaviour of a wing A. The flaps of the wing A are inclined by 8° with respect to a main plane of the wing A.

(144) The panel B in FIG. 10 shows on the one hand, drag force coefficients C.sub.z evaluated according to a method implemented by a flight simulator, according to a method of the prior art using nomograph pointers. These values of the drag force coefficients C.sub.z are shown by the continuous curves, for example the curve (a). Moreover, the panel B in FIG. 10 shows drag force coefficients C.sub.z evaluated according to the invention, and corresponding to the curves marked with stars, for example the curve (b). The evaluations of the C.sub.Z according to the invention are obtained by weighting the scalar mass flow rates by a factor 0.25 and weighting the vector mass flow rates by a factor 0.1 so as to adjust the evaluations to the real behaviour of a wing A. The flaps of the wing A are inclined by 8° with respect to a main plane of the wing A.

(145) The panel C in FIG. 10 shows on the one hand, drag force coefficients C.sub.X evaluated according to a method implemented by a flight simulator, according to a method of the prior art using nomograph pointers. These values of the drag force coefficients C.sub.X are shown by the continuous curves, for example the curve (a). Moreover, the panel C in FIG. 10 shows drag force coefficients C.sub.X evaluated according to the invention, and corresponding to the curves marked with stars, for example the curve (b). The C.sub.X evaluations according to the invention are obtained by weighting the scalar mass flow rates by a factor 0.75 and weighting the vector mass flow rates by a factor 0.1 so as to adjust the evaluations to the real behaviour of a wing A. The flaps of the wing A are inclined by 42° with respect to a main plane of the wing A.

(146) The panel D in FIG. 10 shows on the one hand, drag force coefficients C, evaluated according to a method implemented by a flight simulator, according to a method of the prior art using nomograph pointers. These values of the drag force coefficients C.sub.z are shown by the continuous curves, for example the curve (a). Moreover, the panel D in FIG. 10 shows drag force coefficients C.sub.z evaluated according to the invention, and corresponding to the curves marked with stars, for example the curve (b). The evaluations of the C.sub.Z according to the invention are obtained by weighting the scalar mass flow rates by a factor 0.75 and weighting the vector mass flow rates by a factor 0.1 so as to adjust the evaluations to the real behaviour of a wing A. The flaps of the wing A are inclined by 42° with respect to a main plane of the wing A.

(147) FIG. 1 shows a wing A and an arrangement of sources S according to the invention. In variants of a method for the construction of a model according to the invention, advantageously fewer than one hundred sources S are arranged, preferentially fewer than ten sources S. In this way, it is possible to reduce the calculation time, by a simulator SIM implementing a simulation method according to the invention, to a period less than 6 ms.

(148) FIG. 12 is a photograph of a simulator SIM according to the invention.