Search NASA⌕ Search

SEARCH · Search NASA

Results for “temporal continuity”

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 19 records

Exact enforcement of temporal continuity in sequential physics-informed neural networks

The use of deep learning methods in scientific computing represents a potential paradigm shift in engineering problem solving. One of the most prominent developments is Physics-Informed Neural Networks (PINNs), in which neural networks are trained to satisfy partial differential equations (PDEs). While this method shows promise, the standard version has been shown to struggle in accurately predicting the dynamic behavior of time-dependent problems. To address this challenge, methods have been proposed that decompose the time domain into multiple segments, employing a distinct neural network in each segment and directly incorporating continuity between them in the loss function of the minimization problem. In this work we introduce a method to exactly enforce continuity between successive time segments via a solution ansatz. This hard constrained sequential PINN (HCS-PINN) method is simple to implement and eliminates the need for any loss terms associated with temporal continuity. The method is tested for a number of benchmark problems involving both linear and non-linear PDEs. Examples include various first order time dependent problems in which traditional PINNs struggle, namely advection, Allen–Cahn, and Korteweg–de Vries equations. Furthermore, second and third order time-dependent problems are demonstrated via wave and Jerky dynamics examples, respectively. Notably, the Jerky dynamics problem is chaotic, making the problem especially sensitive to temporal accuracy. Finally, the numerical experiments conducted with the proposed method demonstrated superior convergence and accuracy over both traditional PINNs and the soft-constrained counterparts.

42 ENGINEERING↗

A 3-Channel Algorithm for Retrieving Spatially and Temporally Continuous Cloud Properties Across Different Geostationary Satellite Imagers

Cloud property retrieval algorithms for passive satellite imagers are generally designed to take advantage of all the useful spectral information available for a particular satellite. This strategy optimizes accuracy and reduces misidentification and retrieval biases, particularly for modern satellites with many spectral channels. However, the application of dissimilar algorithms tailored for different satellite sensor scan present a problem within the climate data record (CDR). Algorithm inconsistencies can introduce artificial trends in the CDR that are tied to instrument changes rather than physical changes, especially when older satellites with limited spectral information are included. The NASA CERES (Clouds and the Earth’s Radiant Energy System) data record provides global cloud property retrievals across 23 years and more than 25 satellites. With the goal of producing a spatially and temporally continuous record of cloud properties, the CERES cloud working group has developed algorithms that use only 3 channels that are common to most geostationary satellite imagers: 0.65, 3.9, and 10.8 μm.

Sarah Bedka↗

Temporally continuous thermofluidic–thermomechanical modeling framework for metal additive manufacturing

Additive manufacturing (AM) is known to generate large magnitudes of residual stresses (RS) within builds due to steep and localized thermal gradients. In the current state of commercial AM technology, manufacturers generally perform heat treatments in effort to reduce the generated RS and its detrimental effects on part distortion and in-service failure. Computational models that effectively simulate the deposition process can provide valuable insights to improve RS distributions. Accordingly, it is common to employ Computational fluid dynamics (CFD) models or finite element (FE) models. While CFD can predict geometric and thermal-fluid behavior, it cannot predict the structural response (e.g., stress–strain) behavior. On the other hand, an FE model can predict mechanical behavior, but it lacks the ability to predict geometric and fluid behavior. Thus, an effectively integrated thermofluidic–thermomechanical modeling framework that exploits the benefits of both techniques while avoiding their respective limitations can offer valuable predictive capability for AM processes. In contrast to previously published efforts, the work herein describes a one-way coupled CFD-FEA framework that abandons major simplifying assumptions, such as geometric steady-state conditions, the absence of material plasticity, and the lack of detailed RS evolution/accumulation during deposition, as well as insufficient validation of results. Here, the presented framework is demonstrated for a directed energy deposition (DED) process, and experiments are performed to validate the predicted geometry and RS profile. Both single- and double-layer stainless steel 316L builds are considered. Geometric data is acquired via 3D optical surface scans and X-ray micro-computed tomography, and residual stress is measured using neutron diffraction (ND). Comparisons between the simulations and measurements reveal that the described CFD-FEA framework is effective in capturing the coupled thermomechanical and thermofluidic behaviors of the DED process. The methodology presented is extensible to other metal AM processes, including power bed fusion and wire-feed-based AM.

