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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Computation of diffuse sky irradiance from multidirectional radiance measurements

Accurate determination of the diffuse solar spectral irradiance directly above the land surface is important in characterizing the reflectance properties of these surfaces, especially vegetation canopies. This determination is also needed to infer the net radiation budget of the earth-atmosphere system above these surfaces. An algorithm is developed here for the computation of hemispheric diffuse irradiance using the measurements from an instrument called PARABOLA, which rapidly measures upwelling and downwelling radiances in three selected wavelength bands. The validity of the algorithm is established from simulations. The standard reference data set of diffuse radiances of Dave (1978), obtained by solving the radiative transfer equation numerically for realistic atmospheric models, is used to simulate PARABOLA radiances. Hemispheric diffuse irradiance is estimated from a subset of simulated radiances by using the algorithm described. The algorithm is validated by comparing the estimated diffuse irradiance with the true diffuse irradiance of the standard data set. The validations include sensitivity studies for two wavelength bands (visible, 0.65-0.67 micron; near infrared, 0.81-0.84 micron), different atmospheric conditions, solar elevations, and surface reflectances. In most cases the hemispheric diffuse irradiance computed from simulated PARABOLA radiances and the true irradiance obtained from radiative transfer calculations agree within 1-2 percent. This technique can be applied to other sampling instruments designed to estimate hemispheric diffuse sky irradiance.

Ahmad, Suraiya P.↗
NASA NTRS1992

Effects of buoyancy-driven flow and thermal boundary conditions on physical vapor transport

A 2D numerical model was developed in order to ascertain if reduced gravity conditions are beneficial to physical vapor transport (PVT) and to determine its tolerance limits to residual accelerations. This was solved using the PHOENICS finite-volume code. Reduction of gravitational accelerations to less than 0.1 g0 was found to be sufficient to suppress buoyancy-driven convection to an extent that diffusion was the dominant transport mode, whence a greater uniformity in the growth rate could be obtained. It is shown that a uniform temperature gradient on the ampoule walls causes the vapor to be supersaturated throughout the ampoule, potentially resulting in undesirable nucleation at the walls. A 'hump' in the wall temperature profile can be used to avoid this. The prevailing transport conditions determine the size of the hump needed.

Nadarajah, Arunan↗
NASA NTRS1999

The Effects of Buoyancy and Dilution on the Structure and Lift-off of Coflow Laminar Diffusion Flames

The ability to predict the coupled effects of complex transport phenomena with detailed chemical kinetics in diffusion flames is critical in the modeling of turbulent reacting flows and in understanding the processes by which soot formation and radiative transfer take place. In addition, an understanding of the factors that affect flame extinction in diffusion flames is critical in the suppression of fires and in improving engine efficiency. The goal of our characterizations of coflow laminar diffusion flames is to bring to microgravity the multidimensional diagnostic tools available in normal gravity, and in so doing provide a broader understanding of the successes and limitations of current combustion models. This will lead to a more detailed understanding of the interaction of convection, diffusion and chemistry in both buoyant and nonbuoyant environments. As a sensitive marker of changes in the flame shape, the number densities of excited-state CH (A(exp 2)delta, denoted CH*), and excited-state OH (A(exp 2)Sigma, denoted OH*) are measured in mu-g and normal gravity. Two-dimensional CH* and OH* number densities are deconvoluted from line-of-sight chemiluminescence measurements made on the NASA KC-135 reduced-gravity aircraft. Measured signal levels are calibrated, post-flight, with Rayleigh scattering. Although CH* and OH* kinetics are not well understood, the CH*, OH*, and ground-state CH distributions are spatially coincident in the flame anchoring region. Therefore, the ground-state CH distribution, which is easily computed, and the readily measured CH*/OH* distributions can be used to provide a consistent and convenient way of measuring lift-off height and flame shape in the diffusion flame under investigation. Given that the fuel composition affects flame chemistry and that buoyancy influences the velocity profile of the flow, we have the opportunity to computationally and experimentally study the roles of fluids and chemistry. In performing this microgravity study, improvements to the computational model have been made and new calculations performed for a range of gravity and flow conditions. Furthermore, modifications to the experimental approach were required as a consequence of the constraints imposed by existing microgravity facilities. Results from the computations and experiments are presented in the following sections.

Walsh, Kevin T.↗
NASA NTRS2010

The Role of Kinetic Processes in Magnetospheric Dynamics

