SEARCH · Engineering Papers

Results for “numerical continuation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2023

SCALE 6.2 Lattice Physics Performance Assessment

The US Nuclear Regulatory Commission relies on the lattice physics analysis capabilities of the SCALE code system to perform confirmatory licensing analyses. Either SCALE lattice physics code—TRITON/NEWT or Polaris—can be used to generate cross section data used by the PARCS nodal core simulator for full-core neutronics calculations. This report presents an assessment of the accuracy of SCALE lattice physics codes for preparation of lattice physics data that are used to support simulator codes such as the NRC’s PARCS, for UO 2 -mixed oxide (MOX)/Zr fueled light water reactor (LWR) analyses. Due to the nature of lattice physics calculations, critical reactor experiment benchmarks cannot be modeled in explicit detail in a lattice physics code. However, this limitation does not mean that these measurement data are not usable for lattice physics studies. Therefore, either geometry approximations or axial buckling must be implemented to determine the critical water height. These modeling limitations have led to development of a three-phase assessment strategy. In the first phase, selected critical experiment benchmarks are modeled using the SCALE 3D continuous-energy (CE) Monte Carlo (MC) code KENO, which is the most rigorous neutron transport method available in SCALE, with no approximations in the spatial, angle, or energy treatments. Biases and statistical uncertainties in quantities of interest such as k eff and pin power distributions are determined by comparing CE KENO results to experimental data. This first phase in the accuracy assessment is to demonstrate and establish the use of CE KENO as a reference solution for the second phase. In the second phase, 14 numerical test suites are used to compare SCALE lattice physics calculations with CE KENO as a reference solution. In the third phase, SCALE lattice physics depletion calculations are performed, and the spent fuel isotopic results are compared with available radiochemical assay measurements. Isotopic measurement comparisons provide quantitative assessment of isotopic density distribution predictions with the depletion models in SCALE lattice physics codes. This report documents results for all test suites. The assessment was performed using standard production techniques unless otherwise noted. Both TRITON/NEWT and Polaris exhibited acceptable accuracy for most test cases. For the few test cases in which acceptable accuracy criteria were not met, further code and data development are planned. The computer codes used in this assessment are as follows: SCALE 6.2 rev19189 (pre-release of 6.2.1) was used for CE KENO, TRITON/NEWT and Polaris (PWR only) calculations; SCALE 6.2.2 was used for Polaris calculations to address several updates in support of support boiling water reactor (BWR) lattice geometries. For this work, there is no difference between SCALE 6.2.1 and SCALE 6.2.2 beyond the inclusion of the Polaris BWR analysis capability; The ENDF/B-VII.1 continuous energy and 252 group libraries deployed in SCALE 6.2 were generated with AMPX 6.2 that is distributed with SCALE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Electric field control of magnetism: multiferroics and magnetoelectrics

This article is written on behalf of a large number of colleagues, collaborators, and researchers in the field of complex oxides as well as current and former students and postdocs who continue to enable and undertake cutting-edge research in the field of multiferroics, magnetoelectrics, and the pursuit of electric-field control of magnetism. What we present is something that is extremely exciting from both a fundamental science and applications perspective and has the potential to revolutionize our world. Needless to say, to realize this potential will require numerous new innovations, both in the fundamental science arena as well as translating these scientific discoveries into real applications. Thus, this article will attempt to bridge the gap between fundamental condensed-matter physics and the actual manifestations of the physical concepts into real-life applications. Here, we hope this article will help spur more translational research within the broad materials physics community.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗

DOE OSTI2023

Predictive scale-bridging simulations through active learning

Throughout computational science, there is a growing need to utilize the continual improvements in raw computational horsepower to achieve greater physical fidelity through scale-bridging over brute-force increases in the number of mesh elements. For instance, quantitative predictions of transport in nanoporous media, critical to hydrocarbon extraction from tight shale formations, are impossible without accounting for molecular-level interactions. Similarly, inertial confinement fusion simulations rely on numerical diffusion to simulate molecular effects such as non-local transport and mixing without truly accounting for molecular interactions. With these two disparate applications in mind, we develop a novel capability which uses an active learning approach to optimize the use of local fine-scale simulations for informing coarse-scale hydrodynamics. Our approach addresses three challenges: forecasting continuum coarse-scale trajectory to speculatively execute new fine-scale molecular dynamics calculations, dynamically updating coarse-scale from fine-scale calculations, and quantifying uncertainty in neural network models.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Pulse-by-pulse transient thermal deformation in crystal optics under high-repetition-rate FEL

