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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A hybrid molecular dynamics/fluctuating hydrodynamics method for modelling liquids at multiple scales in space and time

A new 3D implementation of a hybrid model based on the analogy with two-phase hydrodynamics has been developed for the simulation of liquids at microscale. The idea of the method is to smoothly combine the atomistic description in the molecular dynamics zone with the Landau-Lifshitz fluctuating hydrodynamics representation in the rest of the system in the framework of macroscopic conservation laws through the use of a single “zoom-in” user-defined function s that has the meaning of a partial concentration in the two-phase analogy model. In comparison with our previous works, the implementation has been extended to full 3D simulations for a range of atomistic models in GROMACS from argon to water in equilibrium conditions with a constant or a spatially variable function s. Preliminary results of simulating the diffusion of a small peptide in water are also reported.

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

Thermal transport modeling of laser-irradiated spheres

Thermal transport of uniformly laser-irradiated spheres of various materials is investigated computationally. One-dimensional simulations of low- to mid-Z materials (Be, Al, and Cu) are performed to evaluate the impact of nonlocal electron transport on experimental observables under laser intensities of relevance to direct-drive inertial confinement fusion. Here we compare thermal transport models of different levels of fidelity: flux-limited Spitzer–Harm diffusion, the Schurtz–Nicolai–Busquet (SNB) reduced-order nonlocal model, and a Fokker–Planck description. Spitzer–Harm diffusion with different flux-limiter factors are compared with different implementations of the SNB model in the HYDRA radiation hydrodynamics code. Under the conditions of interest, the peak heat flux in the thermal front with the SNB model shows good agreement with Fokker–Planck calculations, with the largest errors below 10% at 10 15 W/cm 2 laser intensity. From HYDRA-SNB simulations, two experimentally relevant effects are observed from nonlocal heat transport when compared to flux-limited Spitzer–Harm modeling: coronal temperatures are cooler due to reduced heat fluxes in the expanding plasma and (for mid-Z materials) x-ray emissions are enhanced due to preheating in the dense plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Three-dimensional dendritic needle network model for alloy solidification

Here, we present a three-dimensional (3D) formulation of the multiscale Dendritic Needle Network (DNN) model for dendritic microstructure growth. This approach is aimed at simulating quantitatively the solidification dynamics of complex hierarchical networks in spatially extended dendritic arrays, hence bridging the scale gap between phase-field simulations at the scale of a few dendrites and coarse-grained simulations on the larger scale of entire polycrystalline structures. In the DNN model, the dendritic network is represented by a network of branches that interact through the solutal diffusion field. The tip velocity V(t) and tip radius ρ(t) of each needle is determined by combining a standard solvability condition that fixes the product ρ 2 V and a solute flux conservation condition that fixes the product ρV 2 in 2D and ρV in 3D as a function of a solutal flux intensity factor ℱ(t). The latter measures the intensity of the solute flux in the dendrite tip region and can be calculated by contour (2D) or surface (3D) integration around the tip of each needle. We first present an extended formulation of the 2D DNN model where needles have a finite thickness and parabolic tips. This formulation remains valid for a larger range of tip Péclet number than the original thin needle formulation and is readily extended to 3D needles with paraboloidal tips. The 3D DNN model based on this thick-needle formulation is developed for both isothermal and directional solidification. Model predictions are validated by comparisons with known analytical solutions that describe the early transient and steady-state growth regimes. We exploit the power of the DNN model to characterize the competitive growth of well-developed secondary branches in 3D on the scale of the diffusion length. The results show that the length of active secondary branches increases as a power law of distance behind the tip with an exponent in good quantitative agreement with experimental measurements. Finally, we apply the model to simulate the three-dimensional directional solidification of an Al-7wt% Si alloy, which we directly compare to observed microstructures from microgravity experiments onboard the International Space Station. The predictions of selected microstructural features, such as dendrite arm spacings, show a good agreement with experiments. The computationally-efficient DNN model opens new avenues for investigating the dynamics of large dendritic arrays at length and time scales relevant to solidification experiments and processes.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A Prediction of Uranium(VI) Diffusion in Montmorillonite at Various Chemical Solution Conditions