The question of whether the micro scale controls the macro scale or vice-versa remains one of the most challenging problems in plasmas. A particular topic of interest within this context is collisionless magnetic reconnection, where both points of views are espoused by different groups of researchers. This presentation will focus on this topic. We will begin by analyzing the global properties of magnetic reconnection. Following that, we discuss properties of electron diffusion region dynamics both for guide field and antiparallel reconnection, and how they can be scaled to different inflow conditions. As a next step, we will study typical temporal variations of the microscopic dynamics with the objective of understanding the potential for secular changes to the macroscopic system. The research will be based on a combination of analytical theory and numerical modeling.

Hesse, Michael↗
DOE OSTI2022

Analytical homogenization techniques applied to the Fickian diffusion: Effective diffusivity coefficient

For multiple applications in nuclear energy, the ability to accurately represent material behavior with a simplified model is important to facilitate practical engineering-scale simulations. In this work, we focus on the homogenized thermal response of a medium containing spherical inclusions, similar to a fuel form (compact or pebble) containing TRISO particles. An extensive survey on effective thermal conductivity modeling was performed in our previous study, considering a random distribution of mono-sized spherical inclusions in a continuous matrix. Using the analogy between heat conduction and the simplified Fickian diffusion (or fission product species conservation), we can use the same analytical homogenization methods to obtain ETC as for the effective diffusivity coefficient (EDC). We performed several numerical experiments at varying conditions to assess the validity of our hypothesis for EDC calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Connecting the failure of K-theory inside and above vegetation canopies and ejection-sweep cycles by a large eddy simulation

