Non-linear resonance in strings under narrow band random excitation, part I: Planar response and stability

Non-linear planar response of a string to planar narrow band random excitation is investigated in this paper. A response equation for the mean square deflection σ2 is obtained under a single mode approximation by using the equivalent linearization technique. It is shown that the response is triple valued, as in the case of harmonic excitation, if the centre frequency of excitation Ω lies in a certain specified range. The triple valued response occurs only if the excitation bandwidth β is smaller than a critical value βcrit which is a monotonically increasing function of the intensity of excitation. An approximate method of investigating the almost sure asymptotic stability of the solution is presented and regions of instability in the Ω-σ2 plane have been charted. It is shown that planar response can become unstable either due to an unbounded growth of the in-plane component of motion or due to a spontaneous appearance of an out-of-plane component.

The accuracy of varietal selection using factor analytic models for multi-environment plant breeding trials

Modeling of cultivar x trial effects for multienvironment trials (METs) within a mixed model framework is now common practice in many plant breeding programs. The factor analytic (FA) model is a parsimonious form used to approximate the fully unstructured form of the genetic variance-covariance matrix in the model for MET data. In this study, we demonstrate that the FA model is generally the model of best fit across a range of data sets taken from early generation trials in a breeding program. In addition, we demonstrate the superiority of the FA model in achieving the most common aim of METs, namely the selection of superior genotypes. Selection is achieved using best linear unbiased predictions (BLUPs) of cultivar effects at each environment, considered either individually or as a weighted average across environments. In practice, empirical BLUPs (E-BLUPs) of cultivar effects must be used instead of BLUPs since variance parameters in the model must be estimated rather than assumed known. While the optimal properties of minimum mean squared error of prediction (MSEP) and maximum correlation between true and predicted effects possessed by BLUPs do not hold for E-BLUPs, a simulation study shows that E-BLUPs perform well in terms of MSEP.

Ghost fishing in the pot fishery for blue swimmer crabs Portunus pelagicus in Queensland, Australia

Blue swimmer crabs (Portunus pelagicus) are an economically important crab caught in baited traps throughout the Indo-west Pacific and Mediterranean. In Australia they are traditionally caught using rigid wire traps (approximate to pots) but there has been a recent increase in the use of collapsible pots constructed from polyethylene trawl mesh. Two experiments were conducted in Moreton Bay, Queensland, to determine the ghost fishing potential of lost crab pots on both target and bycatch species and to evaluate the differences between traditional and contemporary pot designs. A lost contemporary, collapsible trawl mesh pot will catch between 3 and 223 R pelagicus per year after the bait has been exhausted, while a traditional wire mesh pot would catch 11-74 crabs peryear. As most fishers now use the collapsible trawl mesh pots, ghost fishing mortality could be as high as 111,811-670,866 crabs per year. Bycatch retention was also higher in contemporary designs. Periods of strong winds appeared to increase the ghost fishing potential of lost pots. The use of escape gaps, larger mesh sizes and construction options that allow for the deterioration of entrance funnels to minimise ghost fishing are recommended to reduce environmental impacts.

Noble Metal Ionic Catalysts

Because of growing environmental concerns and increasingly stringent regulations governing auto emissions, new more efficient exhaust catalysts are needed to reduce the amount of pollutants released from internal combustion engines. To accomplish this goal, the major pollutants in exhaust-CO, NOx, and unburned hydrocarbons-need to be fully converted to CO2, N-2, and H2O. Most exhaust catalysts contain nanocrystalline noble metals (Pt, Pd, Rh) dispersed on oxide supports such as Al2O3 or SiO2 promoted by CeO2. However, in conventional catalysts, only the surface atoms of the noble metal particles serve as adsorption sites, and even in 4-6 nm metal particles, only 1/4 to 1/5 of the total noble metal atoms are utilized for catalytic conversion. The complete dispersion of noble metals can be achieved only as ions within an oxide support. In this Account, we describe a novel solution to this dispersion problem: a new solution combustion method for synthesizing dispersed noble metal ionic catalysts. We have synthesized nanocrystalline, single-phase Ce1-xMxO2-delta and Ce1-x-yTiyMxO2-delta (M = Pt, Pd, Rh; x = 0,01-0.02, delta approximate to x, y = 0.15-0.25) oxides in fluorite structure, In these oxide catalysts, pt(2+), Pd2+, or Rh3+ ions are substituted only to the extent of 1-2% of Ce4+ ion. Lower-valent noble metal ion substitution in CeO2 creates oxygen vacancies. Reducing molecules (CO, H-2, NH3) are adsorbed onto electron-deficient noble metal ions, while oxidizing (02, NO) molecules are absorbed onto electron-rich oxide ion vacancy sites. The rates of CO and hydrocarbon oxidation and NOx reduction (with >80% N-2 selectivity) are 15-30 times higher in the presence of these ionic catalysts than when the same amount of noble metal loaded on an oxide support is used. Catalysts with palladium ion dispersed in CeO2 or Ce1-xTixO2 were far superior to Pt or Rh ionic catalysts. Therefore, we have demonstrated that the more expensive Pt and Rh metals are not necessary in exhaust catalysts. We have also grown these nanocrystalline ionic catalysts on ceramic cordierite and have reproduced the results we observed in powder material on the honeycomb catalytic converter. Oxygen in a CeO2 lattice is activated by the substitution of Ti ion, as well as noble metal ions. Because this substitution creates longer Ti-O and M-O bonds relative to the average Ce-O bond within the lattice, the materials facilitate high oxygen storage and release. The interaction among M-0/Mn+, Ce4+/Ce3+, and Ti4+/Ti3+ redox couples leads to the promoting action of CeO2, activation of lattice oxygen and high oxygen storage capacity, metal support interaction, and high rates of catalytic activity in exhaust catalysis.

