Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Metric geometry”

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

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Geometric aspects in digital analysis of Multi-Spectral Scanner (MSS) data

Present automated systems of interpretation which apply pattern recognition techniques on MSS data do not fully consider the geometry of the acquisition system. In an effort to improve the usefulness of the MSS data when digitally treated, geometric aspects are analyzed and discussed. Attempts to correct for scanner instabilities in position and orientation by affine and polynomial transformations, as well as by modified collinearity equations are described. Methods of accounting for panoramic and relief effects are also discussed. It is anticipated that reliable area as well as position determinations can be accomplished during the process of automatic interpretation. A concept for a unified approach to the treatment of remote sensing data, both metric and nonmetric is presented.

Mikhail, E. M.↗
DOE OSTI2024

Evaluating the Importance of Conformers for Understanding the Vacuum-Ultraviolet Spectra of Oxiranes: Experiment and Theory

Vacuum-ultraviolet (VUV) absorption spectroscopy enables electronic transitions that offer the unambiguous identification of molecules. As target molecules become more complex, multifunctional species present a great challenge to both experimental and computational spectroscopy. This research reports both experimental and theoretical studies of oxiranes. Computationally, the nuclear ensemble approach has been used to accurately predict experimental spectra for a variety of molecules. However, this approach incurs great computational cost, as ensembles generally consist of thousands of geometries. The present study aims to drastically reduce the ensemble by evaluating the significance of the conformers to the predicted spectra. This approach was applied to 11 substituted oxiranes using the Conformer Rotamer Ensemble Sampling Tool (CREST) of Grimme to generate an ensemble of unique conformers determined by their Boltzmann populations. Five TD-DFT functionals (BMK, CAM-B3LYP, M06-2X, MN15, ωB97X-D) and EOM-CCSD were used to simulate the spectrum of each substituted oxirane ensemble. Computed spectra were then compared to the experiment using both qualitative and quantitative metrics. Based on these metrics, it was observed that certain conformers may not be necessary to characterize this set of oxiranes despite the temperature (323 K) of the experiment. A single conformer can then be used with TD-DFT and EOM-CCSD to replicate the experimental spectra of these medium-sized combustion species.

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

Minimizing Ground Risk in Cellular-Connected Drone Corridors With mmWave Links

Unmanned Aircraft Systems (UASs) have been receiving significant interest and support from academia, industry, and regulatory bodies over the past decade due to their various use cases. To safely integrate UAS operations into the national airspace, particularly overpopulated regions, the risk posed to ground users, buildings, and vehicles due to unmanned aerial vehicle (UAV) flight should be minimized. This risk can be represented by a numerical metric, which we refer to in this article as the “ground risk.” Many UAS applications also depend on the presence of a reliable wireless communication link between the UAV and a control station for the transmission of UAV position, surveillance video, UAV payload commands, and other mission-related data. Such wireless communication requirements also need to be considered in the design of UAS operations. In this article, we consider both these aspects and study the design of nonintersecting trajectories for UAS operations to minimize ground risk, subject to constraints on the wireless signal strength and geometry of the trajectory, specified in terms of: 1) an enclosing cylinder within which the trajectory must lie and 2) an integrated angular change along the UAV's trajectory. The performance of a computationally expensive optimal algorithm is compared with that of a computationally faster heuristic approach within the dense urban environment of Manhattan, NY, USA. Performance evaluation using ray-tracing simulations shows that the heuristic approach performs close to the optimal algorithm at a reduced computation cost. In conclusion, this research can be utilized to make UAS operations safe and reliable and accelerate their adoption.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Soft particles and infinite-dimensional geometry

