Search NASA⌕ Search

SEARCH · Search NASA

Results for “H CODES”

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

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

H-AMR: A New GPU-accelerated GRMHD Code for Exascale Computing with 3D Adaptive Mesh Refinement and Local Adaptive Time Stepping

General relativistic magnetohydrodynamic (GRMHD) simulations have revolutionized our understanding of black hole accretion. Here, we present a GPU-accelerated GRMHD code H-AMR with multifaceted optimizations that, collectively, accelerate computation by 2–5 orders of magnitude for a wide range of applications. First, it introduces a spherical grid with 3D adaptive mesh refinement that operates in each of the three dimensions independently. This allows us to circumvent the Courant condition near the polar singularity, which otherwise cripples high-resolution computational performance. Second, we demonstrate that local adaptive time stepping on a logarithmic spherical-polar grid accelerates computation by a factor of ≲10 compared to traditional hierarchical time-stepping approaches. Jointly, these unique features lead to an effective speed of ~10 9 zone cycles per second per node on 5400 NVIDIA V100 GPUs (i.e., 900 nodes of the OLCF Summit supercomputer). We illustrate H-AMR's computational performance by presenting the first GRMHD simulation of a tilted thin accretion disk threaded by a toroidal magnetic field around a rapidly spinning black hole. With an effective resolution of 13,440 × 4608 × 8092 cells and a total of ≲22 billion cells and ~0.65 × 10 8 time steps, it is among the largest astrophysical simulations ever performed. We find that frame dragging by the black hole tears up the disk into two independently precessing subdisks. The innermost subdisk rotation axis intermittently aligns with the black hole spin, demonstrating for the first time that such long-sought alignment is possible in the absence of large-scale poloidal magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Direct Numerical Simulation of Flow Over a Wall-Mounted Cube with the Nek5000 Spectral Element Code: DNS at Re h = 3900

The direct numerical simulation (DNS) of the canonical wall-mounted unit cube subjected to two distinct incident velocity profiles is performed at a Reynolds number where the bulk flow characteristics are known to become relatively Reynolds number insensitive. The aim of this work is to highlight the sensitivity of such bluff-body flows to mean shear and to provide a representative set of data for such flow scenarios where common turbulence modeling techniques often fail. Simple boundary conditions and a small domain are selected to allow for cost-effective and easy comparisons for model development purposes. In addition to mean velocity profiles, select turbulence statistics are presented in detail. Further, the basic efficacy of eddy viscosity models is evaluated and found to be adequate only for the shear stress components for the bluff-body flows examined here. Other failure mechanisms for RANS are proposed.

42 ENGINEERING↗
DOE OSTI2021

General relativistic MHD simulations of non-thermal flaring in Sagittarius A*

Sgr A* exhibits regular variability in its multiwavelength emission, including daily X-ray flares and roughly continuous near-infrared (NIR) flickering. The origin of this variability is still ambiguous since both inverse Compton and synchrotron emission are possible radiative mechanisms. The underlying particle distributions are also not well constrained, particularly the non-thermal contribution. In this work, we employ the GPU-accelerated general relativistic magnetohydrodynamics code H-AMR to perform a study of flare flux distributions, including the effect of particle acceleration for the first time in high-resolution 3D simulations of Sgr A*. For the particle acceleration, we use the general relativistic ray-tracing code bhoss to perform the radiative transfer, assuming a hybrid thermal+non-thermal electron energy distribution. We extract ~60 h light curves in the sub-millimetre, NIR and X-ray wavebands, and compare the power spectra and the cumulative flux distributions of the light curves to statistical descriptions for Sgr A* flares. Our results indicate that non-thermal populations of electrons arising from turbulence-driven reconnection in weakly magnetized accretion flows lead to moderate NIR and X-ray flares and reasonably describe the X-ray flux distribution while fulfilling multiwavelength flux constraints. These models exhibit high rms per cent amplitudes, $\gtrsim 150{{\ \rm per\ cent}}$ both in the NIR and the X-rays, with changes in the accretion rate driving the 230 GHz flux variability, in agreement with Sgr A* observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Deployment of the HFIRCON transport and depletion tool for plutonium-238 production studies

