Results for “Design Methods Modeling and Validation” — Page 51

DOE OSTI2022

Performance evaluation of a wrapped around condenser for heat pump water heater applications

In this report, accurate performance evaluation of a wrapped-around condenser for heat pump water heaters (WHPWH) is critical since the COP of the system depends heavily on thermal stratification and condenser design. An analytical, quasi-steady state heat balance method has been developed to determine the optimal spacing between adjacent channels, tube diameter, tube shape, and total refrigerant charge amount for the condenser. The heat transfer rate is compared among three approaches: the refrigerant thermodynamic model, the condenser-wrap fin model, and the analytical natural convection model on the inside of the tank for the three regions based on the refrigerant phase (e.g., super-heated, saturated, or subcooled). The heat transfer rate was predicted by a combination of slug-plug and annular flow for a D-shaped helical tube. Experimental data were used as boundary conditions for validating the model. The rate of heat transfer based on tube shape and tank wall temperatures was compared by CFD analysis. Parametric analysis indicates a tradeoff between refrigerant mass, pressure drop, tube diameter, and tube length to maximize the heat transfer rate. The model suggests the saturated region length can be extended by 400%, and the condenser pressure drop can be reduced by 23% with an optimal spacing pattern.

42 ENGINEERING↗

DOE OSTI2018

Flaw characterization and correlation with cracking strength in Engineered Cementitious Composites (ECC)

Engineered Cementitious Composite (ECC) is a class of fiber reinforced composites showing strain hardening behavior. The variation of cracking strength among various sections of a tensile specimen is a key factor governing the cracking behaviors of ECC. This study employs X-ray computed tomography method to investigate the correlation between cracking strength and flaw distribution in ECC, and identifies the dimensions of pre-existing flaws to be the main influencing parameters. Specifically, the largest cross-sectional area of a flaw perpendicular to stress can be well correlated with cracking strength. To validate the observation, the classic model on first crack strength is modified in this study by adopting a refined fiber bridging relation and employing iterated crack profile. The simulated cracking strength vs flaw size relation is in good agreement with test results. This study provides an improved understanding of the multiple cracking mechanism in ECC, which is useful for material design.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Computing resonant modes of accelerator cavities by solving nonlinear eigenvalue problems via rational approximation

Here, we present an efficient and reliable algorithm for solving a class of nonlinear eigenvalue problems arising from the modeling of particle accelerator cavities. The eigenvalue nonlinearity in these problems results from the use of waveguides to couple external power sources or to allow certain excited electromagnetic modes to exit the cavity. We use a rational approximation to reduce the nonlinear eigenvalue problem first to a rational eigenvalue problem. We then apply a special linearization procedure to turn the rational eigenvalue problem into a larger linear eigenvalue problem with the same eigenvalues, which can be solved by existing iterative methods. By using a compact scheme to represent both the linearized operator and the eigenvectors to be computed, we obtain a numerical method that only involves solving linear systems of equations of the same dimension as the original nonlinear eigenvalue problem. We refer to this method as a compact rational Krylov (CORK) method. We implemented the CORK method in the Omega3P module of the Advanced Computational Electromagnetic 3D Parallel (ACE3P) simulation suite and validated it by comparing the computed cavity resonant frequencies and damping Q factors of a small model problem to those obtained from a fitting procedure that uses frequency responses computed by another ACE3P module called S3P. We also used the CORK method to compute trapped modes damped in an ideal eight 9-cell SRF cavity cryomodule. This was the first time it was possible to compute these modes directly. The damping Q factors of the computed modes match well with those measured in experiments and the difference in resonant frequencies is within the range introduced by cavity imperfection. Therefore, the CORK method is an extremely valuable tool for computational cavity design.

43 PARTICLE ACCELERATORS↗

NASA NTRS1992

Study of natural convection and interface shape in directional solidification of succinonitrile

