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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Energetics of boron near tungsten surfaces: A first-principles study

Interest exists in utilizing boron (B) wall conditioning of fusion tokamaks containing tungsten (W) plasma facing components, in order to improve plasma confinement. To understand the interactions of B with W surfaces, first-principles density functional theory calculations have been performed to model the adsorption, diffusion, and solution of B near the W(100), W(110), and W(111) surfaces. The results show that B within a distance of 0.6 nm above the surfaces is adsorbed to the surfaces without activation barriers. B atoms are strongly adsorbed on the W(100) surface with an adsorption energy of 7.80 eV, which is 1.22 and 1.35 eV larger than on the W(110) and W(111) surfaces. B diffusion on the W(100), W(110), and W(111) surfaces has an activation energies of 2.08, 1.12, and 1.47 eV, respectively; while, diffusion from the adsorption sites into the bulk requires 2.2–2.3 eV. The B solution energy below a clean W(100) surface is the lowest, followed by the W(111) and W(110) surfaces. B clustering and B-induced surface deformation as a function of B coverage has been investigated. B on the W(100) surface occupy epitaxial sites at coverages of 0–1.25 ML, but form clusters at higher coverages. B clustering on the W(110) and W(111) surfaces is expected throughout the adsorption process. Compared to a clean surface, B atoms on the W(100) reduce the surface effect on the B solution energy below the surface, while the presence of B on the W(110) and W(111) surfaces generally decreases or increases the B solution energy below surfaces, respectively.

Yang, L. (ORCID:0000000322166071)↗
DOE OSTI2020

Where did the tumor start? An inverse solver with sparse localization for tumor growth models

In this work, we present a numerical scheme for solving an inverse problem for parameter estimation in tumor growth models for glioblastomas, a form of aggressive primary brain tumor. The growth model is a reaction–diffusion partial differential equation (PDE) for the tumor concentration. We use a PDE-constrained optimization formulation for the inverse problem. The unknown parameters are the reaction coefficient (proliferation), the diffusion coefficient (infiltration), and the initial condition field for the tumor PDE. Segmentation of magnetic resonance imaging (MRI) scans drive the inverse problem where segmented tumor regions serve as partial observations of the tumor concentration. Like most cases in clinical practice, we use data from a single time snapshot. Moreover, the precise time relative to the initiation of the tumor is unknown, which poses an additional difficulty for inversion. We perform a frozen-coefficient spectral analysis and show that the inverse problem is severely ill-posed. We introduce a biophysically motivated regularization on the structure and magnitude of the tumor initial condition. In particular, we assume that the tumor starts at a few locations (enforced with a sparsity constraint on the initial condition of the tumor) and that the initial condition magnitude in the maximum norm is equal to one. We solve the resulting optimization problem using an inexact quasi-Newton method combined with a compressive sampling algorithm for the sparsity constraint. Our implementation uses PETSc and AccFFT libraries. We conduct numerical experiments on synthetic and clinical images to highlight the improved performance of our solver over a previously existing solver that uses standard two-norm regularization for the calibration parameters. The existing solver is unable to localize the initial condition. Our new solver can localize the initial condition and recover infiltration and proliferation. In clinical datasets (for which the ground truth is unknown), our solver results in qualitatively different solutions compared to the two-norm regularized solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Study of The Compression of an Imploding Plasma and the Magnetic- Field Flux Using Advanced, Non-Intrusive, Spectroscopic Techniques

The main motivation of the project was to measure the magnetic field in pulsed power driven plasma compression using advanced non-intrusive spectroscopic techniques. The strength of the project is that it encompasses a combination of laboratory experiments with diagnostics for both the magnetic field and plasma properties and computational modeling. The following issues were addressed: (1) Plasma and magnetic field dynamics, (2) Field compression and diffusion, (3) Energy balance in the system, (4) Onset of instabilities, and (5) Correlation of the plasma and field evolutions with the initial conditions. The unique feature of the research was the use of advanced spectroscopy to measure both the magnetic field and various plasma properties in multiple locations using both side-on and end-on views, simultaneously. The experimental system allowed for imploding different gas mixtures (which facilitated local measurements) and for flexibility in modifying the initial conditions. Two current generators were used with significantly different current rise times (200 ns and 1.6 µs) but at the same current level of 400 kA. The slower pulsed power generator (at rise time 1.6 µs) was based at the Weizmann Institute of Science (WIS) and imploded plasma into which an axial magnetic field (up to 1 T) was initially embedded. The faster current rise time generator (at rise time 200 ns) based at University of California San Diego (UCSD) was used for plasma parameter studies initially, with collaborative efforts we successfully implemented the advanced spectroscopy techniques at UC San Diego. One graduate student (Nick Aybar) and a postdoctoral scholar (Maylis Dozeries) were trained. They both are based at the Lawrence Livermore National Laboratory now.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Conditions for diffusion-limited and reaction-limited recombination in nanostructured solar cells

