Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Power Iteration”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Real-space imaging of periodic nanotextures in thin films via phasing of diffraction data

New properties and exotic quantum phenomena can form due to periodic nanotextures, including Moire patterns, ferroic domains, and topologically protected magnetization and polarization textures. Despite the availability of powerful tools to characterize the atomic crystal structure, the visualization of nanoscale strain-modulated structural motifs remains challenging. Here, we develop nondestructive real-space imaging of periodic lattice distortions in thin epitaxial films and report an emergent periodic nanotexture in a Mott insulator. Specifically, we combine iterative phase retrieval with unsupervised machine learning to invert the diffuse scattering pattern from conventional X-ray reciprocal-space maps into real-space images of crystalline displacements. Our imaging in PbTiO 3 /SrTiO 3 superlattices exhibiting checkerboard strain modulation substantiates published phase-field model calculations. Furthermore, the imaging of biaxially strained Mott insulator Ca 2 RuO 4 reveals a strain-induced nanotexture comprised of nanometer-thin metallic-structure wires separated by nanometer-thin Mott-insulating-structure walls, as confirmed by cryogenic scanning transmission electron microscopy (cryo-STEM). The nanotexture in Ca 2 RuO 4 film is induced by the metal-to-insulator transition and has not been reported in bulk crystals. We expect the phasing of diffuse X-ray scattering from thin crystalline films in combination with cryo-STEM to open a powerful avenue for discovering, visualizing, and quantifying the periodic strain-modulated structures in quantum materials.

Science & Technology - Other Topics↗
DOE OSTI2025

Thermal-Fluid and Thermal-Structural Response of the T-Tube Modular Divertor to Spatiotemporally Varying Heat Loads

Tungsten (W) is the leading candidate for divertor target plates because of its high melting point (>3000°C), thermal conductivity, and ultimate tensile stress. While W and its alloys are the only solid materials that can survive the high heat fluxes incident on the divertor, W’s low-ductility high ductile-to-brittle transition temperature of ~600°C and relatively low recrystallization temperature (RT) of ~1300°C pose structural (among other) challenges. The objective of this work is to estimate the thermal-fluid and thermal-structural performance of the helium (He)-cooled T-tube divertor, which was originally developed by the Advanced Reactor Innovation and Evaluation Study (ARIES) using numerical simulations. Here, predictions of temperature distributions across the plasma-facing structural component and surface pressures from computational fluid dynamics simulations are used to determine stress distributions using commercial structural finite element modeling software over a range of fusion-relevant conditions. The maximum allowable incident heat fluxes are determined based on the temperature limits imposed by the ITER elastic Structural Design Criteria for In-vessel Components (SDC-IC) and the maximum RT over a range of He mass flow rates and presented in the form of performance design charts. Our recent work found that thermal- structural criteria accounting for the low ductility of W in a finger-type modular divertor constrain the maximum incident heat fluxes to values well below the ITER specifications, and those based on considering only the RT demonstrate that integrated thermal-fluid and elastic structural performance evaluation are required for accurate assessment of divertor performance. This novel analysis of the T-tube considers how nonuniform and transient incident heat fluxes affect its thermal-fluid and thermal-structural performance, as well as the effect of volumetric heating, which can be as great as 27% of the power incident on the divertor surface. The W tile of the T-tube, with its relatively large plasma-facing area of ~15 cm 2 , will likely experience significant spatial variations in incident heat flux. This work therefore assesses whether steady-state incident heat flux profiles with a peak of 10 MW/m 2 and maximum heat flux gradients of 200 MW/m 2 per m exceed the structural limits imposed by the ITER elastic SDC-IC and the maximum RT over a range of fusion-relevant conditions. The effect of transient heat fluxes typical of plasma detachment and reattachment from the target plate due, for example, to gas injection are also evaluated

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Feasibility of raised inner strike point equilibria scenario in ITER for detritiation from beryllium co-deposits