Irradiation of {sup 237}Np-bearing targets in Oak Ridge National Laboratory's (ORNL) High Flux Isotope Reactor (HFIR) results in the efficient production of {sup 238}Pu, which, in the form of heat source PuO{sub 2}, is used as a reliable power source for deep-space and planetary NASA missions. A technology demonstration subproject was initiated at ORNL in 2011 to develop and implement the technology required to establish a {sup 238}Pu supply chain. A systematic progression of NpO{sub 2}/Al cermet (20 vol.% NpO{sub 2}) activities to date has successfully demonstrated target fabrication, irradiation, and chemical recovery processes. Recent program tasks have included the development of the HFIRCON transport and depletion tool for efficient reactor physics analyses and the evaluation of increased NpO{sub 2} loadings (i.e., beyond 20 vol.%) and NpN-based targets. This paper documents the deployment of the HFIRCON code to assess various Np concentrations in NpO{sub 2}- and NpN-based targets in HFIR's inner small vertical experiment facilities. Results indicate that {sup 238}Pu production and quality can be enhanced with increased Np loadings; however, target conversion rates are reduced. The results recorded in this paper, thermal and material balance evaluations, and testing requirement planning will be used to determine whether increased NpO{sub 2} loadings or NpN-based targets will be further considered. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Linear source method of characteristics in HELIOS2

The Studsvik Scandpower code HELIOS2 is a fast lattice transport code that has gained wide use for reactor fuel cycle analysis in commercial, regulatory, and research environments. It contains two transport solvers - Collision Probabilities (CP) and Method of Characteristics (MoC) - with a general geometry ability and generalized depletion for a variety of fuel types, experimental facilities, and isotope production applications. The HELIOS2 lattice code has been extended with the development of a linear source MoC option. The linear source equations are derived within the scope of the HELIOS2 general geometry and non-cyclic ray tracing. The implementation was tested for fine and coarse spatial mesh models of the 2D C5G7 MOX computational benchmark, which provides a pure and extreme test of a transport solver. The linear source yields at least comparable accuracy with reduced computational resource requirements relative to the flat-source approximation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations

There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell

Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Because flows in stellar convection zones are highly turbulent, multidimensional hydrodynamic simulations are fundamental to accurately capture the physics of mixing processes. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dynamics in the convection zone and at its boundaries. In this study, we used our fully compressible S EVEN -L EAGUE H YDRO code to run detailed and highly resolved three-dimensional magnetohydrodynamic simulations of turbulent convection, dynamo amplification, and convective boundary mixing in a simplified setup whose stratification is similar to that of an oxygen-burning shell in a star with an initial mass of 25 M ⊙ . We find that the random stretching of magnetic field lines by fluid motions in the inertial range of the turbulent spectrum (i.e., a small-scale dynamo) naturally amplifies the seed field by several orders of magnitude in a few convective turnover timescales. During the subsequent saturated regime, the magnetic-to-kinetic energy ratio inside the convective shell reaches values as high as 0.33, and the average magnetic field strength is ∼10 10 G. Such strong fields efficiently suppress shear instabilities, which feed the turbulent cascade of kinetic energy, on a wide range of spatial scales. The resulting convective flows are characterized by thread-like structures that extend over a large fraction of the convective shell. The reduced flow speeds and the presence of magnetic fields with strengths up to 60% of the equipartition value at the upper convective boundary diminish the rate of mass entrainment from the stable layer by ≈20% as compared to the purely hydrodynamic case.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Tidal disruption discs formed and fed by stream–stream and stream–disc interactions in global GRHD simulations

