Simulation supplements field studies to determine no-till dryland corn population recommendations for semiarid western Nebraska

In a 2-yr multiple-site field study conducted in western Nebraska during 1999 and 2000, optimum dryland corn (Zea mays L.) population varied from less than 1.7 to more than 5.6 plants m(-2), depending largely on available water resources. The objective of this study was to use a modeling approach to investigate corn population recommendations for a wide range of seasonal variation. A corn growth simulation model (APSIM-maize) was coupled to long-term sequences of historical climatic data from western Nebraska to provide probabilistic estimates of dryland yield for a range of corn populations. Simulated populations ranged from 2 to 5 plants m(-2). Simulations began with one of three levels of available soil water at planting, either 80, 160, or 240 mm in the surface 1.5 m of a loam soil. Gross margins were maximized at 3 plants m(-2) when starting available water was 160 or 240 mm, and the expected probability of a financial loss at this population was reduced from about 10% at 160 mm to 0% at 240 mm. When starting available water was 80 mm, average gross margins were less than $15 ha(-1), and risk of financial loss exceeded 40%. Median yields were greatest when starting available soil water was 240 mm. However, perhaps the greater benefit of additional soil water at planting was reduction in the risk of making a financial loss. Dryland corn growers in western Nebraska are advised to use a population of 3 plants m(-2) as a base recommendation.

Effect of surface-perturbed intermolecular interaction on adsorption of simple gases on a graphitized carbon surface

In this paper, we investigate the effect of the solid surface on the fluid-fluid intermolecular potential energy. This modified fluid-fluid interaction energy due to the inducement of a solid surface is used in the grand canonical Monte Carlo (GCMC) simulation of various noble gases, nitrogen, and methane on graphitized thermal carbon black. This effect is such that the effective interaction potential energy between two particles close to surface is less than the potential energy if the solid substrate is not present. With this modification the GCMC simulation results agree extremely well with the experimental data over a wide range of pressures while the simulation results with the unmodified potential energy give rise to a shoulder near the neighborhood of monolayer coverage and the significant overprediction of the second and higher layer coverages. In particular the unmodified GCMC results exhibit very sharp change in those higher layers while the experimental data have a much gradual change in the uptake. We will illustrate this theory with adsorption data of argon, xenon, neon, nitrogen, and methane on graphitized thermal carbon black.

Adsorption of quadrupolar, diatomic nitrogen onto graphitized thermal carbon black and in slit-shaped carbon pores. Effects of surface mediation

In this paper, we study the effect of solid surface mediation on the intermolecular potential energy of nitrogen, and its impact on the adsorption of nitrogen on a graphitized carbon black surface and in carbon slit-shaped pores. This effect arises from the lower effective interaction potential energy between two particles close to the surface compared to the potential energy of the same two particles when they are far away from the surface. A simple equation is proposed to calculate the reduction factor and this is used in the Grand Canonical Monte Carlo (GCMC) simulation of nitrogen adsorption on graphitized thermal carbon black. With this modification, the GCMC simulation results agree extremely well with the experimental data over a wide range of pressure; the simulation results with the original potential energy (i.e. no surface mediation) give rise to a shoulder in the neighbourhood of monolayer coverage and a significant over-prediction of the second and higher layer coverages. The influence of this surface mediation on the dependence of the pore-filling pressure on the pore width is also studied. It is shown that such surface mediation has a significant effect on the pore-filling pressure. This implies that the use of the local isotherms obtained from the potential model without surface mediation could give rise to a serious error in the determination of the pore-size distribution.

Simulation of galvanic corrosion of magnesium coupled to a steel fastener in NaCl solution