Abstract In ITER, tritium retention primarily occurs through co-deposition with beryllium. To avoid exceeding the strict tritium inventory limit, efficient tritium recovery techniques are essential. Baking is the ITER baseline for tritium recovery, but its effectiveness in removing tritium from thick beryllium layers is limited. A raised strike point scenario is considered an alternative method for removing tritium from the ITER inner vertical divertor target by heating components via plasma flux. This paper presents SOLPS-ITER code simulations conducted under various conditions, assessing the divertor performance and tritium outgassing of the raised strike point scenario. As the strike point is raised, recycled neutrals are not efficiently baffled by the dome and scrape-off layer, significantly changing the neutral trajectory and ionization source distribution. This improves detachment accessibility but worsens core-edge compatibility compared to the baseline scenario. However, in the partially detached condition, the impact of raising the strike point, perpendicular transport, and q 95 on target heat flux is not significant, as it primarily scales with the input power. Target heat flux is translated to target surface temperature using a simplified heat transfer model that considers the 3D target monoblock geometry and active cooling condition, excluding Be layer thermal properties. For partially detached divertor conditions, the bulk tungsten monoblock surface temperature remains below the baking temperature, which is insufficient for efficient tritium outgassing under the actively cooled ITER divertor condition. However, considering the potential thermal contact resistance between the beryllium and tungsten layers, which may significantly impact temperature distribution, the temperature of the beryllium layer can be raised to a level sufficient for efficient tritium outgassing. Therefore, the raised strike point scenario can be considered as an alternative in-vessel tritium removal technique.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Distributed Model Identification Algorithm for Multi-Agent Systems: Preprint

In this study, we investigate agent-based approach for system model identification with emphasis on power distribution system applications. Departing from conventional practices of relying on historical data for offline model identification, we adopt online update approach utilizing real-time data by employing the latest data points for gradient computation. This methodology offers advantages including a large reduction in the communication network's bandwidth requirements by minimizing the data exchanged at each iteration and enabling the model to adapt in real-time to disturbances. Furthermore, we extend our model identification process from linear frameworks to more complex non-linear convex models. This extension is validated through numerical studies demonstrating improved control performance for a synthetic IEEE test case.

data-driven control↗
DOE OSTI2022

Recent progress in L-H transition studies at JET: Tritium, Helium, Hydrogen and Deuterium

In this report we present an overview of results from a series of L-H transition experiments undertaken at JET since the installation of the ITER-like-wall (JET-ILW), with Beryllium wall tiles and a Tungsten divertor. Tritium, Helium and Deuterium plasmas have been investigated. Initial results in Tritium show ohmic L-H transitions at low density and the power threshold for the L-H transition (P LH ) is lower in Tritium plasmas than in Deuterium ones at low densities, while we still lack contrasted data to provide a scaling at high densities. In Helium plasmas there is a notable shift of the density at which the power threshold is minimum (n e,min ) to higher values relative to Deuterium and Hydrogen references. Above n e,min (He) the L-H power threshold at high densities is similar for D and He plasmas. Transport modelling in slab geometry shows that in Helium neoclassical transport competes with interchange-driven transport, unlike in Hydrogen isotopes. Measurements of the radial electric field in Deuterium plasmas show that E r shear is not a good indicator of proximity to the L-H transition. Transport analysis of ion heat flux in Deuterium plasmas show a non-linearity as density is decreased below n e,min . Lastly, a regression of the JET-ILW Deuterium data is compared to the 2008 ITPA scaling law.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental study on the role of the target electron temperature as a key parameter linking recycling to plasma performance in JET-ILW*

