Reactive processing of polyolefins using antioxidant systems

Various monoacrylic compounds containing a hindered phenol function (e.g.3,5-di-tert.-butyl-4-hydroxy benzyl alcohol, DBBA and vinyl-3-[3',5'-di-tert.-butyl-4-hydroxy phenyl] propionate, VDBP), and a benzophenone function (2-hydroxy-4-[beta hydroxy ethoxy] benzophenone, HAEB) were synthesised and used as reactive antioxidants (AO's) for polypropylene (PP). These compounds were reacted with PP melt in the presence of low concentration of a free radical generator such a peroxide (reactive processing) to produce bound-antioxidant concentrates. The binding reaction of these AO's onto PP was found to be low and this was shown to be mainly due to competing reactions such as homopolymerisation of the antioxidant. At high concentrations of peroxide, higher binding efficiency resulted, but, this was accompanied by melt degradation of the polymer. In a special reactive processing procedure, a di- or a trifunctional reactant (referred to as coagent), e.g.tri-methylol propane tri-acrylate, Tris, and Divinyl benzene, DVB, were used with the antioxidant and this has led to an enhanced efficiency of the grating reaction of antioxidants on the polymer in the melt. The evidence suggests that this is due to copolymerisation of the antioxidants with the coagent as well as grafting of the copolymers onto the polymer backbone. Although the 'bound' AO's containing a UV stabilising function showed lower overall stabilisation effect than the unbound analogues before extraction, they were still much more effective when subjected to exhaustive solvent extraction. Furthermore, a very effective synergistic stabilising activity when two reactive AO's containing thermal and UV stabilising functions e.g. DBBA and HAEB, were reactively processed with PP in the presence of a coagent. The stabilising effectiveness of such a synergist was much higher than that of the unbound analogues both before and after extraction. Analysis using the GPC technique of concentrates containing bound-DBBA processed in the presence of Tris coagent showed higher molecular weight (Mn), compared to that of a polymer processed without the coagent, but was still lower than that of the control processed PP with no additives. This indicates that Tris coagent may inhibit further melt degradation of the polymer. Model reactions of DBBA in liquid hydrocarbon (decalin) and analysis of the products using FTIR and NMR spectroscopy showed the formation of grafted DBBA onto decalin molecules as well as homopolymerisation of the AO. In the presence of Tris coagent, copolymerisation of DBBA with the Tris inevitably occured; which was followed by grafting of the copolymer onto the decalin, FTIR and NMR results of the polymer concentrates containing bound-DBBA processed with and without Tris, showed similar behaviour as the above model reactions. This evidence supports the effect of Tris in enhancing the efficiency of the reaction of DBBA in the polymer melt. Reactive procesing of HAEB in polymer melts exhibited crosslinking formation In the early stages of the reaction, however, in the final stage, the crosslinked structure was 'broken down' or rearranged to give an almost gel free polymer with high antioxidant binding efficiency.

The adsorptive separation of aromatic hydrocarbon mixtures

A detailed study has been made of the feasibility of adsorptive purification of slack waxes from traces of aromatic compounds using type 13X molecular sieves to achieve 0.01% aromatics in the product. The limited literature relating to the adsorption of high molecular weight aromatic compounds by zeolites was reviewed. Equilibrium isotherms were determined for typical individual aromatic compounds. Lower molecular weight, or more compact, molecules were preferentially adsorbed and the number of molecules captured by one unit cell decreased with increasing molecular weight of the adsorbate. An increase in adsorption temperature resulted in a decrease in the adsorption value. The isosteric heat of adsorption of differnt types of aromatic compounds was determined from pairs of isotherms at 303 K to 343 K at specific coverages. The lowest heats of adsorption were for dodecylbenzene and phenanthrene. Kinetics of adsorption were studied for different aromatic compounds. The diffusivity decreased significantly when a long alkyl chain was attached to the benzene ring e.g. in dodecylbenzene; molecules with small cross-sectional diameter e.g. cumene were adsorbed most rapidly. The sorption rate increased with temperature. Apparent activation energies increased with increasing polarity. In a study of the dynamic adsorption of selected aromatic compounds from binary solutions in isooctane or n-alkanes, naphthalene exhibited the best dynamic properties followed by dibenzothiophene and finally dodecylbenzene. The dynamic adsorption of naphthalene from different n-alkane solvents increased with a decrease in solvent molecular weight. A tentative mathematical approach is proposed for the prediction of dynamic breakthrough curves from equilibrium isotherms and kinetic data. The dynamic properties of liquid phase adsorption of aromatics from slack waxes were studied at different temperatures and concentrations. The optimum operating temperature was 543 K. The best dynamic performance was achieved with feeds of low aromatic content. The studies with individual aromatic compounds demonstrated the affinity of type NaX molecular sieves to adsorb aromatics in the concentration range 3% - 5% . Wax purification by adsorption was considered promising and extension of the experimental programme was recommended.

Nuclear magnetic resonance studies of molecular interactions and structures

This thesis is concerned with the investigation, by nuclear magnetic resonance spectroscopy, of the molecular interactions occurring in mixtures of benzene and cyclohexane to which either chloroform or deutero-chloroform has been added. The effect of the added polar molecule on the liquid structure has been studied using spin-lattice relaxation time, 1H chemical shift, and nuclear Overhauser effect measurements. The main purpose of the work has been to validate a model for molecular interaction involving local ordering of benzene around chloroform. A chemical method for removing dissolved oxygen from samples has been developed to encompass a number of types of sample, including quantitative mixtures, and its supremacy over conventional deoxygenation technique is shown. A set of spectrometer conditions, the use of which produces the minimal variation in peak height in the steady state, is presented. To separate the general diluting effects of deutero-chloroform from its effects due to the production of local order a series of mixtures involving carbon tetrachloride, instead of deutero-chloroform, have been used as non-interacting references. The effect of molecular interaction is shown to be explainable using a solvation model, whilst an approach involving 1:1 complex formation is shown not to account for the observations. It is calculated that each solvation shell, based on deutero-chloroform, contains about twelve molecules of benzene or cyclohexane. The equations produced to account for the T1 variations have been adapted to account for the 1H chemical shift variations in the same system. The shift measurements are shown to substantiate the solvent cage model with a cage capacity of twelve molecules around each chloroform molecule. Nuclear Overhauser effect data have been analysed quantitatively in a manner consistent with the solvation model. The results show that discrete shells only exist when the mole fraction of deutero-chloroform is below about 0.08.

Studies of the effects of molecular encounters on N.M.R. spin-lattice relaxation times

This thesis is concerned with investigations of the effects of molecular encounters on nuclear magnetic resonance spin-lattice relaxation times, with particular reference to mesitylene in mixtures with cyclohexane and TMS. The purpose of the work was to establish the best theoretical description of T1 and assess whether a recently identified mechanism (buffeting), that influences n.m.r. chemical shifts, governs Tl also. A set of experimental conditions are presented that allow reliable measurements of Tl and the N. O. E. for 1H and 13C using both C. W. and F.T. n.m.r. spectroscopy. Literature data for benzene, cyclohexane and chlorobenzene diluted by CC14 and CS2 are used to show that the Hill theory affords the best estimation of their correlation times but appears to be mass dependent. Evaluation of the T1 of the mesitylene protons indicates that a combined Hill-Bloembergen-Purcell-Pound model gives an accurate estimation of T1; subsequently this was shown to be due to cancellation of errors in the calculated intra and intemolecular components. Three experimental methods for the separation of the intra and intermolecular relaxation times are described. The relaxation times of the 13C proton satellite of neat bezene, 1,4 dioxane and mesitylene were measured. Theoretical analyses of the data allow the calculation of Tl intra. Studies of intermolecular NOE's were found to afford a general method of separating observed T1's into their intra and intermolecular components. The aryl 1H and corresponding 13C T1 values and the NOE for the ring carbon of mesitylene in CC14 and C6H12-TMS have been used in combination to determine T1intra and T1inter. The Hill and B.P.P. models are shown to predict similarly inaccurate values for T1linter. A buffeting contribution to T1inter is proposed which when applied to the BPP model and to the Gutowsky-Woessner expression for T1inter gives an inaccuracy of 12% and 6% respectively with respect to theexperimentally based T1inter.

Studies of molecular interactions using physical techniques

The investigations described in this thesis concern the molecular interactions between polar solute molecules and various aromatic compounds in solution. Three different physical methods were employed. Nuclear magnetic resonance (n.m.r.) spectroscopy was used to determine the nature and strength of the interactions and the geometry of the transient complexes formed. Cryoscopic studies were used to provide information on the stoichiometry of the complexes. Dielectric constant studies were conducted in an attempt to confirm and supplement the spectroscopic investigations. The systems studied were those between nitromethane, chloroform, acetonitrile (solutes) and various methyl substituted benzenes. In the n.m.r. work the dependence of the solute chemical shift upon the compositions of the solutions was determined. From this the equilibrium quotients (K) for the formation of each complex and the shift induced in the solute proton by the aromatic in the complex were evaluated. The thermodynamic parameters for the interactions were obtained from the determination of K at several temperatures. The stoichiometries of the complexes obtained from cryoscopic studies were found to agree with those deduced from spectroscopic investigations. For most systems it is suggested that only one type of complex, of 1:1 stiochiometry, predominates except that for the acetonitrile-benzene system a 1:2 complex is formed. Two sets of dielectric studies were conducted, the first to show that the nature of the interaction is dipole-induced dipole and the second to calculate K. The equilibrium quotients obtained from spectroscopic and dielectric studies are compared. Time-averaged geometries of the complexes are proposed. The orientation of solute, with respect to the aromatic for the 1:1 complexes, appears to be the one in which the solute lies symmetrically about the aromatic six-fold axis whereas for the 1:2 complex, a sandwich structure is proposed. It is suggested that the complexes are formed through a dipole-induced dipole interaction and steric factors play some part in the complex formation.

Antitumour imidazotetrazines: probes for the major groove of DNA

The antitumour imidazotetrazinones are believed to act as prodrugs for the triazene series of alkylating agents, showing a marked pteference for the alkylation of the middle guanine residue in a run of three or more contiguous guanines. However, the. exact nature of the interactions of imidazotetrazinones within the micro~environment of DNA are; as yet unknown. In order to examine such interactions a three pronged approach involving molecular modelling, synthetic chemistry and biological analysis has been undertaken during the course of this project. . Molecular modelling studies have shown that for the 8-carboxamido substituted imidazotetrazinones antitumour activity is dependent upon the. presence of a free NH group which can be involved in the formation of both intramolecular and intermolecular hydrogen bonds, and the presence of a non-bulky substituent with a small negative potential . volume. Modelling studies involving the docking of .mitozolomide into the major groove of DNA in the region of a triguanine sequence has shown that a number of hydrogen bonding interactions are feasible. A series of 8-substituted carboxamide derivatives of mitozolomide have been synthesised via the 8-acid chloride and 8-carboxylic acid derivatives including a number of peptide analogues. The peptide derivatives were based upon the key structural features of the helix-turn-helix motif of DNA-binding proteins with a view to developing agents that are capable of binding to DNA with greater selectivity. An examination of the importance of intramolecular hydrogen bonding in influencing the antitumour activity:of :the imidazotetrazinones has led to the synthesis of the novel pyrimido[4',5' :4,3]pyrazolo[5,1-d]-1,2,3,5-tetrazine ring system. In general, in vitro cytotoxicity assays showed that the new derivatives were less active against the TLX5 lymphoma cell line. than the parent compound mitozolomide despite an increased potential for hydrogen bonding interactions. Due to the high reactivity of the: tetrazinone ring system it is difficult to study the interactions between the imidazotetrazinones and DNA. Consequently a number of structural analogues that are stable under physiological conditions have been. prepared based upon the 1,2,3 triazin-4(3H)-one ring system fused with both benzene and pyrazole rings. Although the 3-methylbenzotriazinones failed to antagonise the cytotoxic activity of temozolomide encouraging results with a 3-methylpyrazolotriazinone may suggest the existence of an imidazotetrazinone receptor site within DNA. The potential of guanine rich sequences to promote the alkylating selectivity of imidazotetrazinones by acting as a catalyst for ring cleavage and thereby generation of the alkylating agent was examined. Experiments involving the monitoring: of the rate of breakdown of mitozolomide incubated in the presence of synthetic oIigonucleotides did not reveal any catalytic effect resulting from the DNA. However, it was noted that the breakdown of mitozolomide was dependent upon the type of buffer used in the incubations and this may indeed mask any catalysis by the oligonucleotides.

Cationic ring opening polymerisation (CROP) of oxetane and its derivatives using oxonium derived initiators

Several cationic initiator systems were developed and used to polymerise oxetane with two oxonium ion initiator systems being investigated in depth. The first initiator system was generated by the elimination of a chloride group from a chloro methyl ethyl ether. Adding a carbonyl co-catalyst to a carbocationic centre generated the second initiator system. It was found that the anion used to stabilise the initiator was critical to the initial rate of polymerisation of oxetane with hexafluoroantimonate resulting in the fastest polymerisations. Both initiator systems could be used at varying monomer to initiator concentrations to control the molecular number average, Mn, of the resultant polymer. Both initiator systems showed living characteristics and were used to polymerise further monomers and generate higher molecular weight material and block copolymers. Oxetane and 3,3-dimethyl oxetane can both be polymerised using either oxonium ion initiator system in a variety of DCM or DCM/1,4-dioxane solvent mixtures. The level of 1,4-dioxane does have an impact on the initial rate of polymerisation with higher levels resulting in lower initial rates of polymerisation but do tend to result in higher polydispersities. The level of oligomer formation is also reduced as the level of 1,4-dioxane is increased. 3,3-bis-bromomethyl oxetane was also polymerised but a large amount of hyperbranching was seen at the bromide site resulting in a difficult to solvate polymer system. Multifunctional initiator systems were also generated using the halide elimination reactions with some success being achieved with 1,3,5-tris-bromomethyl-2,4,6-tris-methyl-benzene derived initiator system. This offered some control over the molecular number average of the resultant polymer system.

Reactive processing of polymers: structural characterization of products by 1H and 13C NMR spectroscopy for glycidyl methacrylate grafting onto EPR in the absence and presence of a reactive comonomer

Grafted GMA on EPR samples were prepared in a Thermo-Haake internal mixer by free radical melt grafting reactions in the absence (conventional system; EPR-g-GMA(CONV)) and presence of the reactive comonomer divinyl benzene, DVB (EPR-g-GMA(DVB)). The GMA-homopolymer (poly-GMA), a major side reaction product in the conventional system, was almost completely absent in the DVB-containing system, the latter also resulted in a much higher level of GMA grafting. A comprehensive microstructure analysis of the formed poly-GMA was performed based on one-dimensional H-1 and C-13 NMR spectroscopy and the complete spectral assignments were supported by two-dimensional NMR techniques based on long range two and three bond order carbon-proton couplings from HMBC (Heteronuclear Multiple Bond Coherence) and that of one bond carbon-proton couplings from HSQC (Heteronuclear Single Quantum Coherence), as well as the use of Distortionless Enhancement by Polarization Transfer (DEPT) NMR spectroscopy. The unambiguous analysis of the stereochemical configuration of poly-GMA was further used to help understand the microstructures of the GMA-grafts obtained in the two different free radical melt grafting reactions, the conventional and comonomer-containing systems. In the grafted GMA, in the conventional system (EPR-g-GMA(CONV)), the methylene protons of the GMA were found to be sensitive to tetrad configurational sequences and the results showed that 56% of the GMA sequence in the graft is in atactic configuration and 42% is in syndiotactic configuration whereas the poly-GMA was predominantly syndiotactic. The differences in the microstructures of the graft in the conventional EPR-g-GMA(CONV) and the DVB-containing (EPR-g-GMA(DVB)) systems is also reported (C) 2009 Elsevier Ltd. All rights reserved.

Reactive processing of polymers:effect of bifunctional and tri-functional comonomers on melt grafting of glycidyl methacrylate onto polypropylene

Two reactive comonomers, divinyl benzene (DVB) and trimethylolpropane triacrylate (TRIS), were evaluated for their role in effecting the melt free radical grafting reaction of the monomer glycidyl methacrylate (GMA) onto polypropylene (PP). The characteristics of the GMA-grafting systems in the presence and absence of DVB or TRIS were examined and compared in terms of the yield of the grafting reaction and the extent of the main side reactions, namely homopolymerisation of GMA (poly-GMA) and polymer degradation, using different chemical compositions of the reactive systems and processing conditions. In the absence of the comonomers, i.e. in a conventional system, high initiator concentrations of peroxides were typically required to achieve the highest possible GMA grafting levels which were found to be generally low. Concomitantly, both poly-GMA and degradation of the polymer by chain scission takes place with increasing initiator amounts. On the other hand, the presence of a small amount of the comonomers, DVB or Tris, in the GMA-grafting system, was shown to bring about a significant increase in the grafting level paralleled by a large reduction in poly-GMA and PP degradation. In the presence of these highly reactive comonomers, the optimum grafting system requires a much lower concentration of the peroxide initiator and, consequently, would lead to the much lower degree of polymer degradation observed in these systems. The differences in the effects of the presence of DVB and that of TRIS in the grafting systems on the rate of the GMA-grafting and homopolymerisation reactions, and the extent of PP degradation (through melt flow changes), were compared and contrasted with a conventional GMA-grafting system.

Phosphorus catalysis in the pyrolysis behaviour of biomass

Phosphorus is a key plant nutrient and as such, is incorporated into growing biomass in small amounts. This paper examines the influence of phosphorus, present in either acid (HPO) or salt ((NH)PO) form, on the pyrolysis behaviour of both Miscanthus × giganteus, and its cell wall components, cellulose, hemicellulose (xylan) and lignin (Organosolv). Pyrolysis-gas chromatography-mass spectrometry (PY-GC-MS) is used to examine the pyrolysis products during thermal degradation, and thermogravimetric analysis (TGA) is used to examine the distribution of char and volatiles. Phosphorus salts are seen to catalyse the pyrolysis and modify the yields of products, resulting in a large increase in char yield for all samples, but particularly for cellulose and Miscanthus. The thermal degradation processes of cellulose, xylan and Miscanthus samples occur in one step and the main pyrolysis step is shifted to lower temperature in the presence of phosphorus. A small impact of phosphorus was observed in the case of lignin char yields and the types of pyrolysis decomposition products produced. Levoglucosan is a major component produced in fast pyrolysis of cellulose. Furfural and levoglucosenone become more dominant products upon P-impregnation pointing to new rearrangement and dehydration routes. The P-catalysed xylan decomposition route leads to a much simpler mixture of products, which are dominated by furfural, 3-methyl-2-cyclopenten-1-one and one other unconfirmed product, possibly 3,4-dihydro-2-methoxy-2H-pyran or 4-hydroxy-5,6-dihydro-(2H)-pyran-2-one. Phosphorus-catalysed lignin decomposition also leads to a modified mixture of tar components and desaspidinol as well as other higher molecular weight component become more dominant relative to the methoxyphenyl phenols, dimethoxy phenols and triethoxy benzene. Comparison of the results for Miscanthus lead to the conclusion that the understanding of the fast pyrolysis of biomass can, for the most part, be gained through the study of the individual cell wall components, provided consideration is given to the presence of catalytic components such as phosphorus.

Phenol methylation over nanoparticulate Co2FeO4 inverse spinel catalysts:the effect of morphology on catalytic performance

A series of CoFe2O4 nanoparticles have been prepared via co-precipitation and controlled thermal sintering, with tunable diameters spanning 7–50 nm. XRD confirms that the inverse spinel structure is adopted by all samples, while XPS shows their surface compositions depend on calcination temperature and associated particle size. Small (<20 nm) particles expose Fe3+ enriched surfaces, whereas larger (∼50 nm) particles formed at higher temperatures possess Co:Fe surface compositions close to the expected 1:2 bulk ratio. A model is proposed in which smaller crystallites expose predominately (1 1 1) facets, preferentially terminated in tetrahedral Fe3+ surface sites, while sintering favours (1 1 0) and (1 0 0) facets and Co:Fe surface compositions closer to the bulk inverse spinel phase. All materials were active towards the gas-phase methylation of phenol to o-cresol at temperatures as low as 300 °C. Under these conditions, materials calcined at 450 and 750 °C exhibit o-cresol selectivities of ∼90% and 80%, respectively. Increasing either particle size or reaction temperature promotes methanol decomposition and the evolution of gaseous reductants (principally CO and H2), which may play a role in CoFe2O4 reduction and the concomitant respective dehydroxylation of phenol to benzene. The degree of methanol decomposition, and consequent H2 or CO evolution, appears to correlate with surface Co2+ content: larger CoFe2O4 nanoparticles have more Co rich surfaces and are more active towards methanol decomposition than their smaller counterparts. Reduction of the inverse spinel surface thus switches catalysis from the regio- and chemo-selective methylation of phenol to o-cresol, towards methanol decomposition and phenol dehydroxylation to benzene. At 300 °C sub-20 nm CoFe2O4 nanoparticles are less active for methanol decomposition and become less susceptible to reduction than their 50 nm counterparts, favouring a high selectivity towards methylation.

Hydrodebromination of bromobenzene over Pt(111)

The surface chemistry of benzene and bromobenzene over Pt(111) has been studied by temperature-programmed XPS/MS and NEXAFS. Time-resolved XPS shows that benzene adopts a single chemically distinguishable environment during low-temperature adsorption within the monolayer, with a saturation coverage at θC6H6 = 0.2 ML. Around 20% of a benzene monolayer desorbs molecularly, while the remainder dehydrogenates to surface carbon. Bromobenzene likewise adsorbs molecularly at 90 K, giving rise to two C 1s environments at 284.4 and 285.3 eV corresponding to the C−H and C−Br functions, respectively. The saturation C6H5Br monolayer coverage is 0.11 ML. NEXAFS reveals that bromobenzene adopts a tilted geometry, with the ring plane at 60 ± 5° to the surface. Bromobenzene multilayers desorb at ∼180 K, with higher temperatures promoting competitive molecular desorption versus C−Br scission within the monolayer. Approximately 30% of a saturated bromobenzene monolayer either desorbs reversibly or as reactively formed hydrocarbons. Debromination yields a stable (phenyl) surface intermediate and atomic bromine at 300 K. Further heating results in desorption of reactively formed H2, C6H6, and HBr; however, there was no evidence for either biphenyl or Br2 formation. Pt(111) is an efficient surface for low-temperature bromobenzene hydrodebromination to benzene and HBr. © 2007 American Chemical Society.

A fast XPS investigation of NO-promoted acetylene cyclotrimerisation on Pd{111}

The strong promotion by NO of acetylene cyclotrimerisation to benzene on Pd{1 1 1} has been investigated by fast XPS, supported by LEED and Δφ data. Islands of NO act to compress co-existing islands of acetylene. The resulting increase in acetylene density dramatically lowers the nominal threshold acetylene coverage required for reaction. In addition, the threshold reaction temperature is reduced from 150 to <100 K, and increased lateral interactions within the organic domains promote formation of tilted benzene. Co-adsorption of NO with benzene itself accelerates the desorption of both tilted and flat-lying benzene at low temperatures. These factors act to increase the yield of reactively formed benzene to 100%, completely suppressing carbon deposition that would otherwise result from decomposition of reactant and product.

Measurement and correlation of the solubility of telmisartan (form A) in nine different solvents from 277.85 to 338.35 K

The solubility of telmisartan (form A) in nine organic solvents (chloroform, dichloromethane, ethanol, toluene, benzene, 2-propanol, ethyl acetate, methanol and acetone) was determined by a laser monitoring technique at temperatures from 277.85 to 338.35 K. The solubility of telmisartan (form A) in all of the nine solvents increased with temperature as did the rates at which the solubility increased except in chloroform and dichloromethane. The mole fraction solubility in chloroform is higher than that in dichloromethane, which are both one order of magnitude higher than those in the other seven solvents at the experimental temperatures. The solubility data were correlated with the modified Apelblat equation and λh equations. The results show that the λh equation is in better agreement with the experimental data than the Apelblat equation. The relative root mean square deviations (σ) of the λh equation are in the range from 0.004 to 0.45 %. The dissolution enthalpies, entropies and Gibbs energies of telmisartan in these solvents were estimated by the Van’t Hoff equation and the Gibbs equation. The melting point and the fusion enthalpy of telmisartan were determined by differential scanning calorimetry.

Development of a mathematical and computer model to assess water quality impacts of hazardous compounds from minor gasoline spills in the inter-tidal zone of the Miami River, Florida

The objective of this study was to develop a model to predict transport and fate of gasoline components of environmental concern in the Miami River by mathematically simulating the movement of dissolved benzene, toluene, xylene (BTX), and methyl-tertiary-butyl ether (MTBE) occurring from minor gasoline spills in the inter-tidal zone of the river. Computer codes were based on mathematical algorithms that acknowledge the role of advective and dispersive physical phenomena along the river and prevailing phase transformations of BTX and MTBE. Phase transformations included volatilization and settling. ^ The model used a finite-difference scheme of steady-state conditions, with a set of numerical equations that was solved by two numerical methods: Gauss-Seidel and Jacobi iterations. A numerical validation process was conducted by comparing the results from both methods with analytical and numerical reference solutions. Since similar trends were achieved after the numerical validation process, it was concluded that the computer codes algorithmically were correct. The Gauss-Seidel iteration yielded at a faster convergence rate than the Jacobi iteration. Hence, the mathematical code was selected to further develop the computer program and software. The model was then analyzed for its sensitivity. It was found that the model was very sensitive to wind speed but not to sediment settling velocity. ^ A computer software was developed with the model code embedded. The software was provided with two major user-friendly visualized forms, one to interface with the database files and the other to execute and present the graphical and tabulated results. For all predicted concentrations of BTX and MTBE, the maximum concentrations were over an order of magnitude lower than current drinking water standards. It should be pointed out, however, that smaller concentrations than the latter reported standards and values, although not harmful to humans, may be very harmful to organisms of the trophic levels of the Miami River ecosystem and associated waters. This computer model can be used for the rapid assessment and management of the effects of minor gasoline spills on inter-tidal riverine water quality. ^