SEARCH · Search NASA

Results for “input device”

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 1,729 records · Page 96

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Sci-Thur AM: YIS – 04: Stopping power-to-Cherenkov power ratios and beam quality specification for clinical Cherenkov emission dosimetry of electrons: beam-specific effects and experimental validation

Purpose: To advance towards clinical Cherenkov emission (CE)-based dosimetry by investigating beam-specific effects on Monte Carlo-calculated electron-beam stopping power-to-CE power ratios (SCRs), addressing electron beam quality specification in terms of CE, and validating simulations with measurements. Methods: The EGSnrc user code SPRRZnrc, used to calculate Spencer-Attix stopping-power ratios, was modified to instead calculate SCRs. SCRs were calculated for 6- to 22-MeV clinical electron beams from Varian TrueBeam, Clinac 21EX, and Clinac 2100C/D accelerators. Experiments were performed with a 20-MeV electron beam from a Varian TrueBeam accelerator, using a diffraction grating spectrometer with optical fiber input and a cooled back-illuminated CCD. A fluorophore was dissolved in the water to remove CE signal anisotropy. Results: It was found that angular spread of the incident beam has little effect on the SCR (≤ 0.3% at d{sub max}), while both the electron spectrum and photon contamination increase the SCR at shallow depths and decrease it at large depths. A universal data fit of R{sub 50} in terms of C{sub 50} (50% CE depth) revealed a strong linear dependence (R{sup 2} > 0.9999). The SCR was fit with a Burns-type equation (R{sup 2} = 0.9974, NRMSD = 0.5%). Below-threshold incident radiation was found to have minimal effect on beam quality specification (< 0.1%). Experiments and simulations were in good agreement. Conclusions: Our findings confirm the feasibility of the proposed CE dosimetry method, contingent on computation of SCRs from additional accelerators and on further experimental validation. This work constitutes an important step towards clinical high-resolution out-of-beam CE dosimetry.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2017

High temperature ceramic heat exchanger for solid oxide fuel cell

A ceramic heat exchanger with high effectiveness and low-pressure drop has been developed for application as a cathode air preheater for a Solid Oxide Fuel Cell (SOFC). At the operating conditions of SOFCs, typical metallic alloys as those used in commercial heat exchangers may undergo chromium volatilization, which is a known cathode degradation mechanism that reduces SOFC performance and life. Use of ceramics such as alumina or alumina-silicate instead of chromium-containing metal is one approach to eliminate the effects of chromium on the SOFC cathode. This project leverages the geometric design of a heat exchanger previously prototyped and tested [1] to fabricate a novel heat exchanger made from a ceramic material. All calculations of heat exchanger performance requirements are based on state of the art SOFC operating conditions. The alumina-silicate material selected to fabricate the modular prototype heat exchanger was made via a heat transfer based trade-off analysis that shows that for low pressure-drop devices, flow convection is the limiting heat transfer mechanism, and therefore wall conductivity becomes less relevant to heat transfer effectiveness. The modular design presented allows for incremental aggregation of modules to target a broad range of operating conditions typical of present and upcoming SOFC applications (e.g., 25 to 400 kWe). The completed prototype was subjected to a series of static leak/pressure tests, dynamic pressure drop tests both at ambient and realistic service temperature conditions, and finally to a matrix of heat transfer experiments to fully characterize its heat transfer effectiveness (including tests at hot gas temperatures up to 1033 K or 1400oF). The stand-alone system was demonstrated in the laboratory to successfully satisfy the design criteria at realistic operating conditions, showing a total pressure drop at the design flow of 8 g/sec of 6 kPa (vs. 3.45 kPa approximate requirement). It was also shown to have heat transfer effectiveness of 76% (with 73% being the design requirement). Moreover, the unit was shown to operate leak-free, and the ceramic elements and seals were capable of withstanding the thermal cycling resulting from the testing procedure. The heat exchanger was then integrated into a modified fuel cell test stand and tested at steady state operation. Unfortunately, the level of performance achieved was lower than expected based on laboratory performance tests. This appears to have been the result of an unforeseen loosening of a heat exchanger structural component, which resulted into a leak that prevented proper supply of hot gas to the heat exchanger. However, in general terms, the heat exchanger/fuel cell integration process provides confidence that a future application of the unit with a SOFC is feasible, assuming structural issues are addressed. It is considered that the difficulties encountered at the last stage of the program provided invaluable input in locating shortcomings of the mechanical design and integration approach. A portion of this work was presented at the ASME Power & Energy Conference, in Charlotte, NC, on June 28th 2016 in a paper entitled “Development of a Ceramic Heat Exchanger for Application as Solid Oxide Fuel Cell Cathode Air Preheater” and the corresponding paper was published in the Proceedings of the ASME 2016 Power & Energy Conference, Paper No.: PowerEnergy2016-59333 [8].

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2015

Flux-driven simulations of turbulence collapse

Using three-dimensional nonlinear simulations of tokamak turbulence, we show that an edge transport barrier (ETB) forms naturally once input power exceeds a threshold value. Profiles, turbulence-driven flows, and neoclassical coefficients are evolved self-consistently. A slow power ramp-up simulation shows that ETB transition is triggered by the turbulence-driven flows via an intermediate phase which involves coherent oscillation of turbulence intensity and E×B flow shear. A novel observation of the evolution is that the turbulence collapses and the ETB transition begins when R{sub T} > 1 at t = t{sub R} (R{sub T}: normalized Reynolds power), while the conventional transition criterion (ω{sub E×B}>γ{sub lin} where ω{sub E×B} denotes mean flow shear) is satisfied only after t = t{sub C} ( >t{sub R}), when the mean flow shear grows due to positive feedback.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

