Search NASA⌕ Search

SEARCH · Search NASA

Results for “process simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Evaluating Energy Savings Measures for Building Retrofits in an Interactive Workshop Environment

As part of the Better Buildings Summit 2025, a workshop was developed to simulate the decision process for implementing a set of energy savings measures on an existing building. As part of the workshop a "game" was created around energy savings strategies, and workshop participants selected renovation cards that provided data around the impacts of that renovation choice. Participants then were able to calculate overall savings and payback from the selected retrofit opportunities. This publication has been produced to package the resources used in the activity in an organized document so that the materials can be distributed and reused more broadly. It provides a comprehensive framework for replicating the workshop, including step-by-step instructions, retrofit "game cards" with data points, and facilitation guidance. By using these materials, readers can conduct their own version of this interactive workshop with stakeholders in their networks, such as building owners, energy managers, or policy makers. The goal is to empower users to simulate the decision-making processes for energy projects, foster collaboration, and drive meaningful discussions around cost-effective retrofit opportunities and goal-oriented outcomes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

The LSST DESC DC2 Simulated Sky Survey

Here, we describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep area of approximately 300 deg 2 , as well as a deep drilling field of approximately 1 deg 2 . We simulate 5 yr of the planned 10 yr survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the data set to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic test bed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time domain cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A GPU Accelerated Mixed‐Precision Finite Difference Informed Random Walker (FDiRW) Solver for Strongly Inhomogeneous Diffusion Problems

In nature, many complex multi‐physics coupling problems exhibit significant diffusivity inhomogeneity, where one process occurs several orders of magnitude faster than others temporally. Simulating rapid diffusion alongside slower processes demands intensive computational resources due to the necessity for small time steps. To address these computational challenges, we have developed an efficient numerical solver named Finite Difference informed Random Walker (FDiRW). In this study, we propose a GPU‐accelerated, mixed‐precision configuration for the FDiRW solver to maximize efficiency through GPU multi‐threaded parallel computation and lower precision computation. Numerical evaluation results reveal that the proposed GPU‐accelerated mixed‐precision FDiRW solver can achieve a 117× speedup over the CPU baseline, while an additional 1.75× speedup is achieved by employing lower precision GPU computation. Notably, for large model sizes, the GPU‐accelerated mixed‐precision FDiRW solver demonstrates strong scaling with the number of nodes used in simulation. When simulating radionuclide absorption processes by porous wasteform particles with a medium‐sized model of 192 × 192 × 192, this approach reduces the total computational time to 10 min, enabling the simulation of larger systems with strongly inhomogeneous diffusivity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Report on High-Fidelity Dynamic Modeling of a Coal-Fired Steam Power Plant

As a result of the growth of renewables including solar and wind energy with fluctuating production, fossil fuel power plants are being required to cycle between high and low power production. This cycling is both at a greater frequency and over a wider range than in the past. In many cases, power plants are not designed for this type of cycling operation but can nonetheless endure these challenging operation requirements under the right conditions. In this project, optimal solution for enhanced flexible operations are being investigated using model based estimation and control techniques. To support the development of model-based estimator and model-based controls at GE Global Research, GE Steam Power configured a dynamic model using a reference steam plant design including the boiler, turbine, and water/steam conditioning systems as well as the controls needed for plant cycling with stability and reliability. The dynamic model was built using the APROS® software from VTT, and then calibrated to multiple load conditions from full load (100%TMCR) to partial loads (75%TMCR, 50%TMCR and 25%TMCR) based on internally developed steady state heat balance models at the unit level. These internal heat balance models are based on first principles and extensive engineering experiences from GE Steam Power as an OEM and a services provider. This topical report presents the structure of the unit level dynamic model, the tuning process, and representative simulation results from typical load cycling simulations using the dynamic model.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Design-Technology Co-Optimization for NVM-based Neuromorphic Processing Elements