Swarm intelligence-based solver for parameter estimation of laboratory through-diffusion transport of contaminants

Theoretical approaches are of fundamental importance to predict the potential impact of waste disposal facilities on ground water contamination. Appropriate design parameters are generally estimated be fitting theoretical models to data gathered from field monitoring or laboratory experiments. Transient through-diffusion tests are generally conducted in the laboratory to estimate the mass transport parameters of the proposed barrier material. Thes parameters are usually estimated either by approximate eye-fitting calibration or by combining the solution of the direct problem with any available gradient-based techniques. In this work, an automated, gradient-free solver is developed to estimate the mass transport parameters of a transient through-diffusion model. The proposed inverse model uses a particle swarm optimization (PSO) algorithm that is based on the social behavior of animals searching for food sources. The finite difference numerical solution of the forward model is integrated with the PSO algorithm to solve the inverse problem of parameter estimation. The working principle of the new solver is demonstrated and mass transport parameters are estimated from laboratory through-diffusion experimental data. An inverse model based on the standard gradient-based technique is formulated to compare with the proposed solver. A detailed comparative study is carried out between conventional methods and the proposed solver. The present automated technique is found to be very efficient and robust. The mass transport parameters are obtained with great precision.

Delineating environmental envelopes to improve mapping of species distributions, via a hurdle model with CART &/or MaxEnT

Species distribution modelling (SDM) typically analyses species’ presence together with some form of absence information. Ideally absences comprise observations or are inferred from comprehensive sampling. When such information is not available, then pseudo-absences are often generated from the background locations within the study region of interest containing the presences, or else absence is implied through the comparison of presences to the whole study region, e.g. as is the case in Maximum Entropy (MaxEnt) or Poisson point process modelling. However, the choice of which absence information to include can be both challenging and highly influential on SDM predictions (e.g. Oksanen and Minchin, 2002). In practice, the use of pseudo- or implied absences often leads to an imbalance where absences far outnumber presences. This leaves analysis highly susceptible to ‘naughty-noughts’: absences that occur beyond the envelope of the species, which can exert strong influence on the model and its predictions (Austin and Meyers, 1996). Also known as ‘excess zeros’, naughty noughts can be estimated via an overall proportion in simple hurdle or mixture models (Martin et al., 2005). However, absences, especially those that occur beyond the species envelope, can often be more diverse than presences. Here we consider an extension to excess zero models. The two-staged approach first exploits the compartmentalisation provided by classification trees (CTs) (as in O’Leary, 2008) to identify multiple sources of naughty noughts and simultaneously delineate several species envelopes. Then SDMs can be fit separately within each envelope, and for this stage, we examine both CTs (as in Falk et al., 2014) and the popular MaxEnt (Elith et al., 2006). We introduce a wider range of model performance measures to improve treatment of naughty noughts in SDM. We retain an overall measure of model performance, the area under the curve (AUC) of the Receiver-Operating Curve (ROC), but focus on its constituent measures of false negative rate (FNR) and false positive rate (FPR), and how these relate to the threshold in the predicted probability of presence that delimits predicted presence from absence. We also propose error rates more relevant to users of predictions: false omission rate (FOR), the chance that a predicted absence corresponds to (and hence wastes) an observed presence, and the false discovery rate (FDR), reflecting those predicted (or potential) presences that correspond to absence. A high FDR may be desirable since it could help target future search efforts, whereas zero or low FOR is desirable since it indicates none of the (often valuable) presences have been ignored in the SDM. For illustration, we chose Bradypus variegatus, a species that has previously been published as an exemplar species for MaxEnt, proposed by Phillips et al. (2006). We used CTs to increasingly refine the species envelope, starting with the whole study region (E0), eliminating more and more potential naughty noughts (E1–E3). When combined with an SDM fit within the species envelope, the best CT SDM had similar AUC and FPR to the best MaxEnt SDM, but otherwise performed better. The FNR and FOR were greatly reduced, suggesting that CTs handle absences better. Interestingly, MaxEnt predictions showed low discriminatory performance, with the most common predicted probability of presence being in the same range (0.00-0.20) for both true absences and presences. In summary, this example shows that SDMs can be improved by introducing an initial hurdle to identify naughty noughts and partition the envelope before applying SDMs. This improvement was barely detectable via AUC and FPR yet visible in FOR, FNR, and the comparison of predicted probability of presence distribution for pres/absence.