Most nuclear waste disposal options include compacted bentonite, predominately consisting of montmorillonite clay, as part of a barrier system to minimize contaminant mobility. Uranium (U) is the primary element in spent nuclear fuel, and a potential contaminant of water resources. Uranium sorption onto clay and its slow diffusive transport away from waste canisters are expected to limit U(VI) mobility, and possibly also control nuclear fuel degradation rates based on mass action considerations. A prediction of U(VI) diffusion in these systems is complicated by two main factors: (1) the complex mineralogical structure of montmorillonite clay, leading to two types of clay porosities and surface sites, and (2) the dependence of U(VI) sorption on U(VI) solution speciation and the specific chemical system conditions. For instance, depending on pH and ionic strength, U(VI) may sorb onto montmorillonite due to weak ion exchange reactions or the formation of stable surface complexes. Furthermore, in contrast to conventional porous media, U(VI) diffusion in montmorillonite is not only limited by contaminant sorption, but also by a (potential) lack of access to the full clay porosity for some U(VI) solution species. For instance, at the high degrees of clay compaction expected in future barrier systems, the partial or full exclusion of U(VI) from negatively-charged clay interlayer spaces may substantially decrease U(VI) fluxes. In this modeling study, we specifically investigated how the charge of U(VI) solution species affects U(VI) sorption and diffusion behavior. For this purpose, we first simplified complex chemical solution speciation diagrams to charge speciation diagrams. Then, we coupled an existing U(VI) surface complexation model (Tournassat et al., Geochim. Cosmochim. Acta, 2018) with an analytical solution of the diffusion equation to simulate U(VI) mobility in lab-scale through-diffusion experiments for various chemical conditions. Our results suggest that there are two important parameters that govern U(VI) diffusion in these systems. First, U(VI) surface complexation on clay edge surfaces leads to contaminant retardation. On the field scale, this will determine how early a plume may cross site boundaries and pose risks to drinking water resources. Second, the degree of access of U(VI) solution species to diffuse layers in clay interlayer spaces largely controls U(VI) diffusive fluxes. In the field, this parameter will govern the level of U(VI) concentrations that can be expected in contaminated waters.

Pistorino, Jonathan C.↗
DOE OSTI2023

Modeling the Effect of Hydrogen on Tungsten Grain Boundary Properties

The objective of this study is to assess the effect of hydrogen isotopes on mechanical properties of tungsten that are relevant to first-wall applications in a tokamak. We focus on the effects of hydrogen segregation to grain boundaries and its impact on the grain boundary structure and properties. Through molecular dynamics and other atomistic techniques, we are able to make predictions of grain boundary structure and properties related to it, including grain boundary mobility and impurity diffusivity along the grain boundary. The grain boundary structure can change significantly and abruptly as conditions change, including due to heating and radiation damage. These grain boundary phase transformations can have a pronounced effect on mobility and impurity diffusivity, even when no hydrogen is present. How this behavior changes in the presence of hydrogen is not known. The phase may affect the level of hydrogen segregation and hydrogen storage, and the hydrogen may induce significant grain boundary structure changes. Our objective is to predict whether they do. The results inform thermomechanical modeling of tungsten at tokamak conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simultaneous reproduction of experimental profiles, fluxes, transport coefficients, and turbulence characteristics via nonlinear gyrokinetic profile predictions in a DIII-D ITER similar shape plasma

Experimental conditions obtained on the DIII-D tokamak in the ITER Similar Shape (ISS) have been compared extensively with nonlinear gyrokinetic simulation using the CGYRO code [Candy et al., J. Comput. Phys. 324, 73–93 (2016)] with comparisons spanning ion and electron heat fluxes, electron and impurity particle transport, and turbulent fluctuation levels and characteristics. Bayesian optimization techniques [Rodriguez-Fernandez et al., Nucl. Fusion 62(7), 076036 (2022)], combined with nonlinear gyrokinetics, have been used to obtain simultaneously Qi, Qe, and Γe flux-matched profiles that are found to be in good agreement with experimental profile measurements. Synthetic diagnostics were used to compare measured beam emission spectroscopy and correlation electron cyclotron emission turbulent fluctuations with nonlinear simulation. Although some disagreements exist, nonlinear simulations are found to be in generally good agreement with measured fluctuation levels, spectral shapes, and measured radial trends in low-k δne/ne and δTe/Te. Low (Li and C) and mid-Z (Ca) impurity transport was also compared with these flux-matched simulations. Fully stripped, low-Z impurities are well reproduced by the gyrokinetic modeling while clear disagreement exists in comparisons with mid-Z impurities. Nonlinear gyrokinetic investigation into the Z dependence of impurity transport in the ISS conditions is also performed, demonstrating clear trends of impurity diffusion with impurity Z (both D∝Z and D∝1/Z) that vary with the radial location studied. These trends are shown to result from the local dominance of ion temperature gradient or ∇n driven trapped electron mode turbulence and may contribute to the disagreement between simulation and experiment in mid-Z impurity transport. The results of this work represent one of the most complete validation studies of the gyrokinetic model performed to date and provide an example of new capabilities for predicting performance in future fusion devices.

