Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 631 records · Page 35

Transient Kinetic Experiments within the High Conversion Domain: The Case of Ammonia Decomposition

In the development of catalytic materials, a set of standard conditions is needed where the kinetic performance of many samples can be compared. This can be challenging when a sample set covers a broad range of activity. Precise kinetic characterization requires uniformity in the gas and catalyst bed composition. This limits the range of convecting devices to low conversion (generally <20%). While steady-state kinetics offer a snapshot of conversion, yield and apparent rates of the slow reaction steps, transient techniques offer much greater detail of rate processes and hence more information as to why certain catalyst compositions offer better performance. In this work, transient experiments in two transport regimes are compared: an advecting differential plug flow reactor (PFR) and a pure-diffusion temporal analysis of products (TAP) reactor. The decomposition of ammonia was used as a model reaction to test three simple materials: polycrystalline iron, cobalt and a bimetallic preparation of the two. These materials presented a wide range of activity and it was not possible to capture transient information in the advecting device for all samples at the same conditions while ensuring uniformity. We push the boundary for the theoretical estimates of uniformity in the TAP device and find reliable kinetic measurement up to 90% conversion. However, what is more advantageous from this technique is the ability to observe the time-dependence of the reaction rate rather than just singular points of conversion and yield. For example, on the iron sample we observed reversible adsorption of ammonia and on cobalt materials we identify two routes for hydrogen production. From the time-dependence of reactants and product, the dynamic accumulation was calculated. This was used to understand the atomic distribution of H and N species regulated by the surface of different materials. When ammonia was pulsed at 550 °C, the surface hydrogen/nitrogen, (H/N), ratios that evolved for Fe, CoFe and Co were 2.4, 0.25 and 0.3 respectively. This indicates that iron will store a mixture of hydrogenated species while materials with cobalt will predominantly store NH and N. While much is already known about iron, cobalt and ammonia decomposition, the goal of this work was to demonstrate new tools for comparing materials over a wider window of conversion and with much greater kinetic detail. As a result, this provides an approach for detailed kinetic discrimination of more complex industrial samples beyond conversion and yield.

36 MATERIALS SCIENCE↗

Nucleation and structural growth of cluster crystals

We study the nucleation of crystalline cluster phases in the generalized exponential model with exponent n = 4. Due to the finite value of this pair potential for zero separation, at high densities the system forms cluster crystals with multiply occupied lattice sites. Here, we investigate the microscopic mechanisms that lead to the formation of cluster crystals from a supercooled liquid in the low-temperature region of the phase diagram. Using molecular dynamics and umbrella sampling, we calculate the free energy as a function of the size of the largest crystalline nucleus in the system, and compare our results with predictions from classical nucleation theory. Employing bond-order parameters based on a Voronoi tessellation to distinguish different crystal structures, we analyze the average composition of crystalline nuclei. We find that even for conditions where a multiply occupied fcc crystal is the thermodynamically stable phase, the nucleation into bcc cluster crystals is strongly preferred. Furthermore, we study the particle mobility in the supercooled liquid and in the cluster crystal. In the cluster crystal, the motion of individual particles is captured by a simple reaction-diffusion model introduced previously to model the kinetics of hydrogen bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An open source fast fluid dynamics model for data center thermal management

Although computational fluid dynamics (CFD) has been widely adopted to improve data center thermal management, the high computational demand limits its applications, such as multivariate optimal design and operation. Fast fluid dynamics (FFD), which has been applied for fast airflow simulation, shows great potential. However, few research applied FFD for optimal design and operation of data center thermal management. This research improves the FFD model for data centers and conducts a comprehensive evaluation and demonstration. First, the FFD model is improved by solving the advection and diffusion equations together using an upwind scheme instead of a semi-Lagrangian advection solver in the conventional FFD model. Second, new features for data centers are added, such as a pressure correction method to simulate plenum airflow and dynamic boundary conditions for IT racks. The new FFD model is first validated with two indoor environment cases and the results show that the new FFD model has slightly better overall prediction accuracy and faster speed compared to the conventional FFD model. It is also observed that both FFD models achieve acceptable accuracy, except for a few localized disparities with experimental data, which might be due to simplified handling of turbulence viscosity near the boundaries. Furthermore, validation with a real data center shows that the FFD model achieves a similar level of accuracy as CFD when compared to the experimental measurements with some level of uncertainties. It is then demonstrated for data center optimal design and operation, which saves 53.4–58.8% of annual energy while still meeting the thermal requirements. In conclusion, with a much faster speed and comparable accuracy compared to CFD, the FFD model parallelized on a graphics processing unit is promising for practical model-based data center early design and operation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Soil Respiration Response to Rainfall Modulated by Plant Phenology in a Montane Meadow, East River, Colorado, USA