42 ENGINEERING↗

On the Relationship Between Sub-Daily Instantaneous and Daily Total Gross Primary Production: Implications for Interpreting Satellite-Based SIF Retrievals

Spatially and temporally continuous estimation of plant photosynthetic carbon fixation (or gross primary production, GPP) is crucial to our understanding of the global carbon cycle and the impact of climate change. Besides spatial, seasonal and interannual variations, GPP also exhibits strong diurnal variations. Satellite retrieved solar-induced chlorophyll fluorescence (SIF) provides a spatially continuous, but temporally discrete measurement of plant photosynthesis, and has the potential to be used to estimate GPP at global scale. However, it remains unclear whether the seasonal time series of SIF snapshots taken at a fixed time of the day can be used to infer daily total GPP variation at spatial and seasonal scales. In this study, we first used GPP estimates from 135 eddy covariance flux sites, covering a wide range of geographic locations and biome types, to investigate the relationship between the instantaneous GPP (GPP(sub inst)) and daily GPP (GPP(sub daily)) on the seasonal course for different times of the day. Latitudinal and diurnal patterns were found to correspond to variations in photosynthetically active radiation (PAR) and light use efficiency (LUE), respectively. We then used the Soil-Canopy Observation Photosynthesis and Energy Balance (SCOPE) model and the FluxCom GPP product to investigate the instantaneous and daily SIF-GPP relationships at five flux tower sites along a latitudinal gradient and at a global scale for different biome types. The results showed that daily SIF had a stronger linear correlation with daily GPP than instantaneous SIF at the seasonal scale, with an instantaneous to daily SIF conversion factor following the latitudinal and seasonal pattern driven by PAR. Our study highlights the necessity to take the latitudinal and diurnal factors into consideration for SIF-GPP relationship analyses or for physiological phenology analyses based on SIF.

photosynthesis model↗

Cluster Observations of Non-Time Continuous Magnetosonic Waves

Equatorial magnetosonic waves are normally observed as temporally continuous sets of emissions lasting from minutes to hours. Recent observations, however, have shown that this is not always the case. Using Cluster data, this study identifies two distinct forms of these non temporally continuous use missions. The first, referred to as rising tone emissions, are characterized by the systematic onset of wave activity at increasing proton gyroharmonic frequencies. Sets of harmonic emissions (emission elements)are observed to occur periodically in the region +/- 10 off the geomagnetic equator. The sweep rate of these emissions maximizes at the geomagnetic equator. In addition, the ellipticity and propagation direction also change systematically as Cluster crosses the geomagnetic equator. It is shown that the observed frequency sweep rate is unlikely to result from the sideband instability related to nonlinear trapping of suprathermal protons in the wave field. The second form of emissions is characterized by the simultaneous onset of activity across a range of harmonic frequencies. These waves are observed at irregular intervals. Their occurrence correlates with changes in the spacecraft potential, a measurement that is used as a proxy for electron density. Thus, these waves appear to be trapped within regions of localized enhancement of the electron density.

Walker, Simon N.↗

Temporal Resource Continuity Increases Predator Abundance in a Metapopulation Model: Insights for Conservation and Biocontrol

The amount of habitat in a landscape is an important metric for evaluating the effects of land cover on biodiversity, yet it fails to capture complex temporal dimensions of resource availability that could be consequential for species population dynamics. Here, we use a spatially-explicit predator–prey metapopulation model to test the effect of different spatiotemporal resource patterns on insect predators and their prey. We examined population responses in model landscapes that varied in both the amount and temporal variability of basal vegetation. Further, we examined cases where prey comprised either a single generalist species or two specialist species that use different resources available either early or late in the growing season. We found that predators and generalist prey benefitted from lower temporal variance of basal resources, which increased landscape-scale abundances. However, increasing the amount of basal resources also increased the variability of generalist prey populations. Specialist prey, on the other hand, did not benefit from less temporally variable basal resources, as they were restricted by habitat type, while also suffering greater predation. Predators achieved greater prey suppression in landscapes with less temporally variable resources, but the overall effects on prey abundance depended on prey habitat specialization. Our simulations demonstrate the joint importance of both the amount and temporal variability of resources for understanding how landscape heterogeneity influences biodiversity and ecosystem services such as the biological control of agricultural pests.