Abstract In the sigma model, soft insertions of moduli scalars enact parallel transport ofS-matrix elements about the finite-dimensional moduli space of vacua, and the antisymmetric double-soft theorem calculates the curvature of the vacuum manifold. We explore the analogs of these statements in gauge theory and gravity in asymptotically flat spacetimes, where the relevant moduli spaces are infinite-dimensional. These models have spaces of vacua parameterized by (trivial) flat connections on the celestial sphere, and soft insertions of photons, gluons, and gravitons parallel transportS-matrix elements about these infinite-dimensional manifolds. We argue that the antisymmetric double-soft gluon theorem ind + 2 bulk dimensions computes the curvature of a connection on the infinite-dimensional space Map ( S d , G ) / G , whereGis the global part of the gauge group. The analogous metrics in abelian gauge theory and gravity are flat, as indicated by the vanishing of the antisymmetric double-soft theorems in those models. In other words, Feynman diagram calculations not only know about the vacuum manifold of Yang–Mills theory, they can also be used to compute its curvature. The results have interesting implications for flat space holography.

Astronomy & Astrophysics↗
DOE OSTI2021

Ab-Initio Simulation Studies of Cr Solvation in Fluoride Molten Salts

Understanding molten salt chemistry is essential in ongoing research of the molten salt nuclear reactor (MSR). In this context, detailed understanding of the mechanisms underlying selective oxidation of metal species, such as Cr, is required to guide the design of effective corrosion mitigation strategies in molten salts. An important starting point for such mechanistic understanding is knowledge of the solvation structure and its role in controlling metal speciation. In this work, we use ab initio molecular dynamics simulations to study the short-range (on the scale of the nearest-neighbor bond lengths) and medium-range (over length scales of several neighbor spacings) structure in three different fluoride melts with and without Cr addition; namely, 2KF-NaF, 2LiF-BeF 2 , and 3LiF-AlF3. We find that Cr 0 ,Cr 2+ ,Cr 3+ can each be coordinated by different numbers of F-, with the variance in coordination number decreasing as oxidation state increases, and that these coordination geometries are largely independent of solvent. The manner by which Cr changes the medium-range structure, however, is found to be solvent-dependent. While 2KF-NaF melts show short and medium range order that is highly dynamic, 2LiF-BeF 2 and 3LiF-AlF 3 are characterized by molecular associates that are relatively long-lived that organize into oligomer structures on larger length scales. Rather than being solvated by F- ions alone, we find that Cr can incorporate into and be solvated within this oligomer structure. Fluoroacidity, alone, may therefore prove too simple a metric for assessing the corrosivity of molten fluorides. As our work suggests, the ability of Cr to solvate must be understood in the context of the short- and medium-range structure of the solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2006

Aircraft Conceptual Design and Risk Analysis Using Physics-Based Noise Prediction

An approach was developed which allows for design studies of commercial aircraft using physics-based noise analysis methods while retaining the ability to perform the rapid trade-off and risk analysis studies needed at the conceptual design stage. A prototype integrated analysis process was created for computing the total aircraft EPNL at the Federal Aviation Regulations Part 36 certification measurement locations using physics-based methods for fan rotor-stator interaction tones and jet mixing noise. The methodology was then used in combination with design of experiments to create response surface equations (RSEs) for the engine and aircraft performance metrics, geometric constraints and take-off and landing noise levels. In addition, Monte Carlo analysis was used to assess the expected variability of the metrics under the influence of uncertainty, and to determine how the variability is affected by the choice of engine cycle. Finally, the RSEs were used to conduct a series of proof-of-concept conceptual-level design studies demonstrating the utility of the approach. The study found that a key advantage to using physics-based analysis during conceptual design lies in the ability to assess the benefits of new technologies as a function of the design to which they are applied. The greatest difficulty in implementing physics-based analysis proved to be the generation of design geometry at a sufficient level of detail for high-fidelity analysis.

Olson, Erik D.↗
DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗
DOE OSTI2023

Estimate of Gamma Dose Rates from Arrays of Fermi-1 Blanket Elements During the MEDE Process

