An overview of APSIM, a model designed for farming systems simulation

The Agricultural Production Systems Simulator (APSIM) is a modular modelling framework that has been developed by the Agricultural Production Systems Research Unit in Australia. APSIM was developed to simulate biophysical process in farming systems, in particular where there is interest in the economic and ecological outcomes of management practice in the face of climatic risk. The paper outlines APSIM's structure and provides details of the concepts behind the different plant, soil and management modules. These modules include a diverse range of crops, pastures and trees, soil processes including water balance, N and P transformations, soil pH, erosion and a full range of management controls. Reports of APSIM testing in a diverse range of systems and environments are summarised. An example of model performance in a long-term cropping systems trial is provided. APSIM has been used in a broad range of applications, including support for on-farm decision making, farming systems design for production or resource management objectives, assessment of the value of seasonal climate forecasting, analysis of supply chain issues in agribusiness activities, development of waste management guidelines, risk assessment for government policy making and as a guide to research and education activity. An extensive citation list for these model testing and application studies is provided. Crown Copyright (C) 2002 Published by Elsevier Science B.V. All rights reserved.

Model-based clustering in gene expression microarrays: An application to breast cancer data

In microarray studies, the application of clustering techniques is often used to derive meaningful insights into the data. In the past, hierarchical methods have been the primary clustering tool employed to perform this task. The hierarchical algorithms have been mainly applied heuristically to these cluster analysis problems. Further, a major limitation of these methods is their inability to determine the number of clusters. Thus there is a need for a model-based approach to these. clustering problems. To this end, McLachlan et al. [7] developed a mixture model-based algorithm (EMMIX-GENE) for the clustering of tissue samples. To further investigate the EMMIX-GENE procedure as a model-based -approach, we present a case study involving the application of EMMIX-GENE to the breast cancer data as studied recently in van 't Veer et al. [10]. Our analysis considers the problem of clustering the tissue samples on the basis of the genes which is a non-standard problem because the number of genes greatly exceed the number of tissue samples. We demonstrate how EMMIX-GENE can be useful in reducing the initial set of genes down to a more computationally manageable size. The results from this analysis also emphasise the difficulty associated with the task of separating two tissue groups on the basis of a particular subset of genes. These results also shed light on why supervised methods have such a high misallocation error rate for the breast cancer data.

Mouse proteome analysis

A general overview of the protein sequence set for the mouse transcriptome produced during the FANTOM2 sequencing project is presented here. We applied different algorithms to characterize protein sequences derived from a nonredundant representative protein set (RPS) and a variant protein set (VPS) of the mouse transcriptome. The functional characterization and assignment of Gene Ontology terms was done by analysis of the proteome using InterPro. The Superfamily database analyses gave a detailed structural classification according to SCOP and provide additional evidence for the functional characterization of the proteome data. The MDS database analysis revealed new domains which are not presented in existing protein domain databases. Thus the transcriptome gives us a unique source of data for the detection of new functional groups. The data obtained for the RPS and VPS sets facilitated the comparison of different patterns of protein expression. A comparison of other existing mouse and human protein sequence sets (e.g., the International Protein Index) demonstrates the common patterns in mammalian proteornes. The analysis of the membrane organization within the transcriptome of multiple eukaryotes provides valuable statistics about the distribution of secretory and transmembrane proteins

Two-locus Linkage Analysis Applied to Putative Quantitative Trait Loci for Lipoprotein(a) Levels

Plasma levels of lipoprotein(a) _ Lp(a) _ are associated with cardiovascular risk (Danesh et al., 2000) and were long believed to be influenced by the LPA locus on chromosome 6q27 only. However, a recent report of Broeckel et al. (2002) suggested the presence of a second quantitative trait locus on chromosome 1 influencing Lp(a) levels. Using a two-locus model, we found no evidence for an additional Lp(a) locus on chromosome 1 in a linkage study among 483 dizygotic twin pairs.

Zygosity Diagnosis in the Absence of Genotypic Data: An Approach Using Latent Class Analysis

