Breakage model development and application with CFD for predicting breakage of whey protein precipitate particles

Particle breakage due to fluid flow through various geometries can have a major influence on the performance of particle/fluid processes and on the product quality characteristics of particle/fluid products. In this study, whey protein precipitate dispersions were used as a case study to investigate the effect of flow intensity and exposure time on the breakage of these precipitate particles. Computational fluid dynamic (CFD) simulations were performed to evaluate the turbulent eddy dissipation rate (TED) and associated exposure time along various flow geometries. The focus of this work is on the predictive modelling of particle breakage in particle/fluid systems. A number of breakage models were developed to relate TED and exposure time to particle breakage. The suitability of these breakage models was evaluated for their ability to predict the experimentally determined breakage of the whey protein precipitate particles. A "power-law threshold" breakage model was found to provide a satisfactory capability for predicting the breakage of the whey protein precipitate particles. The whey protein precipitate dispersions were propelled through a number of different geometries such as bends, tees and elbows, and the model accurately predicted the mean particle size attained after flow through these geometries. © 2005 Elsevier Ltd. All rights reserved.

Paradox lost? Exploring the role of alpha oscillations during externally vs. internally directed attention and the implications for idling and inhibition hypotheses

Although slow waves of the electroencephalogram (EEG) have been associated with attentional processes, the functional significance of the alpha component in the EEG (8.1–12 Hz) remains uncertain. Conventionally, synchronisation in the alpha frequency range is taken to be a marker of cognitive inactivity, i.e. ‘cortical idling’. However, it has been suggested that alpha may index the active inhibition of sensory information during internally directed attentional tasks such as mental imagery. More recently, this idea has been amended to encompass the notion of alpha synchronisation as a means of inhibition of non-task relevant cortical areas irrespective of the direction of attention. Here we test the adequacy of the one idling and two inhibition hypotheses about alpha. In two experiments we investigated the relation between alpha and internally vs. externally directed attention using mental imagery vs. sensory-intake paradigms. Results from both experiments showed a clear relationship between alpha and both attentional factors and increased task demands. At various scalp sites alpha amplitudes were greater during internally directed attention and during increased load, results incompatible with alpha reflecting cortical idling and more in keeping with suggestions of active inhibition necessary for internally driven mental operations.

Dissipative vector solitons with fast evolving states of polarization

We report on a new vector model of an erbium doped fiber laser mode locked with carbon nanotubes. This model goes beyond the limitations of the previously used models based on either coupled nonlinear Schrödinger or Ginzburg-Landau equations. It results in a new family of vector solitons with fast evolving states of polarization experimentally observed in our previous papers.

Dynamics of vector rogue waves in a fiber laser with a ring cavity

A pulse–pulse interaction that leads to rogue wave (RW) generation in lasers was previously attributed either to soliton–soliton or soliton–dispersive-wave interaction. The beating between polarization modes in the absence of a saturable absorber causes similar effects. Accounting for these polarization modes in a laser resonator is the purpose of the distributed vector model of laser resonators. Furthermore, high pump power, high amplitude, and short pulse duration are not necessary conditions to observe pulse attraction, repulsion, and collisions and the resonance exchange of energy between among them. The regimes of interest can be tuned just by changing the birefringence in the cavity with the pump power slightly higher than the laser threshold. This allows the observation of a wide range of RW patterns in the same experiment, as well as to classify them. The dynamics of the interaction between pulses leads us to the conclusion that all of these effects occur due to nonlinearity induced by the inverse population in the active fiber as well as an intrinsic nonlinearity in the passive part of the cavity. Most of the mechanisms of pulse–pulse interaction were found to be mutually exclusive. This means that all the observed RW patterns, namely, the “lonely,” “twins,” “three sisters,” and “cross,” are probably different cases of the same process.

(Table 1) Annual mean concentrations of atmospheric a- and g-Hexachlorocyclohexane at stations Alert, Canada and Zeppelin, Norway