The Enrico Fermi Atomic Power Plant Unit 1 (Fermi-1) was a sodium-cooled fast breeder reactor located in Monroe County, Michigan. The reactor was powered by a core of enriched uranium metal alloy driver fuel, which was enveloped by an axial and radial blanket material consisting of depleted uranium metal alloyed with 3 wt.% molybdenum. There are 406 axial and 559 radial irradiated sodium-bonded Fermi-1 blanket assemblies in storage at INL, totaling 34 metric tons of heavy metal. Disposal of the Fermi-1 blanket material directly into a geological repository is prohibited due to the reactive characteristic of its bond sodium. A melt drain evaporate (MEDE) process can effectively remove bond sodium from Fermi-1 blanket material. Consequently, planning is underway to apply a MEDE process to treat the 34 metric tons of heavy metal of irradiated Fermi-1 blanket material. Given the irradiated Fermi-1 blanket material’s relatively low power history and 50+ years of decay, it can be treated in a shielded glovebox. To assess the requisite shielding, the dose rates of the Fermi-1 blanket elements in various process configurations are needed. The method to perform the dose rate calculations is to first generate an average source term and associated photon source spectra for the Fermi-1 blanket material in SCALE and then use the associated photon spectra to calculate dose rates using MCNP6.2 in seven representative geometries. Dose rates in rem/h were calculated on contact (1 cm from outer geometry surface), 30 cm away from outer geometry surface, and 1 m from outer geometry surface at axial heights spanning the length of the blanket material within the geometry. The maximum average dose rate for the single Fermi-1 radial blanket element is ~0.60 rem/h, for the Fermi-1 radial blanket assembly ~1.10 rem/h, for the single Fermi-1 axial blanket element ~0.5 rem/h, for the Fermi-1 axial blanket assembly ~1.01 rem/h, for the MEDE can ~ 0.98 rem/h, for the Fermi Storage Canister ~ 0.83 rem/h, and for the MEDE cans in the DOE Standard Canister ~0.51 rem/h. Based on these maximum average values, the bounding dose rates are assumed to be 3× the average dose rates calculated for the axial blanket material and 7.6× the average dose rates calculated for the radial blanket material, based on the distribution of Fermi-1 blanket material assembly burnup.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2010

MCNP Simulations of Measurement of Insulation Compaction in the Cryogenic Rocket Fuel Tanks at Kennedy Space Center by Fast/Thermal Neutron Techniques

MCNP simulations have been run to evaluate the feasibility of using a combination of fast and thermal neutrons as a nondestructive method to measure of the compaction of the perlite insulation in the liquid hydrogen and oxygen cryogenic storage tanks at John F. Kennedy Space Center (KSC). Perlite is a feldspathic volcanic rock made up of the major elements Si, AI, Na, K and 0 along with some water. When heated it expands from four to twenty times its original volume which makes it very useful for thermal insulation. The cryogenic tanks at Kennedy Space Center are spherical with outer diameters of 69-70 feet and lined with a layer of expanded perlite with thicknesses on the order of 120 cm. There is evidence that some of the perlite has compacted over time since the tanks were built 1965, affecting the thermal properties and possibly also the structural integrity of the tanks. With commercially available portable neutron generators it is possible to produce simultaneously fluxes of neutrons in two energy ranges: fast (14 Me V) and thermal (25 me V). The two energy ranges produce complementary information. Fast neutrons produce gamma rays by inelastic scattering, which is sensitive to Fe and O. Thermal neutrons produce gamma rays by prompt gamma neutron activation (PGNA) and this is sensitive to Si, Al, Na, K and H. The compaction of the perlite can be measured by the change in gamma ray signal strength which is proportional to the atomic number densities of the constituent elements. The MCNP simulations were made to determine the magnitude of this change. The tank wall was approximated by a I-dimensional slab geometry with an 11/16" outer carbon steel wall, an inner stainless wall and 120 cm thick perlite zone. Runs were made for cases with expanded perlite, compacted perlite or with various void fractions. Runs were also made to simulate the effect of adding a moderator. Tallies were made for decay-time analysis from t=0 to 10 ms; total detected gamma-rays; detected gamma-rays from thermal neutron reactions d. detected gamma-rays from non-thermal neutron reactions and total detected gamma-rays as a function of depth into the annulus volume. These indicated a number of possible independent metrics of perlite compaction. For example the count rate for perlite elements increased from 3600 to 8500 cps for an increase in perlite density from 6 lbs/lcf to 16.5 lbs/cf. Thus the MCNP simulations have confirmed the feasibility of using neutron methods to map the compaction of perlite in the walls of the cryogenic tanks.