Soil respiration is a primary component of the terrestrial carbon cycle. However, predicting the response of soil respiration to climate change remains a challenge due to the complex interactions between environmental drivers, especially plant phenology, temperature, and soil moisture. In this study, we use a 1-D diffusion-reaction model to calculate depth-resolved CO 2 production rates from soil CO 2 concentrations and surface efflux observations in a subalpine meadow in the East River watershed, CO. Modeled rates are compared to in situ soil temperature and moisture conditions and MODIS satellite enhanced vegetation index (EVI) representing plant phenology across three hydrologically distinct growing seasons from 2016-2018. While soil respiration correlated with temperature on diel timescales (p<0.05), seasonal variability was dominated by soil moisture and plant phenology (p<0.05). We observed significant respiration increases in response to precipitation events; however, magnitude and duration were significantly higher in 2017 than 2016 despite similar wetting characteristics. Based on MODIS EVI, we suggest the respiration response to rainfall is controlled by plant phenology, which in turn reflects the capacity of plants to respond to precipitation via increased photosynthesis and autotrophic respiration, behavior that is not captured in typical soil respiration pulse models. Projected changes in montane climate such as earlier snowmelt and prolonged fore-summer drought may decrease soil respiration fluxes by decreasing the overlap between peak productivity and the summer monsoon. Lastly, we observed significant late-season CO 2 fluxes from the deep sub-soil (>165 cm) that support growing evidence for the importance of sub-soil processes in driving integrated respiration fluxes.

54 ENVIRONMENTAL SCIENCES↗

Spherical and sessile droplet dynamics by fluctuating hydrodynamics

We simulate the mesoscopic dynamics of droplets formed by phase-separated fluids at nanometer scales where thermal fluctuations are significant. Both spherical droplets fully immersed in a second fluid and sessile droplets which are also in contact with a solid surface are studied. Our model combines a Cahn–Hilliard formulation with incompressible fluctuating hydrodynamics; for sessile droplets, the fluid–solid contact angle is specified as a boundary condition. Deterministic simulations with an applied body force are used to measure the droplets' mobility from which a diffusion coefficient is obtained using the Einstein relation. Stochastic simulations are independently used to obtain a diffusion coefficient from a linear fit of the variance of a droplet's position with time. In some scenarios, these two measurements give the same value but not in the case of a spherical droplet initialized near a slip wall or in the case of sessile droplets with large contact angles (≥90°) on both slip and no-slip surfaces.

Bell, John B↗

Online quantitative analysis of multispectral images of human body tissues

A method is developed for online monitoring of structural and morphological parameters of biological tissues (haemoglobin concentration, degree of blood oxygenation, average diameter of capillaries and the parameter characterising the average size of tissue scatterers), which involves multispectral tissue imaging, image normalisation to one of its spectral layers and determination of unknown parameters based on their stable regression relation with the spectral characteristics of the normalised image. Regression is obtained by simulating numerically the diffuse reflectance spectrum of the tissue by the Monte Carlo method at a wide variation of model parameters. The correctness of the model calculations is confirmed by the good agreement with the experimental data. The error of the method is estimated under conditions of general variability of structural and morphological parameters of the tissue. The method developed is compared with the traditional methods of interpretation of multispectral images of biological tissues, based on the solution of the inverse problem for each pixel of the image in the approximation of different analytical models. (biomedical optics)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Q3 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q3: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis. • Benchmark between two first wall heat flux mapping methods, identify importance of various heat flux sources and physics impact of using fully coupled CESOL vs post-analysis evaluation. 2. Generate medium fidelity parametrized CAD. • Generate parametrized CAD components for the CAT example case via either user-defined modules called within the geometry generation or by defeatured/parametrized CAD, including DCLL blanket matched to divertor boundary and magnets. Define materials, labels, and boundary conditions for passing the mesh to CFD tools. 3. Demonstrate multiphysics magnet analysis. • Demonstrate magnet analysis workflow called from the FREDA workflow, and 4. Demonstrate nuclear analysis. • Add model to OpenFOAM and/or other codes possibly including Diablo to account for tritium diffusion in solids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Formation of superlattice structure in irradiated materials: a phase field study

