SEARCH · Search NASA

Results for “embedding model”

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 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Transonic analysis of canted winglets

A computational method developed to provide a transonic analysis for upper/lower surface wing-tip mounted winglets is described. Winglets with arbitrary planform, cant and toe angle, and airfoil section can be modeled. The embedded grid approach provides high flow field resolution and the required geometric flexibility. In particular, coupled Cartesian/cylindrical grid systems are used to model the complex geometry presented by canted upper/lower surface winglets. A new rotated difference scheme is introduced in order to maintain the stability of the small-disturbance formulation in the presence of large spanwise velocities. Wing and winglet viscous effects are modeled using a two-dimensional 'strip' boundary layer analysis. Correlations with wind tunnel and flight test data for three transport configurations are included.

Rosen, B. S.↗

NASA NTRS2013

A Validated Set of MIDAS V5 Task Network Model Scenarios to Evaluate Nextgen Closely Spaced Parallel Operations Concepts

The Closely Spaced Parallel Operations (CSPO) scenario is a complex, human performance model scenario that tested alternate operator roles and responsibilities to a series of off-nominal operations on approach and landing (see Gore, Hooey, Mahlstedt, Foyle, 2013). The model links together the procedures, equipment, crewstation, and external environment to produce predictions of operator performance in response to Next Generation system designs, like those expected in the National Airspaces NextGen concepts. The task analysis that is contained in the present report comes from the task analysis window in the MIDAS software. These tasks link definitions and states for equipment components, environmental features as well as operational contexts. The current task analysis culminated in 3300 tasks that included over 1000 Subject Matter Expert (SME)-vetted, re-usable procedural sets for three critical phases of flight; the Descent, Approach, and Land procedural sets (see Gore et al., 2011 for a description of the development of the tasks included in the model; Gore, Hooey, Mahlstedt, Foyle, 2013 for a description of the model, and its results; Hooey, Gore, Mahlstedt, Foyle, 2013 for a description of the guidelines that were generated from the models results; Gore, Hooey, Foyle, 2012 for a description of the models implementation and its settings). The rollout, after landing checks, taxi to gate and arrive at gate illustrated in Figure 1 were not used in the approach and divert scenarios exercised. The other networks in Figure 1 set up appropriate context settings for the flight deck.The current report presents the models task decomposition from the tophighest level and decomposes it to finer-grained levels. The first task that is completed by the model is to set all of the initial settings for the scenario runs included in the model (network 75 in Figure 1). This initialization process also resets the CAD graphic files contained with MIDAS, as well as the embedded operator models that comprise MIDAS. Following the initial settings, the model progresses to begin the first tasks required of the two flight deck operators, the Captain (CA) and the First Officer (FO). The task sets will initialize operator specific settings prior to loading all of the alerts, probes, and other events that occur in the scenario. As a note, the CA and FO were terms used in developing this model but the CA can also be thought of as the Pilot Flying (PF), while the FO can be considered the Pilot-Not-Flying (PNF)or Pilot Monitoring (PM). As such, the document refers to the operators as PFCA and PNFFO respectively.

task network model↗

DOE OSTI2018

Evaluating the Impact of Ocean Eddies in the Global Marine Carbon Cycle

The ocean plays a crucial role in the Earth’s carbon cycle due to its ability to exchange carbon dioxide (CO2) with the atmosphere, and subsequently store it for long periods of time. Part of this sequestration is due to a mechanism known as the Biological Pump, a process where phytoplankton and autotrophic bacteria convert dissolved CO2 to organic matter via photosynthesis, much of which ends up sinking and can remain in the deep ocean for hundreds or thousands of years. If the strength of the Biological Pump changes due the human-induced increase in atmospheric CO2, it will affect the rate of global warming. Earth System Models (ESMs) are now including ocean biogeochemical processes in order to quantify changes in the global carbon cycle under climate change scenarios by embedding a model that represents ocean chemistry and the dynamics of marine ecosystems within an ocean circulation framework. Due to the increased computational cost of including marine biogeochemistry (BGC), they are typically run with relatively low complexity and coarse horizontal grid resolution. Our project utilizes Institutional Computing power to undertake an ESM simulation using a moderate complexity BGC model (35 constituents) and an increased grid resolution that allows the formation of turbulent mesoscale eddies (30 to 200 km in size) in the ocean. Comparing with low resolution runs will allow us to assess whether standard ESMs are biased in their representation of the carbon cycle due to parameterized effects of sub-grid scale physics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Evaluating the Impact of Ocean Eddies in the Global Marine Carbon Cycle