Physics↗
DOE OSTI2018

Assessment of Error Sources in the Modeling of POA Irradiance Under Clear-Sky Conditions: Preprint

A new physics-based model, the Fast All-sky Radiation Model for Solar applications with Narrowband Irradiances on Tilted surfaces (FARMS-NIT), accurately computes plane-of-array (POA) irradiances in both narrow- and broad wavelength bands within the solar region. The evaluation studies suggest that the new model provides increased accuracy than a model assuming isotropic diffuse solar radiation and an empirical model linking the POA irradiances with long-term surface observations at selected locations. This study anlyzes the surface observations used by FARMS-NIT and is intended to understand the error sources under clear-sky conditions as well as their impact to the performance of FARMS-NIT.

14 SOLAR ENERGY↗
DOE OSTI2015

Submodels of model of nonlinear diffusion in the inhomogeneous medium involving absorption

We study the five-parameter model, describing the process of nonlinear diffusion in an inhomogeneous medium in the presence of absorption, for which the differential equation of the model admits a continuous Lie group of transformations, acting on the set of its solutions. We found six submodels of the original model of nonlinear diffusion, with different symmetry properties. Of these six submodels, the five submodels with transient absorption, for which the absorption coefficient depends on time according to a power law, represent the greatest interest with a mathematical point of view and with the point of view of physical applications. For each of these nonlinear submodels, we obtained formulas for producing new solutions that contain arbitrary constants, and we found all invariant submodels. All essentially distinct invariant solutions describing these invariant submodels are found in an explicit form or are reduced to finding the solution of nonlinear integral equations. The presence of the arbitrary constants in the integral equations that determine these solutions provide new opportunities for analytical and numerical study of boundary value problems for the received submodels and, thus, for the original model of nonlinear diffusion. For the received invariant submodels, we studied diffusion processes for which at the initial moment of the time at a fixed point is specified as a concentration and its gradient or as a concentration and its velocity. Solving of boundary value problems describing these processes is reduced to the solving of nonlinear integral equations. We established the existence and uniqueness of solutions of these boundary value problems under some additional conditions. The obtained results can be used to study the diffusion of substances, diffusion of conduction electrons and other particles, diffusion of physical fields and propagation of heat in inhomogeneous medium, and also to study a turbulence (Leith model, differential approximation to wave turbulence)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

In-Situ Spatial Mapping of Hydrogen in Yttrium Hydrides at LANSCE (FY23 Version, Rev. 1)