Time-domain modeling of the thermal deformation of crystal optics can help define acceptable operational ranges across the pulse-energy repetition-rate phase space. In this paper, we have studied the transient thermal deformation of a water-cooled diamond crystal for a cavity-based X-ray free-electron laser (CBXFEL), either an X-ray free-electron laser oscillator (XFELO) or a regenerative amplifier X-ray free-electron laser (RAFEL), by numerical simulations including finite-element analysis and advanced data processing. Pulse-by-pulse transient thermal deformation of a 50 µm-thick diamond crystal has been performed with X-ray pulse repetition rates between 50 kHz and 1 MHz. Results for temperature and thermal deformation have been compared with the results of transient analysis using a continuous wave (CW) power loading. Temperature and thermal deformation results from pulse-by-pulse transient analysis vary with time about the results for the CW case for the same average power. The variation amplitude increases with pulse energy and decreases with repetition rate. When the repetition rate increases to infinity, both temperature and thermal deformation converge to the results for the CW case. Two critical time scales for the operation of crystal optics in a CBXFEL are (1) first-turn time, i.e. the time for the XFEL pulse to complete the first turn around the cavity so that the crystal sees the recirculated XFEL pulse, and (2) period-end time, i.e. the time that the next electron bunch arrives for the amplification, so that the crystal outcouples the amplified FEL power. For the same average power, simulation results show that the crystal thermal deformation seen by the XFEL beam decreases with repetition rate at the first-turn time of a 300 m-long cavity and increases with repetition rate at the period-end time. For the wavefront preservation requirement of the crystal optics, a pulse-energy versus repetition-rate phase space has been established. The upper bounds of the pulse energy at both first-turn and period-end times decreases with repetition rate, especially at the period-end time. The upper bound of the thermal deformation of the crystal at the period-end time for any repetition frequency can be estimated from the CW case. For a water-cooled diamond crystal of dimension 5 mm × 5 mm × 0.05 mm, the time to reach a quasi steady-state is about 50 ms for temperature and 50 µs for thermal deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Galvanostatic Intermittent Titration Technique Reinvented: Part I. A Critical Review

The galvanostatic intermittent titration technique (GITT), introduced in 1977 by Weppner and Huggins, provided a readily accessible means to measuring the chemical diffusion coefficient of electrochemical electrode materials. The method continues to be widely used today, but the reported diffusivity values are highly inconsistent, ranging as much as four orders of magnitude for some Li layered oxide compositions. Even qualitative trends of diffusivity are inconsistent, suggesting significant flaws in the implementation of the method. Other variants of the GITT method also suffer from similar inconsistency problems. Here we identify numerous sources of significant error including composition-dependent reaction overpotentials, mathematical flaws in the relaxation analysis methods, finite-size and non-planar geometry effects, inter-particle inhomogeneity issues, early transient effects, and surface area uncertainties. We propose a simple relaxation analysis scheme using the time variable t relax + τ - t relax , where t relax is relaxation time and τ is the galvanostatic pulse duration. We also propose to use dense diffusion-limited samples to isolate the bulk-diffusion process in the time domain. Chemical diffusivity can be extracted much more reliably with this improved implementation of the GITT method.

25 ENERGY STORAGE↗

DOE OSTI2026

Roadmap to Advance Heliostat Technologies for High Temperature Solar-Thermal Systems