Twelve year datasets of weekly atmospheric concentrations of alpha- and gamma-HCH were compared between the two Arctic monitoring stations of Alert, Nunavut, Canada, and Zeppelin Mountain, Svalbard, Norway. Time-series analysis was conducted with the use of dynamic harmonic regression (DHR), which provided a very good model fit, to examine both the seasonal behaviour in these isomers and the longer-term, underlying trends. Strong spatial differences were not apparent between the two sites, although subtle differences in seasonal behaviour and composition were identified. For example, the composition of gamma-HCH to total HCH (alpha + gamma) was greater at Zeppelin compared to Alert, probably reflecting this site's proximity to major use regions of lindane. Pronounced seasonality in air concentrations for gamma-HCH was marked by a 'spring maximum event' (SME), confirming earlier studies. For alpha-HCH, the SME was much weaker and only evident at Alert, whereas at Zeppelin, seasonal fluctuations for alpha-HCH were marked by elevated concentrations in summer and lower concentrations during winter, with this pattern most apparent for the years after 2000. We attribute this difference in spatial and temporal patterns to the Arctic oscillation. A similar climatic pattern was not evident at either site in the gamma-HCH data. Seasonally adjusted, long-term trends revealed declining concentrations at both sites for alpha- and gamma-HCH over the entire time-series. Recent legislation affecting lindane use appear to account for this decline in gamma-HCH, with little evidence of a delay or 'lag' between the banning of lindane in Europe (a main source region) or Canada, and a decline in air concentrations observed at both Arctic sites.

Holocene precipitation change in different monsoon sub-regions (time-slices and transient data) simulated by different global climate models, with links to model results

The recently proposed global monsoon hypothesis interprets monsoon systems as part of one global-scale atmospheric overturning circulation, implying a connection between the regional monsoon systems and an in-phase behaviour of all northern hemispheric monsoons on annual timescales (Trenberth et al., 2000). Whether this concept can be applied to past climates and variability on longer timescales is still under debate, because the monsoon systems exhibit different regional characteristics such as different seasonality (i.e. onset, peak, and withdrawal). To investigate the interconnection of different monsoon systems during the pre-industrial Holocene, five transient global climate model simulations have been analysed with respect to the rainfall trend and variability in different sub-domains of the Afro-Asian monsoon region. Our analysis suggests that on millennial timescales with varying orbital forcing, the monsoons do not behave as a tightly connected global system. According to the models, the Indian and North African monsoons are coupled, showing similar rainfall trend and moderate correlation in rainfall variability in all models. The East Asian monsoon changes independently during the Holocene. The dissimilarities in the seasonality of the monsoon sub-systems lead to a stronger response of the North African and Indian monsoon systems to the Holocene insolation forcing than of the East Asian monsoon and affect the seasonal distribution of Holocene rainfall variations. Within the Indian and North African monsoon domain, precipitation solely changes during the summer months, showing a decreasing Holocene precipitation trend. In the East Asian monsoon region, the precipitation signal is determined by an increasing precipitation trend during spring and a decreasing precipitation change during summer, partly balancing each other. A synthesis of reconstructions and the model results do not reveal an impact of the different seasonality on the timing of the Holocene rainfall optimum in the different sub-monsoon systems. They rather indicate locally inhomogeneous rainfall changes and show, that single palaeo-records should not be used to characterise the rainfall change and monsoon evolution for entire monsoon sub-systems.

Model reduction for nonlinear aerodynamics and aeroelasticity using a Discrete Empirical Interpolation Method

A novel surrogate model is proposed in lieu of Computational Fluid Dynamics (CFD) solvers, for fast nonlinear aerodynamic and aeroelastic modeling. A nonlinear function is identified on selected interpolation points by
a discrete empirical interpolation method (DEIM). The flow field is then reconstructed using a least square approximation of the flow modes extracted
by proper orthogonal decomposition (POD). The aeroelastic reduce order
model (ROM) is completed by introducing a nonlinear mapping function
between displacements and the DEIM points. The proposed model is investigated to predict the aerodynamic forces due to forced motions using
a N ACA 0012 airfoil undergoing a prescribed pitching oscillation. To investigate aeroelastic problems at transonic conditions, a pitch/plunge airfoil
and a cropped delta wing aeroelastic models are built using linear structural models. The presence of shock-waves triggers the appearance of limit
cycle oscillations (LCO), which the model is able to predict. For all cases
tested, the new ROM shows the ability to replicate the nonlinear aerodynamic forces, structural displacements and reconstruct the complete flow
field with sufficient accuracy at a fraction of the cost of full order CFD
model.

Monitoring and model updating of an FRP pedestrian truss bridge

