Search NASA⌕ Search

SEARCH · Search NASA

Results for “control simulation testing”

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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Reactive transport modeling and characterization of concrete materials with fly ash replacement under carbonation attack - 15477

There currently lacks sufficient knowledge to accurately characterize structural, hydraulic, and chemical performance of cementitious materials that are required to contain and shield nuclear waste for thousands of years. There are a number of degradation methods that can affect the integrity of these materials over this time frame including carbonation which lowers the alkalinity of the material and can induce cracking in materials with steel reinforcement. The rate and extent of carbonation in these materials is impacted by the addition of fly ash and its composition which can vary depending on the fly ash source. To characterize the carbonation reaction, the changes in the microstructure and elemental distribution as a function of extent of carbonation and fly ash type were evaluated using scanning electron microscopy (SEM) coupled with energy dispersive x-ray spectroscopy (EDS). Micro-concretes were prepared with no fly ash (control case) and with fly ash of various calcium compositions. The progression of the carbonation front was observed in the backscattered electron images by the formation of calcium carbonate and depth was confirmed with a phenolphthalein pH indicator test. The migration and deposition of chemical species were observed in carbonated regions from EDS data and were explained based on the changes in material alkalinity and solubility determined from EPA Method 1313 (Liquid-solid partitioning (LSP) as a function of eluate pH). The changes in mass transport properties from carbonation were observed using EPA Method 1315 (Mass transfer rates in monolithic and compacted granular materials using semi-dynamic tank leaching procedures). The carbonation reaction in each material was simulated using LeachXS/ORCHESTRA, a reactive transport modeling tool, which provided further evidence of localized deposition and migration of chemical species due to changes in speciation and solubility. (authors)

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

Correction of Dual-PRF Doppler Velocity Outliers in the Presence of Aliasing

Abstract In Doppler weather radars, the presence of unfolding errors or outliers is a well-known quality issue for radial velocity fields estimated using the dual–pulse repetition frequency (PRF) technique. Postprocessing methods have been developed to correct dual-PRF outliers, but these need prior application of a dealiasing algorithm for an adequate correction. This paper presents an alternative procedure based on circular statistics that corrects dual-PRF errors in the presence of extended Nyquist aliasing. The correction potential of the proposed method is quantitatively tested by means of velocity field simulations and is exemplified in the application to real cases, including severe storm events. The comparison with two other existing correction methods indicates an improved performance in the correction of clustered outliers. The technique proposed is well suited for real-time applications requiring high-quality Doppler radar velocity fields, such as wind shear and mesocyclone detection algorithms, or assimilation in numerical weather prediction models.

58 GEOSCIENCES↗
DOE OSTI2022

Charge Delocalization and Vibronic Couplings in Quadrupolar Squaraine Dyes

Squaraines are prototypical quadrupolar charge-transfer chromophores that have recently attracted much attention as building blocks for solution-processed photovoltaics, fluorescent probes with large two-photon absorption cross sections, and aggregates with large circular dichroism. Their optical properties are often rationalized in terms of phenomenological essential state models, considering the coupling of two zwitterionic excited states to a neutral ground state. As a result, optical transitions to the lowest S 1 excited state are one-photon allowed, whereas the next higher S 2 state can only be accessed by two-photon transitions. A further implication of these models is a substantial reduction of vibronic coupling to the ubiquitous high-frequency vinyl-stretching modes of organic materials. Here, in this study, we combine time-resolved vibrational spectroscopy, two-dimensional electronic spectroscopy, and quantum-chemical simulations to test and rationalize these predictions for nonaggregated molecules. We find small Huang–Rhys factors below 0.01 for the high-frequency, 1500 cm –1 modes in particular, as well as a noticeable reduction for those of lower frequency modes in general for the electronic S 0 → S 1 transition. The two-photon allowed state S 2 is well separated energetically from S 1 and has weak vibronic signatures as well. Thus, the resulting pronounced concentration of the oscillator strength in a narrow region relevant to the lowest electronic transition makes squaraines and their aggregates exceptionally interesting for strong and ultrastrong coupling of excitons to localized light modes in external resonators with chiral properties that can largely be controlled by the molecular architecture.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Effective restoration strategies of interdependent power system and communication network