Device Scale Modeling of Solvent Absorption using MFIX-TFM

Recent climate change is largely attributed to greenhouse gases (e.g., carbon dioxide, methane) and fossil fuels account for a large majority of global CO 2 emissions. That said, fossil fuels will continue to play a significant role in the generation of power for the foreseeable future. The extent to which CO 2 is emitted needs to be reduced, however, carbon capture and sequestration are also necessary actions to tackle climate change. Different approaches exist for CO 2 capture including both post-combustion and pre-combustion technologies, oxy-fuel combustion and/or chemical looping combustion. The focus of this effort is on post-combustion solvent-absorption technology. To apply CO 2 technologies at commercial scale, the availability and maturity and the potential for scalability of that technology need to be considered. Solvent absorption is a proven technology but not at the scale needed by typical power plant. The scale up and down and design of laboratory and commercial packed bed reactors depends heavily on the specific knowledge of two-phase pressure drop, liquid holdup, the wetting efficiency and mass transfer efficiency as a function of operating conditions. Simple scaling rules often fail to provide proper design. Conventional reactor design modeling approaches will generally characterize complex non-ideal flow and mixing patterns using simplified and/or mechanistic flow assumptions. While there are varying levels of complexity used within these approaches, none of these models resolve the local velocity fields. Consequently, they are unable to account for important design factors such as flow maldistribution and channeling from a fundamental perspective. Ideally design would be aided by development of predictive models based on truer representation of the physical and chemical processes that occur at different scales. Computational fluid dynamic (CFD) models are based on multidimensional flow equations with first principle foundations. CFD models can include a more accurate physical description of flow processes and be modified to include more complex behavior. Wetting performance and spatial liquid distribution inside the absorber are recognized as weak areas of knowledge requiring further investigation. CFD tools offer a possible method to investigating such topics and gaining a better understanding of their influence on reactor performance. This report focuses first on describing a hydrodynamic model for countercurrent gas-liquid flow through a packed column and then on the chemistry, heat and mass transfer specific to CO 2 absorption using monoethanolamine (MEA). The indicated model is implemented in MFIX, a CFD open source software package. The user defined functions needed to build this model are described in detail along with the keywords for the corresponding input file. A test case is outlined along with a few results. The example serves to briefly illustrate the developed CFD tool and its potential capability to investigate solvent absorption.

09 BIOMASS FUELS↗

NASA NTRS2002

Molten-Metal Droplet Deposition on a Moving Substrate in Microgravity: Aiding the Development of Novel Technologies for Microelectronic Assembly

Driven by advancements in microelectronics manufacturing, this research investigates the oblique (non-axisymmetric) impact of liquid-metal droplets on flat substrates. The problem of interest is relevant to the development of the novel technology of on-demand dispension (printing) of microscopic solder deposits for the surface mounting of microelectronic devices. The technology, known as solder jetting, features on-demand deposition of miniature solder droplets (30 to 120 microns in diameter) in very fine, very accurate patterns using techniques analogous to those developed for the ink-jet printing industry. Despite its promise, severe limitations exist currently with regards to the throughput rates of the technology; some of these limitations are largely due to the lack of the capability for reliable prediction of solder bump positioning and shapes, especially under ballistic deposition conditions where the droplet impact phenomena are inherently three-dimensional. The study consists of a theoretical and an experimental component. The theoretical work uses a finite element formulation to simulate numerically the non-axisymmetric (3-D) fluid mechanics and heat transfer phenomena of a liquid solder droplet impacting at an angle alpha on a flat substrate. The work focuses on the pre-solidification regime. The modeling of the most challenging fluid mechanics part of the process has been completed successfully. It is based upon the full laminar Navier-Stokes equations employing a Lagrangian frame of reference. Due to the large droplet deformation, the surface (skin) as well as the volumetric mesh have to be regenerated during the calculations in order to maintain the high accuracy of the numerical scheme. The pressure and velocity fields are then interpolated on the newly created mesh. The numerical predictions are being tested against experiments, for cases where wetting phenomena are not important. For the impact parameters used in the example shown (We = 2.38, Fr = 16300, Re = 157), the droplet rolls along the substrate, but its shape remains practically axisymmetric for all impact angles within the range from 0 to 60 deg. Interestingly, the substrate/droplet contact area during the recoiling phase of the impact is not a monotonically decreasing function of time. The experimental component of the research tests the numerical predictions and provides necessary input data (contact angles) for the theoretical model. The experiments are performed in microgravity (2.2s drop tower of the NASA GRC) in order to allow for the use of mm-size solder droplets, which make feasible the performance of accurate measurements, while maintaining similitude of the relevant fluid dynamic groups (Re, Fr, We, Ste). Preliminary oblique impact experiments have been performed using water droplets in normal gravity.

Megaridis, C. M.↗

DOE OSTI2017

Quantum Mechanical Simulations of Complex Nanostructures for Photovoltaic Applications