This paper presents the results of a real bridge field experiment, carried out on a fiber reinforced polymer (FRP) pedestrian truss bridge of which nodes are reinforced with stainless steel plates. The aim of this paper is to identify the dynamic parameters of this bridge by using both conventional techniques and a model updating algorithm. In the field experiment, the bridge was instrumented with accelerometers at a number of locations on the bridge deck, recording both vertical and transverse vibrations. It was excited via jump tests at particular locations along its span and the resulting acceleration signals are used to identify dynamic parameters, such as the bridge mode shape, natural frequency and damping constant. Pedestrianinduced vibrations are also measured and utilized to identify dynamic parameters of the bridge. For a complete analysis of the bridge, a numerical model of the FRP bridge is created whose properties are calibrated utilizing a model updating algorithm. Comparable frequencies and mode shapes to those from the experiment were obtained by the FE models considering the reinforcement by increasing elastic modulus at every node of the bridge by stainless steel plate. Moreover, considering boundary conditions at both ends as fixed in the model resulted in modal properties comparable/similar to those from the experiment. This study also demonstrated that the effect of reinforcement and boundary conditions must be properly considered in an FE model to analyze real behavior of the FRP bridge.

The FIREDASS (Fire Detection and Suppression Simulation) Model

The FIREDASS (FIRE Detection And Suppression Simulation) project is concerned with the development of fine water mist systems as a possible replacement for the halon fire suppression system currently used in aircraft cargo holds. The project is funded by the European Commission, under the BRITE EURAM programme. The FIREDASS consortium is made up of a combination of Industrial, Academic, Research and Regulatory partners. As part of this programme of work, a computational model has been developed to help engineers optimise the design of the water mist suppression system. This computational model is based on Computational Fluid Dynamics (CFD) and is composed of the following components: fire model; mist model; two-phase radiation model; suppression model and detector/activation model. The fire model - developed by the University of Greenwich - uses prescribed release rates for heat and gaseous combustion products to represent the fire load. Typical release rates have been determined through experimentation conducted by SINTEF. The mist model - developed by the University of Greenwich - is a Lagrangian particle tracking procedure that is fully coupled to both the gas phase and the radiation field. The radiation model - developed by the National Technical University of Athens - is described using a six-flux radiation model. The suppression model - developed by SINTEF and the University of Greenwich - is based on an extinguishment crietrion that relies on oxygen concentration and temperature. The detector/ activation model - developed by Cerberus - allows the configuration of many different detector and mist configurations to be tested within the computational model. These sub-models have been integrated by the University of Greenwich into the FIREDASS software package. The model has been validated using data from the SINTEF/GEC test campaigns and it has been found that the computational model gives good agreement with these experimental results. The best agreement is obtained at the ceiling which is where the detectors and misting nozzles would be located in a real system. In this paper the model is briefly described and some results from the validation of the fire and mist model are presented.

Investigating disease persistence in host vector systems: dengue as a case study

The investigation of pathogen persistence in vector-borne diseases is important in different ecological and epidemiological contexts. In this thesis, I have developed deterministic and stochastic models to help investigating the pathogen persistence in host-vector systems by using efficient modelling paradigms. A general introduction with aims and objectives of the studies conducted in the thesis are provided in Chapter 1. The mathematical treatment of models used in the thesis is provided in Chapter 2 where the models are found locally asymptotically stable. The models used in the rest of the thesis are based on either the same or similar mathematical structure studied in this chapter. After that, there are three different experiments that are conducted in this thesis to study the pathogen persistence. In Chapter 3, I characterize pathogen persistence in terms of the Critical Community Size (CCS) and find its relationship with the model parameters. In this study, the stochastic versions of two epidemiologically different host-vector models are used for estimating CCS. I note that the model parameters and their algebraic combination, in addition to the seroprevalence level of the host population, can be used to quantify CCS. The study undertaken in Chapter 4 is used to estimate pathogen persistence using both deterministic and stochastic versions of a model with seasonal birth rate of the vectors. Through stochastic simulations we investigate the pattern of epidemics after the introduction of an infectious individual at different times of the year. The results show that the disease dynamics are altered by the seasonal variation. The higher levels of pre-existing seroprevalence reduces the probability of invasion of dengue. In Chapter 5, I considered two alternate ways to represent the dynamics of a host-vector model. Both of the approximate models are investigated for the parameter regions where the approximation fails to hold. Moreover, three metrics are used to compare them with the Full model. In addition to the computational benefits, these approximations are used to investigate to what degree the inclusion of the vector population in the dynamics of the system is important. Finally, in Chapter 6, I present the summary of studies undertaken and possible extensions for the future work.

Bumps and rings in a two-dimensional neural field: splitting and rotational instabilities