The formation and self-organization of fission gas bubbles under irradiation is a unique phenomenon, which affects the swelling behavior and mechanical stability of nuclear materials. We develop a phase field model to investigate the formation mechanism of fission gas bubbles, and associated material properties under different irradiation condition. The thermodynamic and kinetic origin of superlattice formation is systematically analyzed. In particular, the symmetry of the superlattice is dictated by the anisotropic diffusion of interstitial atoms, while the lattice parameter is determined by the instability of vacancy and gas concentration waves. The spatial distribution of gas bubbles strongly depends on the nucleation, growth and coarsening process, which leads to the formation of self-organized superlattice. This study provides a new insight into designing target experiments to tailor desired microstructure under irradiation.

36 MATERIALS SCIENCE↗

Sodium Pumping via Condensation within a Non-Wetting Porous Structure

The sodium thermal electrochemical converter (Na-TEC) is a heat engine that generates electricity via the isothermal expansion of sodium ions through a beta”-alumina solid-electrolyte. The transfer of liquid sodium from a low-pressure condenser to a high-pressure evaporator is necessary to complete the thermodynamic cycle for this conversion process. A unique sodium capillary pump for the Na-TEC is proposed, whereby low-pressure sodium vapor is condensed within a non-wetting (i.e. contact angle > 90o) stainless steel porous structure. Due to the curvature of the non-wetting liquid-vapor interface, the liquid adjacent to the interface is at a higher pressure than the vapor. This is in contrast to traditional wicks, where the liquid adjacent to the interface has a lower pressure. Thus, the curvature created at this non-wetting condensation interface supplies a directed force upon the liquid sodium that effectively “pushes” it towards a high pressure region, whereas a traditional capillary wick “pulls” the liquid. Porous structures made of stainless steel 316 are commercially available, and this material is non-wetting to sodium at the relevant temperature range of the Na-TEC condenser (500K – 673K). An experimental set-up is used to measure several interfacial and transport properties of the condensed sodium within this porous structure. These properties include the temperature dependent Laplace pressure generated at the phase change interface and the permeability of the porous structure. Liquid sodium is in contact with the porous structure and is progressively pressurized with argon gas until the sodium “breaks through” and the liquid begins to flow. This break-through pressure is determined when electrical contact is established between the flowing liquid metal and a set of electrodes cemented along the flow path. The velocity of the flowing sodium can also be measured with these electrodes, and the permeability is then determined by fitting velocity data to the classical Darcy law. These interfacial properties are then applied to a conjugate heat transfer model, which treats the coupled momentum and thermal transport processes within the non-wetting porous structure as sodium vapor is condensed. Due to the low mass flowrates in the Na-TEC (< 0.1 mg/s), higher order friction terms are neglected and the classical Darcy’s law is used for non-isothermal, compressible sodium vapor flow. At lower temperatures, the sodium vapor density becomes sufficiently small and there is a transition from classic Poiseuille flow (Kn < 0.1) to Knudsen flow (Kn > 10). Furthermore, transport due to binary diffusion is considered for situations where there is a substantial presence of non-condensable argon inside the sodium vapor plenum. Initial modeling results showing expected mass flowrates, temperature profiles, and liquid-vapor interface locations are presented for various physical conditions. Finally, the main design features of an experimental set-up used to demonstrate this capillary pumping is described in detail, including the method used to measure the low sodium mass flowrates. Initial experimental results demonstrating sodium pumping for a range of temperatures (800K – 1100K) will be presented and compared against the model predictions. A pump curve used to gauge the performance of this pumping mechanism is also created.

14 SOLAR ENERGY↗

Integrated simulation of U-10Mo monolithic fuel swelling behavior