Spiesman, Brian↗

Brief communication: Monitoring snow depth using small, cheap, and easy-to-deploy snow–ground interface temperature sensors

Abstract. Temporally continuous snow depth estimates are vital for understanding changing snow patterns and impacts on permafrost in the Arctic. We trained a random forest machine learning model to predict snow depth from variability in snow–ground interface temperature. The model performed well on Alaska's Seward Peninsula where it was trained and at Arctic evaluation sites (RMSE ≤ 0.15 m). It performed poorly at temperate sites with deeper snowpacks, partially due to training data limitations. Small temperature sensors are cheap and easy to deploy, so this technique enables spatially distributed and temporally continuous snowpack monitoring at high latitudes to an extent previously infeasible.

54 ENVIRONMENTAL SCIENCES↗

Code Description for "Brief Communication: Monitoring snow depth using small, cheap, and easy-to-deploy ground surface temperature sensors"

Temporally continuous snow depth estimates are vital for understanding changing snow patterns and impacts on permafrost in the Arctic. We train a random forest machine learning model to predict snow depth from variability in ground surface temperature. To our knowledge, this is the first time that small ground surface temperature sensors have been used to estimate snow depth. The model performs well at sites where the model was trained and at pan-arctic evaluation sites (RMSE <= 0.15 m). Small temperature sensors are cheap and easy-to-deploy, so this technique enables spatially distributed and temporally continuous snowpack monitoring to an extent previously infeasible. The model is flexible and can be applied to datasets retroactively to retrieve snow depth estimates at additional sites. This code package includes a *.joblib file of the trained random forest model and a *.ipynb file showing how to clean input data, train the random forest model, and apply the model.

Bachand, Claire↗

Potential for Monitoring Snow Cover in Boreal Forests by Combining MODIS Snow Cover and AMSR-E SWE Maps

Monitoring of snow cover extent and snow water equivalent (SWE) in boreal forests is important for determining the amount of potential runoff and beginning date of snowmelt. The great expanse of the boreal forest necessitates the use of satellite measurements to monitor snow cover. Snow cover in the boreal forest can be mapped with either the Moderate Resolution Imaging Spectroradiometer (MODIS) or the Advanced Microwave Scanning Radiometer for EOS (AMSR-E) microwave instrument. The extent of snow cover is estimated from the MODIS data and SWE is estimated from the AMSR-E. Environmental limitations affect both sensors in different ways to limit their ability to detect snow in some situations. Forest density, snow wetness, and snow depth are factors that limit the effectiveness of both sensors for snow detection. Cloud cover is a significant hindrance to monitoring snow cover extent Using MODIS but is not a hindrance to the use of the AMSR-E. These limitations could be mitigated by combining MODIS and AMSR-E data to allow for improved interpretation of snow cover extent and SWE on a daily basis and provide temporal continuity of snow mapping across the boreal forest regions in Canada. The purpose of this study is to investigate if temporal monitoring of snow cover using a combination of MODIS and AMSR-E data could yield a better interpretation of changing snow cover conditions. The MODIS snow mapping algorithm is based on snow detection using the Normalized Difference Snow Index (NDSI) and the Normalized Difference Vegetation Index (NDVI) to enhance snow detection in dense vegetation. (Other spectral threshold tests are also used to map snow using MODIS.) Snow cover under a forest canopy may have an effect on the NDVI thus we use the NDVI in snow detection. A MODIS snow fraction product is also generated but not used in this study. In this study the NDSI and NDVI components of the snow mapping algorithm were calculated and analyzed to determine how they changed through the seasons. A blended snow product, the Air Force Weather Agency and NASA (ANSA) snow algorithm and product has recently been developed. The ANSA algorithm blends the MODIS snow cover and AMSR-E SWE products into a single snow product that has been shown to improve the performance of snow cover mapping. In this study components of the ANSA snow algorithm are used along with additional MODIS data to monitor daily changes in snow cover over the period of 1 February to 30 June 2008.

Riggs, George A.↗

Ground Facility for Vicarious Calibration of Skyborne Sensors