This report summarizes the development of neutron imaging capabilities and experimental activities performed at the Los Alamos Neutron Science Center (LANSCE) with the main goal of measuring temperature-driven hydrogen diffusion within bulk-yttrium hydride (YH x ) materials. Yttrium hydride is the leading candidate to serve as a solid neutron moderator in microreactor cores, owing to its high density of hydrogen atoms as well as its superior thermal stability compared to all other metal hydrides. The experimental results and technique developments reported herein support the U.S. Department of Energy Office of Nuclear Energy’s (DOE-NE) Microreactor Program under Technology Maturation. In particular, it addresses the critical need to experimentally validate and verify hydrogen-diffusion models of metal hydrides used in high-temperature microreactor designs by means of high-spatial-resolution neutron imaging. These capabilities were designed to apply large temperature gradients across centimeter-sized YH x pellets to simulate conditions faced in the microreactor environment. In principle, neutron imaging, combined with in-situ sample heating, enables near real-time tracking of hydrogen diffusion in YH x on the sub-millimeter scale. In this report, an overview of neutron imaging methodology and technologies are given in the context of recent spatial measures of hydrogen concentrations in similar metal hydrides. Additionally, the commissioning and operation of a custom-built compact dual-zone furnace is given along with details on three in-situ heating measurements of YH x performed over the 2020 to 2022 LANSCE operation cycles. The aims of these experiments ranged from furnace commissioning, determining sample quality, i.e., hydrogen uniformity via neutron computed tomography, and studying the effects of applied temperature-gradients on YH x pellets. Analyses and results from these neutron imaging measurements are given along with outlooks and guidelines for optimal future hydrogen diffusion measurements. Our conclusions are as follows. Image analyses indicate that centimeter-sized yttrium hydride cylindrical pellets exhibit uniform, whole-body hydrogen desorption and absorption without clear temperature dependence as reflected in the image attenuation at the opposing ends of each sample. This suggests that despite the large magnitude in temperature gradients applied by the furnace heating elements, the sample equilibrates to an unknown intermediate temperature. The origin of this result is likely the combination of short sample length (∼1cm) and use of a TZM can for containment where the latter created a thermal short across the sample. Nevertheless, the results from the most recent measurements indicate that neither significant concentration gradients of hydrogen were formed in centimeter-sized samples through the entire temperature range (25 °C to 950 °C) nor any formed due to temperature gradients on the order of 50 °C/cm up to 700 °C/cm. Furthermore, images from the FY2021 and FY2022 measurements indicate that samples of YH x , fabricated from either the direct hydride or powder metallurgy methods, are highly uniform in their hydrogen concentration to within the measurements’ spatial resolutions. The following questions arise from these latest results: 1) What is the intermediate temperature of the pellets in the TZM cans? 2) How quickly does the temperature equilibrate within the sample? and, 3) Do the observed changes in image attenuation follow known pressure-composition-temperature relations of yttrium hydride?

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Multiscale, mechanistic modeling of irradiation-enhanced silver diffusion in TRISO particles

Tristructural isotropic (TRISO) particles are under consideration for use in several proposed advanced nuclear reactor concepts. The silicon carbide (SiC) layer in TRISO acts as a barrier to prevent the release of the fission products. However, despite remarkable retention, silver (Ag) release has been observed from intact particles, which requires investigation since the Ag isotope ( 110m Ag) has a long half-life. Previous work focused on developing a multiscale, mechanistic model for Ag diffusion accounting for temperature and microstructure effect and has been successfully validated. In this work, we expand the previous model to account for irradiation-enhanced Ag diffusivity in SiC and improve its accuracy over a wider grain size and temperature ranges relevant for advanced reactor conditions. A temperature, grain size, and flux dependent diffusivity is therefore derived using the mesoscale code MARMOT and implemented in the fuel performance code BISON. The irradiation-enhanced Ag diffusivity in SiC is compared against experimental data and validated using BISON against Ag release measurements from the Advanced Gas Reactor Fuel Development and Qualification Program (AGR-1 and AGR-2). Herein, we quantify the impact of SiC grain size, irradiation, and temperature on Ag release. In agreement with previous studies, we find accounting for SiC grain size improves agreement between BISON predictions and experimental observations for most cases. In conclusion, we also find that accounting for irradiation improves agreement for cases where Ag release was underestimated, but the impact was less significant than accounting for microstructure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Conditions for translation and scaling invariance of the neutron diffusion equation

Lie group methods are applied to the time-dependent, monoenergetic neutron diffusion equation in materials with spatial and time dependence. To accomplish this objective, the underlying 2nd order partial differential equation (PDE) is recast as an exterior differential system so as to leverage the isovector symmetry analysis approach. Some of the advantages of this method as compared to traditional symmetry analysis approaches are revealed through its use in the context of a 2nd order PDE. In this context, various material properties appearing in the mathematical model (e.g., a diffusion coefficient and macroscopic cross section data) are left as arbitrary functions of space and time. The symmetry analysis that follows is restricted to a search for translation and scaling symmetries; consequently the Lie derivative yields specific material conditions that must be satisfied in order to maintain the presence of these important similarity transformations. The principal outcome of this work is thus the determination of analytic material property functions that enable the presence of various translation and scaling symmetries within the time-dependent, monoenergetic neutron diffusion equation. The results of this exercise encapsulate and generalize many existing results already appearing in the literature. While the results contained in this work are primarily useful as phenomenological guides pertaining to the symmetry behavior of the neutron diffusion equation under certain assumptions, they may eventually be useful in the construction of exact solutions to the underlying mathematical model. Furthermore, the results of this work are also useful as a starting point or framework for future symmetry analysis studies pertaining to the neutron transport equation and its many surrogates.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Validation of Wichner Predictive Model for Chronic Oxidation by Moisture of Nuclear Graphite