Livingston, R. A.↗
NASA NTRSDate not supplied

Evaluating Retrieval Algorithm Climate Stability: Estimating 3D Optical Thickness Bias Distributions by Cloud Type

Detecting climate trends on large spatiotemporal scales requires accurate, stable measurements and stable retrieval algorithms. We strive to estimate how time-variant retrieval algorithm biases may impact trend detection. Here we focus on the 3D cloud optical thickness (τc) bias, which is among the largest in passive cloud retrieval algorithms. If this bias is time dependent, a possibility with potential decadal changes in cloud morphology, it may obscure genuine trends in τc. Although previous studies have evaluated the cloud- and sun-view geometry-dependent 3D τc bias on small spatial scales, before our current study none have evaluated the stability of this well-known bias on climate-relevant large spatiotemporal scales. These studies must estimate large scale distributions of the 3D τc bias by cloud type and estimate how cloud type amount may change between two climate states. We employ a novel approach to estimate large scale distributions of 3D τc using a proxy of the bias that quantifies the departure of clouds from satisfying the 1D radiative transfer assumption used in passive τc retrievals. This existing globally-distributed proxy is an angular consistency metric that was developed using fused Moderate-Resolution Imaging Spectroradiometer (MODIS) and Multi-angle Imaging Spectroradiometer (MISR) measurements. Calculating the 3D τc bias and the proxy, for known cloud fields enables us to establish statistical relationships between these two quantities, which can be used to calculate large-scale distributions of the 3D τc bias. This approach limits the number of 3D radiative transfer simulations required to only those needed to estimate a statistical relationship between the 3D τc bias for known cloud fields and a proxy of the bias. It is likely that future studies will be needed to evaluate retrieval algorithm bias stability for other geophysical variables as the community develops climate data records from satellite observations and their retrievals. This must be done in addition to monitoring and correcting measurement errors and uncertainties and understanding their impact on retrieved essential climate variables.

Yolanda Shea↗
NASA NTRSDate not supplied

Aerodynamic Design Optimization for Natural Laminar Flow Airfoils

Natural laminar flow technology is a passive laminar flow control (LFC) strategy that seeks to delay the onset of boundary-layer transition (BLT) through shape optimization to reduce the drag of the aerodynamic vehicle. Adjoint-based design optimization for LFC is proposed in an integrated multidisciplinary framework, which includes the computational fluid dynamics (CFD), geometry and grid deformation, and linear stability analysis (LSA) for transition prediction. In particular, the BLT location is predicted using the dual N-factor method that is based on a linear stability theory (LST) eigenvalue problem. The dual N-factor criterion accounts for the amplification of planar Tollmien-Schlichting (TS) and stationary crossflow (CF) boundary-layer instabilities to predict the transition location in three-dimensional boundary-layer flows. The adjoint-based shape optimization procedure is based on an iteratively coupled CFD and LSA methodology to converge the transition location and flow solutions, as well as to calculate the sensitivities of the aerodynamic metrics of interest with respect to the flow and shape design parameters. The RAE 2822 airfoil at 0 and 30 degrees yaw angles, an angle of attack of 0.72 degrees, and subsonic conditions (M∞ = 0.19, Rec = 5.6 × 106 ) are used as baseline configurations for design optimization. The angle of attack and the vertical displacement of free-form-deformation control points are used as design variables to reduce the drag coefficient while reaching a specified lift coefficient. The optimized unswept airfoil designs achieve a 30% drag reduction accompanied by a downstream shift of the transition locations over both suction and pressure sides of the airfoil. The initial design iterations for the swept case also show a favorable trend in the drag reduction with transition delay over both sides.

Transition↗
DOE OSTI2024

Design and Manufacturing of High Performance, Reduced Charge Heat Exchangers