The performance of Dye-sensitized solar cells (DSC) and related devices made of nanostructured semiconductors relies on a good charge separation, which in turn is achieved by favoring charge transport against recombination. Although both processes occur at very different time scales, hence ensuring good charge separation, in certain cases the kinetics of transport and recombination can be connected, either in a direct or an indirect way. In this work, the connection between electron transport and recombination in nanostructured solar cells is studied both theoretically and by Monte Carlo simulation. Calculations using the Multiple-Trapping model and a realistic trap distribution for nanostructured TiO{sub 2} show that for attempt-to-jump frequencies higher than 10{sup 11}–10{sup 13} Hz, the system adopts a reaction limited (RL) regime, with a lifetime which is effectively independent from the speed of the electrons in the transport level. For frequencies lower than those, and depending on the concentration of recombination centers in the material, the system enters a diffusion-limited regime (DL), where the lifetime increases if the speed of free electrons decreases. In general, the conditions for RL or DL recombination depend critically on the time scale difference between recombination kinetics and free-electron transport. Hence, if the former is too rapid with respect to the latter, the system is in the DL regime and total thermalization of carriers is not possible. In the opposite situation, a RL regime arises. Numerical data available in the literature, and the behavior of the lifetime with respect to (1) density of recombination centers and (2) probability of recombination at a given center, suggest that a typical DSC in operation stays in the RL regime with complete thermalization, although a transition to the DL regime may occur for electrolytes or hole conductors where recombination is especially rapid or where there is a larger dispersion of energies of electron acceptors.

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

Transport Studies and Modeling in PEM Fuel Cells

This project’s aim was to develop fuel cell components (i.e. membranes, gas-diffusion media (GDM), bipolar plates and flow fields) that possess specific properties (i.e. water transport and conductivity). A computational fluid dynamics model was developed to elucidate the effect of certain parameters on these specific properties. Ultimately, the model will be used to determine sensitivity of fuel cell performance to component properties to determine limiting components and to guide research. We have successfully reached our objectives and achieved most of the milestones of this project. We have designed and synthesized a variety of hydrocarbon block polymer membranes with lower equivalent weight, structure, chemistry, phase separation and process conditions. These membranes provide a broad selection with optimized water transport properties. We have also designed and constructed a variety of devices that are capable of accurately measuring the water transport properties (water uptake, water diffusivity and electro-osmatic drag) of these membranes. These transport properties are correlated to the membranes’ structures derived from X-ray and microscopy techniques to determine the structure-property relationship. We successfully integrated hydrocarbon membrane MEAs with a current distribution board (CBD) to study the impact of hydrocarbon membrane on water transport in fuel cells. We have designed and fabricated various GDM with varying substrate, diffusivity and micro-porous layers (MPL) and characterized their pore structure, tortuosity and hydrophobicity. We have derived a universal chart (MacMullin number as function of wet proofing and porosity) that can be used to characterize various GDM. The abovementioned GDMs have been evaluated in operating fuel cells; their performance is correlated to various pore structure, tortuosity and hydrophobicity of the GDM. Unfortunately, determining a universal relationship between the MacMullin number and these properties was not achieved. We have simulated fuel cell performance, current distribution and water distribution at various values of the water uptake, membrane diffusivity, and electro-osmotic drag coefficient (EODC) and compared modeling results with segmented-cell data for both serpentine and parallel flow-fields. We have developed iterations of fuel cell flow fields to achieve specific water transport and thermal management targets. This work demonstrated the importance of membrane diffusivity on fuel cell performance, the necessity of a high membrane diffusion coefficient, and the desirability of a low EODC at low levels of relative humidity.

08 HYDROGEN↗
DOE OSTI2022

Application of BISON to UO 2 MiniFuel fission gas release analysis