The ocean plays a crucial role in the Earth’s carbon cycle due to its ability to exchange carbon dioxide (CO 2 ) with the atmosphere, and subsequently store it for long periods of time. Part of this sequestration is due to a mechanism known as the Biological Pump, a process where phytoplankton and autotrophic bacteria convert dissolved CO 2 to organic matter via photosynthesis, much of which ends up sinking and can remain in the deep ocean for hundreds or thousands of years. If the strength of the Biological Pump changes due the human-induced increase in atmospheric CO 2 , it will affect the rate of global warming. Earth System Models (ESMs) are now including ocean biogeochemical processes in order to quantify changes in the global carbon cycle under climate change scenarios by embedding a model that represents ocean chemistry and the dynamics of marine ecosystems within an ocean circulation framework. Due to the increased computational cost of including marine biogeochemistry (BGC), they are typically run with relatively low complexity and coarse horizontal grid resolution. Our project utilizes Institutional Computing power to undertake an ESM simulation using a moderate complexity BGC model (35 constituents) and an increased grid resolution that allows the formation of turbulent mesoscale eddies (30 to 200 km in size) in the ocean. Comparing with low resolution runs will allow us to assess whether standard ESMs are biased in their representation of the carbon cycle due to parameterized effects of sub-grid scale physics.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1991

Computing Availability And Reliability For A System

Elements contributing most to failure of system identified quickly. Reliability/Availability Analysis (RELAV) computer program is comprehensive analytical software tool to determine reliability or availability of any general system modeled as embedded k-out-of-n groups of items (components) and/or subgroups. Used to assess performance of system during late testing phase of design of system, to model candidate designs and/or architectures, or to validate and form predictions during early phases of design. C-language program.

Bowerman, Paul N.↗

NASA NTRS1999

Control of a Thrust Alignment Table for Modeling the Coning Dynamics of a Spinning Spacecraft with a Follower Force

This document presents a system controlling the motion of a spherical air bearing used in the modeling of spacecraft dynamics and controls in a laboratory environment. The system is part of the Spinning Rocket Simulator (SRS), used to simulate the coning of spacecraft during a thrusting stage. The reaction force at the spherical air bearing supporting the spacecraft model must coincide with the thrust axis of the model for proper simulation. Therefore, the bearing is translated in a circular path to introduce a centrifugal force. This horizontal force along with the gravitational reaction force at the bearing combines to simulate the direction of the spacecraft's thrust force. The control system receives attitude information from the spacecraft model via a laser beam embedded in the model that impinges on a photosensitive array. The non-linear system is controlled using high-speed lookup tables and digital techniques. A vector-controlled motor and a stepper motor are given the necessary signals to accurately control the turntable and platform supporting the air bearing. Preliminary performance data is presented. Mechanical elements of the table and platform are described in detail. A wireless (RF) data path for all devices on the spacecraft model to an off-table command computer is also described.

Halsmer, Dominic↗

DOE OSTI2018

A Computer Model to Study the Response of Energetic Materials to a Range of Dynamic Loads

