Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte-Carlo 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Improving Mars-GRAM: Increasing the Accuracy of Sensitivity Studies at Large Optical Depths

Extensively utilized for numerous mission applications, the Mars Global Reference Atmospheric Model (Mars-GRAM) is an engineering-level atmospheric model. In a Monte-Carlo mode, Mars-GRAM's perturbation modeling capability is used to perform high fidelity engineering end-to-end simulations for entry, descent, and landing (EDL). Mars-GRAM has been found to be inexact when used during the Mars Science Laboratory (MSL) site selection process for sensitivity studies for MapYear=0 and large optical depth values such as tau=3. Mars-GRAM is based on the NASA Ames Mars General Circulation Model (MGCM) from the surface to 80 km altitude. Mars-GRAM with the MapYear parameter set to 0 utilizes results from a MGCM run with a fixed value of tau=3 at all locations for the entire year. Imprecise atmospheric density and pressure at all altitudes is a consequence of this use of MGCM with tau=3. Density factor values have been determined for tau=0.3, 1 and 3 as a preliminary fix to this pressure-density problem. These factors adjust the input values of MGCM MapYear 0 pressure and density to achieve a better match of Mars-GRAM MapYear 0 with Thermal Emission Spectrometer (TES) observations for MapYears 1 and 2 at comparable dust loading. These density factors are fixed values for all latitudes and Ls and are included in Mars-GRAM Release 1.3. Work currently being done, to derive better multipliers by including variations with latitude and/or Ls by comparison of MapYear 0 output directly against TES limb data, will be highlighted in the presentation. The TES limb data utilized in this process has been validated by a comparison study between Mars atmospheric density estimates from Mars-GRAM and measurements by Mars Global Surveyor (MGS). This comparison study was undertaken for locations on Mars of varying latitudes, Ls, and LTST. The more precise density factors will be included in Mars-GRAM 2005 Release 1.4 and thus improve the results of future sensitivity studies done for large optical depths.

Justh, Hilary L.↗
NASA NTRS1973

Monte-Carlo analysis of rarefied-gas diffusion including variance reduction using the theory of Markov random walks

Molecular diffusion through a rarefied gas is analyzed by using the theory of Markov random walks. The Markov walk is simulated on the computer by using random numbers to find the new states from the appropriate transition probabilities. As the sample molecule during its random walk passes a scoring position, which is a location at which the macroscopic diffusing flow variables such as molecular flux and molecular density are desired, an appropriate payoff is scored. The payoff is a function of the sample molecule velocity. For example, in obtaining the molecular flux across a scoring position, the random walk payoff is the net number of times the scoring position has been crossed in the positive direction. Similarly, when the molecular density is required, the payoff is the sum of the inverse velocity of the sample molecule passing the scoring position. The macroscopic diffusing flow variables are then found from the expected payoff of the random walks.

Perlmutter, M.↗
NASA NTRS1982

A study of the feasibility of statistical analysis of airport performance simulation

The feasibility of conducting a statistical analysis of simulation experiments to study airport capacity is investigated. First, the form of the distribution of airport capacity is studied. Since the distribution is non-Gaussian, it is important to determine the effect of this distribution on standard analysis of variance techniques and power calculations. Next, power computations are made in order to determine how economic simulation experiments would be if they are designed to detect capacity changes from condition to condition. Many of the conclusions drawn are results of Monte-Carlo techniques.

Myers, R. H.↗
NASA NTRS1997

CGRO Guest Investigator Program