This document reports the results of multi-year extensive research on chronic oxidation of nuclear graphite by traces of moisture that might be present in the helium coolant of high-temperature gas-cooled reactors (HTGRs). The research focuses on four grades of nuclear graphite that are considered potential candidates for component fabrication in advanced nuclear reactors. The four grades are (1) PCEA (medium-fine grain, extruded, manufactured by GrafTech Inc., USA); (2) NBG-17 (medium-fine grain, vibrationally molded, manufactured by SGL Carbon, Germany/France); (3) IG-110 (microfine grain, isostatically molded, manufactured by Toyo Tanso, Japan); and (4) 2114 (microfine grain, isostatically pressed, manufactured by Mersen, USA). The experimental part of this research comprised three objectives that were consistently pursued for each graphite material: (1) accelerated oxidation tests in well-controlled conditions (temperature, gas composition) and analysis of chemical kinetic models and of their parameters for each graphite grade; (2) measurements of water vapor transport rates from H 2 O/He mixtures through thin graphite slabs for estimation of the water effective diffusivity specific to each material; (3) controlled oxidation of graphite specimens in H 2 O/He mixtures and characterization of the effects of oxidation temperatures on the density profiles of oxidized layers. These three groups of experimental results were interpreted using the model developed by Wichner et al. for predicting the penetration depth of oxidants (moisture or oxygen) as a function of oxidation conditions. Basically, starting from the kinetic and transport properties of each graphite grade, the ultimate goal was predicting the stable density profile in the oxidized layer and its dependence on the oxidation temperature. The profiles obtained this way were then compared against direct observations of oxidized layers in specimens obtained in well-controlled oxidation conditions. The general agreement between predictions and observations validated the general applicability of the Wichner model. The results obtained in this project were presented at annual Department of Energy–Advanced Reactor Technologies program reviews and at specialist meetings and conferences and were published in Oak Ridge National Laboratory (ORNL) reports and in professional journals. The results published in 2018 had an immediate impact in 2019, when a group from three US universities used ORNL kinetic data as input information in a complex simulation of chronic oxidation by moisture using the same graphite grades during long-time operation in a modular high-temperature gas-cooled reactor. Their conclusion, in line with the conclusions obtained at the end of this project, was that chronic graphite oxidation by moisture during normal operation conditions of an HTGR will be limited to a thin surface layer of graphite blocks at the reactor bottom, where the highest temperatures are. That will not affect the integrity of graphite components and reactor safety, provided the impurity concentration in the helium coolant is rigorously controlled and maintained at the design specifications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Laser-material interaction during atom probe tomography of oxides with embedded metal nanoparticles

Here, oxide-supported metal nano-particles are of great interest in catalysis but also in the development of new large-spectrum-absorption materials. The design of such nano materials requires three-dimensional characterization with a high spatial resolution and elemental selectivity. The laser assisted Atom Probe Tomography (La-APT) presents both these capacities if an accurate understanding of laser-material interaction is developed. In this paper, we focus on the fundamental physics of field evaporation as a function of sample geometry, laser power, and DC electric field for Au nanoparticles embedded in MgO. By understanding the laser-material interaction through experiments and a theoretical model of heat diffusion inside the sample after the interaction with laser pulse, we point out the physical origin of the noise and determine the conditions to reduce it by more than one order of magnitude, improving the sensitivity of the La-APT for metal-dielectric composites.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Extending the range of validity of Fourier's law into the kinetic transport regime via asymptotic solution of the phonon Boltzmann transport equation

We derive the continuum equations and boundary conditions governing phonon-mediated heat transfer in the limit of a small but finite mean-free path from the asymptotic solution of the linearized Boltzmann equation in the relaxation time approximation. Our method uses the ratio of the mean-free path to the characteristic system length scale, also known as the Knudsen number, as the expansion parameter to study the effects of boundaries on the breakdown of the Fourier description. We show that, in the bulk, the traditional heat conduction equation using Fourier's law as a constitutive relation is valid at least up to second order in the Knudsen number for steady problems and first order for time-dependent problems. Yet, this description does not hold within distances on the order of a few mean-free paths from the boundary; this breakdown is a result of kinetic effects that are always present in the boundary vicinity and require solution of a Boltzmann boundary layer problem to be determined. Matching the inner, boundary layer solution to the outer, bulk solution yields boundary conditions for the Fourier description as well as additive corrections in the form of universal kinetic boundary layers; both are found to be proportional to the bulk-solution gradients at the boundary and parametrized by the material model and the phonon-boundary interaction model (Boltzmann boundary condition). Our derivation shows that the traditional no-jump boundary condition for prescribed temperature boundaries and the no-flux boundary condition for diffusely reflecting boundaries are appropriate only to zeroth order in the Knudsen number; at higher order, boundary conditions are of the jump type. We illustrate the utility of the asymptotic solution procedure by demonstrating that it can be used to predict the Kapitza resistance (and temperature jump) associated with an interface between two materials. All results are validated via comparisons with low-variance deviational Monte Carlo simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Exploration of a Reacting Jet-in-Crossflow in a High-Pressure Axial Stage Combustor