A quantitative understanding of the electronic excitations in nanostructures, especially complex nanostructures, is crucial for making new-generation photovoltaic (PV) cells based on nanotechnology, which have high efficiency and low cost. Yet current quantum mechanical simulation methods are either computationally too expensive or not accurate and reliable enough, hindering the rational design of the nanoscale PV cells. The PI seeks to develop new methodologies to overcome the challenges in this very difficult and long-lasting problem, pushing the field forward so that electronic excitations can be accurately predicted for systems involving thousands of atoms. The primary objective of this project is to develop new approaches for electronic excitation calculations that are more accurate than traditional density functional theory (DFT) and are applicable to systems larger than what current beyond-DFT methods can treat. In this proposal, the PI will first address the excited-state problem within the DFT framework to obtain quasiparticle energies from both Kohn-Sham (KS) eigenvalues and orbitals; and the electron-hole binding energy will be computed based on screened Coulomb interaction of corresponding DFT orbitals. The accuracy of these approaches will be examined against many-body methods of GW/BSE and quantum Monte Carlo (QMC). The PI will also work on improving the accuracy and efficiency of the GW/BSE and QMC methods in electronic excitation computations by using better KS orbitals obtained from orbital-dependent DFT as inputs. Then an extended QMC database of ground- and excited-state properties will be developed, and this will be spot checked and supplemented with data from GW/BSE calculations. The investigation will subsequently focus on the development of an improved exchange-correlation (XC) density functional beyond the current generalized gradient approximation (GGA) level of parameterization, with parameters fitted to the QMC database. This will allow the ground-state properties of focus systems to be more precisely predicted using DFT. These new developments will then be applied to investigate a chosen set of complex nanostructures that have great potential for opening new routes in designing materials with improved transport, electronic, and optical properties for PV and other optoelectronic usages: (1) Hybrid interfaces between materials with distinct electronic and optical properties, such as organic molecules (conjugated polymers, e.g. P3HT) and inorganic semiconducting materials (Si and ZnO). Complicated interface structures, including interface bonding configurations, compositional and geometrical blending patterns, interfacial defects, and various sizes and shapes of inorganic nanomaterials, will be considered for the purpose of understanding the working mechanisms of present organic/nano PV systems and designing optimum interface structures for fast charge separation and injection. (2) Complex-structured semiconducting nanomaterials that could induce charge separation without pn- or hetero-junctions. The new methodology will allow the PI to investigate the performance of realistic semiconducting nanomaterials of internal (impurities, defects, etc.) and external (uneven surface, mechanical twisting and bending, surface chemistry, etc.) complexities on optical absorption and charge transport against charge trapping and recombination. Of particular interest is whether such structural complexity in a single material could even be beneficial for PV usage, for example, charge separation through morphology control. Successful completion of the proposed DFT methodology would have a far-reaching impact on our ability to study and exploit the nature of electronic excitations in complex materials, advancing the design of next-generation electronic and optoelectronic devices in all facets of renewable energy conversion and storage, including photovoltaics, thermoelectricity, photochemistry, etc.

14 SOLAR ENERGY↗

DOE OSTI2015

Andra's strategy and approach for operational safety of French underground facility Cigeo - 15127

According to the French Act on radioactive waste of 28 June 2006, Andra shall design a reversible repository. In response to this demand, Andra is developing the industrial project namely 'Cigeo', a reversible geological disposal facility for HLW and ILW located in Meuse/Haute-Marne. The previous safety cases of 2005, 2009 and 2012 give a sound basis to assess Cigeo's safety both for the operational and post-closure periods. Regarding operational safety, in this new stage of the project, the challenging issues for the preparation of the safety case are the following: - To adapt the various regulatory frameworks (nuclear and non-nuclear) to the underground facility - To design the underground facility and handling equipment's to reduce at first risks in the underground facility as such radioactive particles release, gas explosion, fire and co-activity and limit their consequence. - To establish waste acceptance criteria to be fulfilled by waste producers in line with operational safety and design options. The objective of the paper is to present the experience of Andra in operation safety of the underground facility of Cigeo. Operational safety for Cigeo is assessed as much as possible according to the existing French regulations and guidance for nuclear facilities. Indeed, compared to existing nuclear facilities, a geological repository is a specific case due to the existence of underground facilities with particularities such as: - ventilation management in a large scale underground facility, - fire risk management. As an input for the industrial design of Cigeo for operation, Andra defined a fire guidance and a specific reference safety guidance, based on the adaptation of the French nuclear regulatory texts for and when feasible to mines and tunnels regulatory texts and on its own experience gained through safety operational reports since 2005 and their reviews by the Safety Authority. - underground co-activity of nuclear operation with construction works. According to the common 'defense-in depth' principle for nuclear facilities, these guidelines for the industrial design of Cigeo for operation consider the successive levels of defense of: - Preventing the risk and its spreading by design and construction; - Detecting incidents and early fighting by design and operation dispositions (detection and control systems, survey procedure); - Controlling incidents and limiting of their consequences by safeguarding systems, specific procedures. At that stage, design options for Cigeo constitute an adequate response to the operational safety functions and requirements: - Two levels of containment are set up: the first one is ensured by the waste package, the second one is ensured from the surface facility to the disposal vaults by the transfer cask during the transfer of waste packages in the disposal vaults by the disposal vault itself and a nuclear ventilation with HEPA filter; - Favorable options for fire risk prevention and mitigation: transport system minimizing heat load (cable car for the transfer from the surface, electric transfer engines), twin access tunnels to ease fire-fighting and evacuation...; - Permanent complete physical separation of nuclear zone and construction zone (including ventilation systems). In practice, operation safety assessment of Cigeo, designed following the aforesaid reference safety guidelines, follows a common risk analysis (internal aggressions such as auxiliary losses, chemical and toxic risks, fall or collision, fire, human factors and external aggressions such as earthquake, air traffic, flooding) with the following steps: - Identification of the major risks and the design basis scenarios that have to be considered for dimensioning of the installations and the safety devices; - On the basis of these different situations, those who will have to be quantified are identified; one or more representative and envelop situations are chosen in order to evaluate robustness of the installation. Quantitative evaluation of these situations permits to: - identify the consequences on man, environment and installation, and in particular to gain confidence in the maintaining of the safety functions; - identify or define the technic and/or organisational dispositions classified as 'important elements for protection' (IEP); - Situations beyond dimensioning permit to verify the robustness of the installation. They permit, on one hand to identify possible complementary material or organisational dispositions; on the other hand, to define the 'hard core' i.e. structures, systems and components that permit to maintain the confinement function of radioactive substance. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2014