An emerging use-case of machine learning (ML) is to train a model on a high-performance system and deploy the trained model on energy-constrained embedded systems. Neuromorphic hardware platforms, which operate on principles of the biological brain, can significantly lower the energy overhead of a machine learning inference task, making these platforms an attractive solution for embedded ML systems. In this paper, we present a design-technology tradeoff analysis to implement such inference tasks on the processing elements (PEs) of a Non-Volatile Memory (NVM)-based neuromorphic hardware. Through detailed circuit-level simulations at scaled process technology nodes, we show the negative impact of technology scaling on the information-processing latency, which impacts the quality-of-service (QoS) of an embedded ML system. At a finer granularity, the latency inside a PE depends on 1) the delay introduced by parasitic components on its current paths, and 2) the varying delay to sense different resistance states of its NVM cells. Based on these two observations, we make the following three contributions. First, on the technology front, we propose an optimization scheme where the NVM resistance state that takes the longest time to sense is set on current paths having the least delay, and vice versa, reducing the average PE latency, which improves the QoS. Second, on the architecture front, we introduce isolation transistors within each PE to partition it into regions that can be individually power-gated, reducing both latency and energy. Finally, on the system-software front, we propose a mechanism to leverage the proposed technological and architectural enhancements when implementing a machine-learning inference task on neuromorphic PEs of the hardware. Evaluations with a recent neuromorphic hardware architecture show that our proposed design-technology co-optimization approach improves both performance and energy efficiency of machine-learning inference tasks without incurring high cost-per-bit.

42 ENGINEERING↗
DOE OSTI2018

Advances in Geologic Disposal Safety Assessment and an Unsaturated Alluvium Reference Case

The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Depat ment of Energy (DOE) Office of Nuclear Energy (NE), Office of Fuel Cycle Technology (OFCT) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high level nuclear waste (HLW). Two high priorities for SFWST disposal R&D are design concept development and disposal system modeling (DOE 2011, Table 6). These priorities are directly addressed in the SFWST Geologic Disposal Safety Assessment (GDSA) work package, which is charged with developing a disposal system modeling and analysis capability for evaluating disposal system performance for nuclear waste in geologic media. This report describes specific GDSA activities in fiscal year 2018 (FY 2018) toward the development of GDSA Framework, an enhanced disposal system modeling and analysis capability for geologic disposal of nuclear waste. GDSA Framework employs the PFLOTRAN thermal-hydrologic-chemical multiphysics code (Hammond et al. 2011a; Lichtner and Hammond 2012) and the Dakota uncertainty sampling and propagation code (Adams et al. 2012; Adams et al. 2013). Each code is designed for massivelyparallel processing in a high-performance computing (HPC) environment. Multi-physics representations in PFLOTRAN are used to simulate various coupled processes including heat flow, fluid flow, waste dissolution, radionuclide release, radionuclide decay and ingrowth, precipitation and dissolution of secondary phases, and radionuclide transport through engineered barriers and natural geologic barriers to the biosphere. Dakota is used to generate sets of representative realizations and to analyze parameter sensitivity.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Stand-Alone Dynamic System Simulation of a Fissile Solution System

Dynamic system simulations (DSS) of fissile solution systems (FSS) have been performed at Los Alamos National Laboratory (LANL) for several years. DSS have been designed and executed using the DESIRE (Direct Executing Simulation in Real Time) system, a platform ideal for expert system modelers, as it allows for direct mathematical specification of a model, and facilities for performing experiments and analyses with the model. Models built in DESIRE have been compared to historical data from several historic aqueous homogeneous reactors (AHR) showing close agreement between experimental data and model estimates. Historical systems modeled include: SUPO (Super Power), an AHR that operated at Los Alamos National Laboratory from 1951 to 1974; SILENE, which was operated by CEA at Valduc, France until 2011; KEWB (Kinetics Experiment Water Boiler) A-2 and B-5 cores, operated by North American Atomics during the 1960’s and 1970’s. SUPO is considered a benchmark for steady-state operations, SILENE for pulse, and KEWB as confirmatory for both modes. However, there are some simulation tasks that are not well-suited to DESIRE model implementation. These tasks include engineering design and analysis and providing a hardware simulator backend. It has been determined that executable simulation models would better suit these purposes, hence a process for converting DESIRE simulation models to executable simulation models was designed and implemented, and several models have been converted. DESIRE simulation models are converted to C++ executable models by an automated process. The conversion process is quick, accurate and repeatable, giving confidence that converted models will perform as well as the original DESIRE models. Extensive testing on the SUPO and accelerator driven subcritical system models has shown that the executable models perform almost identically to the DESIRE models. All major trends and events are simulated with excellent agreement and the steady-state results typically agree within a small fraction of a percent (always better than 0.1% and generally better than 0.01% agreement).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Simulation-driven design optimization of reaction injection molding (RIM) process for polydicyclopentadiene (pDCPD): Minimizing cycle time, defects, and warpage