The present work investigates a reacting methane jet in a vitiated cross flow at five atmospheres. The experiments are run at temperatures and velocities relevant to real world gas turbines. High speed chemiluminescence is utilized to obtain flame behavior and visualize flame trajectories. Experimental results are compared with a corresponding RANS CFD model. The characteristics of a diffusion methane jet injected into a vitiated crossflow are explored at various momentum flux ratios and test section inlet boundary conditions. Temperature and velocity measurements are made at the inlet of the test section to ensure uniform boundary conditions for the axial stage. The effect of different momentum flux ratios on flame stabilization mechanisms and jet trajectories are presented. Momentum flux ratios are varied from 0.5-107; at lower momentum flux ratios the diffusion flame is less stable and burns downstream in the shear layer between the jet and crossflow. At higher momentum flux ratios, the diffusion flame burns more upstream and exhibits improved stability compared to the lower momentum flux ratio cases.

03 NATURAL GAS↗
DOE OSTI2014

A phase-field model coupled with lattice kinetics solver for modeling crystal growth in furnaces

In this study, we present a new numerical model for crystal growth in a vertical solidification system. This model takes into account the buoyancy induced convective flow and its effect on the crystal growth process. The evolution of the crystal growth interface is simulated using the phase-field method. Two novel phase-field models are developed to model the crystal growth interface in vertical gradient furnaces with two temperature profile setups: 1) fixed wall temperature profile setup and 2) time-dependent temperature profile setup. A semi-implicit lattice kinetics solver based on the Boltzmann equation is employed to model the unsteady incompressible flow. This model is used to investigate the effect of furnace operational conditions on crystal growth interface profiles and growth velocities. For a simple case of macroscopic radial growth, the phase-field model is validated against an analytical solution. Crystal growth in vertical gradient furnaces with two temperature profile setups have been also investigated using the developed model. The numerical simulations reveal that for a certain set of temperature boundary conditions, the heat transport in the melt near the phase interface is diffusion dominant and advection is suppressed.

Lin, Guang↗
DOE OSTI2018

Use of Physics Based Management Optimization for Groundwater Remedy Design - 18671

Numerical modeling coupled with formal mathematical optimization techniques was used to design a groundwater remedy strategy to address a carbon tetrachloride (CCl{sub 4}) plume covering several square miles. The optimization approach, methodology, and project costs are described. The site consists of a large grain elevator in Nebraska. CCl{sub 4} was used as a fumigant and subsequently impacted soils and groundwater. Remedy design focused on the more diffuse plume, which extended approximately 4.5 miles down gradient of the source area. An eleven-layer model was developed and calibrated to observed groundwater conditions. This model was linked with formal optimization techniques to identify the best design capable of meeting client specified requirements and well location limitations. The planned system included requirements on the location of extraction and injection wells, total system flow rate, number of management periods, and specific constraints limiting further plume expansion. A separate trinitrotoluene (TNT) contaminant plume from a third party was located vertically above the leading edge of the CCl{sub 4} plume within the same aquifer. A further constraint imposed was that the design could not affect or cause additional migration of the TNT plume. Optimal system design was determined using HydroGeoLogic's Physics Based Management Optimization (PBMO{sup TM}) system. Over 38,000 candidate solutions consisting of different combinations of pumping and injection wells locations and rates, were evaluated through flow modeling coupled with either particle tracking or fully convective-dispersive transport simulations to determine the optimal design. Planning for installation of the optimal design is underway. Simulations were completed on a distributed computing platform consisting of 3 Quad Core PCs. Flow and particle tracking simulations required less than 2 minutes per simulation. Coupled flow and transport simulation run times were generally between 30-40 minutes per simulation. The total costs for developing the groundwater flow model and completing the subsequent optimal design were approximately $400,000. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