The following are highlights from the research supported by this grant: (1) Theory of gamma-ray blazars: We studied the theory of gamma-ray blazars, being among the first investigators to propose that the GeV emission arises from Comptonization of diffuse radiation surrounding the jet, rather than from the synchrotron-self-Compton mechanism. In related work, we uncovered possible connections between the mechanisms of gamma-ray blazars and those of intraday radio variability, and have conducted a general study of the role of Compton radiation drag on the dynamics of relativistic jets. (2) A Nonlinear Monte Carlo code for gamma-ray spectrum formation: We developed, tested, and applied the first Nonlinear Monte Carlo (NLMC) code for simulating gamma-ray production and transfer under much more general (and realistic) conditions than are accessible with other techniques. The present version of the code is designed to simulate conditions thought to be present in active galactic nuclei and certain types of X-ray binaries, and includes the physics needed to model thermal and nonthermal electron-positron pair cascades. Unlike traditional Monte-Carlo techniques, our method can accurately handle highly non-linear systems in which the radiation and particle backgrounds must be determined self-consistently and in which the particle energies span many orders of magnitude. Unlike models based on kinetic equations, our code can handle arbitrary source geometries and relativistic kinematic effects In its first important application following testing, we showed that popular semi-analytic accretion disk corona models for Seyfert spectra are seriously in error, and demonstrated how the spectra can be simulated if the disk is sparsely covered by localized 'flares'.

Begelman, Mitchell C.↗
NASA NTRS1990

Typical motions in multiple systems

In very old times, people counted - one, two, many. The author wants to show that they were right. Consider the motions of isolated bodies: (1) N = 1 - simple motion; (2) N = 2 - Keplerian orbits; and (3) N = 3 - this is the difficult problem. In general, this problem can be studied only by computer simulations. The author studied this problem over many years (see, e.g., Agekian and Anosova, 1967; Anosova, 1986, 1989 a,b). The principal result is that two basic types of dynamics take place in triple systems. The first special type is the stable hierarchical systems with two almost Keplerian orbits. The second general type is the unstable triple systems with complicated motions of the bodies. By random choice of the initial conditions, by the Monte-Carlo method, the stable systems comprised about approx. 10% of the examined cases; the unstable systems comprised the other approx. 90% of cases under consideration. In N greater than 3, the studies of dynamics of such systems by computer simulations show that we have in general also the motions roughly as at the cases 1 - 3 with the relative negative or positive energies of the bodies. In the author's picture, the typical trajectories of the bodies in unstable triple systems of the general type of dynamics are seen. Such systems are disrupted always after close triple approaches of the bodies. These approaches play a role like the gravitational slingshot. Often, the velocities of escapers are very large. On the other hand, the movie also shows the dynamical processes of a formation, dynamical evolution and disruption of the temporary wide binaries in triples and a formation of final hard massive binaries in the final evolution of triples.

Anosova, Joanna P.↗
DOE OSTI2025

Moment-preserving Monte-Carlo Coulomb collision method for particle codes

Binary-pairing Monte-Carlo methods are widely used in particle-in-cell codes to capture effects of small angle Coulomb collisions. These methods preserve momentum and energy exactly when the simulation particles have equal weights. However, when the interacting particles are of varying weight, these physical conservation laws are only preserved on average. Here, we 1) extend these methods to weighted particles such that the scattering physics is correct on average, and 2) describe a new method for adjusting the particle velocities post scatter to restore exact conservation of momentum and energy. In conclusion, the efficacy of the model is illustrated with various test problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

Methodolgy For Evaluation Of Technology Impacts In Space Electric Power Systems

The Analysis and Management branch of the Power and Propulsion Office at NASA Glenn Research Center is responsible for performing complex analyses of the space power and In-Space propulsion products developed by GRC. This work quantifies the benefits of the advanced technologies to support on-going advocacy efforts. The Power and Propulsion Office is committed to understanding how the advancement in space technologies could benefit future NASA missions. They support many diverse projects and missions throughout NASA as well as industry and academia. The area of work that we are concentrating on is space technology investment strategies. Our goal is to develop a Monte-Carlo based tool to investigate technology impacts in space electric power systems. The framework is being developed at this stage, which will be used to set up a computer simulation of a space electric power system (EPS). The outcome is expected to be a probabilistic assessment of critical technologies and potential development issues. We are developing methods for integrating existing spreadsheet-based tools into the simulation tool. Also, work is being done on defining interface protocols to enable rapid integration of future tools. Monte Carlo-based simulation programs for statistical modeling of the EPS Model. I decided to learn and evaluate Palisade's @Risk and Risk Optimizer software, and utilize it's capabilities for the Electric Power System (EPS) model. I also looked at similar software packages (JMP, SPSS, Crystal Ball, VenSim, Analytica) available from other suppliers and evaluated them. The second task was to develop the framework for the tool, in which we had to define technology characteristics using weighing factors and probability distributions. Also we had to define the simulation space and add hard and soft constraints to the model. The third task is to incorporate (preliminary) cost factors into the model. A final task is developing a cross-platform solution of this framework.

