Three-dimensional tracer dispersion model

Patent metadata

Field	Value
Publication number	US-8949096-B2
Application number	US-201213467215-A
Country	US
Kind code	B2
Filing date	May 9, 2012
Priority date	May 9, 2012
Publication date	Feb 3, 2015
Grant date	Feb 3, 2015

Abstract

System for solving for the fully three-dimensional advection diffusion reaction (ADR) of dissolved or particulate tracers (biological or chemical materials) in aquatic environments including an input processor and an ADR tracer field processor. Results from a single execution of an ocean circulation model may be used to drive a separate ADR computer simulation and compute a tracer forecast. The velocity fields are not required to be from an ocean circulation model, they could be analysis fields derived from some other source, such as high frequency RADAR observations or satellite-based surface ocean velocity inversion/detection methods.

First claim

What is claimed is:

1. A computer method for forecasting characteristics that can be approximated by mass concentration of a tracer in an incompressible fluid comprising:
(a) defining, by use of a specially-programmed computer, a 2D spatial grid for an area of interest in the incompressible fluid;
(b) defining, by use of the specially-programmed computer, a 3D spatial grid based on the 2D spatial grid and bathymetric data;
(c) defining, by use of the specially-programmed computer, surfaces for each grid point in the 3D spatial grid, the surfaces in combination creating a plurality of grid cells;
(d) defining, by use of the specially-programmed computer, an initial concentration of mass for each of the grid cells over the 3D spatial grid;
(e) defining, by use of the specially-programmed computer, the velocity at each of the grid points in the 3D spatial grid by interpolating a pre-selected spatial 3D flow field to the 3D spatial grid;
(f) calculating, by use of the specially-programmed computer, and using a first order upwind differencing scheme, fluid flow through each of the surfaces based on the velocities at each of the grid points;
(g) computing, by use of the specially-programmed computer, a total flux through each of the surfaces based on the fluid flow;
(h) computing, by use of the specially-programmed computer, a residual mass flux of layers of a water column based on the total flux and the residual mass flux of adjacent of the layers, the residual mass flux including residual velocities above a continuity threshold;
(i) computing, by use of the specially-programmed computer, adjusted flow fields by removing the residual velocities from the velocities;
(j) performing, by use of the specially-programmed computer, time-dependent advection/diffusion calculation to compute the concentration forecast based on the adjusted flow fields;
(k) determining remaining of the residual velocities in the adjusted flow fields;
(l) if the remaining residual velocities are above the continuity threshold, updating the spatial 3D flow field; and
(m) repeating steps (e) through (l).

2. The method as in claim 1 further comprising:
storing the spacing of each of the grid cells as either a scalar for constant spacing or an array for irregular spacing.

3. The method as in claim 1 wherein the area of interest comprises 50-3000 meters of the water column.

4. The method as in claim 1 wherein the bathymetric data comprise a bathymetric data resolution of a 2D grid resolution of the 2D spatial grid.

5. The method as in claim 1 further comprising:
designating a vertical coordinate as a fraction of the depth of the water column.

6. The method as in claim 1 wherein the defining surfaces for each grid point in the 3D spatial grid comprises:
interpolating a simulated finite volume around each grid point in the 3D spatial grid over the bathymetry data, if the bathymetry data is not flat.

7. The method as in claim 1 wherein the initial concentration of mass for each of the grid cells comprises point concentrations over the 3D spatial grid.

8. A computer system for forecasting characteristics that can be approximated by mass concentration of a tracer comprising:
a computer programmed to execute a grid/surface processor defining a 2D spatial grid for an area of interest in an incompressible fluid, a 3D spatial grid based on the 2D spatial grid and bathymetric data, surfaces for each grid point in the 3D spatial grid, the surfaces in combination creating a plurality of grid cells, and an initial concentration of mass for each of the grid cells over the 3D spatial grid;
a computer programmed to execute a velocity grid point processor defining the velocity at each of the grid points in the 3D spatial grid by interpolating a pre-selected spatial 3D flow field to the 3D spatial grid;
a computer programmed to execute a vertical flux processor computing a total flux through each of the surfaces based on the fluid flow;
a computer programmed to execute a residual processor computing a residual mass flux of the layers of a water column based on the total flux and the residual mass flux of adjacent of the layers, the residual mass flux including residual velocities;
a computer programmed to execute a flow field processor calculating, using a first order fluid flow through each of the surfaces based on the velocities at each of the grid points, and adjusting the flow fields by removing the residual velocities from the velocities, the flow field processor determining remaining of the residual velocities in the adjusted flow fields and, if the remaining residual velocities are above a continuity threshold, updating the spatial 3D flow field, the flow field processor recomputing the fluid flow, the total flux, the residual mass, the residual mass flux, the adjusted flow fields, and the concentration forecast;
a computer programmed to execute a velocity receiver receiving velocity components if the adjusted flow fields have an absolute residual velocity above the continuity threshold; and
a computer programmed to execute a tracer field processor performing time-dependent advection/diffusion calculation to compute the concentration forecast based on the adjusted flow fields.

9. The system as in claim 8 wherein the velocity grid point processor accesses the flow fields from a hydrodynamic circulation model, the bathymetric data, the initial concentrations, and the 2D spatial grid for the area of interest.

Assignees

• Jolliff Jason K
• Us Navy

Inventors

• Jolliff Jason K

Classifications

• G06F2111/10

  Numerical modelling · CPC title

• G06F30/20 Primary

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

• G06F17/5018 Primary

  Physics · mapped topic

• G06F30/25

  using particle-based methods · CPC title

• G06F2113/08

  Fluids · CPC title