Three-dimensional tracer dispersion model
AssigneeThe United States of America as represented by the Secretary of the Navy
InventorJason K. Jolliff
Life of this patent
Estimating the term…
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.
Claims9
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).
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.
Dependencies are inferred from the claim text.