Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conditional distributions”

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 487 records · Page 27

Particle size and velocity measurement in flames by laser anemometer

Simultaneous droplet size and velocity measurements by a particle counting Laser Doppler Anemometer (LDA) in kerosene fuel sprays under burning and non-burning conditions are presented. Particle sizes are derived from pulse height analysis of the mean LDA signals and velocities are simultaneously determined by measuring Doppler shift frequencies. The measurements show that droplet velocity is a function of droplet diameter for burning and non-burning conditions, and spatially averaged size distributions are derived from velocity data. A comparison of results obtained under burning and non-burning conditions show changes in size distribution due to preferential vaporization of small droplets, acceleration due to thermal expansion of gases, and corresponding changes in droplet momentum.

Chigier, N. A.↗

Buried Ducts: A Cost-Effective Technical Solution to High Performance Homes

Traditionally, attics have been insulated at the flat ceiling with deliberate venting to manage heat and moisture. Having all the duct air distribution system in the conditioned space is the best building science approach, however this is not always the preferable solution for builders that like to use the roof space for air distribution. Ducts in attics are exposed to very intense heating and cooling conditions. Not only do we have thermal conduction, convection and radiation exchanges occurring but important duct leakage losses that create significant energy balance losses. In the recent IECC 2018 code, deeply buried attics have been credited for reducing the energy losses. Due to the complex thermal performance in each zone adjacent to the buried ductwork these credits seem to be not consistent with the benefits achieve by buring the ductwork in insulation. The thermal performance of buried ducts was examined in both 2-D and 3-D to determine the actual improvement as a function of duct size, duct insulation and attic insulation. In this paper we will describe the field/laboratory and modeling effort undertaken by to quantify the performances of duct systems. Climate zone 2A is included in the field testing, that encompass both hot/dry climates and hot and humid climates. Attics with deeply buried ducts and low air duct leakage had excellent thermal performance and exhibited the best overall thermal performances. Indeed, based on the thermal simulation analysis of deep buried ducts the thermal performance is better than the credit given to them by IRC 2018 code.

Salonvaara, Mikael↗

A novel conditional formulation of the Vlasov–Ampère equations: a conservative, positivity, asymptotic and Gauss law preserving scheme

We propose a novel reformulation of the Vlasov–Ampère equations for plasmas that reveals discrete symmetries that enables simultaneous conservation of mass, momentum and energy; preservation of Gauss’s law; positivity of the distribution function; and consistency with quasi-neutral asymptotics. The approach employs variable and coordinate transformations to yield a coupled system comprising a modified Vlasov equation and associated moment–field equations. The modified Vlasov equation advances a conditional distribution function that excludes mass, momentum and energy densities, which are instead evolved through moment equations enforcing the relevant symmetries, conservation laws and involution constraints. This reformulation aligns naturally with a recent slow-manifold reduction technique, which separates fast electron time scales and simplifies the treatment of the quasi-neutral limit within the reduced moment–field subsystem. Using this framework, we develop a numerical method for the reduced 1D1V subsystem that, for the first time in the literature, satisfies all key physical constraints while maintaining a quasi-neutral asymptotic behaviour. The advantages of the method are demonstrated on canonical electrostatic test problems, including the multiscale ion acoustic shock wave.

1D1V↗

Particle Curtain Experiments at Sandia

Efforts at Sandia National Laboratories have focused on fundamental experiments to understand the dispersal of dense particle distributions in high-speed compressible flow. The experiments are conducted in shock tube facilities where the flow conditions and the initial conditions of the particle distributions are well controlled and well characterized. An additional advantage of the shock tube is that it is more readily able to accommodate advanced measurement diagnostics in comparison to explosive field tests.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SE2-ANL Subchannel Analysis with New Heating Calculation Based on VARIANT Transport Code