Highly-laminated, high-saturation-flux-density, magnetic cores for on-chip inductors in power converter applications

The objective of this project is development of a manufacturing process to realize thick highly-laminated metallic magnetic cores with suppressed eddy current losses for power applications and demonstration of integrated magnetic components using the process. The continual reduction in size and increase in functionality of portable electronic devices drives the development of ultracompact electrical power converters and therefore their associated inductors, which typically occupy the largest volume within the converter. This inductor size reduction can be achieved by increasing both the converter switching frequency as well as inductor operating flux density by utilizing appropriate magnetic core materials. Magnetically soft metallic alloys (e.g., NiFe and CoNiFe) are good core material candidates due to their superior magnetic properties (e.g., high saturation flux density and low coercivity) compared to conventional ferrite. However, to minimize eddy current losses, the core thickness is limited to the skin depth (typically less than 3 µm for 10 MHz operation); cores of this thickness are likely insufficient for handling high power (> 10W). Therefore, there is great interest in the development of highly-laminated metallic alloy cores with single lamination thickness in the sub-micron range to achieve suppressed eddy-current losses and high power handling capacity simultaneously. During the project, highly-laminated metallic alloy cores comprising up to 2000 layers of permalloy (Ni 80 Fe 20 ) or CoNiFe materials with single lamination thickness less than 1 µm have been fabricated based on automated multilayer electrodeposition technology compatible with CMOS manufacturing. The developed highly-laminated metallic cores preserve the high saturation flux density and low coercivity characteristics of metallic alloys, while suppressing eddy-current losses less than 1 W/cm3 at operation flux density of 0.5 T and frequency over 1 MHz. The developed highly-laminated metallic alloy cores have been further integrated into microfabricated windings to realize compact, on-chip inductors. Utilizing CMOS-compatible integration techniques, circular (toroid) and bar-shape (solenoid) microinductors have been fabricated exhibiting a constant inductance higher than 1 µH up to 50 MHz with peak quality factor as high as 35 at 5 MHz. The microfabricated inductor with highly-laminated core was tested in a power converter and demonstrated a converter efficiency of 93% at 5.5 MHz and 42 W under 100V input and 35V output operation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Wave

The U.S. Department of Energy and the National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence (AI) data centers and other variable loads. This dataset entry describes hydrogen production using a single, simulated wave energy conversion device. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel Hydrogen. While the unit supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack. For the wave energy, NLR used a wave energy converter model from PacWave. These devices can be equipped with accumulators and pressure relief values to smooth the power output by storing and releasing hydraulic energy. Using a peak power output of 10 MW, the model created two 25-minute profiles: one with and one without the accumulators and pressure relief valves. To down select the profile data from the native resolution of 20 Hz to 1 Hz, NLR took the mean of every 20 data points. NLR experimented with two simulated wave energy power plants: one that peaks at 10 MW, and one that peaks at 5 MW. These profiles were scaled for the physical 1.25 MW electrolyzer by multiplying the original profiles by one eighth and one quarter, respectively. The first profile matches the capacity rating of eight of the 1.25 MW electrolyzers, while the second matches four electrolyzers. Finally, NLR experimented with two settings for the electrolyzer power supply minimum and maximum current ramp rates (gain and slew): 200 and 400 amperes per second. The simulated profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single wave electrolysis experiment and is formatted as follows: {technology}-{accumulator?}_{number of 1.25 MW electrolyzers connected}-{electrolyzer ramp rate in amperes/second} For instance, “wavePacWave-Noacc_4-400.zip” represents the 25 minute-long experiment using the PacWave’s wave energy converter model, equipped with no accumulator, connected to four 1.25-MW electrolyzers with their power supplies set to a maximum current ramp rate (gain and slew) of 400 A/s. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. An experiment, labeled “characterization_200.zip”, demonstrates the MC250 electrolyzer steady-state response with 30 minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all wave profiles combined into one dataset labeled "combined_wave_experiments.csv". NLR also built an AI/machine-learning predictive model based on these datasets. The model ingests the electrolyzer current command in amperes, as well as various pressures and temperatures across the system, and predicts hydrogen output in kilograms per hour. The complete model can be found at https://huggingface.co/NatLabRockies/ptmelt-hydrogen-electrolysis.