Here, a separate computational branch has been implemented within the DART (Dispersion Analysis Research Tool) computational code to simulate the swelling behavior of U-10Mo monolithic fuel under the operating conditions of high-power research and test reactors (RTRs). The monolithic branch of the DART code implements a mechanistic rate-theory-based fission-gas-behavior model for the calculation of fission gas swelling, as well as a suite of thermal, physical, and mechanical models to take into account various processes occurring in RTR fuels during irradiation. In order to accurately simulate and eventually predict U-10Mo monolithic fuel irradiation behavior, the code uses materials properties calculated with lower length-scale computational methods, such as gas atom diffusivity and U-Mo surface energy from atomic simulations and grain-morphology-specific recrystallization kinetics (recrystallized fuel volume fractions vs. fission density) predicted using the phase-field method. The remainder of fission gas behavior parameters used in the model were calibrated with measured intergranular bubble size distributions. With this integrated simulation approach, the swelling behavior of U-10Mo monolithic fuel was simulated for various initial grain sizes at different operating conditions and compared with measured data. Furthermore, because limited experimental data exist for parameter calibration detailed sensitivity studies for the important parameters used in the fission gas behavior model were performed in order to examine their impact on both intergranular gas bubble morphology at low fission density, and on total porosity at high fission density.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Influence of Shallow Cloud Populations on Transitions to Deep Convection in the Amazon

Abstract In this study, a pair of convection-permitting (2-km grid spacing), month-long, wet-season Weather Research and Forecasting (WRF) Model simulations with and without the eddy-diffusivity mass-flux (EDMF) scheme are performed for a portion of the Green Ocean Amazon (GoAmazon) 2014/15 field campaign period. EDMF produces an ensemble of subgrid-scale convective plumes that evolve in response to the boundary layer meteorological conditions and can develop into shallow clouds. The objective of this study is to determine how different treatments of shallow cumulus clouds (i.e., with and without EDMF) impact the total cloud population and precipitation across the Amazonian rain forest, with emphasis on impacts on the likelihood of shallow-to-deep convection transitions. Results indicate that the large-scale synoptic conditions in the EDMF and control simulations are nearly identical; however, on the local scale their rainfall patterns diverge drastically and the biases decrease in EDMF. The EDMF scheme significantly increases the frequency of shallow clouds, but the frequencies of deep clouds are similar between the simulations. Deep convective clouds are tracked using a cloud-tracking algorithm to examine the impact of shallow cumulus on the surrounding ambient environment where deep convective clouds initiate. Results suggest that a rapid increase of low-level cloudiness acts to cool and moisten the low to midtroposphere during the day, favoring the transition to deep convection.

54 ENVIRONMENTAL SCIENCES↗

Application of the Mixed-Mode Model for Numerical Simulation of Pearlitic Transformation

The multiscale model of pearlitic transformation, including the finite element solution of the second Fick’s equation with moving boundary in the microscale and Fourier equation in the macroscale, is presented in this paper. According to the mixed-mode approach, both the volume diffusion and the interface mobility were considered. Model describes sidewise and frontal growth of cementite and ferrite plates in a single grain of the austenite. Assessment of the possibility of the developed model application to determine the key parameters in terms of strength properties of pearlitic steel, i.e., pearlite grain size and colony size as well as interlamellar spacing for various cooling conditions, was the main aim of this work. The model was validated and verified on the basis of experimental tests performed for two eutectoid steels. In practice, developed model can support design process for a technology of high-strength rods and long products manufacturing. Rails were selected as a case study in this work, and therefore, numerical simulations of accelerated cooling of the rail head were performed and the relationship between the parameters of heat treatment, parameters of the structure and properties of the finished product was determined on the basis of obtained results.

36 MATERIALS SCIENCE↗

Radioactive Waste Water Treatment for Fukushima Daiichi Nuclear Power Plant - 16081