Flow visualization experiments in a Bridgman furnace at zero growth velocity have been performed. The experiments were intended to investigate if the photochromic dye activation method could be used in a crystal growth study. The results from this work have confirmed that with a carefully designed experimental setup, the photochromic-dye method permits qualitative and quantitative evaluation of the flow field near the interface in crystal growth experiments. For the horizontal orientation of a 6 mm x 6 mm ampoule and an axial temperature gradient of 26.5 C/cm, velocity profiles have been obtained accurately at various positions near the interface. The maximum velocity of 1.29 mm/sec was measured in the central vertical plane and the flow was symmetrical about that plane. A flow inversion point was also noted, above which the flow was towards the interface and below it, away from the interface. The results obtained are useful for validating 3-dimensional numerical models and establishing a link between the macroscopic processing conditions and the formation of crystalline defects.

Kawaji, M.↗

NASA NTRS2016

Landing Gear Noise Prediction and Analysis for Tube-and-Wing and Hybrid-Wing-Body Aircraft

Improvements and extensions to landing gear noise prediction methods are developed. New features include installation effects such as reflection from the aircraft, gear truck angle effect, local flow calculation at the landing gear locations, gear size effect, and directivity for various gear designs. These new features have not only significantly improved the accuracy and robustness of the prediction tools, but also have enabled applications to unconventional aircraft designs and installations. Systematic validations of the improved prediction capability are then presented, including parametric validations in functional trends as well as validations in absolute amplitudes, covering a wide variety of landing gear designs, sizes, and testing conditions. The new method is then applied to selected concept aircraft configurations in the portfolio of the NASA Environmentally Responsible Aviation Project envisioned for the timeframe of 2025. The landing gear noise levels are on the order of 2 to 4 dB higher than previously reported predictions due to increased fidelity in accounting for installation effects and gear design details. With the new method, it is now possible to reveal and assess the unique noise characteristics of landing gear systems for each type of aircraft. To address the inevitable uncertainties in predictions of landing gear noise models for future aircraft, an uncertainty analysis is given, using the method of Monte Carlo simulation. The standard deviation of the uncertainty in predicting the absolute level of landing gear noise is quantified and determined to be 1.4 EPNL dB.

Guo, Yueping↗

NASA NTRSDate not supplied

MBSE Applications for the MSR SRC Mars Ascent Vehicle

The objective of the NASA Mars Sample Return (MSR) Campaign is to collect samples from the surface of Mars and return them to Earth for scientific research. The Mars Ascent Vehicle (MAV) will be integrated into a larger Mars Sample Retrieval Lander (SRL) for transit to and storage on Mars. After all Martian samples have been collected and loaded into the MAV payload assembly, MAV will deliver the samples from the Martian surface to Mars orbit. A separate spacecraft, the Earth Return Orbiter (ERO) will retrieve the samples from Mars orbit and return them to Earth. To address common systems engineering challenges associated with using traditional systems engineering practices on complex projects, the MAV systems engineering team has explored implementation of Model-Based Systems Engineering (MBSE) tools and languages. This paper describes the current state of implementation and development of the MAV MBSE model with the Systems Modeling Language (SysML) within the scope of the MAV Systems Requirement Cycle (SRC) systems engineering workflow. The MAV MBSE model has been developed within Magic Draw – a SysML editor commonly used to implement MBSE. The MAV MBSE model has been used to develop mission phase functional flow diagrams for the Concept of Operations, decompose mission to vehicle subsystem functions, develop a functional decomposition, derive functional requirements, trace requirements up to customer-imposed requirements, trace requirements within MAV requirement space, identify requirements trace gaps, define and map the physical design space architecture, allocate requirements to subsystems, develop validation items, define assembly, integration, and test (AI&T) operations, and trace these items across driving goals to develop an integrated digital thread of systems engineering information used to drive design specifications, decision making, and ultimately design verification and validation. Findings and results associated with implementing MBSE in these ways, alongside traditional methods will be discussed.

MBSE↗

NASA NTRSDate not supplied

MBSE Applications for the MSR SRC Mars Ascent Vehicle