Replacing metal components in trucks, trailers, and buses with lightweight polymer composites is challenging due to high temperatures and complex manufacturing. The Reaction Injection Molding (RIM) process using Dicyclopentadiene (DCPD) resin offers a solution by producing robust parts with excellent stiffness, impact strength, and resistance properties. Simulations are essential for optimizing this process, predicting defects, and improving quality. However, most commercial software is tailored for thermoplastics, requiring thermoset users to generate their own datasets. In this study, a material data card for DCPD was developed to perform RIM simulations. Design of Experiments (DOE) was used to identify key factors affecting filling, curing, and warpage, aiming to minimize cycle time and defects. The simulations explored varying injection gate parameters (size, location, number) and process conditions (mold/resin temperature, injection/curing pressure). Results showed that gate design significantly impacts filling behavior and defects. A single central gate provided balanced flow with fewer defects, while two corner gates led to more defects. Additionally, lower injection pressure increased filling time, while higher mold temperature accelerated curing but led to more warpage. In conclusion, this optimization framework aims to enhance DCPD part performance and promote sustainable manufacturing by reducing waste and energy consumption.

42 ENGINEERING↗
DOE OSTI2017

Biogenic sulfide control by nitrate and (per)chlorate – A monitoring and modeling investigation

Biosouring is commonly encountered during secondary oil recovery when seawater or another high sulfate water source is utilized for flooding; as a result, effective souring control is of great interest to the oil industry. Here we describe a laboratory study to evaluate the relative effectiveness of souring interventions through the injection of nitrate, chlorate and perchlorate, collectively (per)chlorate, solutions and whether in-situ galvanic potential measurements can be used for convenient and quantitative tracking of sulfide dynamics. Nitrate has typically been the chemical of choice for souring treatments, while the efficacy of (per)chlorate as a new candidate inhibitor has only been explored recently. (Per)chlorate is known to inhibit oil reservoir souring via mechanisms such as toxicity, bio-competitive exclusion and sulfur redox cycling. Two sets of experiments under different matrix and inoculation conditions were conducted to evaluate treatment efficiency under variable baseline physical and biogeochemical conditions. Our data demonstrated the sensitivity of the galvanic potential signals to sulfide concentrations where the sulfide-galvanic potential correlation is similar to the theoretical predictions based on the Nernst equation, demonstrating the feasibility of using galvanic potential as a quick and economical method for quantifying in situ sulfide concentrations for tracking reservoir souring processes and subsequent intervention effectiveness. Our results show that all three chemicals were effective at suppressing sulfidogenesis. A reactive transport model for perchlorate treatment was developed to simulate the reaction processes and explore the interactions between the underlying competing mechanisms of this inhibitor. A baseline simulation captured the temporal patterns of the effluent chemical species. Subsequent simulations in which individual inhibition mechanisms were systematically removed elucidated the relative role that each inhibition mechanism played at the different phases of the experiment. The simulation results complement the experimental findings. Our study supports the potential advantages of souring control with (per)chlorate treatments, and the application of galvanic signal as an economic, in-situ monitoring approach for tracking souring dynamics and treatment efficacy.