The Great East Japan Earthquake of 11 March 2011 caused serious damage to the Fukushima Daiichi Nuclear Power Station (NPS) of the Tokyo Electric Power Company, Inc.. Toshiba and Shaw/AVANTech installed a Simplified Active Water Retrieve and Recovery System (SARRY{sup TM}) to remove radioactive cesium (Cs) from the contaminated water. Toshiba, in partnership with EnergySolutions, also developed and installed the Multi Radionuclide Removal System (MRRSTM), which can remove radioactive nuclides being treated, except tritium, to a level less than the concentration limit in water. As contaminated water treatment has progressed, the composition of the contaminated water has changed. Both salinity and radioactive Cs concentration in the contaminated water have been decreasing. Thus, a new adsorbent that is suitable for the present water chemistry conditions is required. In order to provide an appropriate contaminated-water treatment system, Toshiba has been examining the removal of not only radioactive Cs but also radioactive strontium (Sr) by SARRY{sup TM}. SARRY{sup TM} has been operated for removing Cs by using a Cs adsorbent. Instead of a Cs adsorbent, a new media that has high selectivity for both Cs and Sr should be installed. This paper describes the results obtained when using the newly developed adsorbent for removing Cs and Sr at once. As the Cs and Sr adsorbent, we have selected crystalline silico-titanate (CST). The distribution coefficient (K{sub d}) of CST was determined by the batch method. The lifetime of the adsorbent was evaluated by a column examination. The selected CST had a K{sub d} value of over 10{sup 4} ml/g for both Cs and Sr. The highest-selectivity CST had a lifetime (defined as the period for the outlet concentration to decrease to one tenth of the inlet level) of more than 125000 BV (bed volume) for Cs and more than 18000 BV for Sr, when the test was conducted at 10 BV/h. Calculation code was also developed for estimating the relation between the BV and the decontamination factor (DF) using an adsorption rate model for diffusion and adsorption. The calculation results of the Cs and Sr adsorption behavior agreed well with the experimental results. To evaluate the replacement cycle of adsorption tower in actual equipment, the radiation exposure during the replacement operation and the required Cs and Sr removal performance were used as limiting conditions for the estimation. The results of the evaluation showed that the selected CST has sufficient performance for practical usage. Based on this finding, Toshiba proposed loading the CST in SARRY{sup TM}, and SARRY{sup TM} loaded with this CST started removing Cs and Sr in December 2014 resulting in decreased contaminated water. It has been continuously working since then. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Micro-Continuum Modeling of Nuclear Waste Glass Corrosion

The purpose of this contribution is to report on preliminary nano-continuum scale modeling of nuclear waste glass corrosion. The focus of the modeling is an experiment involving a French glass SON68 specimen leached for 25 years in a granitic environment (Guittonneau et al., 2011; Gin et al., 2011; Gin et al., 2013). The 25 year experiment carried out by French scientists offers a unique opportunity to examine the controls on the rate of nuclear glass corrosion under controlled experimental conditions using a borosilicate glass of the kind widely considered for nuclear waste storage. The long term experiment makes it possible to move beyond the initial rate stage, with direct measurements of the longer term, residual rates of corrosion. Important in this regard is whether the rate of corrosion shows a parabolic dependence on time, as would be expected from a conventional diffusion model in which the reaction products (amorphous silica-rich gel) create a diffusion barrier that grows in width over time, or whether the time dependence is linear (Guittonneau et al., 2011). A linear dependence on time would call into question the model of a Passive Reactive Interface (PRI) presented by Cailleteau et al. (2008), although it does not entirely preclude such a model if the interface maintains a constant width over time (forming at the leading edge next to the pristine glass and dissolving at the trailing edge). The leachate concentration measurements collected every 56 days over 25 years by Guittonneau et al. (2011) provide data to constrain the long term corrosion rate and its time dependence.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Crack formation in strained β-(AlxGa 1−x ) 2 O 3 films grown on (010) β-Ga 2 O 3 substrates

The cracking and local strain relaxation in (010) (Al x Ga 1−x ) 2 O 3 films grown on Ga 2 O 3 substrates are assessed in terms of film composition and thickness. We utilize x-ray diffraction and electron microscopy techniques combined with simulation and modeling to investigate that film cracking on curvature (flatness) has a directly proportional relationship with film thickness and/or aluminum content. Cross section transmission electron microscopy reveals that cracks along both the (001) and (100) cleavage planes penetrate into the substrate. The diffuse scattered intensity observed in reciprocal space maps (RSMs) is directly correlated with the tilt that is introduced due to the change in the deformation conditions near the cracks. While asymmetric RSMs show that the layers are fully strained, the diffuse scattering distribution in reciprocal space can be interpreted to show that the cracking relaxes and locally tilts the lattice ∼350 nm from the crack edges, which is consistent with the larger radius of curvature associated with the films with higher crack densities. For example, a 200 nm (Al 0.13 Ga 0.87 ) 2 O 3 thick film has an average inter-crack spacing of 3.3 µm, so most of the epitaxial layer is fully strained except near the cracks where it deforms elastically and is consistent with the gallium oxide Poisson ratio. A reciprocal space model was developed, which imports the strain and tilt distributions (based on finite element modeling) to match the features observed in the experimental maps. We also note that previous studies involving (Al x Ga 1−x ) 2 O 3 films may show evidence of cracking as observed in their symmetric and asymmetric RSMs.

