Computer simulation of physical processes including modeling of laminar-to-turbulent transition

1. A method for simulating a fluid flow on a computer, the method comprising:

for one or more locations on or near a boundary surface:

performing a first calculation where a local boundary layer is taken to be a laminar boundary layer;

performing a second calculation where the local boundary layer is taken to be a turbulent boundary layer;

comparing a result from the first calculation where the local boundary layer is taken to be the laminar boundary layer to a result from the second calculation where the local boundary layer is taken to be the turbulent boundary layer; and

based on the comparing, selecting the result of the first calculation where the local boundary layer is taken to be the laminar boundary layer or the result of the second calculation where the local boundary layer is taken to be the turbulent boundary layer; and

inputting the selected result for the one or more locations into a simulation of activity of a fluid in a volume comprising the boundary surface.

2. The method of claim 1 , wherein:

performing the first calculation where the local boundary layer is taken to be the laminar boundary layer comprises calculating a wall momentum flux tensor property for a laminar flow;

performing the second calculation where the local boundary layer is taken to be the turbulent boundary layer comprises calculating a wall momentum flux tensor property for a turbulent flow; and

selecting comprises selecting the laminar wall momentum flux tensor property or the turbulent wall momentum flux tensor property.

3. The method of claim 1 , further comprising:

determining a laminar to turbulent boundary layer transition by:

determining, for each of multiple facets on the boundary surface, a first measure based on the first calculation and a second measure based on the second calculation; and

classifying the fluid flow for at least some of the multiple facets as laminar or turbulent by comparing at least one of the first and second measures to a criterion.

4. The method of claim 3 , wherein selecting comprises:

for facets classified as laminar, selecting a wall momentum flux tensor property for laminar flow; and

for facets classified as turbulent, selecting a wall momentum flux tensor property for turbulent flow.

5. The method of claim 3 , wherein:

the result of the first calculation comprises a measure of laminar wall momentum flux tensor; and

the result of the second calculation comprises a measure of turbulent wall momentum flux tensor; and

the criterion comprises a measure of turbulence intensity.

6. The method of claim 1 , wherein:

performing the first calculation comprises calculating, for each of multiple facets on the boundary surface, a measure of laminar wall momentum flux tensor and performing the second calculation comprises calculating, for each of multiple facets on the boundary surface, a measure of turbulent wall momentum flux tensor using the second calculation;

comparing comprises comparing, for each of the multiple facets on the boundary surface, a calculated measure of turbulence intensity and the measure of turbulent wall momentum flux tensor; and

selecting comprises selecting, for at least some of the multiple facets on the boundary surface, one of the calculated turbulent wall momentum flux tensor and laminar wall momentum flux tensor based on the comparison of the measure of turbulence intensity and the measure of turbulent wall momentum flux tensor.

7. The method of claim 6 , wherein:

comparing, for each of the multiple facets on the boundary surface, the measure of turbulence intensity and the measure of turbulent wall momentum flux tensor comprises determining if the measure of turbulence intensity is greater than the measure of wall momentum flux tensor; and

selecting, for at least some of the multiple facets on the boundary surface, one of the calculated turbulent wall momentum flux tensor and laminar wall momentum flux tensor comprises, for a particular facet, selecting the turbulent wall momentum flux tensor if the measure of turbulence intensity is greater than the measure of turbulent wall momentum flux tensor and selecting the measure of laminar wall momentum flux tensor if the measure of turbulence intensity is less than the measure of turbulent wall momentum flux tensor.

8. The method of claim 6 , further comprising calculating a value of local turbulent kinetic energy.

9. The method of claim 6 , wherein, for a given near-wall fluid velocity the measure of turbulent wall momentum flux tensor is greater than the measure of laminar wall momentum flux tensor.

10. The method of claim 1 , further comprising simulating activity of the fluid in the volume by:

performing interaction operations on state vectors, the interaction operations modeling interactions between elements of different momentum states according to a model; and

performing first move operations of the state vectors to reflect movement of elements to new voxels in the volume according to the model.

11. The method of claim 1 , wherein the second calculation comprises a calculation to determine a measure of turbulent wall momentum flux tensor based on a velocity profile and a wall distance.

12. The method of claim 1 , further comprising:

selecting, for at least some of multiple facets on the boundary surface, a value that is based on a weighted average of the result of the first calculation and the result of the second calculation.

13. The method of claim 3 , further comprising:

selecting, for at least some of the multiple facets on the boundary surface, a wall momentum flux tensor property that is based on a combination of a turbulent wall momentum flux tensor property and a laminar wall momentum flux tensor property.

14. The method of claim 1 , wherein the second calculation comprises a calculation to determine a measure of turbulent wall momentum flux tensor based on local turbulent kinetic energy and a local fluid velocity.

15. The method of claim 1 , wherein a voxel size in a region adjacent to the boundary surface is similar to a voxel size at regions spaced apart from the boundary surface.

16. The method of claim 1 , wherein a voxel size in a region adjacent to the boundary surface is the same as a voxel size at regions spaced apart from the boundary surface.

17. A computer program product tangibly embodied in a computer readable hardware storage device, the computer program product including instructions that, when executed, simulate a fluid flow, the computer program product configured to cause a computer to:

for one or more locations on or near a boundary surface:

perform a first calculation where a local boundary layer is taken to be a laminar boundary layer;

perform a second calculation where the local boundary layer is taken to be a turbulent boundary layer;