02 PETROLEUM↗
DOE OSTI2015

In-situ tube burst testing and high-temperature deformation behavior of candidate materials for accident tolerant fuel cladding

The high resistance of cladding to plastic deformation and burst failure is one of the most essential properties of accident tolerant fuel (ATF) for maintaining structural integrity during a loss-of-coolant accident (LOCA) since the deformation and burst behavior governs the cooling efficiency of flow channels and process of fission product release. To simulate and evaluate such deformation and burst process of thin-walled cladding, an in-situ testing and evaluation method has been developed on the basis of visual imaging and image analysis techniques. The method uses a specialized optics system consisted of a high-resolution video camera, light filtering unit, and monochromatic light sources, and the in-situ testing is performed using a 50 mm long pressurized thin-walled tubular specimen set in a programmable furnace. In this study eleven (11) candidate cladding materials for ATF, i.e., 6 FeCrAl alloys and 5 nanostructured steels, were tested using the newly developed method, and the time-dependent images were analyzed to produce detailed deformation and burst data such as true hoop stress, strain (creep) rate, and failure stress. Relatively soft FeCrAl alloys deformed and burst below 800°C while negligible strain rates were measured for higher strength alloys and/or for relatively thick wall specimens.

Gussev, Maxim N.↗
DOE OSTI2022

Avoiding excess computation in asynchronous evolutionary algorithms

Abstract Asynchronous evolutionary algorithms are becoming increasingly popular as a means of making full use of many processors while solving computationally expensive search and optimization problems. These algorithms excel at keeping large clusters fully utilized, but may sometimes inefficiently sample an excess of fast‐evaluating solutions at the expense of higher‐quality, slow‐evaluating ones. We have previously introduced a steady‐state parent selection strategy, SWEET (“Selection whilE EvaluaTing”), that sometimes selects individuals that are still being evaluated and allows them to reproduce early. We perform a takeover‐time analysis that confirms that this strategy gives slow‐evaluating individuals that have higher fitnesses an increased ability to multiply in the population. We also find that SWEET appears effective at improving optimization performance on problems in which solution quality is positively correlated with evaluation time. We evaluate our approach on six simulated real‐valued optimization problems and three real‐world applications: an autonomous vehicle controller problem that involves tuning a spiking neural network and two adversarial EA problems. We further evaluate SWEET versus a basic asynchronous process in a simulated setting. We present evidence that SWEET outperforms basic asynchronous processes in a use‐case in which performance is positively correlated with evaluation time, and performs comparably (and often better) than basic asynchronous processes in several use‐cases where performance is negatively correlated with evaluation time. That said, in the cases where performance and evaluation time are negatively correlated the variance of outcomes for SWEET is notably high.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Spectra of clinical CT scanners using a portable Compton spectrometer

Purpose: Spectral information of the output of x-ray tubes in (dual source) computer tomography (CT) scanners can be used to improve the conversion of CT numbers to proton stopping power and can be used to advantage in CT scanner quality assurance. The purpose of this study is to design, validate, and apply a compact portable Compton spectrometer that was constructed to accurately measure x-ray spectra of CT scanners. Methods: In the design of the Compton spectrometer, the shielding materials were carefully chosen and positioned to reduce background by x-ray fluorescence from the materials used. The spectrum of Compton scattered x-rays alters from the original source spectrum due to various physical processes. Reconstruction of the original x-ray spectrum from the Compton scattered spectrum is based on Monte Carlo simulations of the processes involved. This reconstruction is validated by comparing directly and indirectly measured spectra of a mobile x-ray tube. The Compton spectrometer is assessed in a clinical setting by measuring x-ray spectra at various tube voltages of three different medical CT scanner x-ray tubes. Results: The directly and indirectly measured spectra are in good agreement (their ratio being 0.99) thereby validating the reconstruction method. The measured spectra of the medical CT scanners are consistent with theoretical spectra and spectra obtained from the x-ray tube manufacturer. Conclusions: A Compton spectrometer has been successfully designed, constructed, validated, and applied in the measurement of x-ray spectra of CT scanners. These measurements show that our compact Compton spectrometer can be rapidly set-up using the alignment lasers of the CT scanner, thereby enabling its use in commissioning, troubleshooting, and, e.g., annual performance check-ups of CT scanners.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Long-term measurements in a mixed-grass prairie reveal a change in soil organic carbon recalcitrance and its environmental sensitivity under warming