Changes in global and edge plasma parameters (H98(y,2), dimensionless collisionality ν*, core density peaking, separatrix density ne,sep) with variations in the D2 fueling rate and divertor configuration are unified into a single trend when mapped to $\langle$Te,ot$\rangle$, the spatially averaged spectroscopically derived outer target electron temperature. Dedicated JET with the ITER-like wall (JET-ILW) experiments in combination with an extended JET-ILW database of unseeded low-triangularity H-mode plasmas spanning a wide range of D2 fueling rates, Ip, Bt and heating power have demonstrated the importance of $\langle$Te,ot$\rangle$ as a key physics parameter linking the recycling particle source and detachment with plasma performance. The remarkably robust H98(y,2) trend with $\langle$Te,ot$\rangle$ is connected to a strong inverse correlation between $\langle$Te,ot$\rangle$, ne,sep and ν*, thus directly linking changes in the divertor recycling moderated by $\langle$Te,ot$\rangle$ with the previously established relationship between ν*, core density peaking and core pressure resulting in a degradation in core plasma performance with decreasing $\langle$Te,ot$\rangle$ (increasing ν*). Additionally, a strong inverse correlation between the separatrix to pedestal density ratio, ne,sep/ne,ped, and $\langle$Te,ot$\rangle$ is also established, with the rise in ne,sep/ne,ped saturating at $\langle$Te,ot$\rangle$ > 10 eV. A strong reduction in H98(y,2) is observed as $\langle$Te,ot$\rangle$ is driven from 30 to 10 eV via additional D2 gas fueling, while the divertor remains attached. Consequently, the pronounced performance degradation in attached divertor conditions has implications for impurity seeding radiative divertor scenarios, in which H98(y,2) is already low (~0.7) before impurities are injected into the plasma since moderate gas fueling rates are required to promote high divertor neutral pressure. A favorable pedestal pressure, pe,ped, dependence on Ip has also been observed, with an overall increase in pe,ped at Ip = 3.4 MA as $\langle$Te,ot$\rangle$ is driven down from attached to high-recycling divertor conditions. In contrast, pe,ped is reduced with decreasing $\langle$Te,ot$\rangle$ in the lower Ip branches. Further work is needed to (i) clarify the potential role of edge opacity on the observed favorable pedestal pressure Ip scaling; as well as to (ii) project the global and edge plasma performance trends with $\langle$Te,ot$\rangle$ to reactor-scale devices to improve predictive capability of the coupling between recycling and confined plasma fueling in what are foreseen to be more opaque edge plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Gyrokinetic prediction of core tungsten peaking in a WEST plasma with nitrogen impurities

Tungsten peaking is predicted in the core of a WEST plasma with total-f gyrokinetic simulations, including both collisional and turbulent transport. This prediction is validated with a synthetic diagnostic of the bolometry. Although nitrogen impurities are shown to reduce the neoclassical peaking of tungsten on-axis, the overall tungsten peaking increases when nitrogen impurities are present, as they reduce the turbulence screening off-axis. This finding is important for the plasma current ramp-up phase of ITER, where light impurities seeding will be desirable to achieve low temperatures at the plasma-facing components and reduce tungsten sputtering. It provides further argument for applying early ECRH heating to maintain margins on the core power balance. The neoclassical peaking factor is cross-verified between XGC and FACIT. The heat flux at separatrix and the heat load width are modeled by XGC and compared to WEST data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Filtered Backprojection in Compton Imaging using a Spherical Harmonic Wiener Filter with Pixelated CdZnTe

Filtered backprojection is an image reconstruction technique for Compton imaging that provides reasonably high resolution at much lower computational costs when compared to iterative methods. Here, this work applies a Wiener filter that has been derived for spherical harmonics on Compton imaging using the OrionUM pixelated CdZnTe imaging-spectrometer. To regularize the filter, an investigation is made into the power spectral density of the signal and noise to develop an appropriate spectral signal-to-noise ratio model for the restoration process. Experimental measurements were conducted with two 228Th sources placed 30° apart. The resulting filtered image of the two sources have an average full-width-at-half-maximum (FWHM) of 9.8° or 7.5° when using a mean squared error and structural similarity optimization approach respectively; an improvement from the 29.0° FWHM image when using simple backprojection.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Atomic electron tomography in three and four dimensions

Atomic electron tomography (AET) has become a powerful tool for atomic-scale structural characterization in three and four dimensions. It provides the ability to correlate structures and properties of materials at the single-atom level. With recent advances in data acquisition methods, iterative three-dimensional (3D) reconstruction algorithms, and post-processing methods, AET can now determine 3D atomic coordinates and chemical species with sub-Angstrom precision, and reveal their atomic-scale time evolution during dynamical processes. In this work, we review the recent experimental and algorithmic developments of AET and highlight several groundbreaking experiments, which include pinpointing the 3D atom positions and chemical order/disorder in technologically relevant materials and capturing how atoms rearrange during early nucleation at four-dimensional atomic resolution.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS1985