The objective of the NASA Mars Sample Return (MSR) Campaign is to collect samples from the surface of Mars and return them to Earth for scientific research. The Mars Ascent Vehicle (MAV) will be integrated into a larger Mars Sample Retrieval Lander (SRL) for transit to and storage on Mars. After all Martian samples have been collected and loaded into the MAV payload assembly, MAV will deliver the samples from the Martian surface to Mars orbit. A separate spacecraft, the Earth Return Orbiter (ERO) will retrieve the samples from Mars orbit and return them to Earth. To address common systems engineering challenges associated with using traditional systems engineering practices on complex projects, the MAV systems engineering team has explored implementation of Model-Based Systems Engineering (MBSE) tools and languages. This paper describes the current state of implementation and development of the MAV MBSE model with the Systems Modeling Language (SysML) within the scope of the MAV Systems Requirement Cycle (SRC) systems engineering workflow. The MAV MBSE model has been developed within Magic Draw – a SysML editor commonly used to implement MBSE. The MAV MBSE model has been used to develop mission phase functional flow diagrams for the Concept of Operations, decompose mission to vehicle subsystem functions, develop a functional decomposition, derive functional requirements, trace requirements up to customer-imposed requirements, trace requirements within MAV requirement space, identify requirements trace gaps, define and map the physical design space architecture, allocate requirements to subsystems, develop validation items, define assembly, integration, and test (AI&T) operations, and trace these items across driving goals to develop an integrated digital thread of systems engineering information used to drive design specifications, decision making, and ultimately design verification and validation. Findings and results associated with implementing MBSE in these ways, alongside traditional methods will be discussed.

MBSE↗

NASA NTRSDate not supplied

FLEET and PLIF Velocimetry Within A Mach 10 Hypersonic Air Flow

Femtosecond laser electronic excitation tagging (FLEET) and planar laser-induced fluorescence (PLIF) velocity measurements utilizing molecular tagging velocity (MTV) methods from three recent test campaigns conducted at the 31-in Mach 10 Air Tunnel at the NASA Langley Research Center are highlighted within. The FLEET measurements reported here include the first direct measurement of freestream velocity at this hypersonic wind tunnel facility. Measurement challenges were exasperated by the low gas density of the Mach 10 air freestream (~0.4% of standard temperature and pressure conditions) and even lower gas densities within the hypersonic wake of a 70-degree sphere-cone model. In addition, the hypersonic freestream and very low speed velocities in the wake also tested the measurement dynamic range. To complement the FLEET measurements in the wake of the sphere-cone model, PLIF velocimetry using seeded nitric oxide was also performed. While NO-PLIF velocimetry has been performed at this facility several times by previous researchers, the use of a 1D diffractive optical element for NO-PLIF velocimetry is reported here for the first time. The 1D DOE enabled the generation of up to 75 laser lines simultaneously and improved the spatial extent of the measurement three times compared to previous work. This enabled a wide velocity measurement plane of approximately 130 mm x 130 mm. The velocimetry methods demonstrated here are expected to improve wind tunnel characterization, provide critical data to validate CFD codes, and improve the design of flight vehicles for planetary entry.

FLEET↗

NASA NTRSDate not supplied

FLEET and PLIF Velocimetry Within A Mach 10 Hypersonic Air Flow

Femtosecond laser electronic excitation tagging (FLEET) and planar laser-induced fluorescence (PLIF) velocity measurements utilizing molecular tagging velocity (MTV) methods from three recent test campaigns conducted at the 31-in Mach 10 Air Tunnel at the NASA Langley Research Center are highlighted within. The FLEET measurements reported here include the first direct measurement of freestream velocity at this hypersonic wind tunnel facility. Measurement challenges were exasperated by the low gas density of the Mach 10 air freestream (~0.4% of standard temperature and pressure conditions) and even lower gas densities within the hypersonic wake of a 70-degree sphere-cone model. In addition, the hypersonic freestream and very low speed velocities in the wake also tested the measurement dynamic range. To complement the FLEET measurements in the wake of the sphere-cone model, PLIF velocimetry using seeded nitric oxide was also performed. While NO-PLIF velocimetry has been performed at this facility several times by previous researchers, the use of a 1D diffractive optical element for NO-PLIF velocimetry is reported here for the first time. The 1D DOE enabled the generation of up to 75 laser lines simultaneously and improved the spatial extent of the measurement three times compared to previous work. This enabled a wide velocity measurement plane of approximately 130 mm x 130 mm. The velocimetry methods demonstrated here are expected to improve wind tunnel characterization, provide critical data to validate CFD codes, and improve the design of flight vehicles for planetary entry.

FLEET↗