The galvanic corrosion of magnesium alloy AZ91D coupled to a steel fastener was studied using a boundary element method (BEM) model and experimental measurements. The BEM model used the measured polarization curves as boundary conditions. The experimental program involved measuring the total corrosion rate as a function of distance from the interface of the magnesium in the form of a sheet containing a mild steel circular insert (5 to 30 mm in diameter). The measured total corrosion rate was interpreted as due to galvanic corrosion plus self corrosion. For a typical case, the self corrosion was estimated typically to be similar to 230 mm/y for an area surrounding the interface and to a distance of about I cm from the interface. Scanning Kelvin Probe Force Microscopy (SKPFM) revealed microgalvanic cells with potential differences of approximately 100 mV across the AZ91D surface. These microgalvanic cells may influence the relative contributions of galvanic and self corrosion to the total corrosion of AZ91D.

A master equation model for bimolecular reaction via multi-well isomerization intermediates

A reversible linear master equation model is presented for pressure- and temperature-dependent bimolecular reactions proceeding via multiple long-lived intermediates. This kinetic treatment, which applies when the reactions are measured under pseudo-first-order conditions, facilitates accurate and efficient simulation of the time dependence of the populations of reactants, intermediate species and products. Detailed exploratory calculations have been carried out to demonstrate the capabilities of the approach, with applications to the bimolecular association reaction C3H6 + H reversible arrow C3H7 and the bimolecular chemical activation reaction C2H2 +(CH2)-C-1--> C3H3+H. The efficiency of the method can be dramatically enhanced through use of a diffusion approximation to the master equation, and a methodology for exploiting the sparse structure of the resulting rate matrix is established.

Improving wheat simulation capabilities in Australia from a cropping systems perspective - II. Testing simulation capabilities of wheat growth

To simulate cropping systems, crop models must not only give reliable predictions of yield across a wide range of environmental conditions, they must also quantify water and nutrient use well, so that the status of the soil at maturity is a good representation of the starting conditions for the next cropping sequence. To assess the suitability for this task a range of crop models, currently used in Australia, were tested. The models differed in their design objectives, complexity and structure and were (i) tested on diverse, independent data sets from a wide range of environments and (ii) model components were further evaluated with one detailed data set from a semi-arid environment. All models were coded into the cropping systems shell APSIM, which provides a common soil water and nitrogen balance. Crop development was input, thus differences between simulations were caused entirely by difference in simulating crop growth. Under nitrogen non-limiting conditions between 73 and 85% of the observed kernel yield variation across environments was explained by the models. This ranged from 51 to 77% under varying nitrogen supply. Water and nitrogen effects on leaf area index were predicted poorly by all models resulting in erroneous predictions of dry matter accumulation and water use. When measured light interception was used as input, most models improved in their prediction of dry matter and yield. This test highlighted a range of compensating errors in all modelling approaches. Time course and final amount of water extraction was simulated well by two models, while others left up to 25% of potentially available soil water in the profile. Kernel nitrogen percentage was predicted poorly by all models due to its sensitivity to small dry matter changes. Yield and dry matter could be estimated adequately for a range of environmental conditions using the general concepts of radiation use efficiency and transpiration efficiency. However, leaf area and kernel nitrogen dynamics need to be improved to achieve better estimates of water and nitrogen use if such models are to be use to evaluate cropping systems. (C) 1998 Elsevier Science B.V.

Improving wheat simulation capabilities in Australia from a cropping systems perspective - III. The integrated wheat model (I_WHEAT)