Liikkumisen ja kokemusten yhteys virtuaaliympäristöissä

The aim of this study was to build a model and analyze how users move in a virtual environment and to explore the experiential dimensions connected with different ways of moving. Due to the lack of previous research on this subject, this was an explorative study. This study also aimed to identify different ways how users move in virtual environments and the background variables connected to them. It was hypothesized that fluent movement in virtual environments is connected to high presence, skill and challenge assessments. Test participants (n = 68) were mostly highly educated young adults. A virtual environment was built using a CAVE -type virtual reality interface. The task was to search for objects that do not belong into a normal house. The participants movement in the virtual house was recorded on a computer. Movement was modelled using a cluster analysis of information entropy based movement measurements, acceleration, amount of stops and time spent being stationary. The experiential dimensions were measured using the EVEQ -questionnaire. We were able to identify four different ways of moving in virtual environments. In respect of background variables, the four groups differed only in the amount of weekly computer usage. However, fluent movement in virtual environments was connected to a high sense of presence. Furthermore, participants who moved fluently in the environment assessed their skills as being high and regarded the use of virtual environment as challenging. The results indicate that different ways of moving affects how people experience virtual environments. Consequently the participants assessment of their skills and level of challenge have an impact on the affective evaluation of the situation at hand. Entropy measures have not been previously applied when studying movement, and in addition the role of movement on the experiential dimensions of virtual environments is an unexplored subject. The movement analysis method introduced here is applicable to other research problems. Finally, this study expands on our knowledge of the special characteristics connected with the experiential dimensions of virtual environments.

Crystal and molecular structure of the quinoxaline antibiotic analog TANDEM (des-N-tetramethyltriostin A)

The crystal structure of TANDEM (des-N-tetramethyltriostin A), a synthetic analogue of the quinoxaline antibiotic triostin A, has been determined independently at -135 and 7 'C and refined to R values of 0.088 and 0.147, respectively. The molecule has approximate 2-fold symmetry, with the quinoxaline chromophores and the disulfide cross-bridge projecting from opposite sides of the peptide ring. The quinoxaline groups are nearly parallel to each other and separated by about 6.5 A. The peptide backbone resembles a distorted antiparallel 13 ribbon joined by intramolecular hydrogen bonds N-H(LVal)--O(L-Ala). At low temperatures, the TANDEM molecule is surrounded by a regular first- and second-order hydration sphere containing 14 independent water molecules. At room temperature, only the first-order hydration shell is maintained. Calculations of the interplanar separation of the quinoxaline groups as a function of their orientation with respect to the peptide ring support the viability of TANDEM to intercalate bifunctionally into DNA.

Exact solutions of coupled multispecies linear reaction–diffusion equations on a uniformly growing domain