For zygosity diagnosis in the absence of genotypic data, or in the recruitment phase of a twin study where only single twins from same-sex pairs are being screened, or to provide a test for sample duplication leading to the false identification of a dizygotic pair as monozygotic, the appropriate analysis of respondents' answers to questions about zygosity is critical. Using data from a young adult Australian twin cohort (N = 2094 complete pairs and 519 singleton twins from same-sex pairs with complete responses to all zygosity items), we show that application of latent class analysis (LCA), fitting a 2-class model, yields results that show good concordance with traditional methods of zygosity diagnosis, but with certain important advantages. These include the ability, in many cases, to assign zygosity with specified probability on the basis of responses of a single informant (advantageous when one zygosity type is being oversampled); and the ability to quantify the probability of misassignment of zygosity, allowing prioritization of cases for genotyping as well as identification of cases of probable laboratory error. Out of 242 twins (from 121 like-sex pairs) where genotypic data were available for zygosity confirmation, only a single case was identified of incorrect zygosity assignment by the latent class algorithm. Zygosity assignment for that single case was identified by the LCA as uncertain (probability of being a monozygotic twin only 76%), and the co-twin's responses clearly identified the pair as dizygotic (probability of being dizygotic 100%). In the absence of genotypic data, or as a safeguard against sample duplication, application of LCA for zygosity assignment or confirmation is strongly recommended.

Biaxial failure model for fiber reinforced concrete

Steel fiber reinforced concrete (SFRC) is widely applied in the construction industry. Numerical elastoplastic analysis of the macroscopic behavior is complex. This typically involves a piecewise linear failure curve including corner singularities. This paper presents a single smooth biaxial failure curve for SFRC based on a semianalytical approximation. Convexity of the proposed model is guaranteed so that numerical problems are avoided. The model has sufficient flexibility to closely match experimental results. The failure curve is also suitable for modeling plain concrete under biaxial loading. Since this model is capable of simulating the failure states in all stress regimes with a single envelope, the elastoplastic formulation is very concise and simple. The finite element implementation is developed to demonstrate the conciseness and the effectiveness of the model. The computed results display good agreement with published experimental data.

An approach to verifying and debugging simulation models governed by ordinary differential equations: Part 2. Residuals analysis and a case study

In Part 1 of this paper a methodology for back-to-back testing of simulation software was described. Residuals with error-dependent geometric properties were generated. A set of potential coding errors was enumerated, along with a corresponding set of feature matrices, which describe the geometric properties imposed on the residuals by each of the errors. In this part of the paper, an algorithm is developed to isolate the coding errors present by analysing the residuals. A set of errors is isolated when the subspace spanned by their combined feature matrices corresponds to that of the residuals. Individual feature matrices are compared to the residuals and classified as 'definite', 'possible' or 'impossible'. The status of 'possible' errors is resolved using a dynamic subset testing algorithm. To demonstrate and validate the testing methodology presented in Part 1 and the isolation algorithm presented in Part 2, a case study is presented using a model for biological wastewater treatment. Both single and simultaneous errors that are deliberately introduced into the simulation code are correctly detected and isolated. Copyright (C) 2003 John Wiley Sons, Ltd.

On the evaluation of some differential formulations for the pom-pom constitutive model

In recent years, progress has been made in modelling long chain branched polymers by the introduction of the so-called pompom model. Initially developed by McLeish and Larson (1998), the model has undergone several improvements or alterations, leading to the development of new formulations. Some of these formulations however suffer from certain mathematical defects. The purpose of the present paper is to review some of the formulations of the pom-pom constitutive model, and to investigate their possible mathematical defects. Next, an alternative formulation is proposed, which does not appear to exhibit mathematical defects, and we explore its modelling performance by comparing the predictions with experiments in non-trivial rheometric flows of an LDPE melt. The selected rheometric flows are the double step strain, as well as the large amplitude oscillatory shear experiments. For LAOS experiments, the comparison involves the use of Fourier-transform analysis.

Evaluating a model of parental influence on youth physical activity

