Simulation modeling and resource management

Simulation models are becoming increasingly important as tools for synthesizing and applying information in almost all aspects of land management. They are particulary valuable for predicting and comparing outcomes of alternative decisions and assumptions. Models also permit managers to consider and integrate the potential influences of a large number of variables.     

The role of stand history in assessing forest impacts

Air pollution, harvesting practices, and natural disturbances can affect the growth of trees and forest development. To make predictions about anthropogenic impacts on forests, we need to understand how these factors affect tree growth. In this study the effect of disturbance history on tree growth and stand structure was examined by using a computer model of forest development. The model was run under the climatic conditions of east Tennessee, USA, and the results compared to stand structure and tree growth data from a yellow poplar-white oak forest. Basal area growth and forest biomass were more accurately projected when rough approximations of the thinning and fire history typical of the measured plots were included in the simulation model. Stand history can influence tree growth rates and forest structure and should be included in any attempt to assess forest impacts.     

CFD-based assessment and visualization of the failure consequences of LPG tankers

This paper presents a simulation analysis of the explosions following an LPG leak and visualizes the consequences of the accident to reduce the consequences of the LPG leak explosion. Firstly, this paper proposes a CFD numerical simulation-based method for visualizing the consequences of LPG tanker failure. The method combines satellite maps and CFD numerical simulation data to visualize the consequences of LPG leaks and explosions, taking into account the influence of obstacles on the danger range of leaks and explosions; Secondly, this paper applies the method to a liquefied petroleum gas accident that occurred in the Wenling section of the Shenhai Expressway and performs CFD numerical simulation on the accident process and visualizes the consequences of the accident. Therefore, this method can provide a theoretical reference for the prior prevention of LPG accidents and the analysis of the consequences of accidents, as well as certain practical guidance instructive.     

An integrated fuzzy-stochastic modeling approach for risk assessment of groundwater contamination

An integrated fuzzy-stochastic risk assessment (IFSRA) approach was developed in this study to systematically quantify both probabilistic and fuzzy uncertainties associated with site conditions, environmental guidelines, and health impact criteria. The contaminant concentrations in groundwater predicted from a numerical model were associated with probabilistic uncertainties due to the randomness in modeling input parameters, while the consequences of contaminant concentrations violating relevant environmental quality guidelines and health evaluation criteria were linked with fuzzy uncertainties. The contaminant of interest in this study was xylene. The environmental quality guideline was divided into three different strictness categories: "loose", "medium" and "strict". The environmental-guideline-based risk (ER) and health risk (HR) due to xylene ingestion were systematically examined to obtain the general risk levels through a fuzzy rule base. The ER and HR risk levels were divided into five categories of "low", "low-to-medium", "medium", "medium-to-high" and "high", respectively. The general risk levels included six categories ranging from "low" to "very high". The fuzzy membership functions of the related fuzzy events and the fuzzy rule base were established based on a questionnaire survey. Thus the IFSRA integrated fuzzy logic, expert involvement, and stochastic simulation within a general framework. The robustness of the modeling processes was enhanced through the effective reflection of the two types of uncertainties as compared with the conventional risk assessment approaches. The developed IFSRA was applied to a petroleum-contaminated groundwater system in western Canada. Three scenarios with different environmental quality guidelines were analyzed, and reasonable results were obtained. The risk assessment approach developed in this study offers a unique tool for systematically quantifying various uncertainties in contaminated site management, and it also provides more realistic support for remediation-related decisions.     

Development of a model to simulate soil heavy metals lateral migration quantity based on SWAT in Huanjiang watershed, China

