Numerical modeling of the effect of variation of boundary conditions on vadose zone hydraulic properties

An accurate estimation of hydraulic fluxes in the vadose zone is essential for the prediction of water, nutrient and contaminant transport in natural systems. The objective of this study was to simulate the effect of variation of boundary conditions on the estimation of hydraulic properties (i.e. water content, effective unsaturated hydraulic conductivity and hydraulic flux) in a one-dimensional unsaturated flow model domain. Unsaturated one-dimensional vertical water flow was simulated in a pure phase clay loam profile and in clay loam interlayered with silt loam distributed according to the third iteration of the Cantor Bar fractal object Simulations were performed using the numerical model Hydrus 1D. The upper and lower pressure heads were varied around average values of -55 cm for the near-saturation range. This resulted in combinations for the upper and lower constant head boundary conditions, respectively, of -50 and -60 cm, -40 and -70 cm, -30 and -80 cm, -20 and -90 cm, and -10 and -100 cm. For the drier range the average head between the upper and lower boundary conditions was set to -550 cm, resulting in the combinations -500 and -600 cm, -400 and -700 cm, -300 and -800 cm, -200 and -900 cm, and -100 and -1,000 cm, for upper and lower boundary conditions, respectively. There was an increase in water contents, fluxes and hydraulic conductivities with the increase in head difference between boundary conditions. Variation in boundary conditions in the pure phase and interlayered one-dimensional profiles caused significant deviations in fluxes, water contents and hydraulic conductivities compared to the simplest case (a head difference between the upper and lower constant head boundaries of 10 cm in the wetter range and 100 cm in the drier range).

Comparison between molecular dynamics abd Monte Carlo simulations of an ordering process in a binary alloy

Ordering in a binary alloy is studied by means of a molecular-dynamics (MD) algorithm which allows to reach the domain growth regime. Results are compared with Monte Carlo simulations using a realistic vacancy-atom (MC-VA) mechanism. At low temperatures fast growth with a dynamical exponent x>1/2 is found for MD and MC-VA. The study of a nonequilibrium ordering process with the two methods shows the importance of the nonhomogeneity of the excitations in the system for determining its macroscopic kinetics.

Soil phosphorus dynamics and availability and irrigated coffee yield

Research data have demonstrated that the P demand of coffee (Coffea arabica L.) is similar to that of short-cycle crops. In this context, the objective of this study was to evaluate the influence of annual P fertilization on the soil P status by the quantification of labile, moderately labile, low-labile, and total P fractions, associating them to coffee yield. The experiment was installed in a typical dystrophic Red Latosol (Oxisol) cultivated with irrigated coffee annually fertilized with triple superphosphate at rates of 0, 50, 100, 200, and 400 kg ha-1 P2O5. Phosphorus fractions were determined in two soil layers: 0-10 and 10-20 cm. The P leaf contents and coffee yield in 2008 were also evaluated. The irrigated coffee responded to phosphate fertilization in the production phase with gains of up to 138 % in coffee yield by the application of 400 kg ha-1 P2O5. Coffee leaf P contents increased with P applications and stabilized around 1.98 g kg-1, at rates of 270 kg ha-1 P2O5 and higher. Soil P application caused, in general, an increase in bioavailable P fractions, which constitute the main soil P reservoir.

Finite-size scaling analysis of the avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics

A numerical study is presented of the third-dimensional Gaussian random-field Ising model at T=0 driven by an external field. Standard synchronous relaxation dynamics is employed to obtain the magnetization versus field hysteresis loops. The focus is on the analysis of the number and size distribution of the magnetization avalanches. They are classified as being nonspanning, one-dimensional-spanning, two-dimensional-spanning, or three-dimensional-spanning depending on whether or not they span the whole lattice in different space directions. Moreover, finite-size scaling analysis enables identification of two different types of nonspanning avalanches (critical and noncritical) and two different types of three-dimensional-spanning avalanches (critical and subcritical), whose numbers increase with L as a power law with different exponents. We conclude by giving a scenario for avalanche behavior in the thermodynamic limit.

Spanning avalanches in the three-dimensional Gaussian random-field Ising model with metastable dynamics: Field dependence and geometrical properties

Spanning avalanches in the 3D Gaussian Random Field Ising Model (3D-GRFIM) with metastable dynamics at T=0 have been studied. Statistical analysis of the field values for which avalanches occur has enabled a Finite-Size Scaling (FSS) study of the avalanche density to be performed. Furthermore, a direct measurement of the geometrical properties of the avalanches has confirmed an earlier hypothesis that several types of spanning avalanches with two different fractal dimensions coexist at the critical point. We finally compare the phase diagram of the 3D-GRFIM with metastable dynamics with the same model in equilibrium at T=0.

Hysteresis and avalanches in the T=0 random field ising model with 2-spin flip dynamics