Circuit analysis method for thin-film solar cell modules

The design of a thin-film solar cell module is dependent on the probability of occurrence of pinhole shunt defects. Using known or assumed defect density data, dichotomous population statistics can be used to calculate the number of defects expected in a module. Probability theory is then used to assign the defective cells to individual strings in a selected series-parallel circuit design. Iterative numerical calculation is used to calcuate I-V curves using cell test values or assumed defective cell values as inputs. Good and shunted cell I-V curves are added to determine the module output power and I-V curve. Different levels of shunt resistance can be selected to model different defect levels.

Burger, D. R.↗
DOE OSTI2022

Methods and systems for evaluating and improving distribution-grid observability

This document describes systems and techniques for evaluating and improving distribution-grid observability. These systems and techniques allow engineers to quantify the observability of a distribution grid, which represents an ability to combine actual measurements and various types of computations (e.g., analytics, estimators, forecasters), from a system model. Distribution engineers can also identify islands of observability where operating parameters, including voltages, currents, and power flows, can be determined from available sensor readings. By exclusion, distribution engineers can similarly identify areas of the distribution grid with observability deficiencies that may require additional instrumentation to maintain proper operation. Distribution engineers, using an iterative or automated process, can determine the observability of the system model with new or relocated sensors to generate a sensor allocation plan. The sensor allocation plan can indicate the number and location of sensors to either maximize observability for a fixed sensor cost or minimize sensor cost for predetermined observability.

Reiman, Anthony P.↗
NASA NTRS2008

Post-Landing Orion Crew Survival in Warm Ocean Areas: A Case Study in Iterative Environmental Design

The Orion crew module (CM) is being designed to perform survivable land and water landings. There are many issues associated with post-landing crew survival. In general, the most challenging of the realistic Orion landing scenarios from an environmental control standpoint is the off-nominal water landing. Available power and other consumables will be very limited after landing, and it may not be possible to provide full environmental control within the crew cabin for very long after splashdown. Given the bulk and thermal insulation characteristics of the crew-worn pressure suits, landing in a warm tropical ocean area would pose a risk to crew survival from elevated core body temperatures, if for some reason the crewmembers were not able to remove their suits and/or exit the vehicle. This paper summarizes the analyses performed and conclusions reached regarding post-landing crew survival following a water landing, from the standpoint of the crew s core body temperatures.

Rains, George E.↗
DOE OSTI2022

Dynamic energy system modeling using hybrid physics-based and machine learning encoder–decoder models

Three model configurations are presented for multi-step time series predictions of the heat absorbed by the water and steam in a thermal power plant. The models predict over horizons of 2, 4, and 6 steps into the future, where each step is a 5-minute increment. The evaluated models are a pure machine learning model, a novel hybrid machine learning and physics-based model, and the hybrid model with an incomplete dataset. The hybrid model deconstructs the machine learning into individual boiler heat absorption units: economizer, water wall, superheater, and reheater. Each configuration uses a gated recurrent unit (GRU) or a GRU-based encoder–decoder as the deep learning architecture. Mean squared error is used to evaluate the models compared to target values. The encoder–decoder architecture is over 11% more accurate than the GRU only models. The hybrid model with the incomplete dataset highlights the importance of the manipulated variables to the system. The hybrid model, compared to the pure machine learning model, is over 10% more accurate on average over 20 iterations of each model. Automatic differentiation is applied to the hybrid model to perform a local sensitivity analysis to identify the most impactful of the 72 manipulated variables on the heat absorbed in the boiler. The models and sensitivity analyses are used in a discussion about optimizing the thermal power plant.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2004

Research Capabilities for Oil-Free Turbomachinery Expanded by New Rotordynamic Simulator Facility