Over the past decade we developed a model to enable computer simulation of the mechanical and subsequent energetic response of explosives and propellants to mechanical insults such as impacts, perforations, drops, and falls. The model is embedded in computer simulation programs that solve the non-linear, large deformation equations of compressible solid and fluid flow in space and time. It is implemented as a user-defined model, which returns the updated stress tensor and composition that result from the simulation supplied strain tensor change. Although it is multi-phase, in that gas and solid species are present, it is single-velocity, in that the gas does not flow through the porous solid. More than 70 time-dependent variables are made available for additional analyses and plotting. The model encompasses a broad range of possible responses: mechanical damage with no energetic response, and a continuous spectrum of degrees of violence including delayed and prompt detonation. This work describes the basic workings of the model.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Theory, Solution Methods, and Implementation of the HERMES Model

The HERMES (high explosive response to mechanical stimulus) model was developed over the past decade to enable computer simulation of the mechanical and subsequent energetic response of explosives and propellants to mechanical insults such as impacts, perforations, drops, and falls. The model is embedded in computer simulation programs that solve the non-linear, large deformation equations of compressible solid and fluid flow in space and time. It is implemented as a user-defined model, which returns the updated stress tensor and composition that result from the simulation supplied strain tensor change. Although it is multi-phase, in that gas and solid species are present, it is single-velocity, in that the gas does not flow through the porous solid. More than 70 time-dependent variables are made available for additional analyses and plotting. The model encompasses a broad range of possible responses: mechanical damage with no energetic response, and a continuous spectrum of degrees of violence including delayed and prompt detonation. This paper describes the basic workings of the model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Extending theories on muon-specific interactions

The proton radius puzzle, the discrepancy between the proton radius measured in muonic hydrogen and electronic hydrogen, has yet to be resolved. There are suggestions that beyond the standard model (BSM) physics could resolve both this puzzle and the muon anomalous magnetic moment discrepancy. Karshenboim et al. point out that simple, nonrenormalizable, models in this direction involving new vector bosons have serious problems when confronting high energy data. The prime example is radiative corrections to W to μν decay which exceed experimental bounds. We show how embedding the model in a larger and arguably renormalizable theory restores gauge invariance of the vector particle interactions and controls the high energy behavior of decay and scattering amplitudes. Thus BSM explanations of the proton radius puzzle can still be viable.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1995

Coronal roots of solar wind streams: 3-D MHD modeling

Weak (discontinuous) solutions of the 3-D MHD equations look like a promising tool to model the transonic solar wind with structural elements: current sheets, coronal plumes etc. Using the observational information about various coronal emissions one can include these structural elements into the 3-D MHD solar wind model by embedding the discontinuities of given type. Such 3-D MHD structured solar wind is calculated self-consistently: variants are examined via numerical experiments. In particular, the behavior of coronal plumes in the transonic solar wind flow, is modeled. The input information for numerical modeling (for example, the magnetic field map at the very base of the solar corona) can be adjusted so that fast stream arises over the center of the coronal hole, over the coronal hole boundaries and, even, over the region with closed magnetic topology. 3-D MHD equations have the analytical solution which can serve as a model of supersonic trans-alfvenic solar wind in the (5-20) solar radii heliocentric distance interval. The transverse, nonradial total (gas + magnetic field) pressure balance in the flow is the corner-stone of this solution. The solution describes the filamentation (ray-like structure of the solar corona) and streaming (formation of high-speed streams with velocities up to 800 km/sec) as a consequence of the magnetic field spatial inhomogeneous structure and trans-alfvenic character of the flow. The magnetic field works in the model as a 'controller' for the solar wind streaming and filamentation.

Pisanko, Yu. V.↗

DOE OSTI2025

Liquid Phase Modeling in Porous Media: Adsorption of Methanol and Ethanol in H-MFI in Condensed Water