The multi-assembly, steady-state sub-channel analysis code SE2-ANL has been widely used in thermal-hydraulic analysis of sodium cooled fast reactors (SFR). It is used to compute the core-wide temperature profiles for a given flow distribution and to determine the assembly flow rates yielding a desired (assembly-wise) peak temperature distribution. It is a modified version of SUPERENERGY-2 coupled with the DIF3D/GAMSOR/DIF3D procedure for coupled neutron and gamma heating calculations, developed at Argonne National Laboratory (ANL). In the current version of SE2-ANL, the neutron and gamma fluxes are calculated using the finite-difference diffusion theory option of the DIF3D code with six triangular meshes per hexagonal assembly. In addition, the pin power distributions are estimated by assuming that the power distribution within an assembly is separable in the radial and axial directions. Furthermore, the total heat generation rate in an assembly duct wall is simply determined by a user-specified fraction of the total gamma heating within that assembly. In order to determine the pin power distributions more accurately by eliminating these limitations, a new coupled heating calculation procedure has recently developed based on the neutron and gamma transport calculations with the VARIANT code. For this, the GAMSOR code has been extended to generate the intra-nodal gamma source distribution to be consistent with the trial functions of VARIANT. Using the computed neutron and gamma flux distributions, the power distribution in each fuel pin and duct wall is determined using a new utility code CURVE, which has been developed to evaluate the VARIANT flux solution at the center of each fuel pin and along the duct mid-wall. For each pin segment contained in a VARIANT node, the axial power shape is represented as a quadratic profile. The SE2-ANL code has been modified to interface with the new, coupled neutron and gamma heating calculation scheme. In addition, several improvements have been made on the computational methods and models of SE2-ANL. To eliminate the limit on the axial mesh size due to the numerical instability problem of the current explicit scheme, the numerical scheme for discretizing the energy equation in the axial direction has been changed to the {theta}-method of the SLTHEN code. An automated flow allocation scheme has also been implemented in SE2-ANL to determine assembly flow rates such that the peak cladding mid-wall temperatures of individual orifice zones are equalized over the burn cycle. The impacts of the new heating calculation scheme based on VARIANT transport solutions on the cladding and fuel temperatures have been tested against the existing method of SE2-ANL. Using the existing and new heating calculation schemes, full-core subchannel analyses were performed for the metal core design of the 1000 MWt Advanced Burner Reactor (ABR-1000), and the resulting coolant, cladding and fuel temperatures were compared. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Pore-scale investigation on stress-dependent characteristics of granular packs and the impact of pore deformation on fluid distribution

Understanding the effect of changing stress conditions on multiphase flow in porous media is of fundamental importance for many subsurface activities including enhanced oil recovery, water drawdown from aquifers, soil confinement, and geologic carbon storage. Geomechanical properties of complex porous systems are dynamically linked to flow conditions, but their feedback relationship is often oversimplified due to the difficulty of representing pore-scale stress deformation and multiphase flow characteristics in high fidelity. In this work, we performed pore-scale experiments of single- and multiphase flow through bead packs at different confining pressure conditions to elucidate compaction-dependent characteristics of granular packs and their impact on fluid flow. A series of drainage and imbibition cycles were conducted on a water-wet, soda-lime glass bead pack under varying confining stress conditions. Simultaneously, X-ray micro-CT was used to visualize and quantify the degree of deformation and fluid distribution corresponding with each stress condition and injection cycle. Micro-CT images were segmented using a gradient-based method to identify fluids (e.g., oil and water), and solid phase redistribution throughout the different experimental stages. Changes in porosity, tortuosity, and specific surface area were quantified as a function of applied confining pressure. Results demonstrate varying degrees of sensitivity of these properties to confining pressure, which suggests that caution must be taken when considering scalability of these properties for practical modeling purposes. Changes in capillary number with confining pressure are attributed to the increase in pore velocity as a result of pore contraction. Furthermore, this increase in pore velocity was found to have a marginal impact on average phase trapping at different confining pressures.

58 GEOSCIENCES↗

Identification and Characterization of Persistent Cold Pool Events from Temperature and Wind Profilers in the Columbia River Basin

Cold pool events occur when deep layers of stable, cold air remain trapped in a valley or basin for multiple days, without mixing out from daytime heating. With large impacts on air quality, freezing events, and es-pecially on wind energy production, they are often poorly forecast by modern mesoscale numerical weather prediction (NWP) models. Understanding the characteristics of cold pools is, therefore, important to provide more accurate forecasts. This study analyzes cold pool characteristics with data collected during the Second Wind Forecast Improvement Project (WFIP2), which took place in the Columbia River basin and Gorge of Oregon and Washington from fall 2015 until spring 2017. A subset of the instrumentation included three microwave radiometer pro?lers, six radar wind pro?lers with radio acoustic sounding systems, and seven sodars, which together provided seven sites with collocated vertical pro?les of temperature, humidity, wind speed, and wind direction. Using these collocated observations, we developed a set of criteria to determine if a cold pool was present based on stability, wind speed, direction, and temporal continuity, and then developed an automated algorithm based on these criteria to identify all cold pool events over the 18 months of the ?eld project. Characteristics of these events are described, including statistics of the wind speed distributions and pro?les, stability conditions, cold pool depths, and descent rates of the cold pool top. The goal of this study is a better understanding of these characteristics and their processes to ultimately lead to improved physical parameterizations in NWP models, and consequently improve forecasts of cold pool events in the study region as well at other locations that experiences similar events.