Since its establishment, the Heliostat Consortium (HelioCon) has made substantial progress toward closing many of the gaps in concentrating solar power (CSP) research. Numerous techno-economic studies have been performed, investigating topics ranging from the trade-off between size and temperature for industrial process heat applications to optimization of the heliostat design itself for various applications. Significant improvements have been made in optical metrology techniques, with first steps toward in situ measurement of heliostat fields. Several standards have been, and continue to be, developed with the coordination of an international group of CSP industry participants. Training programs have been developed, with universities including CSP in their engineering curricula, and many public webinars have been held to provide broad access to the latest CSP research. Improved CSP components such as mirror facets and wireless communication systems have been developed, and the solar tower at Sandia National Laboratories has been upgraded with a testbed for closed-loop controls research and development. Field deployment challenges involving heliostat foundations and sensitive wildlife habitats have been explored, with progress made toward methods for streamlining project development and permitting. Additional knowledge has been added to the body of work on wind behavior of heliostats and arrays of heliostats, with progress made toward a holistic understanding of wind design methods. Finally, techniques have been developed and demonstrated for assessing soiling conditions at a proposed project site, with predictive models for the soiling rate showing good results. Taking these results together, HelioCon has contributed greatly to the global CSP research and development effort over the past several years.

14 SOLAR ENERGY↗

DOE OSTI2024

A simplified all-frequency stable formulation with an implicit Coulomb gauge

A potential-based finite element formulation has been developed in the past to circumvent the low-frequency breakdown issues commonly encountered in electromagnetic (EM) simulations of low-frequency and multiscale problems. In this formulation, the magnetic vector and electric scalar potentials have been employed to express the electric field and magnetic flux, leading to a set of two equations that represents all four Maxwell’s equations and the current continuity equation. To enforce the Coulomb gauge, an auxiliary potential has been introduced, which results in a total of three equations to be solved simultaneously. To reduce the number of equations and unknowns needed in a simulation, a simplified formulation is proposed in this paper to enforce the Coulomb gauge implicitly without the need for the auxiliary potential. A numerical example is given to demonstrate the accuracy of the proposed formulation.

Mekonnen, Minyechil↗

DOE OSTI2024

Stress-hybrid virtual element method on six-noded triangular meshes for compressible and nearly-incompressible linear elasticity

In this paper, we present a first-order Stress-Hybrid Virtual Element Method (SH-VEM) on six-noded triangular meshes for linear plane elasticity. Here, we adopt the Hellinger–Reissner variational principle to construct a weak equilibrium condition and a stress based projection operator. In each element, the stress projection operator is expressed in terms of the nodal displacements, which leads to a displacement based formulation. This stress-hybrid approach assumes a globally continuous displacement field while the stress field is discontinuous across each element. The stress field is initially represented by divergence-free tensor polynomials based on Airy stress functions, but we also present a formulation that uses a penalty term to enforce the element equilibrium conditions, referred to as the Penalty Stress-Hybrid Virtual Element Method (PSH-VEM). Numerical results are presented for PSH-VEM and SH-VEM, and we compare their convergence to the composite triangle FEM and B-bar VEM on benchmark problems in linear elasticity. The SH-VEM converges optimally in the L 2 norm of the displacement, energy seminorm, and the L 2 norm of hydrostatic stress. Furthermore, the results reveal that PSH-VEM converges in most cases at a faster rate than the expected optimal rate, but it requires the selection of a suitably chosen penalty parameter.

42 ENGINEERING↗

DOE OSTI2021

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗

DOE OSTI2025

What are grana in chloroplasts of vascular plants good for?

All plants and green algae contain stacked grana thylakoid membranes in their chloroplasts, underscoring an evolutionary pressure to maintain this unique structural feature. In addition, numerous studies have demonstrated that particular lateral and vertical dimensions of grana facilitate the function, regulation and repair of the photosynthetic machinery responsible for energy conversion. In this review, we present an updated overview of our understanding concerning the structure of grana thylakoids, the forces that contribute to their formation and their architectural dynamics. After establishing the structural foundation, we continue by reviewing the implications of grana formation on light harvesting, electron transport and protein maintenance in the thylakoid membranes of vascular plants. Regarding light harvesting, we discuss the implications of grana formation on energy spillover, macromolecular crowding, lateral versus vertical excitation energy transfer, and state transition. With respect to electron transport, we summarize recent findings regarding the functional connectivity of spatially separated photosystems facilitated by grana formation through small mobile electron carriers. We also explore how the dynamic responses of grana architecture regulate electron transport. Finally, we address how grana formation contributes to the organization of protein repair and biogenesis within thylakoid membranes.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Finite-Time Teleportation Phase Transition in Random Quantum Circuits