When a star passes close to a supermassive black hole (BH), the BH’s tidal forces rip it apart into a thin stream, leading to a tidal disruption event (TDE). In this work, we study the post-disruption phase of TDEs in general relativistic hydrodynamics (GRHD) using our GPU-accelerated code h-amr. We carry out the first grid-based simulation of a deep-penetration TDE (β = 7) with realistic system parameters: a black hole-to-star mass ratio of 10 6 , a parabolic stellar trajectory, and a non-zero BH spin. We also carry out a simulation of a tilted TDE whose stellar orbit is inclined relative to the BH midplane. We show that for our aligned TDE, an accretion disc forms due to the dissipation of orbital energy with ~20 percent of the infalling material reaching the BH. The dissipation is initially dominated by violent self-intersections and later by stream–disc interactions near the pericentre. The self-intersections completely disrupt the incoming stream, resulting in five distinct self-intersection events separated by approximately 12 h and a flaring in the accretion rate. We also find that the disc is eccentric with mean eccentricity e ≈ 0.88. For our tilted TDE, we find only partial self-intersections due to nodal precession near pericentre. Although these partial intersections eject gas out of the orbital plane, an accretion disc still forms with a similar accreted fraction of the material to the aligned case. These results have important implications for disc formation in realistic tidal disruptions. For instance, the periodicity in accretion rate induced by the complete stream disruption may explain the flaring events from Swift J1644+57.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Preparation & shipping of ten neutron irradiated Eurofer97 steel variants to HFIR for SANS experiments

At ORNL, ten variants of Eurofer97 steel were irradiated in the high flux isotope reactor (HFIR) to ~2.94 – 3.24 dpa at 300 ± 30 °C, as part of the EUROfusion Lot-IV collaboration. The irradiations were performed in non-instrumented rabbit capsules in the flux trap region, which included SS-J3 type tensile samples and M4-CVN multi-notch bend bars. Their full set of microstructure and mechanical properties were previously reported in multiple previous publications. In this project, half-broken irradiated and nonirradiated SS-J3 tensile samples from ten alloys, code-named H, I, J, K, L, M, N, O, P and E, were prepared and shipped for facilitating small angle neutron scattering (SANS) experiments at the HFIR general purpose (GP)-SANS beam line. This report summarizes the completed tasks which included canister moves at the Irradiated Materials Examination and Testing (IMET) hot-cell facility to retrieve the samples from long-term storage, loading of the ten irradiated samples at IMET inside lead (Pb) piglets that were specifically provided by ORNL for SANS experiments and radiological shipment from IMET facility to the HFIR hot cells for performing SANS experiments. In addition to the irradiated samples, ten nonirradiated half-broken pieces from the Eurofer97 variants were also provided to HFIR for SANS.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nozzle Shocks, Disk Tearing, and Streamers Drive Rapid Accretion in 3D GRMHD Simulations of Warped Thin Disks

The angular momentum of gas feeding a black hole (BH) may be misaligned with respect to the BH spin, resulting in a tilted accretion disk. Rotation of the BH drags the surrounding spacetime, manifesting as Lense–Thirring torques that lead to disk precession and warping. We study these processes by simulating a thin (H/r = 0.02), highly tilted ($\tau$ = 65°) accretion disk around a rapidly rotating (a = 0.9375) BH at extremely high resolutions, which we performed using the general-relativistic magnetohydrodynamic code H-AMR. The disk becomes significantly warped and continuously tears into two individually precessing subdisks. We find that mass accretion rates far exceed the standard α-viscosity expectations. We identify two novel dissipation mechanisms specific to warped disks that are the main drivers of accretion, distinct from the local turbulent stresses that are usually thought to drive accretion. In particular, we identify extreme scale height oscillations that occur twice an orbit throughout our disk. When the scale height compresses, "nozzle" shocks form, dissipating orbital energy and driving accretion. Separate from this phenomenon, there is also extreme dissipation at the location of the tear. This leads to the formation of low-angular momentum "streamers" that rain down onto the inner subdisk, shocking it. The addition of low-angular momentum gas to the inner subdisk causes it to rapidly accrete, even when it is transiently aligned with the BH spin and thus unwarped. These mechanisms, if general, significantly modify the standard accretion paradigm. Additionally, they may drive structural changes on much shorter timescales than expected in α-disks, potentially explaining some of the extreme variability observed in active galactic nuclei.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Spatially and temporally resolved plasma parameter estimations of laser heated MagLIF relevant gas pipes at NIF