In this paper we consider instabilities of localised solutions in planar neural field firing rate models of Wilson-Cowan or Amari type. Importantly we show that angular perturbations can destabilise spatially localised solutions. For a scalar model with Heaviside firing rate function we calculate symmetric one-bump and ring solutions explicitly and use an Evans function approach to predict the point of instability and the shapes of the dominant growing modes. Our predictions are shown to be in excellent agreement with direct numerical simulations. Moreover, beyond the instability our simulations demonstrate the emergence of multi-bump and labyrinthine patterns. With the addition of spike-frequency adaptation, numerical simulations of the resulting vector model show that it is possible for structures without rotational symmetry, and in particular multi-bumps, to undergo an instability to a rotating wave. We use a general argument, valid for smooth firing rate functions, to establish the conditions necessary to generate such a rotational instability. Numerical continuation of the rotating wave is used to quantify the emergent angular velocity as a bifurcation parameter is varied. Wave stability is found via the numerical evaluation of an associated eigenvalue problem.

A SCR Model based on Reactor and Engine Experimental Studies for a Cu-zeolite Catalyst

A NOx reduction efficiency higher than 95% with NH3 slip less than 30 ppm is desirable for heavy-duty diesel (HDD) engines using selective catalytic reduction (SCR) systems to meet the US EPA 2010 NOx standard and the 2014-2018 fuel consumption regulation. The SCR performance needs to be improved through experimental and modeling studies. In this research, a high fidelity global kinetic 1-dimensional 2-site SCR model with mass transfer, heat transfer and global reaction mechanisms was developed for a Cu-zeolite catalyst. The model simulates the SCR performance for the engine exhaust conditions with NH3 maldistribution and aging effects, and the details are presented. SCR experimental data were collected for the model development, calibration and validation from a reactor at Oak Ridge National Laboratory (ORNL) and an engine experimental setup at Michigan Technological University (MTU) with a Cummins 2010 ISB engine. The model was calibrated separately to the reactor and engine data. The experimental setup, test procedures including a surrogate HD-FTP cycle developed for transient studies and the model calibration process are described. Differences in the model parameters were determined between the calibrations developed from the reactor and the engine data. It was determined that the SCR inlet NH3 maldistribution is one of the reasons causing the differences. The model calibrated to the engine data served as a basis for developing a reduced order SCR estimator model. The effect of the SCR inlet NO2/NOx ratio on the SCR performance was studied through simulations using the surrogate HD-FTP cycle. The cumulative outlet NOx and the overall NOx conversion efficiency of the cycle are highest with a NO2/NOx ratio of 0.5. The outlet NH3 is lowest for the NO2/NOx ratio greater than 0.6. A combined engine experimental and simulation study was performed to quantify the NH3 maldistribution at the SCR inlet and its effects on the SCR performance and kinetics. The uniformity index (UI) of the SCR inlet NH3 and NH3/NOx ratio (ANR) was determined to be below 0.8 for the production system. The UI was improved to 0.9 after installation of a swirl mixer into the SCR inlet cone. A multi-channel model was developed to simulate the maldistribution effects. The results showed that reducing the UI of the inlet ANR from 1.0 to 0.7 caused a 5-10% decrease in NOx reduction efficiency and 10-20 ppm increase in the NH3 slip. The simulations of the steady-state engine data with the multi-channel model showed that the NH3 maldistribution is a factor causing the differences in the calibrations developed from the engine and the reactor data. The Reactor experiments were performed at ORNL using a Spaci-IR technique to study the thermal aging effects. The test results showed that the thermal aging (at 800°C for 16 hours) caused a 30% reduction in the NH3 stored on the catalyst under NH3 saturation conditions and different axial concentration profiles under SCR reaction conditions. The kinetics analysis showed that the thermal aging caused a reduction in total NH3 storage capacity (94.6 compared to 138 gmol/m3), different NH3 adsorption/desorption properties and a decrease in activation energy and the pre-exponential factor for NH3 oxidation, standard and fast SCR reactions. Both reduction in the storage capability and the change in kinetics of the major reactions contributed to the change in the axial storage and concentration profiles observed from the experiments.

DYNAMIC MODEL BASED STATE ESTIMATION IN A HEAVY DUTY DIESEL AFTERTREATMENT SYSTEM FOR ONBOARD DIAGNOSTICS AND CONTROLS