Zeolites are used in the chemical and separation industries for their exceptional selectivity, adsorption capacity, regenerability, and stability in gas and liquid phase processing. Here, we developed an explicit solvation method for predicting solvent/condensed phase effects on adsorption free energies in microporous media such as zeolites based on the hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) technique. Our explicit solvation method for zeolite systems, called eSZS, aims to capture site-specific interactions during the adsorption process at the Brønsted acid sites of H-MFI zeolite while still considering the diverse configuration space of the solvent molecules. This strategy is ideal for chemical reactions or adsorbates that interact with the microporous medium in few distinct adsorbate/transition state configurations, i.e., the harmonic or similar approximations are acceptable for the adsorbate/transition state while such approximations break down for the solvent molecules that require extensive configuration space sampling. In this way, our approach effectively overcomes the limitations of implicit solvation models and classical force field methods for describing solvation effects on chemical reactions within porous materials such as zeolites. Specifically, in this study, we investigated various aspects of our hybrid QM/MM approach, including QM cluster size dependencies in a periodic electrostatically embedded cluster model (PEECM), rules for link atoms at the QM/MM boundary, and functional and basis set considerations for converged and reasonably accurate gas and aqueous phase methanol and ethanol adsorption free energy predictions in H-MFI. For gas phase adsorption of methanol and ethanol in H-MFI at a Brønsted acid site in T12 position, we compute adsorption free energies at 298 K of −0.61 and −0.75 eV, respectively, using a PEECM containing 50 Si and 1 Al atom with ωB97x-D/def2-TZVP level of theory. For solvent effect calculations, we sample the aqueous phase using grand canonical Monte Carlo (GCMC) simulations to (1) obtain a mean field of electrostatic interactions in the reaction system and (2) perform a rigorous free energy perturbation calculation. Similar to the experimentally and computationally observed endergonic solvation effects observed for hydrocarbon adsorption on metal surfaces, we also observe that a condensed aqueous environment destabilizes methanol and ethanol at these acid sites in H-MFI at 298 K. Specifically, the computed solvation free energies of adsorption (ΔΔG solv ) for methanol and ethanol are +0.44 and +0.54 eV, respectively. From this study, it is evident that adsorbates (methanol and ethanol) are competing with water for adsorption space inside the H-MFI zeolite, leading to an endergonic solvation effect. Here, we expect that the endergonic, aqueous solvent effect during adsorption in microporous zeolites is highly tunable by changing the pore size and hydrophobicity of the microporous material as this will affect the water density inside the pore structure.

Adsorption↗

DOE OSTI2015

Influence of string-like cooperative atomic motion on surface diffusion in the (110) interfacial region of crystalline Ni

Although we often think about crystalline materials in terms of highly organized arrays of atoms, molecules, or even colloidal particles, many of the important properties of this diverse class of materials relating to their catalytic behavior, thermodynamic stability, and mechanical properties derive from the dynamics and thermodynamics of their interfacial regions, which we find they have a dynamics more like glass-forming (GF) liquids than crystals at elevated temperatures. This is a general problem arising in any attempt to model the properties of naturally occurring crystalline materials since many aspects of the dynamics of glass-forming liquids remain mysterious. We examine the nature of this phenomenon in the “simple” case of the (110) interface of crystalline Ni, based on a standard embedded-atom model potential, and we then quantify the collective dynamics in this interfacial region using newly developed methods for characterizing the cooperative dynamics of glass-forming liquids. As in our former studies of the interfacial dynamics of grain-boundaries and the interfacial dynamics of crystalline Ni nanoparticles (NPs), we find that the interface of bulk crystalline Ni exhibits all the characteristics of glass-forming materials, even at temperatures well below the equilibrium crystal melting temperature, T{sub m}. This perspective offers a new approach to modeling and engineering the properties of crystalline materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1987

Flux tube models of solar plages

The spectroscopic properties of model solar plages composed of narrow magnetic flux tubes embedded in undisturbed quiet sun are analyzed. The final models are magnetostatic flux tubes for which the level tau5000 = 1 level is 200 km below the same level in the quiet sun; these models are hotter than the quiet sun at all levels and expand by roughly a factor of three between this level and the quiet sun temperature minimum. These models are tested against a wide range of observational constraints. The assumed flux tube diameter is found to be an important parameter for predictions of the center-to-limb variation of the magnetically insensitive lines and of the visible continuum. A diameter of 300 km at tau5000 = 1 is necessary to give qualitative agreement with observations. A model plage consisting of the derived flux tube model embedded in a quiet sun model derived from the Ca II K line wings reproduces this line in plages at the observed value of mu; one of the models also predicts values of disk center plage contrast at 30, 50, 100, and 200 micron wavelength in good agreement with the observations of Jefferies et al. (1982).

