Search PubMed⌕ Search
Home/Archive/PubMed

PubMed · 12656293

Fast transport simulation with an adaptive grid refinement.

Frieder Haefner·Siegrun Boy
View original source ↗Export metadata ↓

Abstract

One of the main difficulties in transport modeling and calibration is the extraordinarily long computing times necessary for simulation runs. Improved execution time is a prerequisite for calibration in transport modeling. In this paper we investigate the problem of code acceleration using an adaptive grid refinement, neglecting subdomains, and devising a method by which the Courant condition can be ignored while maintaining accurate solutions. Grid refinement is based on dividing selected cells into regular subcells and including the balance equations of subcells in the equation system. The connection of coarse and refined cells satisfies the mass balance with an interpolation scheme that is implicitly included in the equation system. The refined subdomain can move with the average transport velocity of the subdomain. Very small time steps are required on a fine or a refined grid, because of the combined effect of the Courant and Peclet conditions. Therefore, we have developed a special upwind technique in small grid cells with high velocities (velocity suppression). We have neglected grid subdomains with very small concentration gradients (zero suppression). The resulting software, MODCALIF, is a three-dimensional, modularly constructed FORTRAN code. For convenience, the package names used by the well-known MODFLOW and MT3D computer programs are adopted, and the same input file structure and format is used, but the program presented here is separate and independent. Also, MODCALIF includes algorithms for variable density modeling and model calibration. The method is tested by comparison with an analytical solution, and illustrated by means of a two-dimensional theoretical example and three-dimensional simulations of the variable-density Cape Cod and SALTPOOL experiments. Crossing from fine to coarse grid produces numerical dispersion when the whole subdomain of interest is refined; however, we show that accurate solutions can be obtained using a fraction of the execution time required by uniformly fine-grid solutions.

Explore related subjects

CalibrationModels, TheoreticalParticle SizeReproducibility of ResultsWater MovementsWater SupplyJournal Article

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Frieder Haefner, Siegrun Boy. Fast transport simulation with an adaptive grid refinement.. https://doi.org/10.1111/j.1745-6584.2003.tb02590.x

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

PubMed2007

Fermentation process tracking through enhanced spectral calibration modeling.

The FDA process analytical technology (PAT) initiative will materialize in a significant increase in the number of installations of spectroscopic instrumentation. However, to attain the greatest benefit from the data generated, there is a need for calibration procedures that extract the maximum information content. For example, in fermentation processes, the interpretation of the resulting spectra is challenging as a consequence of the large number of wavelengths recorded, the underlying correlation structure that is evident between the wavelengths and the impact of the measurement environment. Approaches to the development of calibration models have been based on the application of partial least squares (PLS) either to the full spectral signature or to a subset of wavelengths. This paper presents a new approach to calibration modeling that combines a wavelength selection procedure, spectral window selection (SWS), where windows of wavelengths are automatically selected which are subsequently used as the basis of the calibration model. However, due to the non-uniqueness of the windows selected when the algorithm is executed repeatedly, multiple models are constructed and these are then combined using stacking thereby increasing the robustness of the final calibration model. The methodology is applied to data generated during the monitoring of broth concentrations in an industrial fermentation process from on-line near-infrared (NIR) and mid-infrared (MIR) spectrometers. It is shown that the proposed calibration modeling procedure outperforms traditional calibration procedures, as well as enabling the identification of the critical regions of the spectra with regard to the fermentation process.

Calibration↗
PubMed2007

Particle sizing with a fast polar nephelometer.

We reported previously the design of a polar nephelometer that uses a rotational confocal imaging setup to enable fast scanning of the scattering phase function within a field of view of 55 degrees . The full dynamic range of the detection system can be used by increasing the signal-to-noise ratio by means of averaging successive scans. The calibration of the angular response of the instrument is achieved by obtaining the transfer function of the optical detection system using Rayleigh scatterers. Accurate particle sizing of individual polystyrene spheres (ranging from 1.5 to 9 micro m in diameter) in aqueous suspension is achieved by maximizing a correlation coefficient between precalculated tables of Mie phase functions and data obtained from the polar nephelometer. Good correlation is achieved between experimental and theoretical data, proving the functioning of the instrument as a fast and convenient particle sizer.

Calibration↗
PubMed2007

Online standard additions calibration of transient signals for inductively coupled plasma mass spectrometry.

An online standard additions calibration method for transient signals in ICPMS is demonstrated in which a small volume of standard is injected as a spike into the sample/carrier stream, overlaying the analyte peak. This technique provides the advantages of conventional standard additions but requires only a single sample run. The method corrects for matrix effects and is suitable for transient signals in which the severity of the matrix effect changes over the analyte peak. The method uses a peak-fitting program to determine the area of the underlying peak and is shown to be effective for the determination of trace metal concentrations in both a high ionic strength matrix and in a biological matrix (urine). Eight analytes with concentrations in the range of 0.82-233.2 mug L-1 in urine were simultaneously determined using a standard spiking solution of 75 mug L-1 injected through a 100-muL loop. The measured concentrations for analytes free of spectral interferences agreed with the certified values, and the precision achieved was comparable to that achieved by the certifying agency. Using a conventional cross-flow nebulizer and Scott-type spray chamber, the accuracy obtained for online standard additions calibration was within 2%, and the precision was within 5%.

Calibration↗