An automated ground facility, for vicarious radiometric calibration of airborne and spaceborne sensors of visible and infrared light has been established. In the term "vicarious calibration," "vicarious" is used in the sense of "in place of another," signifying "in place of laboratory calibration." Vicarious calibration involves the use of ground truth in the form of measurements by ground-viewing radiometers, a Sun-viewing photometer, and meteorological instruments positioned in a ground target area. The instrumentation at the facility includes a light-emitting-diode spectrometer (LSpec), which consists of eight tripod-mounted, ground-viewing radiometer units containing LEDs biased to operate as photodetectors (instead of light emitters) at their respective wavelengths. The LSpec provides an essentially continuous stream of measurements at eight discrete wavelengths. These are merged with spectral surface-reflectance measurements made on occasional site visits to obtain temporally continuous coverage with high spectral resolution.

Bruegge, Carol↗

Precursory Off-Fault Deformation in Restraining and Releasing Step Overs: Insights From Discrete Element Method Models

Accelerating geophysical activity is detected preceding some, but not all, large earthquakes. This observation may indicate that no precursors occur before some earthquakes, or that the instrumentation lacks the required sensitivity. Here, to aid crustal monitoring efforts, we use discrete element method models to identify the locations and styles of deformation that may provide useful information about approaching fault reactivation. We model the reactivation of two healed rough faults in a variety of step over configurations, embedded in a host rock with varying amounts of damage subject to shear velocity loading parallel to the faults. Both the fault geometry and ratio of fault to host rock strength control the amount of off-fault deformation. Consistent with field observations, models with larger steps and more preexisting host rock damage produce higher amounts of off-fault deformation. We assess the size of the continuous regions of high velocities and strains to compare the value of the precursory information of each velocity and strain component. Comparing the three components of the velocity vector suggests that the fault-parallel velocity produces the largest and most temporally continuous regions of elevated velocity. The size of these regions increases toward failure, indicating the usefulness of tracking this component. Comparing the volumetric and shear components of the three-dimensional strain tensor suggests that during most of the interseismic period, the shear strain provides more information about approaching fault slip than the volumetric strain. However, in the days and months preceding fault reactivation, both the shear and volumetric strains provide similarly valuable information.

58 GEOSCIENCES↗

Towards an A-train Aerosol Assimilation System: Assimilation of MODIS Aerosol Optical Thickness Retrievals into a Global Aerosol Transport and Radiation Model

To understand the role of atmospheric aerosols in Earth's climate system we ideally require a global, three-dimensional, and temporally continuous system of measurements. The view of aerosols provided by a single polar-orbiting satellite instrument (e.g., MODIS), while providing the greatest spatial coverage of any measurement system available, can offer only a single relevant observation per day of a particular location. Furthermore, the quantitative interpretation of those observations is subject to considerable uncertainty because of we do not know the aerosol optical properties or vertical profiles well enough. Ground-based and in situ measurements can provide missing detail about aerosol vertical and optical properties, but lack the spatial (and sometimes temporal) coverage desired. Chemical transport models treating the aerosol lifecycle complement this picture by providing three-dimensional, global, and continuous aerosol fields, but suffer from poor constraints on their initial conditions (e.g., the aerosol source, sink, and transformation mechanisms). Here we present a global, three- dimensional aerosol transport model and assimilation system. The model assimilates the MODIS aerosol optical thickness product to constrain its aerosol distributions and burdens. We provide an analysis of the model aerosol distributions as constrained by ground-based sun photometer measurements and aerosol products from other satellite measurements.

Colarco, P.↗

Soil Moisture and Sea Surface Salinity Derived from Satellite-Borne Sensors