How long does it take to entangle two distant qubits in a quantum circuit evolved by generic unitary dynamics? Here, we show that if the time evolution is followed by measurements of all but two infinitely separated test qubits, then the entanglement between them can undergo a phase transition and become nonzero at a finite critical time t c . The fidelity of teleporting a quantum state from an input qubit to an infinitely distant output qubit shows the same critical onset. Specifically, these finite-time transitions occur in short-range interacting two-dimensional random unitary circuits and in sufficiently long-range interacting one-dimensional circuits. The phase transition is understood by mapping the random continuous-time evolution to a finite-temperature thermal state of an effective spin Hamiltonian, where the inverse temperature equals the evolution time in the circuit. In this framework, the entanglement between two distant qubits at times t > t c corresponds to the emergence of long-range ferromagnetic spin correlations below the critical temperature. We verify these predictions using numerical simulation of Clifford circuits and propose potential realizations in existing platforms for quantum simulation.

36 MATERIALS SCIENCE↗

DOE OSTI2022

The COMBLE campaign: a study of marine boundary-layer clouds in Arctic cold-air outbreaks

One of the most intense air mass transformations on Earth happens when cold air flows from frozen surfaces to much warmer open water in cold-air outbreaks (CAOs), a process captured beautifully in satellite imagery. Despite the ubiquity of the CAO cloud regime over high-latitude oceans, we have a rather poor understanding of its properties, its role in energy and water cycles, and its treatment in weather and climate models. The Cold-air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) was conducted to better understand this regime and its representation in models. COMBLE aimed to examine the relations between surface fluxes, boundary-layer structure, aerosol, cloud and precipitation properties, and mesoscale circulations in marine CAOs. Processes affecting these properties largely fall in a range of scales where boundary-layer processes, convection, and precipitation are tightly coupled, which makes accurate representation of the CAO cloud regime in numerical weather prediction and global climate models most challenging. COMBLE deployed an Atmospheric Radiation Measurement Mobile Facility at a coastal site in northern Scandinavia (69°N), with additional instruments on Bear Island (75°N), from December 2019 to May 2020. CAO conditions were experienced 19% (21%) of the time at the main site (on Bear Island). A comprehensive suite of continuous in situ and remote sensing observations of atmospheric conditions, clouds, precipitation, and aerosol were collected. Because of the clouds’ well-defined origin, their shallow depth, and the broad range of observed temperature and aerosol concentrations, the COMBLE dataset provides a powerful modeling test bed for improving the representation of mixed-phase cloud processes in large-eddy simulations and large-scale models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Temperature Effects on Droplet Oscillation Decay with Application to Fuel Property Measurement

Observation of oscillation decay in droplets has been shown to be an effective approach in determining physical properties such as viscosity and surface tension for emerging biofuels using μl quantities. Herein this work extends the approach to higher temperatures relevant to fuel injection conditions for internal combustion engines. Experiments are conducted that use high-resolution strobed imaging of moving, heated, μm-sized, fuel droplets to capture shape oscillation decay through image processing and analysis. Two fuels are investigated, iso-butanol and a primary reference fuel (PRF 84), which is a mixture of iso-octane and n-heptane. A piezoelectric droplet generator is used to generate a continuous train of single droplets, which are given an initial perturbation and subsequently undergo damped oscillations, captured using strobed imaging. Surface tension and viscosity calculated at four different temperatures ranging from 30°C-55°C using the frequency and decay time associated with the fundamental mode are found to be within ~ 10% of reference values obtained from the literature. Complementary numerical simulations are performed that utilize a volume-of-fluid approach to track transient droplet oscillation phenomena along with heat and mass transfer in a 2D axisymmetric domain. Simulations, where droplet size and temperature can be independently varied, capture an expected decrease in oscillation frequency and increase in decay time, with increase in fuel temperature. The simulations are further used to study the relative contribution to deviations in surface tension and viscosity predictions due to heat and mass transfer from the droplet, as well as viscous effects violating the inviscid flow assumption in the droplet oscillation theory. Mass loss effects are found to be negligible. Temperature change due to heat transfer is found to have the next highest sensitivity, particularly for the more volatile PRF 84, which has a lower viscosity and associated Ohnesorge number. For isobutanol, which has a higher viscosity and Ohnesorge number, viscous effects contribute the most to deviation in fuel property predictions.