The power system is continuously evolving towards a smarter grid by adding intelligent electrical devices that leverage communication technologies to provide visualisation, control and automation. The integration of power system and communication network has brought challenges into the resilience of interdependent networks. Following an outage or blackout, system operators cannot only restore one network while simply ignoring the interconnection with the other network. Effective restoration strategies that involve co‐ordination between power system and communication network have become critical in enhancing the resilience. An optimisation‐based restoration strategy is developed to coordinate the restoration actions between two networks. The problem is modelled as a mixed integer linear programming problem, with the objective of activating every node in both networks with the minimum number of activation/energisation of branches. The developed strategies are tested on the modified IEEE 14‐bus system incorporating communication network. Simulation results demonstrate that the coordinated restoration strategies are effective and critical in the successful recovery of interdependent networks. The developed model could also provide insights in network design and planning to enhance the resilience of interdependent power system and communication network.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

The impact of baryonic potentials on the gravothermal evolution of self-interacting dark matter haloes

The presence of a central baryonic potential can have a significant impact on the gravothermal evolution of self-interacting dark matter (SIDM) haloes. We extend a semi-analytical fluid model to incorporate the influence of a static baryonic potential and calibrate it using controlled N-body simulations. We construct benchmark scenarios with varying baryon concentrations and different SIDM models, including constant and velocity-dependent self-interacting cross-sections. The presence of the baryonic potential induces changes in SIDM halo properties, including central density, core size, and velocity dispersion, and it accelerates the halo’s evolution in both expansion and collapse phases. Furthermore, we observe a quasi-universality in the gravothermal evolution of SIDM haloes with the baryonic potential, resembling a previously known feature in the absence of the baryons. By appropriately rescaling the physical quantities that characterize the SIDM haloes, the evolution of all our benchmark cases exhibits remarkable similarity. In conclusion, our findings offer a framework for testing SIDM predictions using observations of galactic systems where baryons play a significant dynamical role.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Oxygenated Fuel Blending Effects in a Mixing-Controlled Compression-Ignition Engine Equipped with Ducted Fuel Injection [Slides]

Continued creation of harmful emissions such as NO x and soot from compression-ignition engines utilizing mixing-controlled combustion systems (i.e., diesel engines) remains a problem and is the subject of on-going research. The inherently high efficiency, relatively low cost, and numerous other desirable attributes of such engines, coupled with a widely supported infrastructure, motivates their continued advancement. Recently, a scientifically distinct and mechanically simple technology called ducted fuel injection (DFI) has shown a robust ability to allow such engines to operate with simultaneously low engine-out soot and NO x emissions when it is employed with simulated exhaust-gas recirculation. To better understand the property ranges of sustainable, oxygenated-fuel blending stocks that will most improve engine performance, two oxygenated blendstocks were separately blended with a commercial diesel base fuel and tested within a heavy-duty diesel optical engine equipped with a four-duct DFI configuration. Conventional and crank-angle-resolved optical diagnostics were used to elucidate the effects of fuel ignition quality, oxygenate molecular structure, and overall oxygen content on engine performance.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2016

Field Testing and Modeling of Supermarket Refrigeration Systems as a Demand Response Resource