08 HYDROGEN↗

DOE OSTI2015

T2-weighted four dimensional magnetic resonance imaging with result-driven phase sorting

Purpose: T2-weighted MRI provides excellent tumor-to-tissue contrast for target volume delineation in radiation therapy treatment planning. This study aims at developing a novel T2-weighted retrospective four dimensional magnetic resonance imaging (4D-MRI) phase sorting technique for imaging organ/tumor respiratory motion. Methods: A 2D fast T2-weighted half-Fourier acquisition single-shot turbo spin-echo MR sequence was used for image acquisition of 4D-MRI, with a frame rate of 2–3 frames/s. Respiratory motion was measured using an external breathing monitoring device. A phase sorting method was developed to sort the images by their corresponding respiratory phases. Besides, a result-driven strategy was applied to effectively utilize redundant images in the case when multiple images were allocated to a bin. This strategy, selecting the image with minimal amplitude error, will generate the most representative 4D-MRI. Since we are using a different image acquisition mode for 4D imaging (the sequential image acquisition scheme) with the conventionally used cine or helical image acquisition scheme, the 4D dataset sufficient condition was not obviously and directly predictable. An important challenge of the proposed technique was to determine the number of repeated scans (N{sub R}) required to obtain sufficient phase information at each slice position. To tackle this challenge, the authors first conducted computer simulations using real-time position management respiratory signals of the 29 cancer patients under an IRB-approved retrospective study to derive the relationships between N{sub R} and the following factors: number of slices (N{sub S}), number of 4D-MRI respiratory bins (N{sub B}), and starting phase at image acquisition (P{sub 0}). To validate the authors’ technique, 4D-MRI acquisition and reconstruction were simulated on a 4D digital extended cardiac-torso (XCAT) human phantom using simulation derived parameters. Twelve healthy volunteers were involved in an IRB-approved study to investigate the feasibility of this technique. Results: 4D data acquisition completeness (C{sub p}) increases as NR increases in an inverse-exponential fashion (C{sub p} = 100 − 99 × exp(−0.18 × N{sub R}), when N{sub B} = 6, fitted using 29 patients’ data). The N{sub R} required for 4D-MRI reconstruction (defined as achieving 95% completeness, C{sub p} = 95%, N{sub R} = N{sub R,95}) is proportional to N{sub B} (N{sub R,95} ∼ 2.86 × N{sub B}, r = 1.0), but independent of N{sub S} and P{sub 0}. Simulated XCAT 4D-MRI showed a clear pattern of respiratory motion. Tumor motion trajectories measured on 4D-MRI were comparable to the average input signal, with a mean relative amplitude error of 2.7% ± 2.9%. Reconstructed 4D-MRI for healthy volunteers illustrated clear respiratory motion on three orthogonal planes, with minimal image artifacts. The artifacts were presumably caused by breathing irregularity and incompleteness of data acquisition (95% acquired only). The mean relative amplitude error between critical structure trajectory and average breathing curve for 12 healthy volunteers is 2.5 ± 0.3 mm in superior–inferior direction. Conclusions: A novel T2-weighted retrospective phase sorting 4D-MRI technique has been developed and successfully applied on digital phantom and healthy volunteers.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

2015 Summary Report on Industrial and Regulatory Engagement Activities

The Advanced Instrumentation, Information, and Control (II&C) Systems Technologies pathway of the Light Water Reactor Sustainability(LWRS) Program conducts a vigorous engagement strategy with the U.S. nuclear power industry, including the nuclear operating companies, major support organizations, the Nuclear Regulatory Commission (NRC), and suppliers. The goal of this engagement strategy is to develop a shared vision and common understanding across the nuclear industry of the need for II&C modernization, the performance improvement that can be obtained, and the opportunities for collaboration to enact this vision. The primary means of engaging the nuclear operating companies is through a Utility Working Group (UWG), composed of utility representatives that participate in formal meetings and bi-monthly phone calls to provide input on nuclear plant needs and priorities for II&C technologies. Two working groups were initiated during FY 2015 to provide a means for UWG members to focus on particular technologies of interest. The Outage Improvement Working Group consists of eight utilities that participate in periodic conference calls and have access to a share-point web page for acccess to project materials developed in the Advanced Outage Control Center pilot project. In the area of computer-based procedures and automated work packages, the II&C Pathway has worked with the Nuclear Information Technology Strategic Leadership (NITSL) to set up a monthly conference call with interested utility members to discuss various aspects of mobile worker technologies. Twenty one technical and project reports were delivered to the UWG during FY 2015, reflecting the work of the II&C Pathway pilot projects during the year. Distribution of these reports is one of the primary means of transferring to the nuclear industry the knowledge and experience gained during the development of advanced II&C technologies in support of LWR sustainability. Site visits to discuss pilot project activities and future plans were made to Arizona Public Service, Exelon, Duke Energy, Pacific Gas & Electric, SCANA, Southern Nuclear, South Texas Project, STARS Alliance, Tennessee Valley Authority, and Xcel. Discussions were also held on the pathway goals and activities with major industry support organizations during FY 2102, including the Institute of Nuclear Power Operations (INPO), the Nuclear Information Technology Strategic Leadership (NITSL), the Nuclear Energy Institute (NEI), and the Electric Power Research Institute. The Advanced II&C Pathway work was presented at five major industry conferences and Informal discussions were held with key NRC managers at industry conferences. In addition, discussions were held with NRC senior managers on digital regulatory issues through participation on the NEI Digital I&C Working Group. Meetings were held with major industry suppliers and consultants, to explore opportunities for collaboration and to provide a means of pilot project technology transfer. In the international area, discussions were held with Electricite’ de France (EdF) concerning possible collaboration in the area NPP configuration control using intelligent wireless devices.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