A new test rig has been developed for simulating high-speed turbomachinery shafting using Oil-Free foil air bearing technology. Foil air journal bearings are self-acting hydrodynamic bearings with a flexible inner sleeve surface using air as the lubricant. These bearings have been used in turbomachinery, primarily air cycle machines, for the past four decades to eliminate the need for oil lubrication. More recently, interest has been growing in applying foil bearings to aircraft gas turbine engines. They offer potential improvements in efficiency and power density, decreased maintenance costs, and other secondary benefits. The goal of applying foil air bearings to aircraft gas turbine engines prompted the fabrication of this test rig. The facility enables bearing designers to test potential bearing designs with shafts that simulate the rotating components of a target engine without the high cost of building actual flight hardware. The data collected from this rig can be used to make changes to the shaft and bearings in subsequent design iterations. The rest of this article describes the new test rig and demonstrates some of its capabilities with an initial simulated shaft system. The test rig has two support structures, each housing a foil air journal bearing. The structures are designed to accept any size foil journal bearing smaller than 63 mm (2.5 in.) in diameter. The bearing support structures are mounted to a 91- by 152-cm (3- by 5-ft) table and can be separated by as much as 122 cm (4 ft) and as little as 20 cm (8 in.) to accommodate a wide range of shaft sizes. In the initial configuration, a 9.5-cm (3.75-in.) impulse air turbine drives the test shaft. The impulse turbine, as well as virtually any number of "dummy" compressor and turbine disks, can be mounted on the shaft inboard or outboard of the bearings. This flexibility allows researchers to simulate various engine shaft configurations. The bearing support structures include a unique bearing mounting fixture that rotates to accommodate a laserbased alignment system. This can measure the misalignment of the bearing centers in each of 2 translational degrees of freedom and 2 rotational degrees of freedom. In the initial configuration, with roughly a 30.5-cm- (12-in.-) long shaft, two simulated aerocomponent disks, and two 50.8-cm (2-in.) foil journal bearings, the rig can operate at 65,000 rpm at room temperature. The test facility can measure shaft displacements in both the vertical and horizontal directions at each bearing location. Horizontal and vertical structural vibrations are monitored using accelerometers mounted on the bearing support structures. This information is used to determine system rotordynamic response, including critical speeds, mode shapes, orbit size and shape, and potentially the onset of instabilities. Bearing torque can be monitored as well to predict the power loss in the foil bearings. All of this information is fed back and forth between NASA and the foil bearing designers in an iterative fashion to converge on a final bearing and shaft design for a given engine application. In addition to its application development capabilities, the test rig offers several unique capabilities for basic bearing research. Using the laser alignment system mentioned earlier, the facility will be used to map foil air journal bearing performance. A known misalignment of increasing severity will be induced to determine the sensitivity of foil bearings to misalignment. Other future plans include oil-free integral starter generator testing and development, and dynamic load testing of foil journal bearings.

Howard, Samuel A.↗
DOE OSTI2020

Iterative and Self-Consistent Estimation of Degradation and Soiling Loss in PV Systems - A Case Study

In many parts of the world, soiling is a significant loss mechanism for PV energy production. Methods for accurately quantifying historical and current soiling losses are needed. Furthermore, state-of-the-art methods to estimate degradation rates based on production data do not handle soiling in a satisfactory manner. In order to accurately estimate either of these quantities, it is in many cases helpful to consider their combined effect on PV system performance. In this work we make use of the Combined Degradation and Soiling (CODS) algorithm, which iteratively decomposes time series of daily normalized energy yield into four components: a soiling ratio, a degradation trend, seasonal variations, and noise. We demonstrate the use of the model on field performance data from a utility-scale power plant in a high-soiling environment in the Middle East. We apply the algorithm to the temperature-corrected performance index of the DC power on inverter level, quantifying how soiling varies over time and through the plant. The irradiation-weighted average soiling loss in this plant is 2.3 %, with a median soiling rate of -0.086 %/day. Between the inverters in the site in study, estimated degradation rates range between -0.86 %/year and -0.15 %/year, with a mean rate of -0.48 %/year.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2024

Endogenizing Probabilistic Resource Adequacy Risks in Deterministic Capacity Expansion Models