Holda, Julie↗
NASA NTRS1960

An Exploratory Statistical Analysis of a Planet Approach-Phase Guidance Scheme Using Angular Measurements with Significant Error

An exploratory analysis of vehicle guidance during the approach to a target planet is presented. The objective of the guidance maneuver is to guide the vehicle to a specific perigee distance with a high degree of accuracy and minimum corrective velocity expenditure. The guidance maneuver is simulated by considering the random sampling of real measurements with significant error and reducing this information to prescribe appropriate corrective action. The instrumentation system assumed includes optical and/or infrared devices to indicate range and a reference angle in the trajectory plane. Statistical results are obtained by Monte-Carlo techniques and are shown as the expectation of guidance accuracy and velocity-increment requirements. Results are nondimensional and applicable to any planet within limits of two-body assumptions. The problem of determining how many corrections to make and when to make them is a consequence of the conflicting requirement of accurate trajectory determination and propulsion. Optimum values were found for a vehicle approaching a planet along a parabolic trajectory with an initial perigee distance of 5 radii and a target perigee of 1.02 radii. In this example measurement errors were less than i minute of arc. Results indicate that four corrections applied in the vicinity of 50, 16, 15, and 1.5 radii, respectively, yield minimum velocity-increment requirements. Thrust devices capable of producing a large variation of velocity-increment size are required. For a vehicle approaching the earth, miss distances within 32 miles are obtained with 90-percent probability. Total velocity increments used in guidance are less than 3300 feet per second with 90-percent probability. It is noted that the above representative results are valid only for the particular guidance scheme hypothesized in this analysis. A parametric study is presented which indicates the effects of measurement error size, initial perigee, and initial energy on the guidance requirements. Measurement error size significantly affects both guidance accuracy and velocity-increment expenditure. The initial trajectory, as given by its perigee and energy, affects the velocity-increment expenditure but not final guidance accuracy.

Friedlander, Alan L.↗
NASA NTRS2010

Radiation Beamline Testbeds for the Simulation of Planetary and Spacecraft Environments for Human and Robotic Mission Risk Assessment

The Center for Radiation Engineering and Science for Space Exploration (CRESSE) at Prairie View A&M University, Prairie View, Texas, USA, is establishing an integrated, multi-disciplinary research program on the scientific and engineering challenges faced by NASA and the international space community caused by space radiation. CRESSE focuses on space radiation research directly applicable to astronaut health and safety during future long term, deep space missions, including Martian, lunar, and other planetary body missions beyond low earth orbit. The research approach will consist of experimental and theoretical radiation modeling studies utilizing particle accelerator facilities including: 1. NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory; 2. Proton Synchrotron at Loma Linda University Medical Center; and 3. Los Alamos Neutron Science Center (LANSCE) at Los Alamos National Laboratory. Specifically, CRESSE investigators are designing, developing, and building experimental test beds that simulate the lunar and Martian radiation environments for experiments focused on risk assessment for astronauts and instrumentation. The testbeds have been designated the Bioastronautics Experimental Research Testbeds for Environmental Radiation Nostrum Investigations and Education (BERT and ERNIE). The designs of BERT and ERNIE will allow for a high degree of flexibility and adaptability to modify experimental configurations to simulate planetary surface environments, planetary habitats, and spacecraft interiors. In the nominal configuration, BERT and ERIE will consist of a set of experimental zones that will simulate the planetary atmosphere (Solid CO2 in the case of the Martian surface.), the planetary surface, and sub-surface regions. These experimental zones can be used for dosimetry, shielding, biological, and electronic effects radiation studies in support of space exploration missions. BERT and ERNIE are designed to be compatible with the experimental areas associated with the above facilities. CRESSE has broad expertise in space radiation in the areas of space radiation environment modeling, Monte-Carlo radiation transport modeling, space radiation instrumentation and dosimetry, radiation effects on electronics, and multi-functional composite shielding materials. The BERT and ERNIE testbeds will be utilized in individual and collaborative research incorporating this expertise. The research goal is to maximize the technical readiness level (TRL) of radiation instrumentation for human and robotic missions, optimizing the return value of CRESSE for NASA exploration and international co-operative missions. Outcomes and knowledge from research utilizing BERT and ERNIE will be applied to a variety of scientific and engineering disciplines vital for safe and reliable execution of future space exploration missions, which can be negatively impacted by the space radiation environment. The testbeds will be central to a variety of university educational activities and educational goals of NASA. Specifically, BERT and ERNIE will enhance educational opportunities in science, technology, engineering and mathematics (STEM) disciplines for engineering and science students at PVAMU, a historically black college/university. Preliminary data on prototype testbed configurations, including simulated lunar regolith (JSC-1A stimulant based on Apollo 11 samples), regolith/polyethylene composites, and dry ice, will be presented to demonstrate the usefulness of BERT and ERNIE in radiation beam line experiments.