Objective: To test a conceptual model linking parental physical activity orientations, parental support for physical activity, and children's self-efficacy perceptions with physical activity participation. Participants and Setting: The sample consisted of 380 students in grades 7 through 12 (mean age, 14.0 +/- 1.6 years) and their parents. Data collection took place during the fall of 1996. Main Outcome Measures: Parents completed a questionnaire assessing their physical activity habits, enjoyment of physical activity, beliefs regarding the importance of physical activity, and supportive behaviors for their child's physical activity. Students completed a 46-item inventory assessing physical activity during the previous 7 days and a 5-item physical activity self-efficacy scale. The model was tested via observed variable path analysis using structural equation modeling techniques (AMOS 4.0). Results: An initial model, in which parent physical activity orientations predicted child physical activity via parental support and child self-efficacy, did not provide an acceptable fit to the data. Inclusion of a direct path from parental support to child physical activity and deletion of a nonsignificant path from parental physical activity to child physical activity significantly improved model fit. Standardized path coefficients for the revised model ranged from 0.17 to 0.24, and all were significant at the p < 0.0001 level. Conclusions: Parental support was an important correlate of youth physical activity, acting directly or indirectly through its influence on self-efficacy. Physical activity interventions targeted at youth should include and evaluate the efficacy of individual-level and community-level strategies to increase parents' capacity to provide instrumental and motivational support for their children's physical activity.

Adsorption of supercritical fluids in non-porous and porous carbons: analysis of adsorbed phase volume and density

In this paper, we revisit the surface mass excess in adsorption studies and investigate the role of the volume of the adsorbed phase and its density in the analysis of supercritical gas adsorption in non-porous as well as microporous solids. For many supercritical fluids tested (krypton, argon, nitrogen, methane) on many different carbonaceous solids, it is found that the volume of the adsorbed phase is confined mostly to a geometrical volume having a thickness of up to a few molecular diameters. At high pressure the adsorbed phase density is also found to be very close to but never equal or greater than the liquid phase density. (C) 2003 Elsevier Science Ltd. All rights reserved.

Analysis of multicomponent adsorption kinetics on activated carbon

An integrated mathematical model for the kinetics of multicomponent adsorption on microporous carbon was developed. Transport in this bidisperse solid is represented by balance equations in the macropore and micropore phases, in which gas-phase diffusion dominates the mass transfer in the macropores, with the phenomenological diffusivities represented by the generalized Maxwell-Stefan (GMS) formulation. Viscous flow also contributes to the macropore fluxes and is included in the MS expressions. Diffusion of the adsorbed phase controls the mass transfer in the micro ore phase, p which is also described in a similar way by the MS method. The adsorption isotherms are represented by a new heterogeneous modified vacancy solution theory formulation of adsorption, which has proved to be a robust method for adsorption on activated carbons. The model is applied to the coadsorption and codesorption of C2H6 and C3H8 on Ajax and Norit carbon, as well as the displacement on Ajax carbon. The effect of the viscous flow in the macropore phase is not significant for the cases studied. The model accurately predicts the overshoot behavior and rollup of C2H6 during coadsorption. The prediction for the heavier compound C3H8 is always satisfactory, though at higher C3H8 mole fraction, the overshoot extent of C2H6 is overpredicted, possibly due to neglect of heat effects.

An analysis of helium primordial nucleosynthesis with a variable cosmological coupling

The synthesis of helium in the early Universe depends on many input parameters, including the value of the gravitational coupling during the period when the nucleosynthesis takes place. We compute the primordial abundance of helium as function of the gravitational coupling, using a semi-analytical method, in order to track the influence of G in the primordial nucleosynthesis. To be specific, we construct a cosmological model with varying G, using the Brans-Dicke theory. The greater the value of G at nucleosynthesis period, the greater the predicted abundance of helium. Using the observational data for the abundance of primordial helium, constraints for the time variation of G are established.

Integrating environmental and economic life cycle analysis in product development: a material selection case study