NASA NTRS2018

Real-Time Adaptive Drag Minimization Wind Tunnel Investigation of a Flexible Wing with Variable Camber Continuous Trailing Edge Flap System

This paper reports the results of a recently completed real-time adaptive drag minimization wind tunnel investigation of a highly flexible wing wind tunnel model equipped with the Variable Camber Continuous Trailing Flap (VCCTEF) technology at the University of Washington Aeronautical Laboratory (UWAL). The wind tunnel investigation is funded by NASA SBIR Phase II contract with Scientific Systems Company, Inc. (SSCI) and University of Washington (UW) as a subcontractor. The wind tunnel model is a sub-scale Common Research Model (CRM) wing constructed of foam core and fiberglass skin and is aeroelastically scaled to achieve a wing tip deflection of 10% of the wing semi-span which represents a typical wing tip deflection for a modern transport such as Boeing 787. The jig-shape twist of the CRM wing is optimized using a CART3D aero-structural model to achieve the minimum induced drag for the design cruise lift coefficient of 0.5. The wing is equipped with two chord wise cambered segments for each of the six span wise flap sections for a total of 12 individual flap segments that comprise the VCCTEF system. Each of the 12 flap segments is actively controlled by an electric servo-actuator. The real-time adaptive drag optimization strategy includes an on-board aerodynamic model identification, a model excitation, and a real-time drag optimization. The on-board aerodynamic model is constructed parametrically as a function of the angle of attack and flap positions to model the lift and drag coefficients of the wing. The lift coefficient models include a linear model and a second-order model. The drag coefficient models include a quadratic model and a higher-order up to 6th-order model to accurately model the drag coefficient at high angles of attack. The onboard aerodynamic model identification includes a recursive least-squares (RLS) algorithm and a batch least-squares (BLS) algorithm designed to estimate the model parameters. The model excitation method is designed to sample the input set that comprises the angle of attack and the flap positions. Three model excitation methods are developed: random excitation method, sweep method, and iterative angle-of-attack seeking method. The real-time drag optimization includes a generic algorithm developed by SSCI and several optimization methods developed by NASA which include a second-order gradient Newton-Raphson optimization method, an iterative gradient optimization method, a pseudo-inverse optimization method, an analytical optimization method, and an iterative refinement optimization method. The first wind tunnel test entry took place in September 2017. This test revealed major hardware issues and required further redesign of the flap servo mechanisms. The second test entry took place in April 2018. However, the test was not successful due to the issues with the onboard aerodynamic model identification RLS algorithm which incorrectly identified model parameters. This test also provides an experimental comparison study between the VCCTEF and a variable camber discrete trailing edge flap (VCDTEF) without the elastomer transition mechanisms. The experimental result confirms the benefit of the VCCTEF which produces lower drag by 5% than the VCDTEF. The third and final test entry took place in June 2018 after the issues with the RLS algorithm have been identified and corrected. Additional improvements were implemented. These include the BLS algorithm, the iterative angle-of-attack seeking method, the iterative gradient optimization method, and the pseudo-inverse optimization method. The test objectives were successfully demonstrated as the real-time drag optimization identifies several optimal solutions at off-design lift coefficients. The iterative gradient optimization method is found to achieve up to 4.7% drag reduction for the off-design lift coefficient of 0.7. The pseudo-inverse optimization method which does not require the drag coefficient model is found to be quite effective in reducing drag. Up to 9.4% drag reduction for the off-design lift coefficient of 0.7 is achieved with the pseudo-inverse optimization method. The wind tunnel investigation demonstrates the potential of real-time drag optimization technology. Several new capabilities are developed that could enable future adaptive wing technologies for flexible wings equipped with drag control devices such as the VCCTEF.

Drag Optimization↗

DOE OSTI2017

Assimilation of Wave Imaging Radar Observations for Real-time Wave-by-Wave Forecasting