The monitoring of soil moisture and sea surface salinity over the Earth has been profoundly enhanced during the last thirteen years owing to new generation of satellite sensors. L-band radiometry is currently the only technology providing direct measurements of soil moisture, insensitive to surface roughness and distribution of elements in the soil, and the only technology that allows to measure sea surface salinity from space. The SMOS (Soil Moisture and Ocean Salinity) and SMAP (Soil Moisture Active Passive) satellite missions resolve global and local variability with a spatial resolution of approximately 43 km, a swath width close to 1000 km, and a sampling time, for each mission, of at least twice every 3 days. These resolutions and samplings can be increased by either merging data from the two sensors, and with complementary informations gathered with other passive or active sensors, or with in situ information at higher spatial resolution. Numerous scientific studies based on the use of this new type of measurement have led to a better understanding and constraint of the processes governing the variability of the water cycle, ocean circulation and the Earth's climate. The continuity of measurements, and the increased spatial and radiometric resolution is critical for fulfilling scientific needs. Future L-band radiometry missions currently being planned in Europe (the Copernicus Imaging Microwave Radiometer), and in China (the Ocean Salinity Mission) should provide better constraint on auxiliary parameters by combining multiple frequencies, but they will not have improved spatial resolution beyond SMOS and SMAP. The temporal continuity with SMOS and SMAP will most likely not be ensured. In parallel, new concepts are being developed to increase spatial resolution of both land and ocean parameters.

J. Boutin↗

Impact of a Merged Precipitation Data on Global Soil Moisture Variability

Accurate soil moisture information has proved to be important to climate simulations and climate and weather forecasts. However, many difficulties exist that limit our understanding of soil moisture distribution and variability. One of them is the lack of accurate precipitation with appropriate spatial and temporal resolution. Precipitation as an input forcing to the land surface greatly influences soil moisture characteristics and variability. To improve precipitation data quality, an algorithm has been developed to generate a spatially and temporally continuous 3-hourly global precipitation data for the period of 1987 to present. This precipitation product is a combination of the precipitation from Special Sensor Microwave Imager (SSMI) with the Goddard Earth Observing System-1 Data Assimilation System (GEOS-1 DAS) employing a Physical-space Statistical Analysis System (PSAS). In this study we investigate the impact of this merged/analyzed precipitation data on the global soil moisture variability using an Off-line Land-surface GEOS Assimilation (OLGA) system. Two OLGA integrations starting from 1987 to 1993 are performed forced with the analyzed and GEOS-1 DAS precipitation respectively. We examine the spatial and temporal characteristics of soil moisture variability in response to the analyzed precipitation. The influence of this merged precipitation on the soil moisture variability and regional hydrological budget is estimated throughout the comparison with the results forced with the GEOS-1 DAS precipitation only. In the OLGA the sut@-grid scale horizontal heterogeneity is explicitly represented on the tile space. This provides a means to assess the role of the surface moisture heterogeneity in the interaction with the surface atmosphere and surface hydrological budget, and to validate OLGA results at tile space with in situ observation. ABRACOS (Anglo-Brazilian Amazonian Climate Observation Study), FIFE (First ISLSCP Field Experiment) I and HAPEX data will be used in the validation.

Yang, Runhua↗

Using the SPoRT POES/GOES Hybrid Product in OCONUS Forecasting

The SPoRT (Short-term Prediction and Research Transition) Program at the NASA/Marshall Space Flight Center has been providing unique NASA and NOAA data and techniques to partner Weather Forecast Offices (WFOs) for ten years. Data are provided in the Decision Support System used by WFO forecasters: AWIPS. For the last couple of years, SPoRT has been producing the POES/GOES Hybrid. This suite of products combines the strength ofl5- minute animations of GOES imagery - providing temporal continuity, with the higher resolution, relatively random availability, of polar orbiting (POES) imagery data. The product was first introduced with only MODIS data from NASA's Terra and Aqua satellites, but recently the VIIRS instrument onboard the Suomi-NPP satellite was added, providing better high-resolution coverage. These products represent SPoRT's efforts to prepare for higher resolution, higher frequency GOES-R imagery - as well as helping to move VIIRS (JPSS) data into the mainstream of weather forecasting. SPoRT generates 5 products for this dataset: Visible, Longwave Infrared (11 micrometers), Shortwave IR (3.7 micrometers), Water Vapor (6.7 micrometers), and Fog (Difference of 11 micrometer and 3.7 micrometer channels). The Water Vapor hybrid product has a Red-Blue-Green image from MODIS inlaid, since it provides even more qualitative information than water vapor alone. Animated examples of the products will be shown in this presentation. While the resolution at nadir of GOES imagery is nominally Han (4km for IR channels), the inlaid polar orbiter imagery has a resolution of 250m (lkm for IR channels). This has tremendous application in the continental US. However, in high latitudes, since the usefulness of GOES degrades poleward rapidly, the contrast of GOES and POES data is stark. The consistent temporal nature of GOES, even though at a reduced resolution at high latitudes, provides basic situational awareness, but the introduction of polar data is very helpful in seeing the big picture with clarity - even if only briefly. This presentation will offer real situations where these products helped forecasters make better informed decisions quickly. Plans to augment the product further with the addition of data from several A VHRR instruments will be described.