Previous work has identified several short-comings in the ability of four spring wheat and one barley model to simulate crop processes and resource utilization. This can have important implications when such models are used within systems models where final soil water and nitrogen conditions of one crop define the starting conditions of the following crop. In an attempt to overcome these limitations and to reconcile a range of modelling approaches, existing model components that worked demonstrably well were combined with new components for aspects where existing capabilities were inadequate. This resulted in the Integrated Wheat Model (I_WHEAT), which was developed as a module of the cropping systems model APSIM. To increase predictive capability of the model, process detail was reduced, where possible, by replacing groups of processes with conservative, biologically meaningful parameters. I_WHEAT does not contain a soil water or soil nitrogen balance. These are present as other modules of APSIM. In I_WHEAT, yield is simulated using a linear increase in harvest index whereby nitrogen or water limitations can lead to early termination of grainfilling and hence cessation of harvest index increase. Dry matter increase is calculated either from the amount of intercepted radiation and radiation conversion efficiency or from the amount of water transpired and transpiration efficiency, depending on the most limiting resource. Leaf area and tiller formation are calculated from thermal time and a cultivar specific phyllochron interval. Nitrogen limitation first reduces leaf area and then affects radiation conversion efficiency as it becomes more severe. Water or nitrogen limitations result in reduced leaf expansion, accelerated leaf senescence or tiller death. This reduces the radiation load on the crop canopy (i.e. demand for water) and can make nitrogen available for translocation to other organs. Sensitive feedbacks between light interception and dry matter accumulation are avoided by having environmental effects acting directly on leaf area development, rather than via biomass production. This makes the model more stable across environments without losing the interactions between the different external influences. When comparing model output with models tested previously using data from a wide range of agro-climatic conditions, yield and biomass predictions were equal to the best of those models, but improvements could be demonstrated for simulating leaf area dynamics in response to water and nitrogen supply, kernel nitrogen content, and total water and nitrogen use. I_WHEAT does not require calibration for any of the environments tested. Further model improvement should concentrate on improving phenology simulations, a more thorough derivation of coefficients to describe leaf area development and a better quantification of some processes related to nitrogen dynamics. (C) 1998 Elsevier Science B.V.

Small molecular probes for G-protein-coupled C5a receptors: Conformationally constrained antagonists derived from the C terminus of the human plasma protein C5a

Activation of the human complement system of plasma proteins in response to infection or injury produces a 4-helix bundle glycoprotein (74 amino acids) known as C5a. C5a binds to G-protein-coupled receptors on cell surfaces triggering receptor-ligand internalization, signal transduction, and powerful inflammatory responses. Since excessive levels of C5a are associated with autoimmune and chronic inflammatory disorders, inhibitors of receptor activation may have therapeutic potential. We now report solution structures and receptor-binding and antagonist activities for some of the first small molecule antagonists of C5a derived from its hexapeptide C terminus. The antagonist NMe-Phe-Lys-Pro-D-Cha-Trp-D-Arg-CO2H (1) surprisingly shows an unusually well-defined solution structure as determined by H-1 NMR spectroscopy. This is one of the smallest acyclic peptides found to possess a defined solution conformation, which can be explained by the constraining role of intramolecular hydrogen bonding. NOE and coupling constant data, slow deuterium exchange, and a low dependence on temperature for the chemical shift of the D-Cha-NH strongly indicate an inverse gamma turn stabilized by a D-Cha-NH ... OC-Lys hydrogen bond. Smaller conformational populations are associated with a hydrogen bond between Trp-NH ... OC-Lys, defining a type II beta turn distorted by the inverse gamma turn incorporated within it. An excellent correlation between receptor-affinity and antagonist activity is indicated for a limited set of synthetic peptides. Conversion of the C-terminal carboxylate of 1 to an amide decreases antagonist potency 5-fold, but potency is increased up to 10-fold over 1 if the amide bond is made between the C-terminal carboxylate and a Lys/Orn side chain to form a cyclic analogue. The solution structure of cycle 6 also shows gamma and beta turns; however, the latter occurs in a different position, and there are clear conformational changes in 6 vs 1 that result in enhanced activity. These results indicate that potent C5a antagonists can be developed by targeting site 2 alone of the C5a receptor and define a novel pharmacophore for developing powerful receptor probes or drug candidates.

Flush air data system calibration using numerical simulation