Wilkins, Richard↗
NASA NTRS2007

A Risk-Based Approach for Aerothermal/TPS Analysis and Testing

The current status of aerothermal and thermal protection system modeling for civilian entry missions is reviewed. For most such missions, the accuracy of our simulations is limited not by the tools and processes currently employed, but rather by reducible deficiencies in the underlying physical models. Improving the accuracy of and reducing the uncertainties in these models will enable a greater understanding of the system level impacts of a particular thermal protection system and of the system operation and risk over the operational life of the system. A strategic plan will be laid out by which key modeling deficiencies can be identified via mission-specific gap analysis. Once these gaps have been identified, the driving component uncertainties are determined via sensitivity analyses. A Monte-Carlo based methodology is presented for physics-based probabilistic uncertainty analysis of aerothermodynamics and thermal protection system material response modeling. These data are then used to advocate for and plan focused testing aimed at reducing key uncertainties. The results of these tests are used to validate or modify existing physical models. Concurrently, a testing methodology is outlined for thermal protection materials. The proposed approach is based on using the results of uncertainty/sensitivity analyses discussed above to tailor ground testing so as to best identify and quantify system performance and risk drivers. A key component of this testing is understanding the relationship between the test and flight environments. No existing ground test facility can simultaneously replicate all aspects of the flight environment, and therefore good models for traceability to flight are critical to ensure a low risk, high reliability thermal protection system design. Finally, the role of flight testing in the overall thermal protection system development strategy is discussed.

Wright, Michael J.↗
NASA NTRS2003

Radiation Risk Projections for Space Travel

Space travelers are exposed to solar and galactic cosmic rays comprised of protons and heavy ions moving with velocities close to the speed of light. Cosmic ray heavy ions are known to produce more severe types of biomolecular damage in comparison to terrestrial forms of radiation, however the relationship between such damage and disease has not been fully elucidated. On Earth, we are protected from cosmic rays by atmospheric and magnetic shielding, and only the remnants of cosmic rays in the form of ground level muons and other secondary radiations are present. Because human epidemiology data is lacking for cosmic rays, risk projection must rely on theoretical understanding and data from experimental models exposed to space radiation using charged particle accelerators to simulate space radiation. Although the risks of cancer and other late effects from cosmic rays are currently believed to present a severe challenge to space travel, this challenge is centered on our lack of confidence in risk projections methodologies. We review biophysics and radiobiology data on the effects of the cosmic ray heavy ions, and the current methods used to project radiation risks . Cancer risk projections are described as a product of many biological and physical factors, each of which has a differential range of uncertainty due to lack of data and knowledge. Risk projections for space travel are described using Monte-Carlo sampling from subjective error di stributions that represent the lack of knowledge in each factor that contributes to the projection model in order to quantify the overall uncertainty in risk projections. This analysis is applied to space mi ssion scenarios including lunar colony, deep space outpost, and a Mars mission. Results suggest that the number of days in space where cancer mortality risks can be assured at a 95% confidence level to be below the maximum acceptable risk for radi ation workers on Earth or the International Space Station is only on the order of 100-200 days. Approaches to reduce these unceI1ainties and mitigate risks are described.