09 BIOMASS FUELS↗

DOE OSTI2025

New particle pusher with hadronic interactions for modeling multimessenger emission from compact objects

We propose novel numerical schemes based on the Boris method in curved spacetime, incorporating both hadronic and radiative interactions for the first time. Once the proton has lost significant energy due to radiative and hadronic losses, and its gyroradius has decreased below typical scales on which the electromagnetic field varies, we apply a guiding center approximation (GCA). We fundamentally simulate collision processes either with a Monte-Carlo method or, where applicable, as a continuous energy loss, contingent on the local optical depth. To test our algorithm for the first time combining the effects of electromagnetic, gravitational, and radiation fields including hadronic interactions, we simulate highly relativistic protons traveling through various electromagnetic fields and proton backgrounds. We provide unit tests in various spatially dependent electromagnetic and gravitational fields and background photon and proton distributions, comparing the trajectory against analytic results. We propose that our method can be used to analyze hadronic interactions in black hole accretion disks, jets, and coronae to study the neutrino abundance from active galactic nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Materials for a Sustainable Microelectronics Future: Electric Field Control of Magnetism with Multiferroics

This article is written on behalf of many colleagues, collaborators, and researchers in the field of complex oxides as well as current and former students and postdocs who continue to enable and undertake cutting-edge research in the field of multiferroics, magnetoelectrics, and the pursuit of electric-field control of magnetism. What I present is something that is extremely exciting from both a fundamental science and applications perspective and has the potential to revolutionize our world, particularly from a sustainability perspective. To realize this potential will require numerous new innovations, both in the fundamental science arena as well as translating these scientific discoveries into real applications. Thus, this article will attempt to bridge the gap between fundamental materials physics and the actual manifestations of the physical concepts into real-life applications. Further, I hope this article will help spur more translational research within the broad materials community.

36 MATERIALS SCIENCE↗

DOE OSTI2025

CSAC observations from 2025 CHEDAR review [Slides]

Overall, the CSAC is impressed by the strong start of CHEDAR, a multi-institutional center with broader scope than its predecessor. The CSAC was impressed and energized by discussions with center students and staff; the enthusiasm for science is universal and the sense of the contribution in context to the NNSA mission is mostly understood. The Center is well positioned to build on a strong start and (1) enhance integration of code-development targeted to, and numerical design in on-going experimental campaigns, (2) focus and better tell the story of the astrophysical connection for on-going experimental campaigns, and (3) expand and deepen ties between Center students (and faculty) and the NNSA National Laboratories. Overall, the CSAC was encouraged and energized by the materials presented. Now, the Center has to recognize the need to continue growing strategically. There is abundant access to experimental facilities – is that a challenge to focused effort? The PIs will need to balance the pace of experimental campaigns and analysis and simulation so that informed progress and focused investigation can develop in each campaign. There was a clear need to connect simulations to the goals of experimental campaigns – what is the best way to cultivate that connection? Are there adequate resources in the near term to develop code capability? What is the balance between close personal mentorship and collaboration across the separate institutions that compose the Center? The CSAC feels confident that the team leading CHEDAR will find creative solution to these challenges.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

NH 4 OH Looping with Membrane CO 2 Absorber and Distributed Stripper for Enhanced Algae Growth