Supermarkets offer a substantial demand response (DR) resource because of their high energy intensity and use patterns; however, refrigeration as the largest load has been challenging to access. Previous work has analyzed supermarket DR using heating, ventilating, and air conditioning; lighting; and anti-sweat heaters. This project evaluated and quantified the DR potential inherent in supermarket refrigeration systems in the Bonneville Power Administration service territory. DR events were carried out and results measured in an operational 45,590-ft2 supermarket located in Hillsboro, Oregon. Key results from the project include the rate of temperature increase in freezer reach-in cases and walk-ins when refrigeration is suspended, the load shed amount for DR tests, and the development of calibrated models to quantify available DR resources. Simulations showed that demand savings of 15 to 20 kilowatts (kW) are available for 1.5 hours for a typical store without precooling and for about 2.5 hours with precooling using only the low-temperature, non-ice cream cases. This represents an aggregated potential of 20 megawatts within BPA's service territory. Inability to shed loads for medium-temperature (MT) products because of the tighter temperature requirements is a significant barrier to realizing larger DR for supermarkets. Store owners are reluctant to allow MT case set point changes, and laboratory tests of MT case DR strategies are needed so that owners become comfortable testing, and implementing, MT case DR. The next-largest barrier is the lack of proper controls in most supermarket displays over ancillary equipment, such as anti-sweat heaters, lights, and fans.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

PI Loop Resonance Control for Dark Photon Experiment at 2 K Using a 2.6 GHz SRF Cavity

Two 2.6 GHz SRF cavities are being used for a dark photon search at the vertical test stand (VTS) in FNAL, for the second phase of the Dark SRF experiment. During testing at 2 K the cavities experience frequency detuning caused by microphonics and slow frequency drifts. The experiment requires that the two cavities have the same frequency within the cavity's bandwidth. These two cavities are equipped with frequency tuners consisting of three piezo actuators. The piezo actuators are used for fine-fast frequency tuning. A proportional-integral (PI) loop utilizing the three piezos on the emitter was used to stabilize the cavity frequency and match the receiver cavity frequency. The results from this implementation will be discussed. The integration time was also calculated via simulation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Controlling and leveraging small-scale information in tomographic galaxy–galaxy lensing

ABSTRACT The tangential shear signal receives contributions from physical scales in the galaxy–matter correlation function well below the transverse scale at which it is measured. Since small scales are difficult to model, this non-locality has generally required stringent scale cuts or new statistics for cosmological analyses. Using the fact that uncertainty in these contributions corresponds to an uncertainty in the enclosed projected mass around the lens, we provide an analytic marginalization scheme to account for this. Our approach enables the inclusion of measurements on smaller scales without requiring numerical sampling over extra free parameters. We extend the analytic marginalization formalism to retain cosmographic (shear-ratio) information from small-scale measurements that would otherwise be removed due to modelling uncertainties, again without requiring the addition of extra sampling parameters. We test the methodology using simulated likelihood analysis of a Dark Energy Survey Year 5-like galaxy–galaxy lensing and galaxy clustering datavector. We demonstrate that we can remove parameter biases due to the presence of an unmodelled one-halo contamination of the galaxy–galaxy lensing signal, and use the shear-ratio information on small scales to improve cosmological parameter constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

DEMO-FTES: Development, Monitoring, and Control of Fracture Thermal Energy Storage in Crystalline Rock Formations (CRADA Final Report)

The DEMO-FTES project investigated the feasibility of Fracture Thermal Energy Storage (FTES) as a seasonal energy storage solution in crystalline rock formations. FTES leverages hydraulically induced fractures to exchange heat between circulating fluids and the surrounding rock mass, enabling long-term thermal energy retention due to the high specific heat and low thermal conductivity of rock. This approach has the potential to reduce heating and cooling energy demands and enhance building energy resilience. The project combined dimensional analysis, numerical modeling, laboratory experiments, and meso-scale field tests to evaluate FTES performance and advance its technology readiness level from 3 to 5. Scaling analysis identified key dimensionless parameters governing heat transfer and fluid flow, ensuring laboratory and field tests were representative of larger-scale systems. Numerical simulations using TOUGH and iTOUGH2 frameworks supported experiment design and interpretation, modeling fracture geometry, thermal-hydraulic behavior, and thermo-mechanical coupling. Laboratory tests at EPFL involved creating single and multiple fractures in 25 cm cubic samples of Gabbro and Granite under true triaxial stress.

25 ENERGY STORAGE↗
DOE OSTI2020

Raspberry colloid-templated approach for the synthesis of palladium-based oxidation catalysts with enhanced hydrothermal stability and low-temperature activity