We study the nonequilibrium behavior of the three-dimensional Gaussian random-field Ising model at T=0 in the presence of a uniform external field using a two-spin-flip dynamics. The deterministic, history-dependent evolution of the system is compared with the one obtained with the standard one-spin-flip dynamics used in previous studies of the model. The change in the dynamics yields a significant suppression of coercivity, but the distribution of avalanches (in number and size) stays remarkably similar, except for the largest ones that are responsible for the jump in the saturation magnetization curve at low disorder in the thermodynamic limit. By performing a finite-size scaling study, we find strong evidence that the change in the dynamics does not modify the universality class of the disorder-induced phase transition.

Diluted three-dimensional random field Ising model at zero temperature with metastable dynamics.

The influence of vacancy concentration on the behavior of the three-dimensional random field Ising model with metastable dynamics is studied. We have focused our analysis on the number of spanning avalanches which allows us a clear determination of the critical line where the hysteresis loops change from continuous to discontinuous. By a detailed finite-size scaling analysis we determine the phase diagram and numerically estimate the critical exponents along the whole critical line. Finally, we discuss the origin of the curvature of the critical line at high vacancy concentration.

Influence of the driving mechanism on the response of systems with athermal dynamics: The example of the random-field Ising model

We investigate the influence of the driving mechanism on the hysteretic response of systems with athermal dynamics. In the framework of local mean-field theory at finite temperature (but neglecting thermally activated processes), we compare the rate-independent hysteresis loops obtained in the random field Ising model when controlling either the external magnetic field H or the extensive magnetization M. Two distinct behaviors are observed, depending on disorder strength. At large disorder, the H-driven and M-driven protocols yield identical hysteresis loops in the thermodynamic limit. At low disorder, when the H-driven magnetization curve is discontinuous (due to the presence of a macroscopic avalanche), the M-driven loop is reentrant while the induced field exhibits strong intermittent fluctuations and is only weakly self-averaging. The relevance of these results to the experimental observations in ferromagnetic materials, shape memory alloys, and other disordered systems is discussed.

Clear Zone – A Synthesis of Practice and an Evaluation of the Benefi ts of Meeting the 10-ft Clear Zone Goal on Urban Streets, November 2008

In urban communities, there are often limited amounts of right-of-way available for establishing a large setback distance from the curb for fixed objects. Urban communities must constantly weigh the cost of purchasing additional right-of-way for clear zones against the risk of fixed object crashes. From 2004 to 2006, this type of crash on curbed roads represented 15% of all fatal crashes and 3% of all crashes in the state of Iowa. Many states have kept the current minimum AASHTO recommendations as their minimum clear zone standards; however, other states have decided that these recommendations are insufficient and have increased the required minimum clear zone distance to better suit the judgment of local designers. This report presents research on the effects of the clear zone on urban curbed streets. The research was conducted in two phases. The first phase involved a synthesis of practice that included a literature review and a survey of practices in jurisdictions that have developmental and historical patterns similar to those of Iowa. The second phase involved investigating the benefits of a 10 ft clear zone, which included examining urban corridors in Iowa that meet or do not meet the 10 ft clear zone goal. The results of this study indicate that a consistent fixed object offset results in a reduction in the number of fixed object crashes, a 5 ft clear zone is most effective when the goal is to minimize the number of fixed object c ashes, and a 3 ft clear zone is most effective when the goal is to minimize the cost of fixed object crashes.

Statistics of transverse mode turn-on dynamics in VCSEL's

The turn-on process of a multimode VCSEL is investigated from a statistical point of view. Special attention is paid to quantities such as time jitter and bit error rate. The single-mode performance of VCSEL¿s during current modulation is compared to that of edge-emitting lasers.

Examining Driver Behavior in Response to Work Zone Interventions: A Driving Simulator Study, 2010

Construction zones pose a significant threat to both workers and drivers causing numerous injuries and deaths each year. Innovations in work zone safety could reduce these numbers. However, implementing work zone interventions before they are validated can undermine rather than enhance safety. The objective of this research is to demonstrate how driving simulators can be used to evaluate the effect of various work zone interventions on driver performance.

Application of Visualization and Simulation Program to Improve Work Zone Safety and Mobility, January 2010

A previous study sponsored by the Smart Work Zone Deployment Initiative, “Feasibility of Visualization and Simulation Applications to Improve Work Zone Safety and Mobility,” demonstrated the feasibility of combining readily available, inexpensive software programs, such as SketchUp and Google Earth, with standard two-dimensional civil engineering design programs, such as MicroStation, to create animations of construction work zones. The animations reflect changes in work zone configurations as the project progresses, representing an opportunity to visually present complex information to drivers, construction workers, agency personnel, and the general public. The purpose of this study is to continue the work from the previous study to determine the added value and resource demands created by including more complex data, specifically traffic volume, movement, and vehicle type. This report describes the changes that were made to the simulation, including incorporating additional data and converting the simulation from a desktop application to a web application.