NASA NTRSDate not supplied

Modeling Self-Pressurization and Spray Bar Pressure Control of A Cryogenic Storage Tank in Normal Gravity

This paper presents computational fluid dynamics (CFD) models for simulating self-pressurization and spray-bar pressure control processes in a large-scale liquid hydrogen storage tank under normal gravity conditions. For self-pressurization, the model employs the kinetics-based Schrage equation alongside the volume-of-fluid (VOF) method to account for interfacial mass transfer. The CFD predictions of pressure and temperature are compared with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) experiment, and the predicted interfacial mass transfer rates are also presented. A CFD model simulating pressure control using a spray bar has also been developed. An Eulerian-Lagrangian approach models the interactions between discrete droplets and the continuous ullage (vapor) phase. The spray model is coupled with the VOF method by tracking droplets in the ullage and removing them when they reach the liquid interface. The T-sat model calculates droplet-ullage heat and mass transfer, where droplets warm up to the saturation temperature corresponding to the ullage vapor pressure before evaporating while remaining at the saturation temperature. The evolution of tank pressure, vapor temperature, and liquid temperature predicted by the CFD model is validated against data from the MHTB spray-bar mixing experiment. Overall, the CFD models agree with experimental data, demonstrating their capability to simulate self-pressurization and pressure control processes in large-scale cryogenic storage tanks. These models can be valuable tools for designing and optimizing cryogenic fluid management systems in future applications.

Computational Fluid Dynamics↗

NASA NTRSDate not supplied

Modeling Self-Pressurization and Spray Bar Pressure Control of A Cryogenic Storage Tank in Normal Gravity

This paper presents computational fluid dynamics (CFD) models for simulating self-pressurization and spray-bar pressure control processes in a large-scale liquid hydrogen storage tank under normal gravity conditions. For self-pressurization, the model employs the kinetics-based Schrage equation alongside the volume-of-fluid (VOF) method to account for interfacial mass transfer. The CFD predictions of pressure and temperature are compared with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) experiment, and the predicted interfacial mass transfer rates are also presented. A CFD model simulating pressure control using a spray bar has also been developed. An Eulerian-Lagrangian approach models the interactions between discrete droplets and the continuous ullage (vapor) phase. The spray model is coupled with the VOF method by tracking droplets in the ullage and removing them when they reach the liquid interface. The T-sat model calculates droplet-ullage heat and mass transfer, wherein droplets warm up to the saturation temperature corresponding to the ullage vapor pressure before evaporating while remaining at the saturation temperature. The evolution of tank pressure, vapor temperature, and liquid temperature predicted by the CFD model is validated against data from the MHTB spray-bar mixing experiment. Overall, the CFD models agree with experimental data, demonstrating their capability to simulate self-pressurization and pressure control processes in large-scale cryogenic storage tanks. These models can be valuable tools for designing and optimizing cryogenic fluid management systems in future applications.

Self-Pressurization↗

DOE OSTI2024

Estimation of mean electron energy in helium surface ionization waves on dielectric substrates

Abstract The determination of basic plasma parameters in atmospheric pressure discharges is critical to advancing their use in applications. Atmospheric pressure plasma jets have found use in the fields of medicine, agriculture, material modification and others. Atmospheric pressure plasma jets often generate plasma surface ionization waves (SIW) which interact with and propagate over surfaces. Electrical diagnostics are challenging in SIW due to high collision frequencies and small scale of the plasma discharge. This work employs a passive optical emission line ratio technique to estimate the mean electron energy in SIW over planar dielectric substrates. The method uses an intensity ratio of two helium triplet lines: He( 3 3 S ) at 706.5 nm and He( 3 3 D ) at 587.56 nm. A collisional-radiative model is used to correlate line ratio to mean electron energy and determine dependencies on electron density and He/air gas mixture. Mean electron energies ranging from 3–8 eV are determined in He/air mixtures and are found to remain constant as the surface wave propagates radially. This work provides a 2D, time-resolved, mean electron energy diagnostic for surface ionization wave propagation and validation of numerical modeling in atmospheric pressure systems with spatially varying He/air gas mixtures. The model in question is designed for use with any He line ratio in the n = 3 excitation level.