It is becoming increasingly urgent to develop and utilize novel, more efficient and stable materials for mobile and stationary emission control applications as the deleterious consequences of anthropogenic air pollution are becoming more evident and pressing. Tightening regulations, particularly related to automotive exhaust treatment, together with continued improvements in engine design, that result in lowering the engine operating temperatures and inadvertently lead to the release of an overwhelming proportion of pollutants during the cold start, present new challenges for materials design, specifically for oxidation catalysts. In particular, improvements in the low-temperature activity while maintaining catalyst stability at high temperatures are required from the next-generation catalyst. Typical catalysts for removal of pollutants from automotive exhaust streams incorporate platinum group metals (PGMs). They tend to be inefficient at low temperatures (below 250 °C), thus accounting for the cold start problem, yet sinter and lose their activity at high temperatures that are frequently encountered during catalyst operation. High PGM loadings are often employed to compensate for catalyst inefficiencies and fast degradation, ultimately resulting in high-cost catalytic converters. Here, we have developed a new approach for the design and formation of catalytic materials that allows for both significantly more efficient PGM incorporation and improved overall catalyst performance at reduced PGM loadings. The method provides control over the composition and geometry of the support through self-assembly of sacrificial composite template — “raspberry” polymeric colloids decorated with catalytic particles — accompanied by infiltration with metal-oxide precursor and subsequent removal of the colloids. This method simultaneously structures the porous network and organizes the catalytic particles within it. Uniquely, the resulting catalytic particles are partially embedded in the support matrix and partially exposed to the pore interior, producing catalytic sites that are both stable and accessible. Herein, the feasibility of this novel and versatile approach for automotive catalytic conversion is demonstrated: the studies include testing alumina-based raspberry-colloid-templated (RCT) catalysts containing Pd nanoparticles (RCT Pd/Al2O3) for oxidation of propane and carbon monoxide under simulated diesel exhaust conditions and hydrothermal aging at 800 °C for 50 h in the simulated stream. The RCT Pd/Al2O3 catalysts exhibit exceptional activity toward CO oxidation, reduced reaction onset temperature, and high stability to elevated temperatures (demonstrated through prolonged exposure to temperatures up to 950 °C) and reactive gas streams, without migration, sintering or loss of the precious metal NPs. Notably, the novel catalyst shows the same or slightly better performance than the commercial catalysts even when the PGM load is reduced by ~80 % compared to the commercial counterparts. These results provide confidence for the utilization of the RCT approach for the fabrication of robust nanostructured catalysts for next-generation, energy-efficient catalytic converters with improved performance at low and high temperatures and reduced costs. The RCT methodology is, in addition, highly generalizable, and can thus be applied for the design of a wide range of catalytic systems in the automotive sector and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Artificial Neural Network-based State Estimation for Low Observable, Unbalanced Microgrids for Microgrid Building Blocks

The microgrid building blocks (MBB) were proposed as microgrid components with combined sub-components with power conversion, communication, and microgrid control capability, or a subset of such sub-components. This work addresses the microgrid controller, present in an MBB, which requires accurate state estimation to perform its tasks, including for monitoring, power flow (dispatch), fault detection, etc. In this paper, an artificial neural network (ANN)-based framework for state estimation is proposed for an MBB, especially for unbalanced and low observable microgrids. To overcome the challenge of low observability in unbalanced systems, a concept of extended adjacent matrix is introduced to reduce the required number of measurements for state estimation. Addressing the challenges, a feed forward neural network (FNN) is utilized to enhance estimation accuracy and reliability with the reduced number of measurements. The proposed state estimation is validated through extensive simulations on a microgrid, which was achieved from the modified IEEE 34-bus distribution test feeder with multiple distributed energy resources (DERs) and demonstrated superior performance in estimation accuracy and low observability.

Choi, Jongchan↗
DOE OSTI2020

A Scalable and Distributed Algorithm for Managing Residential Demand Response Programs using Alternating Direction Method of Multipliers (ADMM)