The University of Kentucky Center for Applied Energy (UK CAER) has devised a unique, integrated CO2 capture and utilization technology. CO2 from coal-fired power generation flue gas is first captured at half the operating cost of a typical aqueous CO2 capture system (CCS), distributed in an aqueous stream and then fixed by algae in bioreactors where the algae production is increased by 50% over that with a typical intermittent nutrient feeding system. Lower CCS operating cost is achieved by eliminating the flue gas pretreatment step for cooling and SO2 removal, eliminating steam extraction from the power generation steam cycle for solvent regeneration, and eliminating CO2 compression. Higher algae production is achieved by continuous, just-in-time nutrient feed to the bioreactors directly from a distributed solvent regenerator, which maintains the bioreactor pH for optimum growth. The process starts with a uniquely configured membrane absorber, where the flue gas is indirectly contacted with an ammonium hydroxide (NH4OH) solvent. Dissolved NH3 is attractive for both CO2 capture and as an algae nutrient. For CO2 capture it is inexpensive, has a low regeneration energy, is thermally- and oxidatively-stable and has a viscosity near that of water, which makes is easy to transport. Numerous studies have shown that the scrubbing capacity of NH3 is approximately 0.9-1.2 kg of CO2/kg of NH3, with a CO2 removal efficiency of ~99% and half the solvent regeneration energy than that of 30 wt% MEA[1, 2, 3]. NH3 is attractive as an algae nutrient due to its low cost. The rich NH4OH solvent is pumped to a set of distributed regenerators which are co-located with the algae bioreactors. Solvent pumping, transport and distribution reduces the balance of plant (BOP) cost compared to a typical aqueous CCS related to the flue gas duct and boost fan required to transport the flue gas. The energy required for the distributed solvent regeneration is supplied by solar-thermal panels eliminating the need for steam extraction from the power generation steam cycle. After solvent regeneration, the product stream contains both the CO2 captured from the flue gas and volatized NH3 from the solvent. This product stream is fed directly to the bioreactors, eliminating the need for compression of the CO2 stream. The relative amounts of CO2 and NH3 in the product stream are adjusted and controlled by a controlling the regeneration conditions (pressure and temperature). The continuous feed of the right ratio of nutrients overcomes the typical inhibition of algae growth resulting from frequent pH swings in the bioreactor due to unbalanced (intermittent) feeding systems for CO2 and N. Also, because the regenerators will operate at pressure and be located in close proximity to the bioreactors, there is no worry about pressure drop when sparging the gas into the algae. Sparging produces small bubbles which is beneficial for mass transfer efficiency. One known challenge when using an NH4OH solvent is high NH3 emission. Hydrophobic membranes are used for CO2 capture using an aqueous NH3 solution[4, 5] without the direct contact between flue gas and aqueous solution. Additionally, UK CAER CO2 capture and utilization process manages NH3 slip in three extra measures. First, NH3 slip is minimized by working with minimal species partial pressure, which is proportional to the concentration in the liquid. Hence, lowering the capture solvent concentration will lower the NH3 partial pressure. Second, UK CAER’s previous work has demonstrated that the addition of Zn2+ into NH3 solutions to chelate the NH3 can reduce NH3 volatility. Third, the configuration of the membrane CO2 absorber utilizes condensed water from the flue gas to continually wash the gas-side of the membrane to reduce fouling and recapture NH3 slip. Additional details about the UK CAER unique, integrated CO2 capture and utilization technology will be presented along with technology development plans. Diao, N., Q. Li, and Z. Fang. 2004. Heat transfer in ground heat exchangers with groundwater advection. International Journal of Thermal Sciences. 43: 1203-1211, He, Q., M. Chen, L. Meng, K. Liu, and W. Pan. 2004. Study on Carbon Dioxide Removal from Flue Gas by Absorption of Aqueous Ammonia. Western Kentucky University. Yeh, A.C., and H. Bai. 1999. Comparison of ammonia and monoethanolamine solvents to reduce CO2 greenhouse gas emissions. The Science of the Total Environment. 228: 121-133, Villeneuve, K., D. Roizard, J.C. Remigy, M. Iacono, and S. Rode. 2018. CO2 capture by aqueous ammonia with hollow fiber membrane contactors: Gas phase reactions and performance stability. Separation and Purification Technology, 199: 189-197, Toro Molina, C., and C. Bouallou. 2016. Carbon dioxide absorption by ammonia intensified with membrane contactors. Clean Techn Environ Policy 18, 2133–2146 (2016)

20 FOSSIL-FUELED POWER PLANTS↗