Air-to-fluid heat exchangers (HXs) play a critical role as the main heat transfer component in Heating, Ventilation, Air-Conditioning, and Refrigeration (HVAC&R) systems. However, their airside thermal resistance significantly inhibits their overall performance. Furthermore, these HXs must be continually more compact to meet the latest refrigerant charge limits to reduce emissions. Recent literature suggests that traditional HX geometries (e.g., round or flat tubes with fins), have reached their limits, and more sophisticated shape- and topology-optimized designs are required to achieve the next jump in performance. This research sheds light on the next generation of air-to-refrigerant HXs and aims to address several practical issues to commercialization such as novelty challenges (improved performance for significant charge reduction; modeling expertise & time investment), manufacturing challenges (non-round tube manufacturing; tube-header integration; product qualification, e.g., burst pressure testing, extreme operational environment, etc.), and operational challenges (flow maldistribution, fouling & wetting, noise & vibration). For example, a >20% improvement on one (or more) HX-level performance metrics (e.g., envelope volume, airside pressure drop, face area, capacity, refrigerant charge, weight, cost, etc.) must be achieved before a HX design is considered for commercialization. We present a new, comprehensive and experimentally validated air-to-refrigerant HX optimization framework with simultaneous thermal-hydraulic performance and mechanical strength considerations for novel, non-round, shape- and topology-optimized tubes capable of optimizing single and two-phase HX designs for any refrigerant choice and performance requirement with significant engineering time savings compared to conventional design practices. The framework was exercised for a wide range of applications and refrigerants, resulting in HXs which achieved greater than 20% improved performance, 20% reductions in size, and 25% reductions in refrigerant charge. To enable non-round tube bundle use in next generation HVAC&R equipment, novel manufacturing techniques were investigated, including the development of conventional manufacturing methods for small diameter, non-round tubes and novel tube-header integration strategies. In total, ten HX prototypes were manufactured, nine using conventional methods directly attributed to this project and one using advanced additive manufacturing methods. The five-year manufacturing feasibility of the proposed HXs was found to have a good outlook. The non-round tube HX simulated performance was validated through comprehensive experimental testing, including nine in-house component-level tests, one independent component-level test at an industry partner laboratory, and in-house system-level tests of using a commercially-available, residential packaged A/C unit which was retrofitted with a non-round tube prototype HX. It was found that HX designs proposed by the new framework can successfully predict experimental thermal-hydraulic performance within ±10-20% the first time with no manual design changes, eliminating the need for time-consuming and expensive prototyping efforts. This work will accelerate design and time to market for next generation HXs while simultaneously facilitating industry transition to new refrigerants at lower charge.

42 ENGINEERING↗
DOE OSTI2024

Interlaced Characterization and Calibration (ICC) for Improved Computational Simulation Credibility

Accurate material characterization and model calibration are pivotal for simulations used for high-consequence engineering decisions. Current characterization and calibration methods (1) use simplified test specimen geometries and global data, (2) cannot guarantee that sufficient characterization data is collected for a specific model of interest, (3) provide only mean parameter values with no uncertainty quantification, and (4) are sequential, inflexible, and time-consuming. This work developed a new paradigm—coined Interlaced Characterization and Calibration (ICC)—which drives forward the state-of-the-art in model calibration by bringing together recent advancements into one improved workflow. The ICC paradigm (1) employs tools to efficiently use full-field data to calibrate high-fidelity material models, (2) aligns the data needed with the data collected by adopting an optimal experimental design protocol, (3) provides uncertainty metrics on the calibrated model parameters, and (4) incorporates these advances into a quasi real-time feedback loop. The ICC framework was validated synthetically with both low-fidelity and high-fidelity simulations paired with several different elastoplastic material models, and was also demonstrated experimentally with an aluminum 6061 cruciform exemplar specimen. Results showed that the ICC framework—in which Bayesian optimal experimental design actively guided the experiment— resulted in calibrations with similar or better accuracy than predetermined experiments based on subject matter expertise. Moreover, the ICC framework produced a complete model calibration— with quantified uncertainties on model parameters—in 1 week, a 5 - 10× increase in efficiency over traditional approaches. Thus, the ICC paradigm improves both the calibration process and quality, by (1) improving efficiency, which increases agility of solid mechanics modeling and enables utilization of computational simulation (CompSim) at earlier stages of the design cycle and (2) providing quantified, and in some cases reduced, parameter uncertainties, which increases confidence in model predictions and supports credible decision making.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Triple interference, non-linear Talbot effect and gravitization of the quantum