This project addressed Topic 3: “Wave Measurement Instrumentation for Feed Forward Controls” under the FOA number DE-FOA-0000971. The overall goal of the program was to develop a phase-resolving wave forecasting technique for application to the active control of Wave Energy Conversion (WEC) devices. We have developed an approach that couples a wave imaging marine radar with a phase-resolving linear wave model for real-time wave field reconstruction and forward propagation of the wave field in space and time. The scope of the project was to develop and assess the performance of this novel forecasting system. Specific project goals were as follows: Develop and verify a fast, GPU-based (Graphical Processing Unit) wave propagation model suitable for phase-resolved computation of nearshore wave transformation over variable bathymetry; Compare the accuracy and speed of performance of the wave model against a deep water model in their ability to predict wave field transformation in the intermediate water depths (50 to 70 m) typical of planned WEC sites; Develop and implement a variational assimilation algorithm that can ingest wave imaging radar observations and estimate the time-varying wave conditions offshore of the domain of interest such that the observed wave field is best reconstructed throughout the domain and then use this to produce model forecasts for a given WEC location; Collect wave-resolving marine radar data, along with relevant in situ wave data, at a suitable wave energy test site, apply the algorithm to the field data, assess performance, and identify any necessary improvements; and Develop a production cost estimate that addresses the affordability of the wave forecasting technology and include in the Final Report. The developed forecasting algorithm (“Wavecast”) was evaluated for both speed and accuracy against a substantial synthetic dataset. Early in the project, performance tests definitively demonstrated that the system was capable of forecasting in real-time, as the GPU-based wave model backbone was very computationally efficient. The data assimilation algorithm was developed on a polar grid domain in order to match the sampling characteristics of the observation system (wave imaging marine radar). For verification purposes, a substantial set of synthetic wave data (i.e. forward runs of the wave model) were generated to be used as ground truth for comparison to the reconstructions and forecasts produced by Wavecast. For these synthetic cases, Wavecast demonstrated very good accuracy, for example, typical forecast correlation coefficients were between 0.84-0.95 when compared to the input data. Dependencies on shadowing, observational noise, and forecast horizon were also identified. During the second year of the project, a short field deployment was conducted in order to assess forecast accuracy under field conditions. For this, a radar was installed on a fishing vessel and observations were collected at the South Energy Test Site (SETS) off the coast of Newport, OR. At the SETS site, simultaneous in situ wave observations were also available owing to an ongoing field project funded separately. Unfortunately, the position and heading information that was available for the fishing vessel were not of sufficient accuracy in order to validate the forecast in a phase-resolving sense. Instead, a spectral comparison was made between the Wavecast forecast and the data from the in situ wave buoy. Although the wave and wind conditions during the field test were complex, the comparison showed a promising reconstruction of the wave spectral shape, where both peaks in the bimodal spectrum were represented. However, the total reconstructed spectral energy (across all directions and frequencies) was limited to 44% of the observed spectrum. Overall, wave-by-wave forecasting using a data assimilation approach based on wave imaging radar observations and a physics-based wave model shows promise for short-term phase-resolved predictions. Two recommendations for future work are as follows: first, we would recommend additional focused field campaigns for algorithm validation. The field campaign should be long enough to capture a range of wave conditions relevant to the target application and WEC site. In addition, it will be crucial to make sure the vessel of choice has high accuracy position and heading instrumentation (this instrumentation is commercially available but not standard on commercial fishing vessels). The second recommendation is to expand the model physics in the wave model backbone to include some nonlinear effects. Specifically, the third-order correction to the wave speed due to amplitude dispersion would be the next step in order to more accurately represent the phase speeds of large amplitude waves.

16 TIDAL AND WAVE POWER↗

DOE OSTI2025

Neutron generation dynamics inside a MA-class dense plasma focus Z-pinch

Dense plasma focii (DPFs) are appealing as energy efficient sources of short pulses of ions, neutrons, and x rays. The output of these sources is expected to scale with input current (I4), but has been shown to drop at the MA level [S. K. H. Auluck, “On the failure of neutron yield scaling in the dense plasma focus,” Phys. Plasmas 30, 080701 (2023)]. New results on the MegaJOuLe Neutron Imaging Radiography DPF showed neutron yield production in agreement with the input current scaling beyond the previously observed drop. This work provides insight into the pinch formation on a DPF and reports on the two different mechanisms leading to neutron generation inside a DPF using a combination of kinetic simulations and experimental data. A combination of particle-in-cell (PIC) and 1D shock theory results are used to describe the pinch formation and disassembly and the corresponding thermonuclear and beam-target mechanisms. The temporal evolution of the pinch column predicted by the PIC simulations shows qualitative agreement with the experimental data from plasma photon emission as well as temporal neutron pulse shapes. In MJ-class DPFs, both thermonuclear and beam-target mechanisms can occur over the course of the implosion and contribute to the total neutron production. Hence the neutron source size of a DPF will change throughout the implosion. Experimental neutron radiographs show the increase in source size as the pinch breaks apart, in agreement with simulation's prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Modeling, measurement, and 3-D equilibrium reconstruction of the bootstrap current in the Helically Symmetric Experiment

