Simulation of granular soil behaviour using the bullet physics library

A physics engine is computer software which provides a simulation of certain physical systems, such as rigid body dynamics, soft body dynamics and fluid dynamics. Physics engines were firstly developed for using in animation and gaming industry ; nevertheless, due to fast calculation speed they are attracting more and more attetion from researchers of the engineering fields. Since physics engines are capable of performing fast calculations on multibody rigid dynamic systems, soil particles can be modeled as distinct rigid bodies. However, up to date, it is not clear to what extent they perform accurately in modeling soil behaviour from a geotechnical viewpoint. To investigate this, examples of pluviation and vibration-induced desification were simulated using the physics engine called Bullet physics library. In order to create soil samples, first, randomly shaped polyhedrons, representing gravels, were generated using the Voronoi tessellation approach. Then, particles were pluviated through a funnel into a cylinder. Once the soil particles settled in a static state, the cylinder was subjected to horizontal sinusoidal vibration for a period of 20 seconds. The same procedure for sample perparation was performed in the laboratory. The results of pluviation and vibration tests weere recorded and compared to those of simulations. A good agreement has been found between the results of simulations and laboratory tests. The findings in this study reinforce the idea that physics engines can be employed as a geotechnical engineering simulation tool

Simulating Tsunami Inundation and Soil Response in a Large Centrifuge.

Tsunamis are rare, extreme events and cause significant damage to coastal infrastructure, which is often exacerbated by soil instability surrounding the structures. Simulating tsunamis in a laboratory setting is important to further understand soil instability induced by tsunami inundation processes. Laboratory simulations are difficult because the scale of such processes is very large, hence dynamic similitude cannot be achieved for small-scale models in traditional water-wave-tank facilities. The ability to control the body force in a centrifuge environment considerably reduces the mismatch in dynamic similitude. We review dynamic similitudes under a centrifuge condition for a fluid domain and a soil domain. A novel centrifuge apparatus specifically designed for exploring the physics of a tsunami-like flow on a soil bed is used to perform experiments. The present 1:40 model represents the equivalent geometric scale of a prototype soil field of 9.6 m deep, 21 m long, and 14.6 m wide. A laboratory facility capable of creating such conditions under the normal gravitational condition does not exist. With the use of a centrifuge, we are now able to simulate and measure tsunami-like loading with sufficiently high water pressure and flow velocities. The pressures and flow velocities in the model are identical to those of the prototype yielding realistic conditions of flow-soil interaction

Survey of methods for soil moisture determination

Existing and proposed methods for soil moisture determination are discussed. These include: (1) in situ investigations including gravimetric, nuclear, and electromagnetic techniques; (2) remote sensing approaches that use the reflected solar, thermal infrared, and microwave portions of the electromagnetic spectrum; and (3) soil physics models that track the behavior of water in the soil in response to meteorological inputs (precipitation) and demands (evapotranspiration). The capacities of these approaches to satisfy various user needs for soil moisture information vary from application to application, but a conceptual scheme for merging these approaches into integrated systems to provide soil moisture information is proposed that has the potential for meeting various application requirements

Some aspects of soil physics applicable to trickle irrigation : a thesis presented in partial fulfilment of the requirements for the degree of Master of Science in Soil Science at Massey University

Irrigation of crops is one of the more widely used techniques to increase yeilds. Trickle irrigation is one such method and is more suited to horticultural crops. In New Zealand, with horticulture assuming more importance, appropriate methods of design and operation of trickle irrigation systems are required. In this study a simple approximation to Wooding's solution for steady infiltration from a shallow ponded source, much like that found under trickle emitters is examined. This may aid in irrigation design and practice. The approximation also allowed for the development of a method to concurrently measure the saturated hydraulic conductivity and sorptivity from simple field infiltration measurements with a minimum of soil disturbance. Saturated hydraulic conductivities and sorptivities are of great use in soil water studies in general. A commercial trickle irrigation system was also examined to determine the suitability of such irrigation systems to particular soils, and to examine the present irrigation scheduling. The approximation to Wooding's solution was found to perform well in the field in many respects, particularly in determining steady ponded zone sizes. Ponded zone sizes are important in that they control the volume of soil wetted by irrigation to a large degree. Much of this agreement is due to the use of parameters determined by the simple field method developed from this theory. Sorptivities and saturated hydraulic conductivities obtained by this method were found to be more realistic for trickle irrigation than those determined by other existing methods. Systematic errors in these other methods, mainly soil disturbance and the concomitant creation of continuous flow paths for water, as well as soil smearing, are thought to be the main cause of this difference. Temporal and spatial variation in soil physical properties are however, found to hinder the use of soil physics theory in the field. Macropores (due to soil biological activity) were found to profoundly influence infiltration processes and soil-water distribution. These effects were particularly marked for the site with a commercial trickle irrigation system. Here the efficiency of the present system is thought to be low, and evidence indicates that irrigation was in excess of plant requirements. The utility of Wooding's solution, and the method to measure soil physical parameters developed from this, is further demonstrated in this orchard