There has been a recent push to accelerate fuel qualification by developing revolutionary capabilities to reduce irradiation periods, and thereby, reduce the time required to qualify a new fuel system. One such capability is the MiniFuel irradiation capsule designed to miniaturize fuel samples and irradiate “mini” fuel samples under isothermal temperature conditions. MiniFuel allows steady-state irradiations to decouple the traditionally coupled fission rate (i.e., power) and temperature parameters to understand and generate microstructures observed in fuel operated in a commercial reactor. Furthermore, this process offers the possibility to gather in situ data as well as postirradiation or transient data such as thermal conductivity, specific heat, fission gas diffusion and release, etc. However, accelerating fuel qualification is not solely reliant on generating large amounts of data but also on developing an informed test matrix designed to rapidly generate impactful data. Additionally, this process is reliant on fuel performance codes, such as BISON, to evaluate MiniFuel irradiations using existing material models. This process pinpoints model/data gaps, identifies desired irradiation conditions, and subsequently supports model validation and development. This work describes the use of BISON to perform a number of sensitivity studies designed to understand conditions that lead to fission gas release (FGR) under steady-state isothermal irradiation conditions and temperature transient conditions. The model is applied to a UO 2 MiniFuel example and shows an overall good qualitative agreement with experimental FGR annealing tests under different temperature conditions. It also accounts well for microstructural effects on FGR. When quantitatively compared with FGR data from previously irradiated 103 MWd/kgU UO 2 discs under thermal annealing, the model shows a less satisfactory agreement with the experimental data. Finally, a UO 2 MiniFuel test matrix is proposed to help to extend the model's operational range and validate the new FGR model capabilities to higher burnups and transient conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Flow visualisation in real-size optical injectors of conventional, additised, and renewable gasoline blends

Research on renewable and alternative fuels is crucial for improving the energy and environmental efficiency of modern gasoline internal combustion engines. To highlight the influence of fuel rheological and thermodynamic properties on phase change and atomisation processes, three types of gasoline blends were tested. More specifically, the campaign comprised a reference gasoline, an ethanol/gasoline blend (10% v/v) representative of renewable fuels, and an additised gasoline sample treated with viscoelasticity-inducing agents. High-speed imaging of the transient two-phase flow field arising in the internal geometry and the near-nozzle spray region of gasoline injectors was performed employing Diffuse Backlight Illumination. The metallic body of a commercial injector was modified to fit transparent tips realising two nozzle layouts, namely a two-hole real size model resembling the Engine Combustion Network spray G injector and an enraged replica with an offset hole. Experiments were conducted at realistic operating conditions comprising an injection pressure of 100 bar and ambient pressures in the range of 0.1–6.0 bar to cover the entire range of chamber pressures prevailing in Gasoline Direct Injection engines. The action of viscoelastic additives was verified to have a suppressive effect on in-nozzle cavitation (6% reduction in cavitation extent) , while also enhancing spray atomisation at flash-boing conditions, in a manner resembling the more volatile gasoline/ethanol blends. Finally, persisting liquid ligaments were found to form after the end of injection for the additised sample, owing to the surfactant nature of the additives.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2015

An experimental and kinetic modeling study on dimethyl carbonate (DMC) pyrolysis and combustion

Because of the absence of C–C bonds and the large oxygen content in its molecular structure, dimethyl carbonate (DMC) is a promising oxygenated additive or substitute for hydrocarbon fuels. In order to understand its chemical oxidation and combustion kinetics, flow reactor pyrolysis at different pressures (40, 200 and 1040 mbar) and low-pressure laminar premixed flames with different equivalence ratios (1.0 and 1.5) were investigated. Mole fraction profiles of many reaction intermediates and products were obtained within estimated experimental uncertainties. From theoretical calculations and estimations, a detailed kinetic model for DMC pyrolysis and high-temperature combustion consisting of 257 species and 1563 reactions was developed. The performance of the kinetic model was then analyzed using detailed chemical composition information, primarily from the present measurements. In addition, it was examined against the chemical structure of an opposed-flow diffusion flame, relying on global combustion properties such as the ignition delay times and laminar burning velocities. Furthermore, these extended comparisons yielded overall satisfactory agreement, demonstrating the applicability of the present model over a wide range of high-temperature conditions.

09 BIOMASS FUELS↗
DOE OSTI2019

Multiscale dendritic needle network model of alloy solidification with fluid flow