Walton, Stephen R.↗

NASA NTRS1999

Effects of Precipitation on Ocean Mixed-Layer Temperature and Salinity as Simulated in a 2-D Coupled Ocean-Cloud Resolving Atmosphere Model

A two-dimensional coupled ocean-cloud resolving atmosphere model is used to investigate possible roles of convective scale ocean disturbances induced by atmospheric precipitation on ocean mixed-layer heat and salt budgets. The model couples a cloud resolving model with an embedded mixed layer-ocean circulation model. Five experiment are performed under imposed large-scale atmospheric forcing in terms of vertical velocity derived from the TOGA COARE observations during a selected seven-day period. The dominant variability of mixed-layer temperature and salinity are simulated by the coupled model with imposed large-scale forcing. The mixed-layer temperatures in the coupled experiments with 1-D and 2-D ocean models show similar variations when salinity effects are not included. When salinity effects are included, however, differences in the domain-mean mixed-layer salinity and temperature between coupled experiments with 1-D and 2-D ocean models could be as large as 0.3 PSU and 0.4 C respectively. Without fresh water effects, the nocturnal heat loss over ocean surface causes deep mixed layers and weak cooling rates so that the nocturnal mixed-layer temperatures tend to be horizontally-uniform. The fresh water flux, however, causes shallow mixed layers over convective areas while the nocturnal heat loss causes deep mixed layer over convection-free areas so that the mixed-layer temperatures have large horizontal fluctuations. Furthermore, fresh water flux exhibits larger spatial fluctuations than surface heat flux because heavy rainfall occurs over convective areas embedded in broad non-convective or clear areas, whereas diurnal signals over whole model areas yield high spatial correlation of surface heat flux. As a result, mixed-layer salinities contribute more to the density differences than do mixed-layer temperatures.

Li, Xiaofan↗

DOE OSTI2019

Comparative Assessment of Experimental Testing of Instrument with an Embedded Digital Device Using Model-Based and Conventional Methods