Purpose Achieving sustainability by rethinking products, services and strategies is an enormous challenge currently laid upon the economic sector, in which materials selection plays a critical role. In this context, the present work describes an environmental and economic life cycle analysis of a structural product, comparing two possible material alternatives. The product chosen is a storage tank, presently manufactured in stainless steel (SST) or in a glass fibre reinforced polymer composite (CST). The overall goal of the study is to identify environmental and economic strong and weak points related to the life cycle of the two material alternatives. The consequential win-win or trade-off situations will be identified via a Life Cycle Assessment/Life Cycle Costing (LCA/LCC) integrated model. Methods The LCA/LCC integrated model used consists in applying the LCA methodology to the product system, incorporating, in parallel, its results into the LCC study, namely those of the Life Cycle Inventory (LCI) and the Life Cycle Impact Assessment (LCIA). Results In both the SST and CST systems the most significant life cycle phase is the raw materials production, in which the most significant environmental burdens correspond to the Fossil fuels and Respiratory inorganics categories. The LCA/LCC integrated analysis shows that the CST has globally a preferable environmental and economic profile, as its impacts are lower than those of the SST in all life cycle stages. Both the internal and external costs are lower, the former resulting mainly from the composite material being significantly less expensive than stainless steel. This therefore represents a full win-win situation. As a consequence, the study clearly indicates that using a thermoset composite material to manufacture storage tanks is environmentally and economically desirable. However, it was also evident that the environmental performance of the CST could be improved by altering its End-of-Life stage. Conclusions The results of the present work provide enlightening insights into the synergies between the environmental and the economic performance of a structural product made with alternative materials. Further, they provide conclusive evidence to support the integration of environmental and economic life cycle analysis in the product development processes of a manufacturing company, or in some cases even in its procurement practices.

Travel preferences of the residents of Guimarães: a revised model of push and pull motives

The historic center of the Portuguese city of Guimarães is a world heritage site (UNESCO) since 2001, having hosted the European Capital of Culture (ECOC) in 2012. In this sense, Guimarães has made a major effort in promoting tourism, positioning itself as an urban and cultural tourism destination. The present paper has two objectives. The first, to examine if an existing push and pull motivation model finds statistical support with regard to the population of the municipality of Guimarães, a cultural tourism destination. The second, to study the role that important socio-demographic variables, such as gender, age, and education, play in determining travel motivations of residents from this municipality. Insight on tourism motivation may be an important policy tool for tourism planners and managers in the development of products and marketing strategies. The empirical analysis is undertaken based on questionnaires administered in 2012 to residents of Guimarães. The present study shows that gender, age and education make a difference with regard to travel motivations.

Estimating Local Part Thickness in Midplane Meshes for Finite Element Analysis

Within the development of motor vehicles, crash safety (e.g. occupant protection, pedestrian protection, low speed damageability), is one of the most important attributes. In order to be able to fulfill the increased requirements in the framework of shorter cycle times and rising pressure to reduce costs, car manufacturers keep intensifying the use of virtual development tools such as those in the domain of Computer Aided Engineering (CAE). For crash simulations, the explicit finite element method (FEM) is applied. The accuracy of the simulation process is highly dependent on the accuracy of the simulation model, including the midplane mesh. One of the roughest approximations typically made is the actual part thickness which, in reality, can vary locally. However, almost always a constant thickness value is defined throughout the entire part due to complexity reasons. On the other hand, for precise fracture analysis within FEM, the correct thickness consideration is one key enabler. Thus, availability of per element thickness information, which does not exist explicitly in the FEM model, can significantly contribute to an improved crash simulation quality, especially regarding fracture prediction. Even though the thickness is not explicitly available from the FEM model, it can be inferred from the original CAD geometric model through geometric calculations. This paper proposes and compares two thickness estimation algorithms based on ray tracing and nearest neighbour 3D range searches. A systematic quantitative analysis of the accuracy of both algorithms is presented, as well as a thorough identification of particular geometric arrangements under which their accuracy can be compared. These results enable the identification of each technique’s weaknesses and hint towards a new, integrated, approach to the problem that linearly combines the estimates produced by each algorithm.