In this work, we demonstrate how power system capacity expansion models can understate the stochastic effects of thermal outages when considering resource availabilities on an hourly expected value basis, yielding system designs with multiple orders of magnitude more shortfall risk than stated adequacy targets. We develop a novel approximation approach to efficiently endogenize awareness of this risk in a deterministic, linear capacity expansion framework. We compare this approach to exogenous tuning of an energy reserve margin, the leading alternative method to compensate for unmodeled probabilistic shortfall risk. Empirical results from a test system show that the new endogenous method cost-effectively meets all regional reliability targets with a single optimization solve, and produces a near-identical system design as the incumbent method without the need for repeated re-optimizations to find an appropriate reserve level. The endogenous method may also use iterative re-optimizations to further improve solution quality, although these incremental benefits were modest in the system studied.

capacity expansion modeling↗
DOE OSTI2023

Physics of the thermal quench in MST tokamak plasmas (Final Technical Report)

This project at the University of Wisconsin-Madison targeted measurements in tokamak plasmas in the Madison Symmetric Torus (MST) to help elucidate the physics of the tokamak disruption. We intended to focus on the thermal quench, during which energy is lost from the plasma in a very short time. Understanding and ultimately mitigating the thermal quench are among the highest priority research thrusts for ITER, a large international fusion experiment presently under construction. The measurements in MST exploited an advanced set of diagnostics for internal measurements with high time resolution, targeting linear and nonlinear magnetohydrodynamic (MHD) activity as well as the electron temperature and density profiles. We also exploited other features of the MST device, including very flexible power supplies for the toroidal magnetic field and toroidal plasma current that are capable of producing almost arbitrary waveforms. While the work funded by this project was purely experimental, we continued and expanded an existing collaboration with UW-Madison experts in nonlinear MHD computational modeling of tokamak plasmas. That modeling incorporated realistic, experimentally measured parameters. The results of the modeling were anticipated to deepen considerably our understanding of the experimental measurements, and in turn, the experimental measurements were anticipated to help validate the modeling results. The science of transients such as disruptions is one of the highest priority research areas in the U.S. and international fusion programs. The work carried out in this project was targeted to help advance the science of tokamak disruptions which will in turn contribute to the success of ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Performance and Energy Study of GPU-Resident Preconditioners for Conjugate Gradient Solvers: In the Context of Existing and Novel Approaches

Optimizing a particular subprogram out of the set of Basic (sparse) Linear Algebra Subprograms (BLAS) for a given architecture is a common topic of research. In applications, however, these BLAS functions rarely appear in isolation; usually, many of them are used together, in various combinations and with varying inputs. As the need to solve a large, sparse linear system is ubiquitous throughout HPC applications, linear solvers constitute a realistic, sufficiently complex and well-defined representative use case for composite BLAS routines. To this end, based on a representative set of matrices drawn from a diverse set of fields, we present a framework to study, from the performance and energy perspective, the efficacy of GPU- resident parallel Conjugate Gradient (CG) linear solver with different preconditioner options, including Gauss-Seidel, Jacobi, and incomplete Cholesky. We also propose a novel GPU-based preconditioner, in which the triangular solves are approximated by an iterative process. The development of this preconditioner was motivated by solving large graph Laplacian linear systems, for which the existing preconditioners either perform slow on GPU-based platforms or are not applicable. We compare the performance of these preconditioners on different hardware accelerator architectures, i.e., AMD MI250X, MI100, Nvidia A100, V100, and Jetson. Our experiments reveal performance trade-offs and provide information on how to select the best strategy for the given linear system, dictated by its properties, and the platform of interest. We demonstrate the application of our novel preconditioner for solving CG and graph Laplacian systems. Overall, the framework can be utilized as a benchmark to guide informed decisions in choosing a specific preconditioner, i.e., whether it is better to rely on the performance of a triangular solver or on the performance of sparse matrix-vector product. Finally, by considering power consumption to solve the linear systems, we report the energy footprint for the solvers.

Preconditioned Conjugate Gradient, GPUs, iterative↗