The use of digital equipment has not been widespread in U.S. nuclear facilities. A major contributing factor to this situation arises from regulatory concerns about common cause failures (CCFs) stemming from the use of identical software in redundant systems, questions as to how to quantify the reliability of software-based products, validity of tools that are used for software verification and validation (V&V), and the effectiveness of software QA procedures. Compounding the concerns are sporadic incidents of digital equipment malfunctions in nuclear power plants, including lockups resulting from software memory leaks, feedwater control malfunctions due to excessive traffic on Ethernet networks, software logic errors, and other problems leading to plant trips. Additionally, the limited amount of demonstrated research and a lack of industry consensus on quantitative software reliability methodologies has hindered the widespread adoption of digital I&C in the nuclear industry. In a recent Regulatory Issue Summary (RIS), the U.S. Nuclear Regulatory Commission (NRC) identified concerns about the impact of greater use of instrumentation and control (I&C) equipment with an embedded digital device (EDD) [1]. Specifically, the NRC staff states that increased use of such devices “may increase a facility’s vulnerability to a CCF” (i.e., common-cause failure). The prevailing NRC guidance regarding software CCF (SCCF) identifies two design attributes that are acceptable for eliminating CCF concerns: (1) diversity or (2) testability (specifically, 100% testability). Either solution can result in high costs and remaining licensing uncertainty (i.e., how much diversity is enough? how to ensure test coverage of every possible sequence of device states?). Consequently, many utilities and reactor designers have limited or avoided more extensive use of digital technology to minimize licensing, scheduling, and financial risk. Without development of cost effective qualification methods to satisfy regulatory requirements and address the potential for CCF vulnerability associated with EDDs, the nuclear power industry may not be able to realize the benefits of digital technology achieved by other industries. In order to address these challenges, the U.S. Department of Energy (DOE) is sponsoring this research to assess the prospect of CCF in digital I&C equipment. The purpose of the current research is to develop an effective approach employing science-based methods to resolve concerns about CCF vulnerability that serve to inhibit deployment of advanced instrumentation (e.g., sensors, actuators, microcontrollers) with EDDs in nuclear power applications. The research objectives address the challenge of establishing high levels of safety and reliability assurance needed for the qualification of EDDs (e.g., microprocessors, programmable logic devices) that are subject to software design faults, complex failure modes, and CCF vulnerability. Specific objectives are: (1) assess the regulatory context for treatment of CCF vulnerability a graded approach to their qualification, (3) develop and extend model-based testing methods to enable effective demonstration of whether devices are subject to CCF, which may arise from vulnerabilities introduced at any stage of the design lifecycle, (4) establish a cost-effective testing framework that incorporates automation and test scenario prioritization, and (5) demonstrate the qualification approach through selection and testing of candidate digital device(s). This report documents an assessment of the results from experimental testing of a representative instrument with an embedded digital device. The experiment serves as a proof-of-concept demonstration of a model-based testing (MBT) framework. The objective of MBT is to ensure that an EDD undergoes suitably comprehensive testing to establish that it is not subject to SCCF. The approach by which MBT accomplishes this objective is through development and application of an effective test suite that provides the appropriate coverage to detect the full range of postulated faults that may arise at the requirements, design, and code (implementation) level of the EDD software. The findings of the experiment involve the application of the MBT framework by Ohio State University (OSU) through testing of the selected demonstration instrument (i.e., the prototype smart pressure sensor developed by Virginia Commonwealth University, VCU). In addition, the experiment findings also involve application of a more conventional black-box testing approach that serves as a baseline against which the MBT methodology can be compared. The primary outcome documented in this report is the comparative analysis of the capabilities demonstrated for each method. The results of the analysis demonstrate clear benefits that arise from the MBT methodology. The comparison indicates greater effectiveness from the more comprehensive evidence available through MBT and enhanced efficiency achieved through automation within the MBT framework. The structure of the report addresses the approach and application of each test method as well as the comparative analysis of the experimental results. Chapter 2 gives a summary overview of the demonstration approach. Chapter 3 covers the baseline testing and Chapter 4 documents the model-based testing. Chapter 5 describes the findings of the comparative analysis and presents the conclusions that can be drawn about the effectiveness and efficiency of the MBT methodology.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Development and Validation of SAM Multi-dimensional Flow Model for Thermal Mixing and Stratification Modeling