36 MATERIALS SCIENCE↗

Air separation via two-step solar thermochemical cycles based on SrFeO 3-δ and (Ba,La) 0.15 Sr 0.85 FeO 3-δ perovskite reduction/oxidation reactions to produce N 2 : rate limiting mechanism(s) determination

Two-step solar thermochemical cycles based on reversible reactions of SrFeO 3–δ and (Ba,La) 0.15 Sr 0.85 FeO 3–δ perovskites were considered for air separation. The cycle steps encompass (1) the thermal reduction of SrFeO 3–δ or (Ba,La) 0.15 Sr 0.85 FeO 3–δ perovskites driven by concentrated solar irradiation and (2) oxidation in air to remove O 2 and produce N 2 . Rate limiting mechanisms were examined for both reactions using a combination of isothermal and non-isothermal thermogravimetry for temperature-swings between 673 and 1373 K, heating rates of 10, 20, and 50 K min –1 , and O 2 pressure-swings between 20% O 2 /Ar and 100% Ar at atmospheric pressure. Evolved O 2 and associated lag due to transport behavior were measured with gas chromatography and used with measured sample temperatures to predict equilibrium compositions from a compound energy formalism thermodynamic model. Measured and predicted chemical equilibrium changes in deviation from stoichiometry were compared. Rapid chemical kinetics were observed as the samples equilibrated rapidly for all conditions, indicative that heat and mass transfer were the rate limiting mechanisms. Furthermore, the effects of bulk diffusion (or gas diffusion through the bed or pellet) were examined using pelletized and loose powdered samples and determined to have no discernable impact.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Data-Driven RANS Turbulence Closures for Forced Convection Flow in Reactor Downcomer Geometry

Recent progress in data-driven turbulence modeling has shown its potential to enhance or replace traditional equation-based Reynolds-averaged Navier-Stokes (RANS) turbulence models. Here, this work utilizes invariant neural network (NN) architectures to model Reynolds stresses and turbulent heat fluxes in forced convection flows (when the models can be decoupled). As the considered flow is statistically one dimensional, the invariant NN architecture for the Reynolds stress model reduces to the linear eddy viscosity model. To develop the data-driven models, direct numerical and RANS simulations in vertical planar channel geometry mimicking a part of the reactor downcomer are performed. Different conditions and fluids relevant to advanced reactors (sodium, lead, unitary-Prandtl-number fluid, and molten salt) constitute the training database. The models enabled accurate predictions of velocity and temperature, and compared to the baseline k–τ turbulence model with the simple gradient diffusion hypothesis, do not require tuning of the turbulent Prandtl number. The data-driven framework is implemented in the open-source graphics processing unit–accelerated spectral element solver nekRS and has shown the potential for future developments and consideration of more complex mixed convection flows.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Electromagnetic Extended Finite Elements for High-Fidelity Multimaterial Problems LDRD Final Report

Surface effects are critical to the accurate simulation of electromagnetics (EM) as current tends to concentrate near material surfaces. Sandia EM applications, which include exploding bridge wires for detonator design, electromagnetic launch of flyer plates for material testing and gun design, lightning blast-through for weapon safety, electromagnetic armor, and magnetic flux compression generators, all require accurate resolution of surface effects. These applications operate in a large deformation regime, where body-fitted meshes are impractical and multimaterial elements are the only feasible option. State-of-the-art methods use various mixture models to approximate the multi-physics of these elements. The empirical nature of these models can significantly compromise the accuracy of the simulation in this very important surface region. We propose to substantially improve the predictive capability of electromagnetic simulations by removing the need for empirical mixture models at material surfaces. We do this by developing an eXtended Finite Element Method (XFEM) and an associated Conformal Decomposition Finite Element Method (CDFEM) which satisfy the physically required compatibility conditions at material interfaces. We demonstrate the effectiveness of these methods for diffusion and diffusion-like problems on node, edge and face elements in 2D and 3D. We also present preliminary work on h -hierarchical elements and remap algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