Cucinotta, Francis↗
DOE OSTI2025

Data-Efficient Strategies for Probabilistic Voltage Envelopes under Network Contingencies

This work presents an efficient data-driven method to construct probabilistic voltage envelopes (PVE) using power flow learning in grids with network contingencies. First, a network-aware Gaussian process (GP) termed Vertex-Degree Kernel (VDK-GP), developed in prior work, is used to estimate voltage–power functions for a few network configurations. The paper introduces a novel multi-task vertex degree kernel (MT-VDK) that amalgamates the learned VDK-GPs to determine power flows for unseen networks, with a significant reduction in the computational complexity and hyperparameter requirements compared to alternate approaches. Simulations on the IEEE 30-Bus network demonstrate the retention and transfer of power flow knowledge in both N-1 and N-2 contingency scenarios. The MT-VDK-GP approach achieves over 50 % reduction in mean prediction error for novel N-1 contingency network configurations in low training data regimes (50–250 samples) over VDK-GP. Additionally, MT-VDK-GP outperforms a hyper-parameter based transfer learning approach in over 75 % of N-2 contingency network structures, even without historical N-2 outage data. Furthermore, the proposed method demonstrates the ability to achieve PVEs using sixteen times fewer power flow solutions compared to Monte-Carlo sampling-based methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2002

Partially Decentralized Control Architectures for Satellite Formations

In a partially decentralized control architecture, more than one but less than all nodes have supervisory capability. This paper describes an approach to choosing the number of supervisors in such au architecture, based on a reliability vs. cost trade. It also considers the implications of these results for the design of navigation systems for satellite formations that could be controlled with a partially decentralized architecture. Using an assumed cost model, analytic and simulation-based results indicate that it may be cheaper to achieve a given overall system reliability with a partially decentralized architecture containing only a few supervisors, than with either fully decentralized or purely centralized architectures. Nominally, the subset of supervisors may act as centralized estimation and control nodes for corresponding subsets of the remaining subordinate nodes, and act as decentralized estimation and control peers with respect to each other. However, in the context of partially decentralized satellite formation control, the absolute positions and velocities of each spacecraft are unique, so that correlations which make estimates using only local information suboptimal only occur through common biases and process noise. Covariance and monte-carlo analysis of a simplified system show that this lack of correlation may allow simplification of the local estimators while preserving the global optimality of the maneuvers commanded by the supervisors.

Carpenter, J. Russell↗
DOE OSTI2025

Enhancement of PyARC for Westinghouse Electric Company’s Lead Fast Reactor Design and Modeling (Final TCF Report)

Westinghouse Electric Company is a nuclear reactor vendor headquartered in the U.S. that is developing advanced reactor technology for the U.S. and global markets. Westinghouse has been relying on the neutronics Argonne Reactor Codes (ARC) executed through the NEAMS Workbench and its PyARC module that are developed under the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) and Advanced Reactor Technology (ART) – Fast Reactor programs. Through this user experience, Westinghouse identified several enhancements that would benefit the ARC codes’ usability by the US industry and therefore its commercialization potential. The enhancements were proposed to deliver both improvements in workflow and analysis capabilities to better support effective fast reactor core design and analysis to the nuclear industry. The PyARC workflow was extended in this project by integrating non-neutronic ARC codes DASSH and NUBOW-3D. The Ducted Assembly Steady-State Heat equation (DASSH) code is developed at ANL to perform steady-state thermal hydraulic sub-channel analysis in liquid metal fast reactor assemblies to determine optimized coolant flow and temperature distributions, which in this project was updated and validated for lead fast reactor (LFR) applications. The interface between REBUS and NUBOW-3D were improved in this project to assess the impact of the core restraint design and thermal induced expansion effects on the reactivity of the core, and to model the deformations of the fuel assemblies induced by temperature and irradiation. Finally, the ARC models that were extensively verified and validated through various SFR-based modeling benchmarks are extended in this project through code-to-code comparison on relevant LFR-specific neutronics benchmarks against Monte-Carlo neutronic solutions. Overall, this work enables verification of the capability of the ARC codes for a wide range of Generation-IV reactor designs. The outcome of this project is the release of a comprehensive modeling toolkit of validated, robust and efficient codes, as well as their user interface, that enables industry to perform a wide range of fast reactor analyses for design and licensing of their concepts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Decoupling thermal and irradiation effects on grain boundary segregation