For effective engagement of residential demand-side resources and to ensure efficient operation of distribution networks, we must overcome the challenges of controlling and coordinating residential components and devices at scale. To overcome this challenge, we present a distributed and scalable algorithm with a three-level hierarchical information exchange architecture for managing the residential demand response programs. First, a centralized optimization model is formulated to maximize community social welfare. Then, this centralized model is solved in a distributed manner with alternating direction method of multipliers (ADMM) by decomposing the original problem to utility-level and house-level problems. The information exchange between the different layers is limited to the primary residual (i.e., supply-demand mismatch), Lagrangian multipliers, and the total load of each house to protect each customer’s privacy. Simulation studies are performed on the IEEE 33 bus test system with 605 residential customers. The results demonstrate that the proposed approach can save customers’ electricity bills and reduce the peak load at the utility level without much affecting customers’ comfort and privacy. Finally, a quantitative comparison of the distributed and centralized algorithms shows the scalability advantage of the proposed ADMM-based approach, and it gives benchmarking results with achievable value for future research works.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

What are the Principles Controlling Biomimetic Heteropolymer Secondary Structure? (Final Technical Report)

The goal of the project was to develop improved theories to understand how nonbiological oligomers could be designed to cooperatively fold into 3D structures. These studies would lay the groundwork for materials made of such molecules, making it possible to create controlled and ordered materials for electron transport, efficient protein-like catalysts that work under extreme conditions, and sensors with highly-specific chemical responsiveness. Two different simulation thrusts were investigated, one focused on programs to identify stable low energy folded structures at a coarse-level of description of oligomers, and another to calculate thermodynamics of such oligomers. We used these theories to answer several specific questions about what properties of oligomers lead to cooperative transitions, and to identify how oligomer knots could serve as secondary structure elements. We also carried out significant collaborative investigation with Dr. Samuel Gellman (UW-Madison, National Academy of Sciences member) on stability for foldamers of interest to them. Only one of the experimentally tested foldamers stably folded, which was indicated by simulations as being the most likely to fold. Finally, we developed new theoretical descriptions of foldamers, showing how cooperativity was determined primarily by the entropy difference between the folded and unfolded state. The research did not answer all questions laid out in the original proposal but laid the groundwork for later efforts to design folded oligomers materials with high switchability.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Sensitivity Analysis of an Experimental Breeder Reactor II Fuel Assembly

The Experimental Breeder Reactor II (EBR-II) was selected as a candidate for a reactor physics benchmark evaluation due to its unique features and capabilities. In connection with the benchmark evaluation, the goal of this work was to perform a sensitivity analysis. The EBR-II was a sodium cooled, pool type, metal fueled, fast spectrum reactor. EBR-II operated from 1964 through 1994. It had a maximum heat output of 62.5 MW which could produce 20 MW of electricity. The design of EBR-II consisted of 637 hexagonal shaped, removable assemblies. The assemblies were divided into three regions: the core, an inner blanket, and an outer blanket. The core region contained assemblies with 91 fuel rods in each. The fuel was enriched uranium metal and was clad within stainless steel. The core region also included two safety assemblies and eight control assemblies with the difference being that these assemblies only contained 61 fuel rods. Additional assemblies in the core region included stainless steel dummies, half worth drivers, and experimental/instrumentation assemblies. The inner blanket region initially consisted of depleted uranium for fuel breeding, but was replaced, after proving the concept of breeding, with stainless steel reflectors for the purpose of irradiation. The outer blanket region consisted of depleted uranium. Figure 1 shows a cross section view of the EBR-II core layout. Each assembly was approximately 234 cm long with the uranium fuel section being 33 cm long. Above and below the fuel were stainless steel neutron reflectors. Additionally, orifices at the bottom of the assembly allowed for sodium coolant to flow upward, with larger holes in the subassemblies at the center where the greatest heat was produced. The core of EBR-II was enclosed by a stainless steel reactor vessel which was comprised of a radial shell, a lower grid plenum structure, and an upper reactor vessel cover that served as a neutron shield. Surrounding the reactor vessel radially were layers of graphite and borated graphite shielding. The entire assembly is submerged 3 meters in a pool of sodium. The pool contained up to 337,000 liters of sodium which was forced through the core and blanket regions using two centrifugal pumps that were also submerged within the sodium pool. The design of EBR-II included many safety features. The compatibility between the sodium and stainless steel prevented corrosion and therefore the risk of a radioactive release due to corrosion. The sodium coolant was kept at near atmospheric pressure which made any leak more controllable. The metal fuel, cladding, and coolant allowed for heat to be removed easily. The large volume of sodium allowed for heat to be absorbed without problems occurring. Lastly, the reactor would expand as temperature increased which would result in greater neutron leakage. This extra leakage caused a large negative temperature reactivity feedback which automatically shut the reactor down in case of an accident. EBR-II had unique capabilities which were proven during the Integral Fast Reactor initiative. This program lasted from 1984 to 1994. The goal of this program was to solve many of the perceived problems with nuclear energy through scientific means. The safety features of EBR-II were put to the test in April 1986. The tests simulated loss of reactor coolant flow and loss of a heat sink. For both of these tests, the reactor was at 100% power and was not scrammed. During both tests, the reactor shut itself down and achieved a safe steady state temperature all due to the design of the reactor. It was because of these tests that EBR-II was chosen for a reactor physics benchmark evaluation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Theoretical methods to design and test quantum simulators for the compact Abelian Higgs model