Recently we have discussed a new approach to the problem of quantum gravity in which the quantum mechanical structures that are traditionally fixed, such as the Fubini–Study metric in the Hilbert space of states, become dynamical and so implement the idea of gravitizing the quantum. Here, in this paper we elaborate on a specific test of this new approach to quantum gravity using triple interference in a varying gravitational field. Our discussion is driven by a profound analogy with recent triple-path interference experiments performed in the context of non-linear optics. We emphasize that the triple interference experiment in a varying gravitational field would deeply influence the present understanding of the kinematics of quantum gravity and quantum gravity phenomenology. We also discuss the non-linear Talbot effect as another striking phenomenological probe of gravitization of the geometry of quantum theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Solar Field Layout and Aimpoint Strategy Optimization

The existing methods that determine heliostat aiming strategies for concentrating solar power (CSP) central receiver plants typically use heuristics and/or are computationally expensive, and they lack flexibility for different desired flux profiles and receiver geometries. Because of the interaction between layout and aimpoint strategy, considering the former without accounting for the latter may yield solutions with superfluous heliostats that cannot be used efficiently without compromising receiver flux constraints. To that end, we develop a software decision tool that uses innovative optimization methods to both optimize aimpoint strategies and improve candidate layouts for the solar collection field of a CSP central receiver plant. A CSP plant’s effectiveness relies on the optical efficiency of the solar field, which may be limited by losses due to (i) the cosine effect, (ii) atmospheric attenuation, (iii) interference (i.e., shading and blocking) between heliostats, (iv) spillage as a result of heliostat positioning and geometry, and (iv) some heliostats’ inability to direct irradiance to the receiver without damage due to excessive thermal flux. The goal of this work is to obtain optimized aiming strategies and improved solar field layouts that reduce capital cost and increase field optical efficiency and utilization, while meeting the power requirements of a given CSP receiver design. We formulate the aimpoint optimization problem as a mixed-integer linear programming model, which we then decompose into submodels that we solve in parallel. The decomposition subdivides the solar field into sections, and aimpoint strategies for each section are obtained independently of the others. To improve existing layouts, we develop a utilization-weighted efficiency metric that we use to relocate heliostats to sections of the solar field with similar efficiency and higher utilization. Finally, to connect our software to high-fidelity flux models, we develop a Python application programming interface for SolarPILOT, a mature software package that characterizes solar field performance and generates the heliostat layouts and flux maps that serve as input to our models.

14 SOLAR ENERGY↗
DOE OSTI2025

Challenge Problem 1: Preliminary Results of the Direct Numerical Simulation of Transient Flows

This report presents the first direct numerical simulations (DNS) of transient mixed convection in an idealized downcomer-like channel (Challenge Problem 1, Phase II). Using the GPU-accelerated NekRS solver, we modeled a sudden decay in driving pressure, mimicking loss-of-flow events, and tracked the resulting evolution of Reynolds number, boundary-layer structure, turbulence statistics, and heat-transfer metrics. Key findings include the systematic thickening and eventual asymmetry of velocity and thermal boundary layers under buoyant deceleration; minimal “memory” lag in Reynolds shear stress and TKE profiles when sampled at matching Re, yet clear shifts of peak locations toward the cooled wall; overshoots in transient eddy-viscosity and eddy-diffusivity (and corresponding sub-unity turbulent Prandtl numbers) on the cooled side; and a pronounced transient Nusselt-number enhancement driven by wall-temperature inertia and residual eddy mixing. These effects combined to offer a temporary cooling margin above steady-state predictions during reactor LOF transients. Future work will extend this work to a more complex “Case II” geometry (90° turn + lower plenum) and generate multi-Re/Pr datasets for data-driven turbulence closures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