compare a result from the first calculation where the local boundary layer is taken to be the laminar boundary layer to a result from the second calculation where the local boundary layer is taken to be the turbulent boundary layer; and

based on the compare, select the result of the first calculation where the local boundary layer is taken to be the laminar boundary layer or the result of the second calculation where the local boundary layer is taken to be the turbulent boundary layer; and

input the selected result for the one or more locations into a simulation of activity of a fluid in a volume comprising the boundary surface.

18. The computer program product of claim 17 , wherein:

the instructions to perform the first calculation comprise instructions to calculate a wall momentum flux tensor property for a laminar flow;

the instructions to perform the second calculation comprise instructions to calculate a wall momentum flux tensor property for a turbulent flow; and

the instructions to select comprise instructions to select the laminar wall momentum flux tensor property or the turbulent wall momentum flux tensor property.

19. The computer program product of claim 17 , further comprising instructions to determine a laminar to turbulent boundary layer transition by comprising instructions to:

determine, for each of multiple facets on the boundary surface, a first measure based on the first calculation and a second measure based on the second calculation; and

classify the fluid flow for at least some of the multiple facets as laminar or turbulent by comparing at least one of the first and second measures to a criterion.

20. The computer program product of claim 19 , wherein the instructions to select comprise:

for facets classified as laminar, instructions to select a wall momentum flux tensor property for laminar flow; and

for facets classified as turbulent, instructions to select a wall momentum flux tensor property for turbulent flow.

21. The computer program product of claim 19 , wherein:

the result of the first calculation comprises a measure of laminar wall momentum flux tensor; and

the result of the second calculation comprises a measure of turbulent wall momentum flux tensor; and

the criterion comprises a measure of turbulence intensity.

22. The computer program product of claim 17 , wherein:

the instructions to perform the first calculation comprise instructions to calculate, for each of multiple facets on the boundary surface, a measure of laminar wall momentum flux tensor and perform the second calculation comprises instructions to calculate, for each of multiple facets on the boundary surface, a measure of turbulent wall momentum flux tensor using the second boundary layer calculation; and

the instructions to compare comprise instructions to compare, for each of the multiple facets on the boundary surface, a calculated measure of turbulence intensity and the measure of turbulent wall momentum flux tensor; and

the instructions to select comprise instructions to select, for at least some of the multiple facets on the boundary surface, one of the calculated turbulent wall momentum flux tensor and laminar wall momentum flux tensor based on the comparison of the measure of turbulence intensity and the measure of turbulent wall momentum flux tensor.

23. A computer system for simulating a fluid flow, the system being configured to:

for one or more locations on or near a boundary surface:

perform a first calculation where a local boundary layer is taken to be a laminar boundary layer;

perform a second calculation where the local boundary layer is taken to be a turbulent boundary layer;

compare a result from the first calculation where the local boundary layer is taken to be the laminar boundary layer to a result from the second calculation where the local boundary layer is taken to be the turbulent boundary layer; and

based on the compare, select the result of the first calculation where the local boundary layer is taken to be the laminar boundary layer or the result of the second calculation where the local boundary layer is taken to be the turbulent boundary layer; and

input the selected result for the one or more locations into a simulation of activity of a fluid in a volume comprising the boundary surface.

24. The system of claim 23 , wherein:

the configurations to perform the first calculation for the laminar boundary layer flow comprise configurations to calculate a wall momentum flux tensor property for the laminar flow;

the configurations to perform the second calculation for the turbulent boundary layer flow comprise configurations to calculate a wall momentum flux tensor property for the turbulent flow; and

the configurations to select the results of the first boundary layer calculation or the results of the second calculation for a turbulent boundary layer flow comprise configurations to select the laminar wall momentum flux tensor property or the turbulent wall momentum flux tensor property.

25. The system of claim 23 , further comprising configurations to determine a laminar to turbulent boundary layer transition for the boundary layer that comprise configurations to:

determine, for each of multiple facets on the boundary surface, a first measure based on the first calculation and a second measure based on the second calculation; and

classify the fluid flow for at least some of the multiple facets as laminar or turbulent by comparing at least one of the first and second measures to a criterion.

26. The system of claim 25 , wherein the configurations to select comprise configurations to:

for facets classified as laminar, configurations to select a wall momentum flux tensor property value for laminar flow; and

for facets classified as turbulent, configurations to select a wall momentum flux tensor property for turbulent flow.

27. The system of claim 25 , wherein:

the result of the first calculation comprises a measure of laminar wall momentum flux tensor;

the result of the second calculation comprises a measure of wall momentum flux tensor; and

the criteria comprises a measure of turbulence intensity.

28. The system of claim 23 , wherein:

the configurations to perform the first calculation comprise configurations to calculate, for each of multiple facets on the boundary surface, a measure of laminar wall momentum flux tensor and the configurations to perform the second calculation comprise configurations to calculate, for each of multiple facets on the boundary surface, a measure of turbulent wall momentum flux tensor using the second calculation; and

the configurations to compare comprise configurations to compare, for each of the multiple facets on the boundary surface, a calculated measure of turbulence intensity and the measure of turbulent wall momentum flux tensor; and

the configurations to select comprise configurations to:

select, for at least some of the multiple facets on the boundary surface, one of the calculated turbulent wall momentum flux tensor and laminar wall momentum flux tensor based on the comparison of the measure of turbulence intensity and the measure of wall momentum flux tensor.