Lateral transportation of soil heavy metals in rainfall events could significantly increase the scope of pollution. Therefore, it is necessary to develop a model with high accuracy to simulate the migration quantity of heavy metals. A model for heavy metal migration simulation was developed based on the SWAT(Soil and Water Assessment Tool) model. This model took into consideration the influence of soil p H value, soil particle size, runoff volume, sediment amount,concentration of water-soluble heavy metals dissolved in runoff and insoluble absorbed to the soil particles. This model was reasonable in Huanjiang watershed, Guangxi Zhuang Autonomous Region, south China, covering an area of 273 km~2. The optimal drainage area threshold was determined by analyzing the effects of watershed subdivision on the simulation results to ensure the simulation accuracy. The main conclusions of this paper were:(1) watershed subdivision could affect simulation migration quantity of heavy metals;(2) the quantity of heavy metals transported by sediment accounted for 97%–99% of the total migration quantity in the study watershed. Therefore, sediment played the most important role in heavy metal migration;(3) the optimal drainage area threshold percentage to ensure high simulation accuracy was determined to be 2.01% of the total watershed;(4) with the optimal threshold percentage, this model could simulate the migration quantity of As, Pb and Cd accurately at the total watershed and subwatershed level. The results of this paper were useful for identifying the key regions with heavy metal migration.     

Case study: Assessment on large scale LPG BLEVEs in the 2011 Tohoku earthquakes

After the 2011 Tohoku earthquakes, several chemical and oil complexes on the Pacific Ocean shoreline of northeast Japan experienced massive losses. In Chiba, a refinery operated by Cosmo Oil lost 17 LPG storage vessels which were either heavily damaged or totally destroyed by fires and explosions in the refinery. These large vessels ranged in size from 1000 to 5000 m³. The estimated volume of LPG at the time of the incident was between 400 and 5000 m³ for each vessel. Five boiling liquid expanding vapor explosions (BLEVEs) of LPG occurred, resulting in huge fire balls measuring about 500 m in diameter.A BLEVE is defined as the explosive release of expanding vapor and boiling liquid when a container holding a pressure-liquefied gas fails catastrophically. It is thus important to estimate the physical properties of superheated liquids: the thermodynamic and transport properties, the intrinsic limits to superheating and depressurization, and the nature of thermodynamic paths. Also it is hoped to provide better understanding of the vessels designed, manufactured, installed, and operated to reduce or eliminate the probability that a sequence of events will result in BLEVE or loss of primary containment. Knowledge of these matters is still incomplete. The objective of this research is to estimate the significant BLEVE phenomenon in very large scale spherical vessels based on published information in Japan. There are some models predicting BLEVEs. However, it is essential to know if this is true for very large scales such as spheres since validation is usually rare to provide confidence in estimating the superheated liquids behaviors. To this end, comparing with the information on this event, the conditions in the five LPG vessels at the time of the BLEVE were determined in terms of: duration of vessel failure (time to BLEVE); mass fraction in the vessel with time; temperature distribution in the liquid and vapor region and pressure within the vessel (e.g. initial pressure and internal high-speed transient pressure during failure), by means of a computer program AFFTAC Analysis of Fire Effects on Tank Cars, which solves heat conduction, stress and a failure model of the tank, a thermodynamic model of its fluid contents, and a flow model for the lading flowing through the safety relief device. Subsequently, the consequences from the sphere BLEVE, such as the expected fireball diameter and duration and the expected blast overpressure produced by the BLEVE failures, are also subjects of active research. Here the blast using the methods of PHAST and SFPE Handbook of Fire Protection Engineering was calculated.Results suggest that methodologies here used gave reasonable estimations for such real and huge BLEVEs in a validated way, which may provide valuable guidance for risk mitigation strategy with regard to LPG facility in design, emergency planning, resiliency, operations, and risk management.     

A simulation model of rhizome networks for Fallopia japonica (Japanese knotweed) in the United Kingdom

Fallopia japonica (Japanese knotweed) is an aggressively invasive herbaceous perennial that causes substantial economic and environmental damage in the United Kingdom (UK). As such, it is of considerable concern to councils, environmental groups, private landowners and property developers. We construct a 3D correlated random walk model of the development of the subterranean rhizome network for a single stand of F. japonica. The formulation of this model uses detailed knowledge of the morphology and physiology of the plant, both of which differ in the UK to that of its native habitat due to factors including a lack of predation and competition, longer growth seasons and favourable environmental conditions in the UK. Field data obtained as a part of this study are discussed and used in the model for parameterisation and validation. The simulation captures the field data well and predicts, for example, quadratic growth in time for the stand area. Furthermore, the role of a selection of parameters on long-term stand development are discussed, highlighting some key factors affecting vegetative spread rates.     