Soil respiration, the major pathway for ecosystem carbon (C) loss, has the potential to enter a positive feedback loop with the atmospheric CO2 due to climate warming. For reliable projections of climate-carbon feedbacks, accurate quantification of soil respiration and identification of mechanisms that control its variability are essential. Process-based models simulate soil respiration as functions of belowground C input, organic matter quality, and sensitivity to environmental conditions. However, evaluation and calibration of process-based models against the long-term in situ measurements are rare. Here, we evaluate the performance of the Terrestrial ECOsystem (TECO) model in simulating total and heterotrophic soil respiration measured during a 16-year warming experiment in a mixed-grass prairie; calibrate model parameters against these and other measurements collected during the experiment; and explore whether the mechanisms of C dynamics have changed over the years. Calibrating model parameters against observations of individual years substantially improved model performance in comparison to pre-calibration simulations, explaining 79–86% of variability in observed soil respiration. Interannual variation of the calibrated model parameters indicated increasing recalcitrance of soil C and changing environmental sensitivity of microbes. Overall, we found that (1) soil organic C became more recalcitrant in intact soil compared to root-free soil; (2) warming offset the effects of increasing C recalcitrance in intact soil and changed microbial sensitivity to moisture conditions. Furthermore, these findings indicate that soil respiration may decrease in the future due to C quality, but this decrease may be offset by warming-induced changes in C cycling mechanisms and their responses to moisture conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Benchmark Problems of the Geothermal Technologies Office Code Comparison Study

A diverse suite of numerical simulators is currently being applied to predict or understand the performance of enhanced geothermal systems (EGS). To build confidence and identify critical development needs for these analytical tools, the United States Department of Energy, Geothermal Technologies Office has sponsored a Code Comparison Study (GTO-CCS), with participants from universities, industry, and national laboratories. A principal objective for the study was to create a community forum for improvement and verification of numerical simulators for EGS modeling. Teams participating in the study were those representing U.S. national laboratories, universities, and industries, and each team brought unique numerical simulation capabilities to bear on the problems. Two classes of problems were developed during the study, benchmark problems and challenge problems. The benchmark problems were structured to test the ability of the collection of numerical simulators to solve various combinations of coupled thermal, hydrologic, geomechanical, and geochemical processes. This class of problems was strictly defined in terms of properties, driving forces, initial conditions, and boundary conditions. Study participants submitted solutions to problems for which their simulation tools were deemed capable or nearly capable. Some participating codes were originally developed for EGS applications whereas some others were designed for different applications but can simulate processes similar to those in EGS. Solution submissions from both were encouraged. In some cases, participants made small incremental changes to their numerical simulation codes to address specific elements of the problem, and in other cases participants submitted solutions with existing simulation tools, acknowledging the limitations of the code. The challenge problems were based on the enhanced geothermal systems research conducted at Fenton Hill, near Los Alamos, New Mexico, between 1974 and 1995. The problems involved two phases of research, stimulation, development, and circulation in two separate reservoirs. The challenge problems had specific questions to be answered via numerical simulation in three topical areas: 1) reservoir creation/stimulation, 2) reactive and passive transport, and 3) thermal recovery. Whereas the benchmark class of problems were designed to test capabilities for modeling coupled processes under strictly specified conditions, the stated objective for the challenge class of problems was to demonstrate what new understanding of the Fenton Hill experiments could be realized via the application of modern numerical simulation tools by recognized expert practitioners.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Long-time integration of parametric evolution equations with physics-informed DeepONets