The bootstrap current for three electron cyclotron resonance heated plasma scenarios in a quasihelically symmetric stellarator (the Helically Symmetric Experiment) are analyzed and compared to a neoclassical transport code PENTA. The three conditions correspond to 50 kW input power with a resonance that is off-axis, 50 kW on-axis heating and 100 kW on-axis heating. When the heating location was moved from off-axis to on-axis with 50 kW heating power, the stored energy and the extrapolated steady-state current were both observed to increase. When the on-axis heating power was increased from 50 kW to 100 kW, the stored energy continued to increase while the bootstrap current slightly decreased. This trend is qualitatively in agreement with the calculations which indicate that a large positive electric field for the 100 kW case was driving the current negative in a small region close to the magnetic axis and accounting for the decrease in the total integrated current. This trend in the calculations is only observed to occur when momentum conservation between particle species is included. Without momentum conservation, the calculated bootstrap current increases monotonically. We show that the magnitude of the bootstrap current as calculated by PENTA agrees better with the experiment when momentum conservation between plasma species is included in the calculation. The total current was observed in all cases to flow in a direction to unwind the transform, unlike in a tokamak in which the bootstrap current adds to the transform. The 3-D inductive response of the plasma is simulated to predict the evolution of the current profile during the discharge. The 3-D equilibrium reconstruction code V3FIT is used to reconstruct profiles of the plasma pressure and current constrained by measurements with a set of magnetic diagnostics. The reconstructed profiles are consistent with the measured plasma pressure profile and the simulated current profile when the reconstruction is constrained by the measured data from a diagnostic array that is internal to the vacuum chamber.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Programmable Phase Selection between Altermagnetic and Noncentrosymmetric Polymorphs of MnTe on InP via Molecular Beam Epitaxy

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (ACS Appl. Mater. Interfaces 2026, 18, 15654-15664). The calculations model the interfacial energetics of two MnTe polymorphs — NiAs-MnTe (hexagonal, alpha phase) and ZnS-MnTe (cubic, gamma phase) — on InP(111) substrates with two surface terminations: In-terminated InP(111)A and P-terminated InP(111)B. This gives four interface configurations: NiAs on In-terminated (experimentally observed), NiAs on P-terminated (computed for comparison), ZnS on In-terminated (computed for comparison), and ZnS on P-terminated (experimentally observed). The dataset is organized into four calculation types, each covering all four polymorph/termination combinations: (i) Slabs: Pristine MnTe/InP heterostructure slabs used to compute total energies and interface energy densities (Eint) for all four configurations, as reported in Fig. 6 of the main text. (ii) Disorder: Same slab geometries with a P_Te + Te_P antisite defect pair introduced near the interface, used to assess chemical intermixing effects on interface stability (Fig. S8, SI). (iii) Strain: Pristine slab calculations with in-plane lattice parameters strained by -1% and +1% relative to the InP lattice constant, used to evaluate strain-dependent interface energetics (Fig. S9, SI). (iv) Charge_Density: Single-point calculations on the full heterostructure, the isolated InP slab, and the isolated MnTe slab at fixed geometry, used to compute differential charge density plots showing interfacial charge accumulation and depletion as a function of surface termination (Fig. S10, SI). Each calculation folder contains INCAR, KPOINTS, POSCAR, CONTCAR, OUTCAR, and POTCAR_info.txt (PAW potential information, excluding the full POTCAR due to VASP licensing restrictions). The calculations were performed using VASP 6.4.3 with PBE exchange-correlation, PAW potentials, a Hubbard correction of Ueff = 5 eV on Mn d-states, and A-type AFM spin initialization.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL Radiological Response Training Range

This assessment uses screening-level models to calculate potential environmental impacts from proposed tests at the Idaho National Laboratory (INL) Radiological Response Training Range (RRTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K2O, LaBr3, KBr, Cu, Zr, F, Ga, Ga2O3, NaNO2, Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface-soil impacts were calculated with the computer program mixing-cell model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated assuming 12 tests are conducted annually using all 11 material types for a period of 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer-code input, output, executable files, batch files, scripts, and spreadsheet files are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS). It is highly unlikely the test scenarios evaluated in this ECAR will adversely impact human health based on comparisons of calculated dose and concentration against regulatory standards and risk-based screening levels. Conservative estimates of dose to workers and the public from atmospheric transport of possible radionuclide releases are far below federal radiation protection standards. Conservative estimates of potential contaminant concentrations in groundwater are less than federal drinking water standards or screening levels. Predicted radionuclide concentrations in surface soils are below risk-based screening levels, except for Ge-68 (material Ga-68) for the worker. The Ge-68 soil concentration can be made less than the worker PRG, if the number of annual tests using Ga-68 is reduced from 12 to 6. However, the sum of ratios still exceeds one because of the high K-40 ratio. If the EF of the worker (number of days the worker is in the contaminated testing area) is reduced from 225 days/yr (default value for full time worker) to 112 days/yr, the Ge-68 ratio is less than one and the sum of ratios is less than one. Actual radiation doses and groundwater and surface-soil concentrations are likely to be much less than those calculated because of the conservative assumptions and parameters employed in the modeling. For example, atmospheric-transport calculations assume the entire inventory of each material type is readily released to the atmosphere and no plume deposition, depletion, or radioactive decay occurs during transport. The calculations also assume the same meteorological conditions (e.g., wind velocity, wind direction, stability class) that produce the maximum 95th percentile concentration (i.e., concentration representing the 95th percentile of a distribution of concentrations derived from 3 years of hourly meteorological data) at each receptor location are the same for all 12 tests during the year, and each receptor is assumed to be present during all 12 tests. The surface-soil assessment assumes the entire inventory of each test is deposited in the top 5 cm of soil. No atmospheric dispersal is assumed, and the radionuclides are subject only to leaching and radioactive decay. The groundwater-pathway modeling is conservative in that it is one-dimensional in the unsaturated zone (no lateral spreading/dilution) and assumes the entire inventory of contaminants infiltrates into the ground at the same location for every test. This is especially conservative for particulate radionuclides because they would have to dissolve or corrode first and some would be dispersed into the atmosphere. The groundwater receptor is also assumed to consume water directly from a hypothetical well positioned in the location of maximum concentration. In addition, conservative degradation rates were used, and volatilization was not considered for the nonradioactive chemical

54 ENVIRONMENTAL SCIENCES↗

NASA NTRSDate not supplied

Modeling and Simulation of Tank Pressure Control using Zero-Boil Off Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

Cameron J. Hines↗