Estimating un-measurable states is an important component for onboard diagnostics (OBD) and control strategy development in diesel exhaust aftertreatment systems. This research focuses on the development of an Extended Kalman Filter (EKF) based state estimator for two of the main components in a diesel engine aftertreatment system: the Diesel Oxidation Catalyst (DOC) and the Selective Catalytic Reduction (SCR) catalyst. One of the key areas of interest is the performance of these estimators when the catalyzed particulate filter (CPF) is being actively regenerated. In this study, model reduction techniques were developed and used to develop reduced order models from the 1D models used to simulate the DOC and SCR. As a result of order reduction, the number of states in the estimator is reduced from 12 to 1 per element for the DOC and 12 to 2 per element for the SCR. The reduced order models were simulated on the experimental data and compared to the high fidelity model and the experimental data. The results show that the effect of eliminating the heat transfer and mass transfer coefficients are not significant on the performance of the reduced order models. This is shown by an insignificant change in the kinetic parameters between the reduced order and 1D model for simulating the experimental data. An EKF based estimator to estimate the internal states of the DOC and SCR was developed. The DOC and SCR estimators were simulated on the experimental data to show that the estimator provides improved estimation of states compared to a reduced order model. The results showed that using the temperature measurement at the DOC outlet improved the estimates of the CO , NO , NO2 and HC concentrations from the DOC. The SCR estimator was used to evaluate the effect of NH3 and NOX sensors on state estimation quality. Three sensor combinations of NOX sensor only, NH3 sensor only and both NOX and NH3 sensors were evaluated. The NOX only configuration had the worst performance, the NH3 sensor only configuration was in the middle and both the NOX and NH3 sensor combination provided the best performance.

Vector self-pulsing in erbium-doped fiber lasers

Insight into instabilities of fiber laser regimes leading to complex self-pulsing operations is an opportunity to unlock the high power and dynamic operation tunability of lasers. Though many models have been suggested, there is no complete covering of self-pulsing complexity observed experimentally. Here, I further generalized our previous vector model of erbium-doped fiber laser and, for the first time, to the best of my knowledge, map tunability of complex vector self-pulsing on Poincare sphere (limit cycles and double scroll polarization attractors) for laser parameters, e.g., power, ellipticity of the pump wave, and in-cavity birefringence. Analysis validated by extensive numerical simulations demonstrates good correspondence to the experimental results on complex self-pulsing regimes obtained by many authors during the last 20 years.

Toxicity of polychlorinated biphenyls and N-phenyl-1-naphthylamine in the juvenile turtle, Chelydra serpentina

Polychlorinated biphenyls (PCBs) and substituted phenylamine antioxidants (SPAs) are two chemical groups that have been used in multiple Canadian industrial processes. Despite the production ban of PCBs in North America in 1977, they are still ubiquitous in the environment and in wildlife tissues. Previous studies of fish, amphibians, birds, and mammals have shown that PCBs are toxic and act as endocrine disruptors. In contrast, SPAs, specifically N-phenyl-1-naphthylamine (PANA), have received very little attention despite their current use in Canada and their expected environmental releases. The effects of PCB and PANA exposures in reptiles remain unknown thus, juvenile Chelydra serpentina were used in this thesis as a model vertebrate to fill in missing toxicity research gaps due to their importance as an environmental indicator. First, food pellets were spiked at an environmentally relevant concentration of the PCB mixture Aroclor 1254 (A1254) to model hepatic bioaccumulation (0.45 μg/g A1254 for 31 days) and depuration (clean food for 50 days) of PCBs in turtles. No significant differences in PCB concentrations were observed between the control and treated animals, suggesting that juvenile turtles exposed to environmentally relevant concentrations of PCBs can likely detoxify low concentrations of PCBs. Additionally, two dose-response experiments were performed using A1254 or PANA spiked food (0-12.7 μg/g and 0-3,446 μg/g, respectively) to determine hepatic toxicity and bioaccumulation in juvenile C. serpentina. An increase in hepatic cyp1a was observed when exposed to the highest dose of both chemicals: 1) for A1254, induction correlated to the significant increase in hepatic PCB congeners that are known to be metabolized by CYP1A; and 2) for PANA, induction suggested that CYP1A has a potential role in its detoxification. PCBs are known endocrine disruptors, but no significant changes were observed for both thyroid receptors (alpha and beta) or by estrogen and androgen receptors. This lack of response, also noted in the PANA exposure, suggests that C. serpentina is less sensitive to endocrine disruption than other vertebrates. Furthermore, the expression of genes involved in cellular stress was not altered in PCB and PANA exposed animals, supporting the resilience of turtles to oxidative stress. This is the first study to demonstrate the toxicity of PCBs and PANA in C. serpentina, demonstrating the turtle’s high tolerance to contamination.