A computational modeling framework for pre-clinical evaluation of cardiac mapping systems

There are a variety of difficulties in evaluating clinical cardiac mapping systems, most notably the inability to record the transmembrane potential throughout the entire heart during patient procedures which prevents the comparison to a relevant “gold standard”. Cardiac mapping systems are comprised of hardware and software elements including sophisticated mathematical algorithms, both of which continue to undergo rapid innovation. The purpose of this study is to develop a computational modeling framework to evaluate the performance of cardiac mapping systems. The framework enables rigorous evaluation of a mapping system’s ability to localize and characterize (i.e., focal or reentrant) arrhythmogenic sources in the heart. The main component of our tool is a library of computer simulations of various dynamic patterns throughout the entire heart in which the type and location of the arrhythmogenic sources are known. Our framework allows for performance evaluation for various electrode configurations, heart geometries, arrhythmias, and electrogram noise levels and involves blind comparison of mapping systems against a “silver standard” comprised of computer simulations in which the precise transmembrane potential patterns throughout the heart are known. A feasibility study was performed using simulations of patterns in the human left atria and three hypothetical virtual catheter electrode arrays. Activation times (AcT) and patterns (AcP) were computed for three virtual electrode arrays: two basket arrays with good and poor contact and one high-resolution grid with uniform spacing. The average root mean squared difference of AcTs of electrograms and those of the nearest endocardial action potential was less than 1 ms and therefore appears to be a poor performance metric. In an effort to standardize performance evaluation of mapping systems a novel performance metric is introduced based on the number of AcPs identified correctly and those considered spurious as well as misclassifications of arrhythmia type; spatial and temporal localization accuracy of correctly identified patterns was also quantified. This approach provides a rigorous quantitative analysis of cardiac mapping system performance. Proof of concept of this computational evaluation framework suggests that it could help safeguard that mapping systems perform as expected as well as provide estimates of system accuracy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Methods for Estimating Hydrogen Fuel Tank Characteristics

The pressure vessels needed to store hydrogen for next-generation hydrogen fuel cell vehicles are expected to be a substantial portion of the total system mass, volume, and cost. Gravimetric capacity, volumetric capacity, and cost per kilogram of usable hydrogen are key performance metrics that the U.S. Department of Energy (DOE) uses to determine the viability of hydrogen fuel cell systems. Research and development related to hydrogen storage systems covers a wide range of potential operating conditions, from cryogenic temperatures to high temperatures (above ambient) and low pressure to high pressure. Researchers at PNNL have developed methods for estimating these key pressure vessel characteristics to support on-board hydrogen storage system design and performance evaluation and to support decision-making about DOE hydrogen storage system research investments. This article describes the pressure tank estimation methodology that has been used as a stand-alone calculation and has been incorporated into larger system evaluation tools. The methodology estimates the geometry, mass, and material cost of type I, type III, and type IV pressure vessels based on operating pressure and material strength at the system's operating temperature, using classical thin-wall and thick-wall pressure vessel stress calculations. The geometry, mass, and material cost requirements of the pressure vessel have significant impacts on the total system performance. For example, hydrogen storage materials that can separately achieve a very high hydrogen density can be deemed impractical for use in fuel cell vehicle hydrogen storage systems because the pressure tank containing them is too large, heavy, or expensive. This article describes the design philosophy and analytical process of the tank characteristic estimation methodology, which has been implemented in spreadsheet calculation tools and system-level analysis tools used by DOE researchers. Each of the three tank types (type I, type III, and type IV) uses a different analysis methodology with some common elements. This article also provides examples of implementing the methodology to perform parametric studies of all three pressure vessel types. The goal of this article is to present the methodology in sufficient detail so it can be implemented in other hydrogen fuel cell vehicle design and analysis tools.

42 ENGINEERING↗