Embryonic development involves diffusion and proliferation of cells, as well as diffusion and reaction of molecules, within growing tissues. Mathematical models of these processes often involve reaction–diffusion equations on growing domains that have been primarily studied using approximate numerical solutions. Recently, we have shown how to obtain an exact solution to a single, uncoupled, linear reaction–diffusion equation on a growing domain, 0 < x < L(t), where L(t) is the domain length. The present work is an extension of our previous study, and we illustrate how to solve a system of coupled reaction–diffusion equations on a growing domain. This system of equations can be used to study the spatial and temporal distributions of different generations of cells within a population that diffuses and proliferates within a growing tissue. The exact solution is obtained by applying an uncoupling transformation, and the uncoupled equations are solved separately before applying the inverse uncoupling transformation to give the coupled solution. We present several example calculations to illustrate different types of behaviour. The first example calculation corresponds to a situation where the initially–confined population diffuses sufficiently slowly that it is unable to reach the moving boundary at x = L(t). In contrast, the second example calculation corresponds to a situation where the initially–confined population is able to overcome the domain growth and reach the moving boundary at x = L(t). In its basic format, the uncoupling transformation at first appears to be restricted to deal only with the case where each generation of cells has a distinct proliferation rate. However, we also demonstrate how the uncoupling transformation can be used when each generation has the same proliferation rate by evaluating the exact solutions as an appropriate limit.

Estimation of red blood cell lifespan from alveolar carbon monoxide measurements