Morsell, Joshua (ORCID:0000000331073002)↗

DOE OSTI2020

Deep Reinforcement Learning Based Volt-VAR Optimization in Smart Distribution Systems

This paper develops a model-free volt-VAR optimization (VVO) algorithm via multi-agent deep reinforcement learning (DRL) in unbalanced distribution systems. This method is novel since we cast the VVO problem in distribution networks to an intelligent deep Q-network (DQN) framework, which avoids solving a specific optimization model directly when facing time-varying operating conditions in the systems. We consider statuses/ratios of switchable capacitors, voltage regulators, and smart inverters installed at distributed generators as the action variables of the agents. A delicately designed reward function guides these agents to interact with the distribution system, in the direction of reinforcing voltage regulation and power loss reduction simultaneously. The forward-backward sweep method for radial three-phase distribution systems provides accurate power flow results within a few iterations to the DRL environment. The proposed method realizes the dual goals for VVO. We test this algorithm on the unbalanced IEEE 13-bus and 123-bus systems. Numerical simulations validate the excellent performance of this method in voltage regulation and power loss reduction.

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

DOE OSTI2020

Development and Validation of a Simulation Tool to Predict the Combined Structural, Electrical, Electrochemical and Thermal Responses of Automotive Batteries

Battery safety is one of the most important design factors of electrified vehicles (EVs) since battery failure may lead to a catastrophic consequence. To improve the safety of battery systems in EVs, battery behaviors under various abuse conditions should be well understood. One way to examine battery response under extreme conditions is to conduct abuse tests in different scenarios. They can provide first-hand pass/fail information on battery safety, but they are usually expensive and time-consuming, which is inconvenient during design iterations and optimization. On the other hand, computational modeling can predict the battery behaviors in a more efficient and cost-effective way, which becomes an important tool to evaluate battery safety. Modeling battery responses and failure is challenging since it involves multiple coupled physical processes, needs to consider complex material properties and is computationally intensive. This project aims to develop a practical simulation tool to predict the combined structural, electrical, electrochemical, and thermal responses of automotive batteries to crash-induced crush and short circuit and validate it for conditions relevant to automotive crash. Advanced material constitutive models will be utilized to capture the mechanical response of cell components. Moreover, methods will be developed to reduce computational complexity of the model and allow battery safety simulations to extend from cell-level to module/pack level with affordable computational resources.

25 ENERGY STORAGE↗

NASA NTRS2010

Establishing Consensus Turbulence Statistics for Hot Subsonic Jets

Many tasks in fluids engineering require knowledge of the turbulence in jets. There is a strong, although fragmented, literature base for low order statistics, such as jet spread and other meanvelocity field characteristics. Some sources, particularly for low speed cold jets, also provide turbulence intensities that are required for validating Reynolds-averaged Navier-Stokes (RANS) Computational Fluid Dynamics (CFD) codes. There are far fewer sources for jet spectra and for space-time correlations of turbulent velocity required for aeroacoustics applications, although there have been many singular publications with various unique statistics, such as Proper Orthogonal Decomposition, designed to uncover an underlying low-order dynamical description of turbulent jet flow. As the complexity of the statistic increases, the number of flows for which the data has been categorized and assembled decreases, making it difficult to systematically validate prediction codes that require high-level statistics over a broad range of jet flow conditions. For several years, researchers at NASA have worked on developing and validating jet noise prediction codes. One such class of codes, loosely called CFD-based or statistical methods, uses RANS CFD to predict jet mean and turbulent intensities in velocity and temperature. These flow quantities serve as the input to the acoustic source models and flow-sound interaction calculations that yield predictions of far-field jet noise. To develop this capability, a catalog of turbulent jet flows has been created with statistics ranging from mean velocity to space-time correlations of Reynolds stresses. The present document aims to document this catalog and to assess the accuracies of the data, e.g. establish uncertainties for the data. This paper covers the following five tasks: Document acquisition and processing procedures used to create the particle image velocimetry (PIV) datasets. Compare PIV data with hotwire and laser Doppler velocimetry (LDV) data published in the open literature. Compare different datasets acquired at roughly the same flow conditions to establish uncertainties. Create a consensus dataset for a range of hot jet flows, including uncertainty bands. Analyze this consensus dataset for self-consistency and compare jet characteristics to those of the open literature. One final objective fulfilled by this work was the demonstration of a universal scaling for the jet flow fields, at least within the region of interest to aeroacoustics. The potential core length and the spread rate of the half-velocity radius were used to collapse of the mean and turbulent velocity fields over the first 20 jet diameters in a highly satisfying manner.