Improving Simulated Soil Moisture Fields Through Assimilation of AMSR-E Soil Moisture Retrievals with an Ensemble Kalman Filter and a Mass Conservation Constraint

Model simulated soil moisture fields are often biased due to errors in input parameters and deficiencies in model physics. Satellite derived soil moisture estimates, if retrieved appropriately, represent the spatial mean of soil moisture in a footprint area, and can be used to reduce model bias (at locations near the surface) through data assimilation techniques. While assimilating the retrievals can reduce model bias, it can also destroy the mass balance enforced by the model governing equation because water is removed from or added to the soil by the assimilation algorithm. In addition, studies have shown that assimilation of surface observations can adversely impact soil moisture estimates in the lower soil layers due to imperfect model physics, even though the bias near the surface is decreased. In this study, an ensemble Kalman filter (EnKF) with a mass conservation updating scheme was developed to assimilate the actual value of Advanced Microwave Scanning Radiometer (AMSR-E) soil moisture retrievals to improve the mean of simulated soil moisture fields by the Noah land surface model. Assimilation results using the conventional and the mass conservation updating scheme in the Little Washita watershed of Oklahoma showed that, while both updating schemes reduced the bias in the shallow root zone, the mass conservation scheme provided better estimates in the deeper profile. The mass conservation scheme also yielded physically consistent estimates of fluxes and maintained the water budget. Impacts of model physics on the assimilation results are discussed

Stellar structure and compact objects before 1940: Towards relativistic astrophysics

Since the mid-1920s, different strands of research used stars as "physics laboratories" for investigating the nature of matter under extreme densities and pressures, impossible to realize on Earth. To trace this process this paper is following the evolution of the concept of a dense core in stars, which was important both for an understanding of stellar evolution and as a testing ground for the fast-evolving field of nuclear physics. In spite of the divide between physicists and astrophysicists, some key actors working in the cross-fertilized soil of overlapping but different scientific cultures formulated models and tentative theories that gradually evolved into more realistic and structured astrophysical objects. These investigations culminated in the first contact with general relativity in 1939, when J. Robert Oppenheimer and his students George Volkoff and Hartland Snyder systematically applied the theory to the dense core of a collapsing neutron star. This pioneering application of Einstein's theory to an astrophysical compact object can be regarded as a milestone in the path eventually leading to the emergence of relativistic astrophysics in the early 1960s.Comment: 83 pages, 4 figures, submitted to the European Physical Journal

Some effects of topography, soil moisture, and sea-surface temperature on continental precipitation as computed with the GISS coarse mesh climate model

The effects of terrain elevation, soil moisture, and zonal variations in sea/surface temperature on the mean daily precipitation rates over Australia, Africa, and South America in January were evaluated. It is suggested that evaporation of soil moisture may either increase or decrease the model generated precipitation, depending on the surface albedo. It was found that a flat, dry continent model best simulates the January rainfall over Australia and South America, while over Africa the simulation is improved by the inclusion of surface physics, specifically soil moisture and albedo variations

Measuring very negative water potentials with polymer tensiometers: principles, performance and applications

In recent years, a polymer tensiometer (POT) was developed and tested to directly measure matric potentials in dry soils. By extending the measurement range to wilting point (a 20-fold increase compared to conventional, water-filled tensiometers), a myriad of previously unapproachable research questions are now open to experimental exploration. Furthermore, the instrument may well allow the development of more water-efficient irrigation strategies by recording water potential rather than soil water content. The principle of the sensor is to fill it with a polymer solution instead of water, thereby building up osmotic pressure inside the sensor. A high-quality ceramic allows the exchange of water with the soil while retaining the polymer. The ceramic has pores sufficiently small to remain saturated even under very negative matric potentials. Installing the sensor in an unsaturated soil causes the high pressure of the polymer solution to drop as the water potentials in the soil and in the POT equilibrate. As long as the pressure inside the polymer chamber remains sufficiently large to prevent cavitation, the sensor will function properly. If the osmotic potential in the polymer chamber can produce a pressure of approximately 2.0 MPa when the sensor is placed in water, proper readings down to wilting point are secured. Various tests in disturbed soil, including an experiment with root water uptake, demonstrate the operation and performance of the new polymer tensiometer and illustrate how processes such as root water uptake can be studied in more detail than before. The paper discusses the available data and explores the long term perspectives offered by the instrument