The ability to control laser pre-heat is an integral part of the inertial confinement fusion concept known as Magnetized Liner Inertial Fusion. This process is studied at the National Ignition Facility (NIF) where 4 of the 192 laser beams are propagated through a 1-cm long gas cell where they deposit >20 kJ of energy into the gaseous fuel via inverse bremsstrahlung absorption. This process ionizes the gas, producing a plasma that follows behind the laser front and expands over the radius of the cell. Emission from this plasma, as viewed by a gated x-ray detector, can be used to build spatially and temporally resolved estimations of the pre-heat plasma's density and temperature profiles. This can then be used to estimate the plasma pressure, internal energy, and radiation losses. Estimations show the evolution of the plasma in magnetized and unmagnetized gas cells filled with ambient temperature neopentane (C5H12) +1% Ar, as well as unmagnetized cryogenically cooled (32 K) deuterium +1% Ne filled targets. This analysis shows the effects of initial gas-fill density, composition, and axial magnetization on the time-dependent plasma parameters. Previously, these parameters at the NIF had not been experimentally characterized, and these estimations provided a potential new means of testing radiation magneto-hydrodynamic predictive capability models. Results in unmagnetized targets have strong agreement with simulations. However, in targets with a 19 T applied axial magnetic field, this method yields electron temperatures up to 100% hotter than those predicted by HYDRA codes.

Bremsstrahlung↗
DOE OSTI2026

GEOSH: Ideal Gas Chemical Equation of State

We present the framework and methodology for the new Los Alamos National Laboratory (LANL) G as chemical E quation O f S tate at H igher temperatures code (GEOSH) which aims to accurately model the behavior of chemically complex gaseous mixtures in equilibrium. Assuming the ideal gas approximation, GEOSH leverages the recursive nature of the Saha ionization and molecular equations in order to eliminate the molecular and ionic degrees of freedom, thereby reducing the problem size to the number of atomic species plus one for the free electrons if ions are included. This approach allows the chemical species, both molecular and ionic, of the mixture to be expressed in terms of the abundances of the elemental species. As a result, the GEOSH framework achieves a reduction in computational expense, increased processing speed, and the capability to efficiently model large-scale chemical networks. This report provides the necessary physical background and theoretical foundations for the GEOSH code, accompanied by benchmarking studies.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Linear simulations of a wide pedestal quiescent H-mode plasma with the extended-MHD code NIMROD

We present linear simulations of a Wide Pedestal QH (WPQH)-mode DIII-D plasma with the NIMROD code to address the role of ion multispecies collisionality on WPQH-modes. We use a full extended-MHD model that captures Alfvénic and magnetoacoustic waves, including two-fluid, ion gyroviscosity and cross heat flux. We show that two-fluid and ion gyroviscous effects destabilize electron-directed peeling–ballooning modes in the pedestal. The linear growth rates of the modes decrease when carbon (C) is added to a deuterium plasma because the multispecies collisional effects increase the plasma resistivity. When replacing C with tungsten (W), the multispecies collisional effects increase the resistivity further, as well as the stabilizing effect on the instabilities. The resistivity impact on mode stability is reversed when both two-fluid and ion gyroviscous effects are ignored, which is consistent with previous works, proving that multispecies collisional effects, together with two-fluid and ion gyroviscosity, should be considered to address the stability of plasmas with significant impurity content.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Determining Unabated Airborne Radionuclide Emissions Monitoring Requirements Using Inventory-Based Methods

Compliance with the airborne radionuclide emission monitoring requirements in the National Emission Standards for Hazardous Air Pollutants (NESHAP; Title 40 of the U.S. Code of Federal Regulations Part 61, Subpart H) and State requirements in Washington Administrative Code 246–247: Radiation Protection – Air Emissions and 173-480: Ambient Air Quality Standards and Emission Limits for Radionuclides were evaluated for Pacific Northwest National Laboratory (PNNL) operations. Additional guidance may be found in the U.S. Department of Energy Handbook, Environmental Radiological Effluent Monitoring and Environmental Surveillance. To meet regulatory requirements, reviews of planned or proposed projects and activities provide the basis for implementing necessary monitoring adjustments or for implementing changes to projects and activities in a timely manner. Potential unabated off-site doses were evaluated for emission locations managed by PNNL and licensed to the Department of Energy. These locations were at facilities in Richland, Washington (i.e., the Hanford Site and PNNL–Richland campus) and in Sequim, Washington, (PNNL–Sequim campus). This report describes the inventory-based methods and provides the results for the NESHAP assessment performed in January 2024 for calendar year ending 2023.

