Theory manual to OctVCE – a cartesian cell CFD code with special application to blast wave problems, Department of Mechanical Engineering Report 2007/12

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries, in particular, from explosions. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel for its simplicity and sufficiency for practical engineering design problems. The code uses a finite-volume formulation of the unsteady Euler equations with a second order explicit Runge-Kutta Godonov (MUSCL) scheme. Gradients are calculated using a least-squares method with a minmod limiter. Flux solvers used are AUSM, AUSMDV and EFM. No fluid-structure coupling or chemical reactions are allowed, but gas models can be perfect gas and JWL or JWLB for the explosive products. This report also describes the code’s ‘octree’ mesh adaptive capability and point-inclusion query procedures for the VCE geometry engine. Finally, some space will also be devoted to describing code parallelization using the shared-memory OpenMP paradigm. The user manual to the code is to be found in the companion report 2007/13.

Finite element analysis of steady-state natural convection problems in fluid-saturated porous media heated from below

In this paper, a progressive asymptotic approach procedure is presented for solving the steady-state Horton-Rogers-Lapwood problem in a fluid-saturated porous medium. The Horton-Rogers-Lapwood problem possesses a bifurcation and, therefore, makes the direct use of conventional finite element methods difficult. Even if the Rayleigh number is high enough to drive the occurrence of natural convection in a fluid-saturated porous medium, the conventional methods will often produce a trivial non-convective solution. This difficulty can be overcome using the progressive asymptotic approach procedure associated with the finite element method. The method considers a series of modified Horton-Rogers-Lapwood problems in which gravity is assumed to tilt a small angle away from vertical. The main idea behind the progressive asymptotic approach procedure is that through solving a sequence of such modified problems with decreasing tilt, an accurate non-zero velocity solution to the Horton-Rogers-Lapwood problem can be obtained. This solution provides a very good initial prediction for the solution to the original Horton-Rogers-Lapwood problem so that the non-zero velocity solution can be successfully obtained when the tilted angle is set to zero. Comparison of numerical solutions with analytical ones to a benchmark problem of any rectangular geometry has demonstrated the usefulness of the present progressive asymptotic approach procedure. Finally, the procedure has been used to investigate the effect of basin shapes on natural convection of pore-fluid in a porous medium. (C) 1997 by John Wiley & Sons, Ltd.

Analysis and application of an equilibrium model for in vitro bioassay systems with three components: Receptor, hormone and hormone-binding-protein

A simple theoretical framework is presented for bioassay studies using three component in vitro systems. An equilibrium model is used to derive equations useful for predicting changes in biological response after addition of hormone-binding-protein or as a consequence of increased hormone affinity. Sets of possible solutions for receptor occupancy and binding protein occupancy are found for typical values of receptor and binding protein affinity constants. Unique equilibrium solutions are dictated by the initial condition of total hormone concentration. According to the occupancy theory of drug action, increasing the affinity of a hormone for its receptor will result in a proportional increase in biological potency. However, the three component model predicts that the magnitude of increase in biological potency will be a small fraction of the proportional increase in affinity. With typical initial conditions a two-fold increase in hormone affinity for its receptor is predicted to result in only a 33% increase in biological response. Under the same conditions an Ii-fold increase in hormone affinity for receptor would be needed to produce a two-fold increase in biological potency. Some currently used bioassay systems may be unrecognized three component systems and gross errors in biopotency estimates will result if the effect of binding protein is not calculated. An algorithm derived from the three component model is used to predict changes in biological response after addition of binding protein to in vitro systems. The algorithm is tested by application to a published data set from an experimental study in an in vitro system (Lim et al., 1990, Endocrinology 127, 1287-1291). Predicted changes show good agreement (within 8%) with experimental observations. (C) 1998 Academic Press Limited.

Sensitive detection of nitric oxide using seeded parametric four-wave mixing