We present a mathematical formulation of a multiscale model for solidification with convective flow in the liquid phase. The model is an extension of the dendritic needle network approach for crystal growth in a binary alloy. We propose a simple numerical implementation based on finite differences and step-wise approximations of parabolic dendritic branches of arbitrary orientation. Results of the two-dimensional model are verified against reference benchmark solutions for steady, unsteady, and buoyant flow, as well as steady-state dendritic growth in the diffusive regime. Simulations of equiaxed growth under forced flow yield dendrite tip velocities within 10% of quantitative phase-field results from the literature. Finally, we perform illustrative simulations of polycrystalline solidification using physical parameters for an aluminum-10wt% copper alloy. Resulting microstructures show notable differences when taking into account natural buoyancy in comparison to a purely diffusive transport regime. The resulting model opens new avenues for computationally and quantitatively investigating the influence of fluid flow and gravity-induced buoyancy upon the selection of dendritic microstructures. Further ongoing developments include an equivalent formulation for directional solidification conditions and the implementation of the model in three dimensions, which is critical for quantitative comparison to experimental measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Numerical evaluation of AGR-2 fission product release

The AGR-2 experiment produced normal operation and accident condition fuel performance data for tri-structural isotropic (TRISO) particles with UCO and UO 2 kernels. Data for compacts with no failed particles and for compacts with one or more failed particles are available. To model the fission product diffusion and release from these compacts, it is important to consider the computational mesh as well as the heat conduction and mass diffusion properties of the materials in a TRISO particle. Code and solution verification studies with the Bison fuel performance code were performed to show that the code is computing correct solutions. Bison comparisons to post-irradiation examination data and PARFUME predictions were made for several sets of compacts. Furthermore, these comparisons considered silver, cesium, strontium, and krypton with and without failed particles. Bison predictions closely match those of PARFUME, while both codes generally overpredict post-irradiation examination data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Non-ideal magnetohydrodynamic simulations of the two-stage fragmentation model for cluster formation

We model molecular cloud fragmentation with thin-disk, non-ideal magnetohydrodynamic simulations that include ambipolar diffusion and partial ionization that transitions from primarily ultraviolet-dominated to cosmic-ray-dominated regimes. These simulations are used to determine the conditions required for star clusters to form through a two-stage fragmentation scenario. Recent linear analyses have shown that the fragmentation length scales and timescales can undergo a dramatic drop across the column density boundary that separates the ultraviolet- and cosmic-ray-dominated ionization regimes. As found in earlier studies, the absence of an ionization drop and regular perturbations leads to a single-stage fragmentation on pc scales in transcritical clouds, so that the nonlinear evolution yields the same fragment sizes as predicted by linear theory. However, we find that a combination of initial transcritical mass-to-flux ratio, evolution through a column density regime in which the ionization drop takes place, and regular small perturbations to the mass-to-flux ratio is sufficient to cause a second stage of fragmentation during the nonlinear evolution. Cores of size ∼0.1 pc are formed within an initial fragment of ∼pc size. Regular perturbations to the mass-to-flux ratio also accelerate the onset of runaway collapse.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SIMULATIONS AND THEORY OF ION INJECTION AT NON-RELATIVISTIC COLLISIONLESS SHOCKS

We use kinetic hybrid simulations (kinetic ions-fluid electrons) to characterize the fraction of ions that are accelerated to non-thermal energies at non-relativistic collisionless shocks. We investigate the properties of the shock discontinuity and show that shocks propagating almost along the background magnetic field (quasi-parallel shocks) reform quasi-periodically on ion cyclotron scales. Ions that impinge on the shock when the discontinuity is the steepest are specularly reflected. This is a necessary condition for being injected, but it is not sufficient. Also, by following the trajectories of reflected ions, we calculate the minimum energy needed for injection into diffusive shock acceleration, as a function of the shock inclination. We construct a minimal model that accounts for the ion reflection from quasi-periodic shock barrier, for the fraction of injected ions, and for the ion spectrum throughout the transition from thermal to non-thermal energies. This model captures the physics relevant for ion injection at non-relativistic astrophysical shocks with arbitrary strengths and magnetic inclinations, and represents a crucial ingredient for understanding the diffusive shock acceleration of cosmic rays.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Implications of Microstructure in Helium-Implanted Nanocrystalline Metals