McCaffrey, Katherine L.↗

Scaling of Floods With Geomorphologic Characteristics and Precipitation Variability Across the Conterminous United States

Abstract Accurate flood risk assessment requires a comprehensive understanding of flood sensitivity to regional drivers and climate factors. This paper presents the scaling of floods (duration, peak, volume) with geomorphologic characteristics of the basin (i.e., drainage area, slope, elevation) and precipitation patterns (rainfall accumulation, variability). Long‐term daily streamflow observations over the 20th and early 21st centuries from Hydro‐Climatic Data Network streamgages across the conterminous United States are used to create a flood event database based on their flood stage information. Antecedent daily rainfall accumulation and variability corresponding to these floods are computed using Global Historical Climatology Network daily data set. Two Bayesian scaling models are developed, and the spatial organization of scaling exponents is investigated. The baseline model quantifies the scaling of floods to geomorphologic characteristics. The dynamic model quantifies the scaling of floods to antecedent precipitation distribution which is further conditioned on geomorphologic characteristics. Results show that small and low‐elevation basins have a stronger response to antecedent rainfall distribution in amplifying flood peaks, while high‐elevation steeper basins have a lower response for flood duration and volume. The dynamic models demonstrate that there are significant variations in the flood scaling rates, with the largest rates up to 40% and 4.5% for flood duration, 64% and 44% for peak, and 98% and 40% for volume found across the Northeast, Coastal Southeast, and Northwest with intensifying rainfall accumulation and variability, respectively. This study advances flood predictions by better informing the flood attributes in the context of dynamical land‐atmosphere perturbations.

54 ENVIRONMENTAL SCIENCES↗

Numerical analysis of flow through oscillating cascade sections

The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and nonzero interblade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the periodic boundary conditions. Approximate nonreflecting boundary conditions were coded for the inlet and exit boundary conditions. Sample unsteady solutions were performed for an oscillating cascade and compared to experimental data. Also, test cases were run for a flat plate cascade to compare with the unsteady, small-perturbation, subsonic analyis. The predictions for oscillating cascades with nonzero interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. The zero degree interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus nonreflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distributions for certain flow conditions.

Huff, Dennis L.↗

Numerical analysis of flow through oscillating cascade sections

The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and nonzero interblade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the perioidic boundary conditions. Approximate nonreflecting boundary conditons were coded for the inlet and exit boundary conditions. Sample unsteady solutions were performed for an oscillating cascade and compared to experimental data. Also, test cases were run for a flat plate cascade to compare with the unsteady, small-perturbation, subsonic analysis. The predictions for oscillating cascades with nonzero interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. The zero degree interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus nonreflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distribution for certain flow conditions.

Huff, Dennis L.↗

Numerical analysis of flow through oscillating cascade sections

The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and non-zero inter-blade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the periodic boundary conditions. Approximate nonreflecting boundary conditions have been coded for the inlet and exit boundary conditions. Sample unsteady solutions have been performed for an oscillating cascade and compared to experimental data. Also, test cases were fun for a flat plate cascade to compare with an unsteady, small-perturbation, subsonic analysis. The predictions for oscillating cascades with non-zero inter-blade phase angles are in good agreement with experimental data and small-perturbation theory. The zero degree inter-blade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus non-reflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distributions for certain flow conditions. This code is expected to be used as a tool for reviewing simpler models that do not include the full nonlinear aerodynamics or as a final check for designs against flutter in turbomachinery.

Huff, Dennis L.↗

Bison Modeling of SiC/SiC Cladding Including Fuel-Pellet Interaction