Dynamics of soil organic matter in primary and secondary forest succession on sandy soils in The Netherlands: An application of the ROMUL model

We applied the simulation model ROMUL of soil organic matter dynamics in order to analyse and predict forest soil organic matter (SOM) changes following stand growth and also to identify gaps of data and modelling problems. SOM build-up was analysed (a) from bare sand to forest soil during a primary succession in Scots pine forest and (b) on mature forest soil under Douglas fir plantations as an example of secondary succession in The Netherlands. As some of the experimental data were unreliable we compiled a set of various scenarios with different soil moisture regime, initial SOM pools and amount and quality of above and below ground litter input. This allowed us to find the scenarios that reflect the SOM dynamics more realistically. In the Scots pine forest, total litter input was estimated as 0.50 kg m⁻² year⁻¹. Two scenarios were defined for the test runs: (a) forest floor moisture regimes—‘dry, mesic and hydric’ and (b) augmenting a root litter pool with three ratios of needles and branches to roots: 1:1, 1:1.5 and 1:2.0. The scenario finally compiled had the following characteristics: (a) climate for dry site with summer drought and high winter moisture of forest floor; (b) a litter input of 0.25 kg m⁻² year⁻¹ above ground and 0.50 kg m⁻² year⁻¹ below ground; (c) a low nitrogen and ash content in all litter fall fractions. The test runs for the estimation of the initial SOM pools and the amount and proportion of above and below ground litter fall were also performed in the Douglas fir plantation. The inputs of above ground litter tested in various combinations were 0.30 and 0.60 kg m⁻² year⁻¹, and below ground litter 0.30, 0.60 and 0.90 kg m⁻² year⁻¹. The scenario that fitted the experimental data had an SOM pool of 20–25 kg m⁻², an aboveground litter input of 0.6 kg m⁻² year⁻¹and a below ground litter input of 0.9 kg m⁻² year⁻¹. The long-term simulation corresponded well with the observed patterns of soil organic matter accumulation associated with the forest soil development in primary and secondary succession. During primary succession in Scots pine forest on dry sand there is a consistent accumulation of a raw humus forest floor. The soil dynamics in the Douglas fir plantation also coincide with the observed patterns of SOM changes during the secondary succession, with SOM decreasing significantly under young forest, and SOM being restored in the older stands.     

板料拉延成形中应力状态的研究

在采用Hill关于各向异性材料的修正屈服函数,并运用板块分析方法对板料轴对称成形进行数值模拟的基础上,着重分析了拉延成形中悬空区域材料的应力状态,提出了对该区域应力状态过渡点即周向应力为零的应力分界圆进行数值求解的公式及方法,同时对影响该分界圆位置的诸因素进行了细致分析.     

A simplified approach to modelling underground migration of radionuclides from contaminated river sediments

Accidental releases of waste water from the first Czechoslovak nuclear power plant, A1, caused contamination of sediments of the Dudváh river, flowing into the Vah river, in Slovakia. Rather high concentrations of ¹³⁷Cs and ⁹⁰Sr (2150Bq dm⁻³ and 215Bq dm⁻³, respectively) were found in bottom sediments of a former channel of the re-engineered river body at a distance of about 250 m from a village, Siladice, with water-supply wells. In order to assess the possibility of contamination of the wells, underground migration of both radionuclides from the contaminated area was simulated using an original layered convection-diffusion model. K_d values determined in laboratory experiments were used. The analysis of the hydrological situation in the area reveals that the critical condition is a dominant horizontal groundwater flow near the water table in the direction from the Váh bank to Siladice, in the periods when the contaminated body lies under the water table. The simulation calculated under conservative conditions showed that the contamination of water-supply wells would not exceed permissible concentration limits.