Ordinary and partial differential equations (ODEs/PDEs) play a paramount role in analyzing and simulating complex dynamic processes across all corners of science and engineering. In recent years machine learning tools are aspiring to introduce new effective ways of simulating such equations, however existing approaches are not able to reliably return stable and accurate predictions across long temporal horizons. We aim to address this challenge by introducing an effective framework for learning evolution operators that map random initial conditions to associated ODE/PDE solutions within a short time interval. Such operators can be parametrized by deep neural networks that are trained in an entirely self-supervised manner without requiring one to generate any paired input-output observations. Global long-time predictions across a range of initial conditions can be then obtained by iteratively evaluating the trained model using each prediction as the initial condition for the next evaluation step. Here, this introduces a new approach to temporal domain decomposition that is shown to be effective in performing accurate long-time simulations for a wide range of parametric ODE and PDE systems, from wave propagation, to reaction-diffusion dynamics and stiff chemical kinetics, introducing a new way of rapidly emulating non-equilibrium processes in science and engineering.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Predicting the structural basis of targeted protein degradation by integrating molecular dynamics simulations with structural mass spectrometry

Targeted protein degradation (TPD) is a promising approach in drug discovery for degrading proteins implicated in diseases. A key step in this process is the formation of a ternary complex where a heterobifunctional molecule induces proximity of an E3 ligase to a protein of interest (POI), thus facilitating ubiquitin transfer to the POI. In this work, we characterize 3 steps in the TPD process. (1) We simulate the ternary complex formation of SMARCA2 bromodomain and VHL E3 ligase by combining hydrogen-deuterium exchange mass spectrometry with weighted ensemble molecular dynamics (MD). (2) We characterize the conformational heterogeneity of the ternary complex using Hamiltonian replica exchange simulations and small-angle X-ray scattering. (3) We assess the ubiquitination of the POI in the context of the full Cullin-RING Ligase, confirming experimental ubiquitinomics results. Differences in degradation efficiency can be explained by the proximity of lysine residues on the POI relative to ubiquitin.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Modeling the Capacitive Deionization Process in Dual-Porosity Electrodes

In many areas of the world, there is a need to increase water availability. Capacitive deionization (CDI) is an electrochemical water treatment process that can be a viable alternative for treating water and for saving energy. A model is presented to simulate the CDI process in heterogeneous porous media comprising two different pore sizes. It is based on a theory for capacitive charging by ideally polarizable porous electrodes without Faradaic reactions or specific adsorption of ions. A two steps volume averaging technique is used to derive the averaged transport equations in the limit of thin electrical double layers. A one-equation model based on the principle of local equilibrium is derived. The constraints determining the range of application of the one-equation model are presented. The effective transport parameters for isotropic porous media are calculated solving the corresponding closure problems. The source terms that appear in the average equations are calculated using theoretical derivations. The global diffusivity is calculated by solving the closure problem.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

When managers forage for pests: Implementing the functional response in pest management

In this work, we explore how the functional response framework can be implemented in pest management. Here, managers take the role of predators foraging on pests and facing monetary costs for survey and control in a spatial domain where the pest distribution and control strategy do not have to be random. To investigate this framework quantitatively, we simulated various management processes on different pest spatial distributions using a spatially-explicit individual-based model and Monte-Carlo simulations, and also confirmed some of the results analytically. By graphing the number of pests controlled versus pest density, we obtained management functional response curves. Whether the management functional response was shaped like a type I, type II or type III functional response depended on the management costs and the search area. Yet, the management spatial strategy and the pest spatial distribution had little effect on the functional response. We applied our model to the management of mountain pine beetle epidemic in Cypress Hills, Saskatchewan, Canada, with simulations matching the real number of attacked trees controlled by managers. We showed how to make an analogy between functional responses in predator–prey interactions and in human–pest interactions and thereby, apply insights from the functional response framework to pest management.

54 ENVIRONMENTAL SCIENCES↗