40 CFR 61 Subpart H↗
DOE OSTI2022

Modeling study of divertor particle and heat flux asymmetries for EAST H-mode discharges

The BOUT++ transport code is run to study the effects of plasma drifts on the divertor out-in asymmetries (DOIAs) of particle and heat fluxes and their decay widths for EAST lower single null H-mode discharges. The diamagnetic drift seems to have no effects on the DOIAs of total particle and heat fluxes due to its divergence-free nature. However, it could significantly increase the DOIAs of peak particle and heat fluxes and the flux decay widths. The E × B drift is found to induce a large plasma flow to the divertor region, enhancing the DOIAs of both total and peak particle and heat fluxes and the flux decay widths. Both the radial and poloidal components of the E × B drift are necessary in increasing the DOIAs, however, the radial E × B drift seems to play a more important role. The effects on the DOIAs caused by both diamagnetic and E × B drifts are reversed with the reverse of toroidal magnetic field. The heat flux decay width λ q and spreading width S q are important physical and engineering parameters for the divertors and could be obtained by fitting the heat flux profiles at divertor targets. The λ q at the outer target from the simulation case with all drifts could well match with the multi-machine scaling proposed by Eich and the DOIA of λ q is in reasonable agreement with the scaling proposed by Goldston.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Correcting for Selection Biases in the Determination of the Hubble Constant from Time-Delay Cosmography

The time delay between multiple images of strongly lensed quasars has been used to infer the Hubble constant. The primary systematic uncertainty for time-delay cosmography is the mass-sheet transform (MST), which preserves the lensing observables while altering the inferred ⁠H 0 . The TDCOSMO collaboration used velocity dispersion measurements of lensed quasars and lensed galaxies to infer that mass sheets are present, which decrease the inferred H 0 by 8 per cent. Here, we test the assumption that the density profiles of galaxy–galaxy and galaxy–quasar lenses are the same. We use a composite star-plus-dark-matter mass profile for the parent deflector population and model the selection function for galaxy–galaxy and galaxy–quasar lenses. We find that a power-law density profile with an MST is a good approximation to a two-component mass profile around the Einstein radius, but we find that galaxy–galaxy lenses have systematically higher mass-sheet components than galaxy–quasar lenses. For individual systems, λ int correlates with the ratio of the half-light radius and Einstein radius of the lens. By propagating these results through the TDCOSMO hierarchical inference code, we find that H 0 is lowered by a further 3 per cent. Using a more recent measurement of velocity dispersions and our fiducial model for selection biases, we infer H 0 = 66 ± 4 (stat) ± 1 (model sys) ± 2 (measurement sys) km s -1 Mpc -1 for the TDCOSMO plus SLACS data set. The first residual systematic error is due to plausible alternative choices in modelling the selection function, and the second is an estimate of the remaining systematic error in the measurement of velocity dispersions for SLACS lenses. Accurate time-delay cosmography requires precise velocity dispersion measurements and accurate calibration of selection biases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Field–potential finite-difference time-domain (FiPo FDTD) technique for computational electromagnetics

Modeling light–matter interactions at the nanoscale requires accurate handling of coupled quantum and electromagnetic systems. This coupling requires information about the electric scalar potential Φ and the magnetic vector potential A, which are not typically calculated in standard computational electromagnetics implementations. To that end, we have developed a field–potential finite-difference time-domain (FiPo FDTD) algorithm, which solves a set of first-order equations for Φ and A alongside equations for the electric and magnetic fields E and H. The FiPo Basic code is essentially conventional FDTD, but with an added module that calculates the potentials. The FiPo Hybrid code self-consistently calculates both fields and potentials and is particularly suitable for coupling with quantum electronic transport solvers because it can be sourced by the potentials themselves. To terminate the domain and mimic infinite space, we have derived and implemented a convolutional perfectly matched layer (CPML) absorbing boundary condition for FiPo FDTD whose performance is on par with state-of-the-art CPMLs for standard FDTD. We present FiPo simulation results on several example systems.

Avazpour, L. [University of Wisconsin-Madison, WI ↗