Helium bubbles are known to form in nuclear reactor structural components when displacement damage occurs in conjunction with helium exposure and/or transmutation. If left unchecked, bubble production can cause swelling, blistering, and embrittlement, all of which substantially degrade materials and—moreover—diminish mechanical properties. On the mission to produce more robust materials, nanocrystalline (NC) metals show great potential and are postulated to exhibit superior radiation resistance due to their high defect and particle sink densities; however, much is still unknown about the mechanisms of defect evolution in these systems under extreme conditions. Here, the performances of NC nickel (Ni) and iron (Fe) are investigated under helium bombardment via transmission electron microscopy (TEM). Bubble density statistics are measured as a function of grain size in specimens implanted under similar conditions. While the overall trends revealed an increase in bubble density up to saturation in both samples, bubble density in Fe was over 300% greater than in Ni. To interrogate the kinetics of helium diffusion and trapping, a rate theory model is developed that substantiates that helium is more readily captured within grains in helium-vacancy complexes in NC Fe, whereas helium is more prone to traversing the grain matrices and migrating to GBs in NC Ni. Our results suggest that (1) grain boundaries can affect bubble swelling in grain matrices significantly and can have a dominant effect over crystal structure, and (2) an NC-Ni-based material can yield superior resistance to irradiation-induced bubble growth compared to an NC-Fe-based material and exhibits high potential for use in extreme environments where swelling due to He bubble formation is of significant concern.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Model for radiation damage-induced grain subdivision and its influence in U 3 Si 2 fuel swelling

Here, the swelling mechanisms of U 3 Si 2 from neutron irradiation under reactor conditions are not unequivocally known. The limited experimental evidence that is available suggests that the main driver of the swelling in this material would be fission gases accumulation at crystalline grain boundaries. The steps that lead to the accumulation of fission gases at these locations are multiple and complex. However, gradually, the gaseous fission products migrate by diffusion. Upon reaching a grain boundary, which acts as a trap, the gaseous fission products start to accumulate, thus leading to formation of bubbles and hence to induced swelling. Therefore, a quantitative model of swelling requires the incorporation of phenomena that increase the presence of grain boundaries and decrease grain sizes, thus creating sites for bubble formation and growth while concurrently deceasing the distance from gas species formation to gas accumulation sites. It is assumed that grain boundary formation results from the conversion of stored energy from accumulated dislocations into energy for the formation of new grain boundaries. This work attempts to develop a quantitative model for grain subdivision in U 3 Si 2 based on the above mentioned phenomena to quantitatively corroborate the presence of such a mechanism. Then the model is used in conjunction with swelling codes to evaluate the total swelling of a fuel pellet during its lifetime in the I 2 S-LWR reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Evaluation of the toroidal torque driven by external non-resonant non-axisymmetric magnetic field perturbations in a tokamak

The toroidal torque driven by external non-resonant magnetic perturbations (neoclassical toroidal viscosity) is an important momentum source affecting the toroidal plasma rotation in tokamaks. The well-known force-flux relation directly links this torque to the non-ambipolar neoclassical particle fluxes arising due to the violation of the toroidal symmetry of the magnetic field. Here, a quasilinear approach for the numerical computation of these fluxes is described, which reduces the dimension of a standard neoclassical transport problem by one without model simplifications of the linearized drift kinetic equation. The only limiting condition is that the non-axisymmetric perturbation field is small enough such that the effect of the perturbation field on particle motion within the flux surface is negligible. Therefore, in addition to most of the transport regimes described by the banana (bounce averaged) kinetic equation also such regimes as, e.g., ripple-plateau and resonant diffusion regimes are naturally included in this approach. Based on this approach, a quasilinear version of the code NEO-2 [W. Kernbichler et al., Plasma Fusion Res. 3, S1061 (2008).] has been developed and benchmarked against a few analytical and numerical models. Results from NEO-2 stay in good agreement with results from these models in their pertinent range of validity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Simulations of Cyclic Voltammetry for Electric Double Layers in Asymmetric Electrolytes: A Generalized Modified Poisson-Nernst-Planck Model