A sensitive near-resonant four-wave mixing technique based on two-photon parametric four-wave mixing has been developed. Seeded parametric four-wave mixing requires only a single laser as an additional phase matched seeder field is generated via parametric four-wave mixing of the pump beam in a high gain cell. The seeder field travels collinearly with the pump beam providing efficient nondegenerate four-wave mixing in a second medium. This simple arrangement facilitates the detection of complex molecular spectra by simply scanning the pump laser. Seeded parametric four-wave mixing is demonstrated in both a low pressure cell and an air/acetylene flame with detection of the two-photon C (2) Pi(upsilon'=0)<--X (2) Pi(upsilon =0) spectrum of nitric oxide. From the cell data a detection limit of 10(12) molecules/cm(3) is established. A theoretical model of seeded parametric four-wave mixing is developed from existing parametric four-wave mixing theory. The addition of the seeder field significantly modifies the parametric four-wave mixing behaviour such that in the small signal regime, the signal intensity can readily be made to scale as the cube of the laser pump power while the density dependence follows a more familiar square law dependence, In general, we find excellent agreement between theory and experiment. Limitations to the process result from an ac Stark shift of the two-photon resonance in the high pressure seeder cell caused by the generation of a strong seeder field, as well as a reduction in phase matching efficiency due to the presence of certain buffer species. Various optimizations are suggested which should overcome these limitations, providing even greater detection sensitivity. (C) 1998 American Institute of Physics, [S0021-9606(98)01014-9].

Sensitive detection of sodium in a flame using parametric four-wave mixing and seeded parametric four-wave mixing

Two-photon resonant parametric four-wave mixing and a newly developed variant called seeded parametric four-wave mixing are used to detect trace quantities of sodium in a flame. Both techniques are simple, requiring only a single laser to generate a signal beam at a different wavelength which propagates collinearly with the pump beam, allowing efficient signal recovery. A comparison of the two techniques reveals that seeded parametric four-wave mixing is more than two orders of magnitude more sensitive than parametric four-wave mixing, with an estimated detection sensitivity of 5 x 10(9) atoms/cm(3). Seeded parametric four-wave mixing is achieved by cascading two parametric four-wave mixing media such that one of the parametric fields generated in the first high-density medium is then used to seed the same four-wave mixing process in a second medium in order to increase the four-wave mixing gain. The behavior of this seeded parametric four-wave mixing is described using semiclassical perturbation theory. A simplified small-signal theory is found to model most of the data satisfactorily. However, an anomalous saturationlike behavior is observed in the large signal regime. The full perturbation treatment, which includes the competition between two different four-wave mixing processes coupled via the signal field, accounts for this apparently anomalous behavior.

Optimising fed-batch production of recombinant proteins using the baculovirus expression vector system

Fed-batch culture can offer significant improvement in recombinant protein production compared to batch culture in the baculovirus expression vector system (BEVS), as shown by Nguyen et al. (1993) and Bedard et al. (1994) among others. However, a thorough analysis of fed-batch culture to determine its limits in improving recombinant protein production over batch culture has yet to be performed. In this work, this issue is addressed by the optimisation of single-addition fed-batch culture. This type of fed-batch culture involves the manual addition of a multi-component nutrient feed to batch culture before infection with the baculovirus. The nutrient feed consists of yeastolate ultrafiltrate, lipids, amino acids, vitamins, trace elements, and glucose, which were added to batch cultures of Spodoptera frugiperda (Sf9) cells before infection with a recombinant Autographa californica nuclear polyhedrosis virus (Ac-NPV) expressing beta-galactosidase (beta-Gal). The fed-batch production of beta-Gal was optimised using response surface methods (RSM). The optimisation was performed in two stages, starting with a screening procedure to determine the most important variables and ending with a central-composite experiment to obtain a response surface model of volumetric beta-Gal production. The predicted optimum volumetric yield of beta-Gal in fed-batch culture was 2.4-fold that of the best yields in batch culture. This result was confirmed by a statistical analysis of the best fed-batch and batch data (with average beta-Gal yields of 1.2 and 0.5 g/L, respectively) obtained from this laboratory. The response surface model generated can be used to design a more economical fed-batch operation, in which nutrient feed volumes are minimised while maintaining acceptable improvements in beta-Gal yield. (C) 1998 John Wiley & Sons, Inc.

Subcycling first- and second-order generalizations of the trapezoidal rule

Subcycling algorithms which employ multiple timesteps have been previously proposed for explicit direct integration of first- and second-order systems of equations arising in finite element analysis, as well as for integration using explicit/implicit partitions of a model. The author has recently extended this work to implicit/implicit multi-timestep partitions of both first- and second-order systems. In this paper, improved algorithms for multi-timestep implicit integration are introduced, that overcome some weaknesses of those proposed previously. In particular, in the second-order case, improved stability is obtained. Some of the energy conservation properties of the Newmark family of algorithms are shown to be preserved in the new multi-timestep extensions of the Newmark method. In the first-order case, the generalized trapezoidal rule is extended to multiple timesteps, in a simple way that permits an implicit/implicit partition. Explicit special cases of the present algorithms exist. These are compared to algorithms proposed previously. (C) 1998 John Wiley & Sons, Ltd.

Solution of adsorption problems involving steep moving profiles

The moving finite element collocation method proposed by Kill et al. (1995) Chem. Engng Sci. 51 (4), 2793-2799 for solution of problems with steep gradients is further developed to solve transient problems arising in the field of adsorption. The technique is applied to a model of adsorption in solids with bidisperse pore structures. Numerical solutions were found to match the analytical solution when it exists (i.e. when the adsorption isotherm is linear). The method is simple yet sufficiently accurate for use in adsorption problems, where global collocation methods fail. (C) 1998 Elsevier Science Ltd. All rights reserved.

Three-dimensional quantum solitons with parametric coupling

We consider the quantum field theory of two bosonic fields interacting via both parametric (cubic) and quartic couplings. In the case of photonic fields in a nonlinear optical medium, this corresponds to the process of second-harmonic generation (via chi((2)) nonlinearity) modified by the chi((3)) nonlinearity. The quantum solitons or energy eigenstates (bound-state solutions) are obtained exactly in the simplest case of two-particle binding, in one, two, and three space dimensions. We also investigate three-particle binding in one space dimension. The results indicate that the exact quantum solitons of this field theory have a singular, pointlike structure in two and three dimensions-even though the corresponding classical theory is nonsingular. To estimate the physically accessible radii and binding energies of the bound states, we impose a momentum cutoff on the nonlinear couplings. In the case of nonlinear optical interactions, the resulting radii and binding energies of these photonic particlelike excitations in highly nonlinear parametric media appear to be close to physically observable values.

Multidimensional parametric quantum solitons

We consider the parametric quantum field theory involving cubic and quartic couplings of two bosonic fields. This is exactly soluble for the two-particle energy eigenstates (or quantum solitons) in one, two, and three space dimensions. We estimate the binding energies and corresponding radii in the case of photonic fields in nonlinear optical materials, and Bose-Einstein condensates. [S1050-2947(98)51110-9].

Theory of multidimensional parametric band-gap simultons

Multidimensional spatiotemporal parametric simultons (simultaneous solitary waves) are possible in a nonlinear chi((2)) medium with a Bragg grating structure, where large effective dispersion occurs near two resonant band gaps for the carrier and second-harmonic field, respectively. The enhanced dispersion allows much reduced interaction lengths, as compared to bulk medium parametric simultons. The nonlinear parametric band-gap medium permits higher-dimensional stationary waves to form. In addition, solitons can occur with lower input powers than conventional nonlinear Schrodinger equation gap solitons. In this paper, the equations for electromagnetic propagation in a grating structure with a parametric nonlinearity are derived from Maxwell's equation using a coupled mode Hamiltonian analysis in one, two, and three spatial dimensions. Simultaneous solitary wave solutions are proved to exist by reducing the equations to the coupled equations describing a nonlinear parametric waveguide, using the effective-mass approximation (EMA). Exact one-dimensional numerical solutions in agreement with the EMA solutions are also given. Direct numerical simulations show that the solutions have similar types of stability properties to the bulk case, providing the carrier waves are tuned to the two Bragg resonances, and the pulses have a width in frequency space less than the band gap. In summary, these equations describe a physically accessible localized nonlinear wave that is stable in up to 3 + 1 dimensions. Possible applications include photonic logic and switching devices. [S1063-651X(98)06109-1].

Characterizing Greenberger-Horne-Zeilinger correlations in nondegenerate parametric oscillation via phase measurements

We present a potential realization of the Greenberger-Horne-Zeilinger all or nothing contradiction of quantum mechanics with local realism using phase measurement techniques in a simple photon number triplet. Such a triplet could be generated using nondegenerate parametric oscillation. [S0031-9007(98)07671-6].

A state-contingent production approach to principal-agent problems with an application to point-source pollution control

Two basic representations of principal-agent relationships, the 'state-space' and 'parameterized distribution' formulations, have emerged. Although the state-space formulation appears more natural, analytical studies using this formulation have had limited success. This paper develops a state-space formulation of the moral-hazard problem using a general representation of production under uncertainty. A closed-form solution for the agency-cost problem is derived. Comparative-static results are deduced. Next we solve the principal's problem of selecting the optimal output given the agency-cost function. The analysis is applied to the problem of point-source pollution control. (C) 1998 Published by Elsevier Science S.A. All rights reserved.