Measurement of alveolar carbon monoxide (CO) presents a facile technique to estimate the lifespan, L, of red blood cells (RBCs) in vivo. Several recent studies employ this technique and calculate L (in days) using the expression, L = 13.8 (Hb)/P-CO(end), where (Hb) is the concentration (in g/dL) of hemoglobin in blood, and P-CO(end) is the endogenous production of CO (in ppm). Implicit in this calculation is the assumption that the fraction, f, of endogenous CO production due to RBC turnover is a constant equal to 0.7, which yields the expected RBC lifespan, L approximate to 120 days, in normal controls. In anemic patients, however, enhanced RBC turnover may increase f substantially above 0.7. The above expression then overestimates L. Here, we deriv an alternative tive expression, L = 3390[Hb]/322P(CO (end)-110, that accounts explicitly for the dependence of f on the rate of RBC turnover and thereby provides more accurate estimates of L without requiring additional measurements. Using the latter expression, we recalculate L from recent measurements on hepatitis C virus infected patients undergoing treatment with ribavirin. We find that our estimates of L in these patients (39 +/- 13 days) are significantly lower than current estimates (46 +/- 14 days), indicating that ribavirin affects RBC survival more severely than expected from current studies. Our expression for L is simple to employ in a clinical setting and would render the broadly applicable technique of alveolar CO measurement for the estimation of RBC lifespan more accurate.

Studies on characteristics of lightning generated currents in an interconnected lightning protection system

An in-depth knowledge about the characteristics of lightning generated currents will facilitate evaluation of the interception efficacy of lightning protection systems. In addition, it would aid in extraction of valuable statistics (from measured current data) on local lightning parameters. Incidentally, present day knowledge on characteristics of lightning induced current in typical lightning protection systems is rather limited. This is particularly true with closely interconnected protection systems, like the one employed in Indian Satellite Launch Pad-II. This system is taken as a specific example in the present study. Various aspects suggest that theoretical modelling would be the best possible approach for the intended work. From the survey of pertinent literature, it is concluded that electromagnetic modelling of lightning return-stroke with current source at the channel base is best suited for this study. Numerical electromagnetic code was used for the required electromagnetic field solution and Fourier transform techniques were employed for computing time-domain results. A validation for the numerical modelling is provided by laborator experiments on a reduced scale model of the system. Apart from ascertaining the influence of various parameters, salient characteristics of tower base currents for different kinds of events are deduced. This knowledge can be used in identifying the type of event, as well as its approximate location. A method for estimation of injected stroke current has also been proposed.

Application of Haar functions for transmission line and transformer differential protection

The development of algorithms, based on Haar functions, for extracting the desired frequency components from transient power-system relaying signals is presented. The applications of these algorithms to impedance detection in transmission line protection and to harmonic restraint in transformer differential protection are discussed. For transmission line protection, three modes of application of the Haar algorithms are described: a full-cycle window algorithm, an approximate full-cycle window algorithm, and a half-cycle window algorithm. For power transformer differential protection, the combined second and fifth harmonic magnitude of the differential current is compared with that of fundamental to arrive at a trip decision. The proposed line protection algorithms are evaluated, under different fault conditions, using realistic relaying signals obtained from transient analysis conducted on a model 400 kV, 3-phase system. The transformer differential protection algorithms are also evaluated using a variety of simulated inrush and internal fault signals.

Fire resistance of light gauge steel frame wall systems lined with gypsum plasterboards

Light gauge steel frame (LSF) wall systems are increasingly used in residential and commercial buildings as load bearing and non-load bearing elements. Conventionally, the fire resistance ratings of such building elements are determined using approximate prescriptive methods based on limited standard fire tests. However, recent studies have shown that in some instances real building fire time-temperature curves could be more severe than the standard fire curve, in terms of maximum temperature and rate of temperature rise. This has caused problems for safe evacuation and rescue activities, and in some instances has also lead to the collapse of buildings earlier than the prescribed fire resistance. Therefore a detailed research study into the performance of LSF wall systems under both standard fire and realistic fire conditions was undertaken using full scale fire tests to understand the fire performance of different LSF wall configurations. Both load bearing and non-load bearing full scale fire tests were performed on LSF walls configurations which included single layer, double layer, externally insulated wall panels made up of different steel sections and thicknesses of gypsum plasterboards. The non-load bearing fire test results were utilized to understand the factors affecting the fire resistance of LSF walls, while loading bearing fire test results led to development of simplified methods to predict the fire resistance ratings of load bearing LSF walls exposed to both standard and realistic design fires. This paper presents the results of full scale experimental study and highlights the effects of standard and realistic fire conditions on fire performance of LSF walls.

Flow over a cylinder with a hinged-splitter plate

Previous work on rigid splitter plates in the wake of a bluff body has shown that the primary vortex shedding can be suppressed for sufficiently long splitter plates. In the present work, we study the problem of a hinged-splitter plate in the wake of a circular cylinder. The splitter plate can rotate about the hinge at the base of the cylinder due to the unsteady fluid forces acting on it, and hence the communication between the two sides of the wake is not totally disrupted as in the rigid splitter plate case. In our study, we investigate this problem in the limit where the stiffness and internal damping associated with the hinge are negligible, and the mass ratio of the splitter plate is small. The experiments show that the splitter plate oscillations increase with Reynolds numbers at low values of Re, and are found to reach a saturation amplitude level at higher Re, Re>4000. This type of saturation amplitude level that appears to continue indefinitely with Re, appears to be related to the fact that there is no structural restoring force, and has been observed previously for transversely oscillating cylinders with no restorin force. In the present case, the saturation tip amplitude level can be tip to 0.45D, where D is the cylinder diameter. For this hinged-rigid splitter plate case, it is found that the splitter plate length to cylinder diameter ratio (L/D) is crucial in determining the character and magnitude of the oscillations. For small splitter plate length (L/D <= 3.0), the oscillations appear to be nearly periodic with tip amplitudes of about 0.45D nearly independent of L/D. The nondiinensional oscillation frequencies (fD/U) on the other hand are found to continuously vary with L/D from fD/U approximate to 0.2 at L/D = 1 to fD/U approximate to 0.1 at L/D = 3. As the splitter plate length is further increased beyond L/D >= 4.0, the character of the splitter plate oscillations suddenly changes. The oscillations become aperiodic with much smaller amplitudes. In this long splitter plate regime, the spectra of the oscillations become broadband, and are reminiscent of the change in character of the wake oscillations seen in the earlier fixed-rigid splitter plate case for L/D >= 5.0. In the present case of the hinged-splitter plate, the sudden transition seen as the splitter plate length (L/D) is increased from 3 to 4 may be attributed to the fact that the wake vortices are no longer able to synchronize with the plate motions for larger splitter plate lengths. Hence, as observed in other vortex-induced vibration problems, the oscillations becomeaperiodic and the amplitude reduces dramatically.

Dynamic Bayesian network inferencing for non-homogeneous complex systems

Dynamic Bayesian Networks (DBNs) provide a versatile platform for predicting and analysing the behaviour of complex systems. As such, they are well suited to the prediction of complex ecosystem population trajectories under anthropogenic disturbances such as the dredging of marine seagrass ecosystems. However, DBNs assume a homogeneous Markov chain whereas a key characteristics of complex ecosystems is the presence of feedback loops, path dependencies and regime changes whereby the behaviour of the system can vary based on past states. This paper develops a method based on the small world structure of complex systems networks to modularise a non-homogeneous DBN and enable the computation of posterior marginal probabilities given evidence in forwards inference. It also provides an approach for an approximate solution for backwards inference as convergence is not guaranteed for a path dependent system. When applied to the seagrass dredging problem, the incorporation of path dependency can implement conditional absorption and allows release from the zero state in line with environmental and ecological observations. As dredging has a marked global impact on seagrass and other marine ecosystems of high environmental and economic value, using such a complex systems model to develop practical ways to meet the needs of conservation and industry through enhancing resistance and/or recovery is of paramount importance.