This paper presents a generalized modified Poisson–Nernst–Planck (MPNP) model derived from first principles based on excess chemical potential and Langmuir activity coefficient to simulate electric double-layer dynamics in asymmetric electrolytes. The model accounts simultaneously for (1) asymmetric electrolytes with (2) multiple ion species, (3) finite ion sizes, and (4) Stern and diffuse layers along with Ohmic potential drop in the electrode. It was used to simulate cyclic voltammetry (CV) measurements for binary asymmetric electrolytes. The results demonstrated that the current density increased significantly with decreasing ion diameter and/or increasing valency |z i | of either ion species. By contrast, the ion diffusion coefficients affected the CV curves and capacitance only at large scan rates. Dimensional analysis was also performed, and 11 dimensionless numbers were identified to govern the CV measurements of the electric double layer in binary asymmetric electrolytes between two identical planar electrodes of finite thickness. A self-similar behavior was identified for the electric double-layer integral capacitance estimated from CV measurement simulations. Two regimes were identified by comparing the half cycle period τ CV and the “RC time scale” τ RC corresponding to the characteristic time of ions’ electrodiffusion. For τ RC ← τ CV , quasi-equilibrium conditions prevailed and the capacitance was diffusion-independent while for τ RC → τ CV , the capacitance was diffusion-limited. The effect of the electrode was captured by the dimensionless electrode electrical conductivity representing the ratio of characteristic times associated with charge transport in the electrolyte and that in the electrode. The model developed here will be useful for simulating and designing various practical electrochemical, colloidal, and biological systems for a wide range of applications.

Wang, Hainan↗
DOE OSTI2021

NRC Multiphysics Analysis Capability Deployment FY21: Part 3

This report details the progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” The deliverables completed for this report are: Deliverable 1c: the capability to model gas mixtures was added to Pronghorn. A test problem mimicking the conditions achieved in a depressurized loss of forced cooling (DLOFC) event was solved with both RELAP-5 and Pronghorn. Pronghorn employed a finite vol ume method with the Kurganov-Tadmor discretization. The comparison between the mass fraction spatial profiles computed with RELAP-5 and Pronghorn clearly shows the presence of numerical artifacts (i.e., overly diffusive behavior at low Mach numbers). We confirmed that the problem disappears at higher Mach numbers. We recommend future work on the implementation of a low Mach finite volume formulation to better treat low Mach number problems. Deliverable 2a: we demonstrated two approaches to model the radiation/conduction/natural convection heat transfer across a stagnant gas for the PBMR-400 design using Pronghorn. The first approach is based on the net radiation method, which relies on the computation of view factors with the Multiphysics Object-Oriented Simulation Environment (MOOSE) ray tracing capability. The second method is a traditional thermal resistance approach. The test problems include both 2D and 3D geometries. In all cases, the results show very good agreement during a DLOFC transient. This confirms that the faster thermal resistance method produces solutions that are equivalent to the net radiation method for this geometry. Deliverable 3d: we demonstrated the use of the advection kernel for the delayed neutron precursor equation in Griffin with a 2D MSFR model. The results appear physical but further verification is recommended. We also recommend the addition of conjugate heat transfer to compute the temperatures and model the thermomechanic behavior of the reflectors and other structures. Significant memory and performance issues were encountered in the 3D axisymmetric model. Future work is recommended in this area. Task 8g: this task allows multidimensional MOOSE applications to be coupled to system codes (RELAP-7 and SAM). We implemented a faster multiphysics iteration coupling algorithm, which provides an overall 6× acceleration of the 3D-1D coupling of the core multidi- mensional fluid flow solver and the 1D primary and secondary loop model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

Characterizing the structure of diffuse emission in Hi-GAL maps

We present a study of the structure of the Galactic interstellar medium (ISM) through the Δ-variance technique, related to the power spectrum and the fractal properties of infrared/submillimeter maps. Through this method, it is possible to provide quantitative parameters, which are useful for characterizing different morphological and physical conditions, and better constraining the theoretical models. In this respect, the Herschel Infrared Galactic Plane Survey, carried out at five photometric bands from 70 to 500 μm, constitutes a unique database for applying statistical tools to a variety of regions across the Milky Way. In this paper, we derive a robust estimate of the power-law portion of the power spectrum of four contiguous 2° × 2° Hi-GAL tiles located in the third Galactic quadrant (217° ≲ ℓ ≲ 225°, –2° ≲ b ≲ 0°). The low level of confusion along the line of sight, testified by CO observations, makes this region an ideal case. We find very different values for the power spectrum slope from tile to tile but also from wavelength to wavelength (2 ≲ β ≲ 3), with similarities between fields attributable to components located at the same distance. Thanks to comparisons with models of turbulence, an explanation of the determined slopes in terms of the fractal geometry is also provided, and possible relations with the underlying physics are investigated. In particular, an anti-correlation between ISM fractal dimension and star formation efficiency is found for the two main distance components observed in these fields. A possible link between the fractal properties of the diffuse emission and the resulting clump mass function is discussed.

79 ASTRONOMY AND ASTROPHYSICS↗