The lattice compact Abelian Higgs model is a nonperturbative regularized formulation of low-energy scalar quantum electrodynamics. In 1+1 dimensions, this model can be quantum simulated using a ladder-shaped optical lattice with Rydberg-dressed atoms. In this setup, one spatial dimension is used to carry the angular momentum of the quantum rotors. One can use truncations corresponding to spin-2 and spin-1 to build local Hilbert spaces associated with the links of the lattice. We argue that ladder-shaped configurable arrays of Rydberg atoms can be used for the same purpose. We make concrete proposals involving two and three Rydberg atoms to build one local spin-1 space (a qutrit). We show that the building blocks of the Hamiltonian calculations are models with one and two spins. We compare target and simulators using perturbative and numerical methods. The two-atom setup provides an easily controllable simulator of the one-spin model while the three-atom setup involves solving nonlinear equations. We discuss approximate methods to couple two spin-1 spaces. The article provides analytical and numerical tools necessary to design and build the proposed simulators with current technology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Stochastic finite volume method for uncertainty quantification of transient flow in gas pipeline networks

We develop a weakly intrusive framework to simulate the propagation of uncertainty in solutions of generic hyperbolic partial differential equation systems on graph-connected domains with nodal coupling and boundary conditions. The method is based on the Stochastic Finite Volume (SFV) approach and can be applied for uncertainty quantification (UQ) of the dynamical state of fluid flow over actuated transport networks. The numerical scheme has specific advantages for modeling intertemporal uncertainty in time-varying boundary parameters, which cannot be characterized by strict upper and lower (interval) bounds. We describe the scheme for a single pipe, and then formulate the controlled junction Riemann problem (JRP) that enables the extension to general network structures. In conclusion, we demonstrate the method's capabilities and performance characteristics using a standard benchmark test network.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Operator splitting method for simulation of dynamic flows in natural gas pipeline networks

Here, we develop an operator splitting method to simulate flows of isothermal compressible natural gas over transmission pipelines. The method solves a system of nonlinear hyperbolic partial differential equations (PDEs) of hydrodynamic type for mass flow and pressure on a metric graph, where turbulent losses of momentum are modeled by phenomenological Darcy-Weisbach friction. Mass flow balance is maintained through the boundary conditions at the network nodes, where natural gas is injected or withdrawn from the system. Gas flow through the network is controlled by compressors boosting pressure at the inlet of the adjoint pipe. Our operator splitting numerical scheme is unconditionally stable and it is second order accurate in space and time. The scheme is explicit, and it is formulated to work with general networks with loops. We test the scheme over range of regimes and network configurations, also comparing its performance with performance of two other state of the art implicit schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