Thermal mixing and stratification in large pools or enclosures are very important phenomena that are critical to nuclear reactor safety. Because of the wide ranges of time and length scales associated with such phenomena, accurate modeling and simulation of thermal mixing and stratification remain as the key unresolved, challenging problems for reactor transient analyses. In traditional system analysis codes, simplified zero-dimensional (0-D) models are widely used for their high numerical efficiency, but they generally suffer from very limited prediction accuracies or range of applicability. Like traditional system analysis codes, the current version of the SAM code has implemented such simplified 0-D and one-dimensional mixing models. On the other hand, high-resolution Computational Fluid Dynamics (CFD) tools are often used to model complex thermal mixing and stratification phenomena. They are, however, generally numerically expensive, and they require large amounts of computational resources. It is therefore desirable to implement advanced and efficient thermal mixing and stratification modeling capabilities embedded in a system analysis code. This approach will improve the accuracy of reactor safety analyses when thermal mixing and stratification are involved, and also avoid using the large computational resources needed for high-resolution CFD analysis. Currently, with the support of the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation program, an effort has been launched to develop and implement a multi-dimensional flow model in the system analysis code SAM, and demonstrate its applications to model thermal mixing and stratification phenomena in large enclosures. The main outcomes of this research and development activity are summarized in this report, which presents an attempt to include a built-in advanced multi-dimensional flow model in a system analysis code with the focus on overcoming the simulation challenges of thermal mixing and stratification phenomena. In this report, we start with the introduction of existing SAM code capabilities to simulate thermal mixing and stratification phenomena, which is followed by a short summary of the multi-dimensional model implemented in the SAM code, including both the physical model and the Finite Element Method code implementation. In this study, two options were implemented in the SAM code to model turbulent flows: a relatively simple built-in turbulence model and an interface to accept externally computed turbulent viscosities (e.g., from a high- resolution CFD simulation). Code validation studies on this newly added capability were then carried out to compare SAM simulation results with experimental data from the SUPERCAVNA facility, which was designed to study the complex flow recirculation and thermal stratification phenomena relevant to sodium fast reactor designs. In this study, one transient and two steady-state test cases were used for code validation. Different approaches have been used to model the complex turbulence flow fields in the SUPERCAVNA facility. A highly simplified zero-equation turbulence model was first used, but it was determined that it is too simple to capture the complex turbulence flow fields in these test cases. Subsequently, the code validation continued with the use of turbulent viscosity data from high-resolution STAR-CCM+ CFD simulations to improve the accuracy of the results. Using this approach, the SAM simulation results showed very good agreement with both the SUPERCAVNA experimental data and STAR-CCM+ simulation results. In this report, we demonstrate the development, implementation, and successful validation of a multi-dimensional flow model in the SAM code, which aims to improve the simulation accuracy for complex thermal mixing and stratification phenomena. Lessons have also been learned, including that in cases where the flow fields are not well predicted by the zero-equation model, the iinput of turbulent viscosities from an external source can enhance the overall predictive capabilities needed to accurately capture complex thermal-hydraulic phenomena. Therefore, future research will be needed to further improve the code’s capabilities, such as by developing a more efficient and robust approach to capture the turbulence effects in the SAM code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Development and Validation of Resistance-Capacitance Model for Phase Change Material Embedded in Porous Media

A simple and effective way to enhance PCM thermal conductivity is by embedding porous materials with high thermal conductivities. In this work, CFD was used to evaluate the melting of PCM embedded in metal foam supplied with constant heat flux. Good agreement between experimental and numerical results has been found with a mean PCM local temperature deviation of 1.62K and 2.97K in two locations. The CFD simulation time was 56 hours for 6.5 hours of flow time. A computationally efficient resistance capacitance-based model (RCM) was developed to improve simulation speed. The proposed RCM exhibited mean local temperature deviations of 1.86K and 2.96K when compared against the experimental data for two thermocouple positions, however, it was approximately 2.5x10 5 times faster than CFD. The charging time deviation was less than 1% compared to CFD, and no significant change in thermal energy stored was observed since the temperature differences were small enough to result in insignificant sensible heat deviation. A simulation cost-effectiveness index (CEI) considering accuracy and computational affordability was introduced; the CEI of the RCM solver is 10 3 times better compared to CFD for the prediction of local PCM temperature.

25 ENERGY STORAGE↗

DOE OSTI2019

Efficient discontinuous finite difference meshes for 3-D Laplace–Fourier domain seismic wavefield modelling in acoustic media with embedded boundaries

Simulation of acoustic wave propagation in the Laplace–Fourier (LF) domain, with a spatially uniform mesh, can be computationally demanding especially in areas with large velocity contrasts. To improve efficiency and convergence, we use 3-D second- and fourth-order velocity-pressure finite difference (FD) discontinuous meshes (DM). Our DM algorithm can use any spatial discretization ratio between meshes. We evaluate direct and iterative parallel solvers for computational speed, memory requirements and convergence. Finally, benchmarks in realistic 3-D models and topographies show more efficient and stable results for DM with direct solvers than uniform mesh results with iterative solvers.

58 GEOSCIENCES↗