Parameterizations of biosphere-atmosphere interaction processes in climate models and other hydrological applications require characterization of turbulent transport of momentum and scalars between vegetation canopies and the atmosphere, which is often modeled using a turbulent analogy to molecular diffusion processes. However, simple flux-gradient approaches (K-theory) fail for canopy turbulence. One cause is turbulent transport by large coherent eddies at the canopy scale, which can be linked to sweep-ejection events, and bear signatures of non-local organized eddy motions. K-theory, that parameterizes the turbulent flux or stress proportional to the local concentration or velocity gradient, fails to account for these non-local organized motions. The connection to sweep-ejection cycles and the local turbulent flux can be traced back to the turbulence triple moment $\overline{C'W'W'}$. In this work, we use large-eddy simulation to investigate the diagnostic connection between the failure of K-theory and sweep-ejection motions. Analyzed schemes are quadrant analysis (QA) and a complete and incomplete cumulant expansion (CEM and ICEM) method. The latter approaches introduce a turbulence timescale in the modeling. Furthermore, we find that the momentum flux needs a different formulation for the turbulence timescale than the sensible heat flux. In conclusion, accounting for buoyancy in stratified conditions is also deemed to be important in addition to accounting for non-local events to predict the correct momentum or scalar fluxes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

Modeling of a Non-Dilute Alloy Solidification Under Terrestrial and Microgravity Conditions

A general finite element model of non-dilute alloy solidification was used to simulate growth of HgCdTe from the melt. 2-D and 3-D models with and without Oberbeck-Boussinesq (O-B) approximation were used. Verification of the 3-D model was based on experimental data. The location and the shape of the interface were determined using melting temperatures obtained from the HgCdTe liquidus curve. Double-diffusive convection in the melt is caused by radial temperature gradients and by rejection at the solid-melt interface of the heavier HgTe-rich solute. Because of these complicating factors, dimensional rather than non-dimensional modeling was performed. Models based on complete set of Navier-Stokes equations for compressible media are superior to the O-B approximation in microgravity, because they provide more accurate estimation of convection.

Bune, Andris V.↗
NASA NTRS1982

Airfoil deposition model

The methodology to predict deposit evolution (deposition rate and subsequent flow of liquid deposits) as a function of fuel and air impurity content and relevant aerodynamic parameters for turbine airfoils is developed in this research. The spectrum of deposition conditions encountered in gas turbine operations includes the mechanisms of vapor deposition, small particle deposition with thermophoresis, and larger particle deposition with inertial effects. The focus is on using a simplified version of the comprehensive multicomponent vapor diffusion formalism to make deposition predictions for: (1) simple geometry collectors; and (2) gas turbine blade shapes, including both developing laminar and turbulent boundary layers. For the gas turbine blade the insights developed in previous programs are being combined with heat and mass transfer coefficient calculations using the STAN 5 boundary layer code to predict vapor deposition rates and corresponding liquid layer thicknesses on turbine blades. A computer program is being written which utilizes the local values of the calculated deposition rate and skin friction to calculate the increment in liquid condensate layer growth along a collector surface.

Kohl, F. J.↗
NASA NTRS2007

Simulation of SEU Cross-sections using MRED under Conditions of Limited Device Information

This viewgraph presentation reviews the simulation of Single Event Upset (SEU) cross sections using the membrane electrode assembly (MEA) resistance and electrode diffusion (MRED) tool using "Best guess" assumptions about the process and geometry, and direct ionization, low-energy beam test results. This work will also simulate SEU cross-sections including angular and high energy responses and compare the simulated results with beam test data for the validation of the model. Using MRED, we produced a reasonably accurate upset response model of a low-critical charge SRAM without detailed information about the circuit, device geometry, or fabrication process

Lauenstein, J. M.↗
DOE OSTI2022

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗
DOE OSTI2020

First-principles equations of state and structures of liquid metals in multi-megabar conditions

Liquid metals at extreme pressures and temperatures are widely interested in the high-pressure community. Based on density functional theory molecular dynamics, we conduct first-principles investigations on the equation of state ( EOS) and structures of four metals (Cu, Fe, Pb, and Sn) at 1.5–5 megabar conditions and 5 × 10 3 –4 × 10 4 K. Our first-principles EOS data enable evaluating the performance of four EOS models in predicting Hugoniot densities and temperatures of the four systems. We find the melting temperature of Cu is 1000–2000 K higher and shows a similar Clapeyron slope, in comparison to those of Fe. Our structure, coordination number, and diffusivity analysis indicates all the four liquid metals form similar simple close-packed structures. Our results set theoretical benchmarks for EOS development and structures of metals in their liquid states and under dynamic compression.

36 MATERIALS SCIENCE↗
NASA NTRS1977

Two phase choke flow in tubes with very large L/D

Data were obtained for two phase and gaseous choked flow nitrogen in a long constant area duct of 16200 L/D with a diverging diffuser attached to the exit. Flow rate data were taken along five isotherms (reduced temperature of 0.81, 0.96, 1.06, 1.12, and 2.34) for reduced pressures to 3. The flow rate data were mapped in the usual manner using stagnation conditions at the inlet mixing chamber upstream of the entrance length. The results are predictable by a two phase homogeneous equilibrium choking flow model which includes wall friction. A simplified theory which in essence decouples the long tube region from the high acceleration choking region also appears to predict the data resonably well, but about 15 percent low.

Hendricks, R. C.↗
NASA NTRS1986

Simulation of the thermospheric tides by use of the NCAR thermospheric general circulation model

Numerical calculations of the thermospheric tidal winds and temperatures at equinox are presented. The calculations were made using the National Center for Atmospheric Research Thermospheric General Circulation Model (TGCM), which includes the effects of viscosity, conductivity, diffusion, ion drag, winds, and temperature gradients. The thermospheric diurnal and semidiurnal tides are excited in situ by solar heating and by ion-neutral momentum coupling. The semidiurnal tidal calculations also include the effects of upward propagating waves generated by heating in the lower atmosphere. This semidiurnal propagating component is modeled by use of the classical tidal perturbations as lower boundary conditions. The model is tuned by adjusting the propagating tidal forcing term until calculated semidiurnal wind and temperature fields best approximate incoherent scatter observations. The tidal TGCM results are consistent with previous theoretical work and successfully reproduce high altitude temperature and meridional velocity data, but they give significantly lower magnitudes for velocities and temperatures near 160 km than are seen by observations.

Fensen, C. G.↗
NASA NTRS1993

Modeling of outgassing and matrix decomposition in carbon-phenolic composites

A new release rate equation to model the phase change of water to steam in composite materials was derived from the theory of molecular diffusion and equilibrium moisture concentration. The new model is dependent on internal pressure, the microstructure of the voids and channels in the composite materials, and the diffusion properties of the matrix material. Hence, it is more fundamental and accurate than the empirical Arrhenius rate equation currently in use. The model was mathematically formalized and integrated into the thermostructural analysis code CHAR. Parametric studies on variation of several parameters have been done. Comparisons to Arrhenius and straight-line models show that the new model produces physically realistic results under all conditions.

Mcmanus, Hugh L.↗
DOE OSTI2020

Challenges in simulating ground interacting nuclear explosions

This paper summarizes recent above-ground nuclear explosion simulations as part of a broader effort to better characterize conditions within a fireball that may influence the chemical evolution of bomb materials and other materials entrained from the local explosion environment. A critical component of this work is validation against historic footage of atmospheric testing, requiring that we understand how the frequency-dependent sensitivity of the utilized film footage influences data captured in such images. We focus first on the early physics of a nuclear explosion in the atmosphere before discussing some of the technical challenges we seek to capture in late-time models that include more complex emplacement conditions and subsurface features. We discuss required physics packages (compressible hydrodynamics, radiation transport, as well as necessary ancillary tables such as equations of state (EOS) and opacities). Additionally, we note reasonable “shortcuts” one may make and their limitations, e.g., using ideal gas EOS, replacing spectrally resolved radiation with spectrally averaged radiation, and exchanging deterministic transport with diffusion. We then discuss an approach to achieving an equilibrated initial stress state for problems where buoyancy and subsurface lithostatic stress are important. Our methodology is presented in the context of LLNL’s ALE3D multiphysics code but may readily be implemented in other codes. In this paper, we start with a description of the challenges of NUDET simulations, followed by a presentation of the simulated intensity (flux) as it would appear on an analysis of the Dixie test. We then progressively introduce additional complexity in subsequent sections (near-surface burst and gravity initialization) before concluding.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

A new method for operating a continuous flow diffusion chamber to investigate immersion freezing: assessment and performance study

Glaciation in mixed-phase clouds predominately occurs through the immersion freezing mode where ice nucleating particles (INPs) immersed within supercooled droplets induce nucleation of ice. Currently, model representations of this process are a large source of uncertainty in simulating cloud radiative properties, and to constrain these estimates, continuous flow diffusion chamber (CFDC)-style INP devices are commonly used to assess the immersion freezing efficiencies of INPs. In this study, a new approach was explored to operating such an ice chamber that provides maximum activation of particles without droplet breakthrough and correction factor ambiguity to obtain high-quality INP measurements in a manner that has not been demonstrated as possible previously. The conditioning section of the chamber was maintained at ~ –20 °C and water relative humidity (RH w ) ~ 113 % conditions to maximize the droplet activation, and the droplets were supercooled with an independently temperature-controlled nucleation section at a steady cooling rate (0.5 °C min –1 ) to induce the freezing of droplets and evaporation of unfrozen droplets. The performance of the modified ice chamber was evaluated using four INP species: K-feldspar, illite-NX, Argentinian soil dust, and airborne arable dust that had shown ice nucleation over a wide span of supercooled temperatures. Dry dispersed and size-selected K-feldspar particles were generated in the laboratory. Illite-NX and soil dust particles were sampled during the second phase of the Fifth International Ice Nucleation Workshop (FIN-02) campaign, and airborne arable dust particles were sampled from the aerosol inlet located on the rooftop of the laboratory. The measured ice nucleation efficiencies of model aerosols with a surface active site density (n s ) metric were higher, but mostly agreed within one order of magnitude compared to literature results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Modeling the influence of nozzle-generated turbulence on diesel sprays

The physical mechanisms governing spray breakup in direct injection engines, such as aerodynamic induced instabilities and nozzle-generated cavitation and turbulence, are not well understood due to the experimental and computational limitations in resolving these processes. Recent x-ray and visible extinction measurements have been con-ducted with a targeted interest in the spray formation region in order to characterize the distribution of droplet sizes throughout the spray. Detailed analysis of these measurements shows promise of yielding insight into likely mechanisms governing atomization, which can inform the improvement of spray models for engine computational fluid dynamic (CFD) codes. In order to investigate potential atomization mechanisms, we employ a joint experimental and computational approach to characterize the structure of the spray formation region using the Engine Combustion Network Spray D injector. X-ray tomography, radiography and ultra-small angle x-ray scattering measurements conducted at the Advanced Photon Source at Argonne National Laboratory quantify the injector geometry, liquid fuel mass and Sauter mean diameter (SMD) distributions under non-vaporizing conditions. Diffused back-illumination imaging measurements, conducted at the Georgia Institute of Technology, characterize the asymmetry of the spray structure. The selected range of injection pressures (50 – 150 MPa) and ambient densities (1.2 – 22.8 kg/m3) allow for the influence of aerodynamic forces on the spray to be studied in a controlled and systematic manner, while isolating the atomization process from the effects of vaporization. In comparison to high ambient density conditions, the spray is observed to be more asymmetric at low ambient density conditions. Although several mechanisms may cause asymmetries in the nozzle exit flow conditions and ultimately the spray distribution, irregularities in the internal nozzle geometry were identified, suggesting an increased sensitivity of the spray structure to internal nozzle surface finish imperfections at such conditions. The presence of these asymmetries may influence the ability to interpret line-of-sight measurements and their derived SMD values and trends from a single viewing angle of the spray. With this consideration in mind, the measured local sensitivities to ambient density suggest that for ambient densities less than 2.4 kg/m3, aerodynamic effects are likely suppressed, allowing the influence of turbulent-induced breakup to be isolated. In concert with the experimental measurements, we utilize three-dimensional, CFD Lagrangian-Eulerian spray simulations in CONVERGE to evaluate the details of the predicted spray structure. In particular, we compare measured and predicted sensitivities of the SMD distribution to changes in injection and ambient conditions from three different atomization models, namely Kelvin Helmholtz (KH), KH Aerodynamics Cavitation Turbulence (KH-ACT), and the newly developed KH-Faeth hybrid model. While none of the existing hybrid spray models were able to replicate the experimentally observed sensitivities, it was found that the scales characterizing the KH-Faeth model show promise of capturing the experimentally observed trends if the effects of secondary droplet breakup are neglected. These results inform recommendations for future experiments and computational studies that can guide the development of an improved spray breakup model.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Calibration of the Diffusivity Predictions of Centipede Using Approximate Bayesian Computation and Applications in Nyx (Engineering Scale) and Xolotl-MARMOT (Meso-Scale) Simulations

Fission gas evolution and release in UO 2 nuclear fuel are important fuel performance metrics and occur in several distinct stages: 1) nucleation, growth and resolution of intra-granular bubbles, 2) diffusion to grain boundaries and 3) nucleation and growth of bubbles at grain boundaries, which eventually form a connected network (percolation) enabling release of gas from grain boundaries through connections to triple junctions, grain edges or free surfaces. The NE-SciDAC project is developing several computational tools to model this problem, which are connected in a hierarchical multi-scale framework. The information transfer in the multi-scale framework is a critical step that, in addition to best-estimates, should include uncertainty quantification. Despite taking a first-principles multi-scale approach, there is a need to perform parameter calibration to ensure consistency with available experimental data. In the present study, uncertainty quantification (UQ) and parameter calibration is demonstrated for one of the lower length scale codes in the multi-scale framework (Centipede) and then the results, including instances of the propagated uncertainties, are used in other codes within the framework, specifically Nyx and Xolotl-MARMOT. We calibrated the model parameters in Centipede, a computer code used to predict diffusivities of uranium (U) and xenon (Xe) in the context of the simulation of fission gas in uranium oxide (UO 2 ) nuclear fuel. The Centipede code depends on 183 parameters, all of which are subject to uncertainty. The three data sets used in our calibration effort are taken from the literature. This data is available as a set of measurements, including measurement errors. Our goal is to calibrate a statistical model that predicts both the value of the measurement and the uncertainty associated with the measurement. We perform a Bayesian calibration of the model parameters using a dedicated approximate Bayesian computation (ABC) likelihood function. To avoid excessive computational costs, we replace the expensive Centipede simulation code by a higher-order surrogate model, constructed using only the 9 most important parameters. These important parameters are identified by a preliminary global sensitivity analysis (GSA) study. Among the important parameters are T0 (the temperature at which UO 2 is perfectly stoichiometric) and Hf_pO2 (the temperature dependence of the oxygen (O) partial pressure) that should be considered as operating conditions to be estimated along with the other parameters. We consider two different cases: one where we define one set of these operating conditions for all data sets, and one where we define distinct operating condition parameters for each data set. The Xe diffusivities predicted by the latter case show distinct features that could not be observed in the former. Next, we use the diffusivity predictions by Centipede as input to Nyx, a reduced order fuel performance code focused on gas behavior alone, in order to estimate quantities associated with inter-granular bubble formation at conditions specified by the experiments. Finally, the diffusivities obtained from the calibrated Centipede runs were used in coupled Xolotl-MARMOT simulations of intra- and inter-granular gas evolution. The results are compared to simulations using the baseline diffusivities from Turnbull et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