SiC/SiC composites are considered leading candidates for accident tolerant fuel cladding in light water reactors. However, when SiC-based materials are exposed to neutron irradiation, they experience significant changes in dimensions and physical properties. In the extreme nuclear reactor environment, SiC-based fuel cladding will see neutron damage, significant heat flux, as well as a corrosive environment. In this report, we evaluate the effect of non-uniform changes in the dimensions caused by neutron irradiation with variable temperatures, along with the closing of the fuel – clad gap, on the stress development in the clad over the course of time. To ensure reliable and safe operation of accident tolerant fuel cladding concepts such as SiC-based materials, it is important to assess its thermo-mechanical performance under in-reactor conditions of irradiation and realistic temperature distributions. In this work, we have used the BISON fuel performance modeling code to perform the thermo-mechanical analysis of SiC/SiC cladding. A constitutive model is constructed and solved numerically to predict the stress distribution and variation in the cladding under normal operating conditions. The dependence of dimensions and physical properties variation with irradiation and temperature has been incorporated. These initial results from the analyses provides insights about the stress distribution and variation with time in the cladding as well as the interaction of fuel pellet with the cladding under different conditions of initial fuel pellet-cladding gap and steady state power levels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Analytic solutions for Asay foil trajectories with implications for ejecta source models and mass measurements

We consider the trajectory of an Asay foil ejecta diagnostic for scenarios where ejecta are produced at a singly shocked planar surface and fly ballistically through a perfect vacuum to the sensor. We do so by building upon a previously established mathematical framework derived for the analytic study of stationary sensors. First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source model. The result is an integrodifferential equation of motion for the foil trajectory. This equation yields an easily calculable closed-form implicit solution for the foil trajectory in instant-production scenarios. From there, we derive a boundary condition that particle velocity distributions must satisfy if their associated foil trajectories are to exhibit a smooth initial acceleration, as occurs in some experiments. This condition is identical to one derived previously from a consideration of piezoelectric voltage data obtained in similar experiments. We also compare techniques for inferring accumulated ejecta masses from foil trajectories, first by deriving the exact solution, and then by quantifying the error imposed by a frequently used approximate solution (both subject to the assumption of instantaneous ejecta production). Finally, we examine the common practice of presenting inferred cumulative ejecta masses as a function of implied ejecta velocity, establishing the conditions under which this methodology is most meaningful.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characteristics of Coupled Nongray Radiating Gas Flows with Ablation Product Effects About Blunt Bodies During Planetary Entries

A computational method was developed for the fully-coupled solution of nongray, radiating gas flows with ablation product effects about blunt bodies during planetary entries. The treatment of radiation accounts for molecular band, continuum, and atomic line transitions with a detailed frequency dependence of the absorption coefficient. The ablation of the entry body was solved as part of the solution for a steady-state ablation process. The method was applied by results at typical conditions during entry to Venus. The radiative heating rates along the downstream region of the body can exceed the stagnation point value. The radiative heating to the body is attenuated in the boundary layer at the downstream region of the body and at the stagnation point of the body. A study of the radiation, inviscid flow about spherically capped, conical bodies during planetary entries shows that the nondimensional, radiative heating distributions are nonsimilar with entry conditions. Caution should be exercised in attempting to extrapolate results from known distributions to other entry conditions for which solutions have not yet been obtained.

Sutton, K.↗

Model-centric distribution automation: Capacity, reliability, and efficiency

A series of analyses along with field validations that evaluate efficiency, reliability, and capacity improvements of model-centric distribution automation are presented. With model-centric distribution automation, the same model is used from design to real-time control calculations. A 14-feeder system with 7 substations is considered. The analyses involve hourly time-varying loads and annual load growth factors. Phase balancing and capacitor redesign modifications are used to better prepare the system for distribution automation, where the designs are performed considering time-varying loads. Coordinated control of load tap changing transformers, line regulators, and switched capacitor banks is considered. In evaluating distribution automation versus traditional system design and operation, quasi-steady-state power flow analysis is used. In evaluating distribution automation performance for substation transformer failures, reconfiguration for restoration analysis is performed. In evaluating distribution automation for storm conditions, Monte Carlo simulations coupled with reconfiguration for restoration calculations are used. As a result, the evaluations demonstrate that model-centric distribution automation has positive effects on system efficiency, capacity, and reliability.

25 ENERGY STORAGE↗

A two-and-a-half dimensional symplectic space-charge solver

The nonlinear space-charge effect plays a significant role in high-intensity accelerators and has been extensively studied using multi-particle tracking methods. In this paper, we present a novel 2.5- dimensional symplectic space-charge solver specifically designed for long beam bunches. We begin by detailing its application to a transverse Gaussian density distribution under open boundary conditions in a straight system, where a semi-analytical expression is derived. We then demonstrate the solver’s adaptation to arbitrary distributions in open space, as well as within rectangular and round conducting pipes. Finally, we discuss the extension of this solver to circular accelerator systems. This study shows that the fast 2.5-dimensional solver can be a good approximation to the fully three-dimensional solver for long bunches in large circular accelerators.

Beam code development & simulation techniques↗