The use of computational fluid dynamics simulations for calibrating a flush air data system is described, In particular, the flush air data system of the HYFLEX hypersonic vehicle is used as a case study. The HYFLEX air data system consists of nine pressure ports located flush with the vehicle nose surface, connected to onboard pressure transducers, After appropriate processing, surface pressure measurements can he converted into useful air data parameters. The processing algorithm requires an accurate pressure model, which relates air data parameters to the measured pressures. In the past, such pressure models have been calibrated using combinations of flight data, ground-based experimental results, and numerical simulation. We perform a calibration of the HYFLEX flush air data system using computational fluid dynamics simulations exclusively, The simulations are used to build an empirical pressure model that accurately describes the HYFLEX nose pressure distribution ol cr a range of flight conditions. We believe that computational fluid dynamics provides a quick and inexpensive way to calibrate the air data system and is applicable to a broad range of flight conditions, When tested with HYFLEX flight data, the calibrated system is found to work well. It predicts vehicle angle of attack and angle of sideslip to accuracy levels that generally satisfy flight control requirements. Dynamic pressure is predicted to within the resolution of the onboard inertial measurement unit. We find that wind-tunnel experiments and flight data are not necessary to accurately calibrate the HYFLEX flush air data system for hypersonic flight.

RWMODEL II: Computer simulation of the Rescorla-Wagner model of Pavlovian conditioning

RWMODEL II simulates the Rescorla-Wagner model of Pavlovian conditioning. It is written in Delphi and runs under Windows 3.1 and Windows 95. The program was designed for novice and expert users and can be employed in teaching, as well as in research. It is user friendly and requires a minimal level of computer literacy but is sufficiently flexible to permit a wide range of simulations. It allows the display of empirical data, against which predictions from the model can be validated.

Coupled experimental and thermodynamic modeling studies for metallurgical smelting and coal combustion slag systems

An extensive research program focused on the characterization of various metallurgical complex smelting and coal combustion slags is being undertaken. The research combines both experimental and thermodynamic modeling studies. The approach is illustrated by work on the PbO-ZnO-Al2O3-FeO-Fe2O3-CaO-SiO2 system. Experimental measurements of the liquidus and solidus have been undertaken under oxidizing and reducing conditions using equilibration, quenching, and electron probe X-ray microanalysis. The experimental program has been planned so as to obtain data for thermodynamic model development as well as for pseudo-ternary Liquidus diagrams that can be used directly by process operators. Thermodynamic modeling has been carried out using the computer system FACT, which contains thermodynamic databases with over 5000 compounds and evaluated solution models. The FACT package is used for the calculation of multiphase equilibria in multicomponent systems of industrial interest. A modified quasi-chemical solution model is used for the liquid slag phase. New optimizations have been carried out, which significantly improve the accuracy of the thermodynamic models for lead/zinc smelting and coal combustion processes. Examples of experimentally determined and calculated liquidus diagrams are presented. These examples provide information of direct relevance to various metallurgical smelting and coal combustion processes.

Hybrid modelling of coupled pore fluid-solid deformation problems

A hybrid formulation for coupled pore fluid-solid deformation problems is proposed. The scheme is a hybrid in the sense that we use a vertex centered finite volume formulation for the analysis of the pore fluid and a particle method for the solid in our model. The pore fluid formally occupies the same space as the solid particles. The size of the particles is not necessarily equal to the physical size of materials. A finite volume mesh for the pore fluid flow is generated by Delaunay triangulation. Each triangle possesses an initial porosity. Changes of the porosity are specified by the translations of the mass centers of particles. Net pore pressure gradients are applied to the particle centers and are considered in the particle momentum balance. The potential of our model is illustrated by means of a simulation of coupled fracture and fluid flow developed in porous rock under biaxial compression condition.

Computer simulations of coupled problems in geological and geochemical systems

The finite element method is used to simulate coupled problems, which describe the related physical and chemical processes of ore body formation and mineralization, in geological and geochemical systems. The main purpose of this paper is to illustrate some simulation results for different types of modelling problems in pore-fluid saturated rock masses. The aims of the simulation results presented in this paper are: (1) getting a better understanding of the processes and mechanisms of ore body formation and mineralization in the upper crust of the Earth; (2) demonstrating the usefulness and applicability of the finite element method in dealing with a wide range of coupled problems in geological and geochemical systems; (3) qualitatively establishing a set of showcase problems, against which any numerical method and computer package can be reasonably validated. (C) 2002 Published by Elsevier Science B.V.