Bridges, James↗

DOE OSTI2023

A fast Monte Carlo cell-by-cell simulation for radiobiological effects in targeted radionuclide therapy using pre-calculated single-particle track standard DNA damage data

Introduction: We developed a new method that drastically speeds up radiobiological Monte Carlo radiation-track-structure (MC-RTS) calculations on a cell-by-cell basis. Methods: The technique is based on random sampling and superposition of single-particle track (SPT) standard DNA damage (SDD) files from a “pre-calculated” data library, constructed using the RTS code TOPAS-nBio, with “time stamps” manually added to incorporate dose-rate effects. This time-stamped SDD file can then be input into MEDRAS, a mechanistic kinetic model that calculates various radiation-induced biological endpoints, such as DNA double-strand breaks (DSBs), misrepairs and chromosomal aberrations, and cell death. As a benchmark validation of the approach, we calculated the predicted energy-dependent DSB yield and the ratio of direct-to-total DNA damage, both of which agreed with published in vitro experimental data. We subsequently applied the method to perform a superfast cell-by-cell simulation of an experimental in vitro system consisting of neuroendocrine tumor cells uniformly incubated with 177 Lu. Results and discussion: The results for residual DSBs, both at 24 and 48 h post-irradiation, are in line with the published literature values. Our work serves as a proof-of-concept demonstration of the feasibility of a cost-effective “in silico clonogenic cell survival assay” for the computational design and development of radiopharmaceuticals and novel radiotherapy treatments more generally.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2020

Bounded motion design in the Earth zonal problem using differential algebra based normal form methods

Establishing long-term relative bounded motion between orbits in perturbed dynamics is a key challenge in astrodynamics to enable cluster flight with minimum propellant expenditure. Herein, we present an approach that allows for the design of long-term relative bounded motion considering a zonal gravitational model. Entire sets of orbits are obtained via high-order Taylor expansions of Poincarè return maps about reference fixed points. The high-order normal form algorithm is used to determine a change in expansion variables of the map into normal form space, in which the phase space behavior is circular and can be easily parameterized by action–angle coordinates. The action–angle representation of the normal form coordinates is then used to parameterize the original Poincarè return map and average it over a full phase space revolution by a path integral along the angle parameterization. As a result, the averaged nodal period and drift in the ascending node are obtained, for which the bounded motion conditions are straightforwardly imposed. Sets of highly accurate bounded orbits are obtained, extending over several thousand kilometers, and valid for decades.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Three-dimensional cooperative inversion of airborne magnetic and gravity gradient data using deep-learning techniques

Using multiple geophysical methods has become a prevailing approach in numerous geophysical applications to investigate subsurface structures and parameters. These multimethod-based exploration strategies have the potential to greatly diminish uncertainties and ambiguities encountered during geophysical data analysis and interpretation. One of the applications is the cooperative inversion of airborne magnetic and gravity gradient data for the interpretation of data obtained in mineral, oil and gas, and geothermal explorations. In this paper, a unified cooperative inversion framework is designed by combining the standard separate inversions with a deep neural network (DNN), which serves as the link between different types of data. A well-trained DNN takes the separately inverted susceptibility and density models as the inputs and provides improved models that will be used as the initial models of deterministic inversions. A two-round iteration strategy is adopted to guarantee the reasonability of the recovered models and overall efficiency of the inversion. In addition, this deep-learning (DL)-based framework demonstrates excellent generalization abilities when tested on models that are entirely distinct from the training data sets. The framework can easily incorporate multiphysics without necessitating any structural changes to the network. Synthetic experiments validate that our DL-based method outperforms conventional separate inversions and cross-gradient-based joint inversion in view of the accuracy of the recovered models and inversion efficiency. Successful application to field data further verifies the effectiveness of our DL-based method.

Geochemistry & Geophysics↗