Smith, Matt↗

The Global LAnd Surface Satellite (GLASS) evapotranspiration product Version 5.0: Algorithm development and preliminary validation

An accurate estimation of spatially and temporally continuous global terrestrial evapotranspiration (ET) is essential in the assessment of surface energy, water and carbon cycles. The Global LAnd Surface Satellite (GLASS) ET product Version 4.0 (v4.0) based on the Bayesian model averaging (BMA) method was generated to estimate global terrestrial ET. However, certain uncertainty for the GLASS ET product v4.0 limits its application. In this study, we introduced the deep neural networks (DNN) merging framework to improve terrestrial ET estimation for GLASS ET product Version 5.0 (v5.0) generation by integrating five satellite-derived ET products [Moderate Resolution Imaging Spectroradiometer (MODIS) ET product (MOD16), Shuttleworth–Wallace dual-source ET product (SW), Priestley–Taylor-based ET product (PT-JPL), modified satellite-based Priestley–Taylor ET product (MS-PT) and simple hybrid ET product (SIM)]. We compared the performance of DNN method against other merging methods, including GLASS ET algorithm v4.0 (BMA), the gradient boosting regression tree (GBRT) method and the random forest (RF) method, based on 195 global eddy covariance (EC) flux towers covering observations from 2000 through 2015. Validations indicated that the DNN had the highest accuracy among four merging methods across different land cover types, yielding the highest average determination coefficients (R 2 , 0.62), root-mean-squared-error (RMSE, 24.1 W/m 2 ) and Kling–Gupta efficiency (KGE, 0.77) with a of 99% confidence interval. Compared with GLASS ET algorithm v4.0, the DNN improved on the R 2 by approximately 7% (p < 0.01) and the KGE by 10%. Based on the DNN, we then generated 8-day GLASS ET product v5.0 globally with a 1 km spatial resolution from 2001 to 2015 driven by GLASS vegetation and surface net radiation (R n ) datasets and Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA2) datasets. Finally, this global terrestrial ET product provides a valuable dataset for monitoring regional and global water resources and environmental changes.

54 ENVIRONMENTAL SCIENCES↗

deadtrees.earth — An open-access and interactive database for centimeter-scale aerial imagery to uncover global tree mortality dynamics

Excessive tree mortality is a global concern and remains poorly understood as it is a complex phenomenon. We lack global and temporally continuous coverage on tree mortality data. Ground-based observations on tree mortality, e.g., derived from national inventories, are very sparse, and may not be standardized or spatially explicit. Earth observation data, combined with supervised machine learning, offer a promising approach to map overstory tree mortality in a consistent manner over space and time. However, global-scale machine learning requires broad training data covering a wide range of environmental settings and forest types. Low altitude observation platforms (e.g., drones or airplanes) provide a cost-effective source of training data by capturing high-resolution orthophotos of overstory tree mortality events at centimeter-scale resolution. Here, we introduce deadtrees.earth, an open-access platform hosting more than two thousand centimeter-resolution orthophotos, covering more than 1,000,000 ha, of which more than 58,000 ha are manually annotated with live/dead tree classifications. This community-sourced and rigorously curated dataset can serve as a comprehensive reference dataset to uncover tree mortality patterns from local to global scales using space-based Earth observation data and machine learning models. This will provide the basis to attribute tree mortality patterns to environmental changes or project tree mortality dynamics to the future. The open nature of deadtrees.earth, together with its curation of high-quality, spatially representative, and ecologically diverse data will continuously increase our capacity to uncover and understand tree mortality dynamics.

Citizen science↗

Dispersion and the speed-limited particle-in-cell algorithm

This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. Here, these results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