Radiation-induced segregation (RIS) is most often measured by peak solute concentration at a boundary. However, this may give an incomplete picture of segregation quantity and phenomena. Radiation-induced and thermal segregation at grain boundaries was investigated in Fe-9.6 at.% Cr after 9 MeV Fe 3+ ion irradiation at 400 °C. The experimental results were compared to kinetic Monte Carlo (KMC) simulations. The study revealed that Cr enrichment (peak segregation) at the grain boundaries was comparable in both the irradiated and non-irradiated conditions, although irradiation resulted in broader segregation profiles in both experiments and simulations, indicating an overall increase in grain boundary segregation due to irradiation. This broadening is attributed to back diffusion into the grain interior. While it is an established phenomenon, this study offers a quantitative evaluation using experimental data and KMC modeling. Further, these results emphasize the importance of analyzing the entire segregation profile and decoupling the thermal and irradiation contributions to solute segregation.

Grain boundary segregation↗
DOE OSTI2026

Modeling and Experimental Demonstration of Flux Spreading in Light Trapping Planar-Cavity Solar Enclosed Particle Receivers

This study experimentally validates and numerically models the flux-spreading effect in a light-trapping planar-cavity solar receiver) for particle-based concentrating solar power systems. The receiver's shallow cavity with vertical planar walls redistributes concentrated solar flux, reducing peak intensity and achieving uniform heat flux. On-sun tests at National Renewable Energy Laboratory's High-Flux Solar Furnace under flux up to 1500 kW/m2 measured cavity wall temperatures, which were compared with Monte Carlo ray-tracing (SolTrace) and computational numerical simulations. Three angular absorptance models were evaluated: constant absorptance, a Pyromark-based directional model, and a Fresnel-based Cr2O3 model. The Fresnel-derived model showed the best agreement with experiments, achieving high correlation (PC > 0.85), structural similarity (SSIM > 0.98), and signal-to-noise ratios (PSNR > 40 dB), with temperature prediction errors of 1-11%. Results confirm that flux spreading mitigates local overheating and validate the integrated modeling approach, supporting the solar receiver scalability for high-efficiency, high-temperature concentrating solar power applications.

14 SOLAR ENERGY↗
DOE OSTI2025

Particle-based modelling of axisymmetric tandem mirror devices

In this work, we describe the use of a 1D-2V quasi-neutral hybrid electrostatic PIC with Monte-Carlo Coulomb collisions and non-uniform magnetic field to model the parallel transport and confinement in an axisymmetric tandem mirror device. End-plugs, based on simple-mirrors, are positioned at each end of the device and fueled with neutral beams (25 and 100 keV) to produce a sloshing ion population and increase the density of the end-plugs relative to the central cell. Results show the formation of a potential difference barrier between the central cell and the end-plugs. This potential confines a large fraction of the low energy thermal ions in the central cell which would otherwise be lost in a simple mirror, demonstrating the advantage of the beam-driven tandem mirror configuration relative to simple mirrors. In addition, we explore the effect of end-plug electron temperature on the confinement time of the device and compare it with theoretical estimates. Finally, we discuss the limitations of the code in its present form and describe the next logical steps to improve its predictive capability such as a fully nonlinear Fokker–Planck collision operator, multiply nested flux surface solutions and modeling the exhaust region up to the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Monte-Carlo model of pitch-angle scattering in solar cosmic ray events

Monte Carlo simulations of the propagation of solar cosmic-rays in interplanetary space are reported, including the effects of pitch-angle scattering and adiabatic focusing. Time intensity profiles agree well with the corresponding spatial diffusion approximation for models close to those used by Gombosi and Owens (1980). Monte Carlo simulations of the same model problem yield results that disagree with those of Gombosi and Owens, and the most likely possible source of error in the Monte Carlo simulation is considered to be in the region near the 90 deg pitch angle where the finite step size introduces certain inconsistencies.

Palmer, I. D.↗