Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution parameters”

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,423 records · Page 79

Forward modeling approach to nuclear reaction cross sections: Applications in neutron inelastic scattering

The development of nuclear reaction models for the production of evaluated nuclear data has traditionally been performed by comparing measured cross sections with predictions from reaction model codes whose physical input parameters are adjusted to obtain the best agreement between measured and modeled results. To more directly probe reaction model inputs, this work introduces a forward modeling approach to experimental reaction cross-section determination, where the most important physical input parameters to reaction model calculations are obtained via 𝜒 2 minimization between measured and calculated observables. This was demonstrated using data collected by the Gamma Energy Neutron Energy Spectrometer for Inelastic Scattering (GENESIS) at the 88-inch cyclotron at Lawrence Berkeley National Laboratory, a detection array consisting of organic liquid scintillators and high-purity germanium (HPGe) detectors. Using a broad-spectrum neutron beam and a 99.98%-enriched 56 Fe target, GENESIS was used to perform a simultaneous measurement of 56 Fe 𝛾-ray production cross sections and secondary neutron energy and angle distributions. The results of the forward modeling approach to the determination of energy-differential 𝛾-ray production cross sections for the yrast 4 + → 2 + and 6 + → 4 + transitions, as well as eight other off-yrast transitions, were compared against those obtained using conventional techniques, and the results are in good agreement. In addition to discrete 𝛾-ray yield total scattered neutron energy-angular distributions as a function of incident neutron energy were also obtained using forward modeling and found to agree with evaluated data, with the exception of elastic scattering at small angles. The fitted reaction model parameters obtained through forward modeling were also used to calculate the cross section for the unobserved (𝑛, 2⁢𝑛) reaction; excellent agreement with the current evaluation was obtained, providing a validation of the predictive capabilities of the forward model approach. This work bridges the gap between nuclear data experiment and evaluation by providing a new means for extracting inelastic neutron-scattering cross sections and neutron-induced 𝛾-ray production data while directly probing reaction model physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hardness Prediction by Incorporating Heat Transfer and Molten Pool Fluid Flow in a Multi-pass, Multi-layer Weld for Onsite Repair of Grade 91 steel

Introduction: In the current fleet of fossil-fired power plants, creep strength enhanced ferritic steels (CSEF) are used to sustain the harsh service conditions. Enhanced properties of Grade 91 steel result from tempered martensite with a fine distribution of MX and M23C6 carbides. Grade 91 steel is subjected to onsite welding repair to remedy their degradation due to extreme service condition. Knowledge of weld repairability of these steels, such as as-welded hardness distribution, is essential to establishing sound repair procedures. Experimental trial and error tests can consume a lot of time as many welding variables need to be studied. For numerical modelling, most of the multi-pass multi-layer models are based on finite element method, which are limited to solve the heat conduction equation and ignore convective heat transfer due to melt flow. Moreover, the mesh has to be pre-built based on a known or assumed weld cross-section geometry. These finite element based models thus have limited predictive capability as defects are not considered and nugget size are pre-assumed. This research aims at developing a thermal and microstructure evolution model incorporating molten pool dynamics in a multi-pass multi-layer material deposition to predict the as-welded hardness distribution. Technical Approach: All the thermal, physical, and metallurgical properties of Grade 91 as a function of temperature are collected from the literature and inputted into the thermo-fluid model based on Flow-3D, a computational fluid dynamics software. A multi-pass, multi-layer material deposition is simulated where the melting of filler wire into the molten pool is directly considered based on the volume of fluid (VOF) method. The flow behaviour of the molten pool is used to understand the formation of deposition geometry and defects. The temperature profiles during the multi-pass, multi-layer welding are calculated. The results computed using the new model are compared against the experimental data of fusion zone geometry and thermal cycles. Hardness prediction in the heat-affected zone (HAZ) are made using Johnson-Mehl-Avrami (JMA) equation for solid-state phase transformation kinetics. The JMA parameters are extracted from the experimental data available in the literature. For comparison, a standard finite element heat conduction model is also developed to predict the thermal cycles and hardness distribution in the multi-pass, multi-layer weld. Expected Result: Results obtained using the molten pool dynamic simulation versus the finite element heat conduction model are compared. Specifically, the effects of convective heat transfer on the accuracy of the calculated thermal history, bead shape and size, and HAZ hardness distribution are examined.

Modeling, Grade 91 steel, Multi-pass multi-layer d↗

Improved limits on the tensor-to-scalar ratio using BICEP and P l a n c k data

In this study, we present constraints on the tensor-to-scalar ratio r using a combination of BICEP/Keck 2018 (BK18) and Planck PR4 data allowing us to fit for r consistently with the six parameters of the Λ CDM model. We discuss the sensitivity of constraints on r to uncertainties in the Λ CDM parameters as defined by the Planck data. In particular, we are able to derive a constraint on the reionization optical depth τ and thus propagate its uncertainty into the posterior distribution for r . While Planck sensitivity to r is slightly lower than the current ground-based measurements, the combination of Planck with BK18 and baryon-acoustic-oscillation data yields results consistent with r = 0 and tightens the constraint to r < 0.032 at 95% confidence.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraining the particle-scale diversity of black carbon light absorption using a unified framework

Abstract. Atmospheric black carbon (BC), the strongest absorber of visible solar radiation in the atmosphere, manifests across a wide spectrum of morphologies and compositional heterogeneity. Phenomenologically, the distribution of BC among diverse particles of varied composition gives rise to enhancement of its light absorption capabilities by over twofold in comparison to that of nascent or unmixed homogeneous BC. This situation has challenged the modeling community to consider the full complexity and diversity of BC on a per-particle basis for accurate estimation of its light absorption. The conventionally adopted core–shell approximation, although computationally inexpensive, is inadequate not only in estimating but also capturing absorption trends for ambient BC. Here we develop a unified framework that encompasses the complex diversity in BC morphology and composition using a single metric, the phase shift parameter (ρBC), which quantifies how much phase shift the incoming light waves encounter across a particle compared to that in its absence. We systematically investigate variations in ρBC across the multi-space distribution of BC morphology, mixing state, mass, and composition as reported by field and laboratory observations. We find that ρBC>1 leads to decreased absorption by BC, which explains the weaker absorption enhancements observed in certain regional BC compared to laboratory results of similar mixing state. We formulate universal scaling laws centered on ρBC and provide physics-based insights regarding core–shell approximation overestimating BC light absorption. We conclude by packaging our framework in an open-source Python application to facilitate community-level use in future BC-related research. The package has two main functionalities. The first functionality is for forward problems, wherein experimentally measured BC mixing state and assumed BC morphology are input, and the aerosol absorption properties are output. The second functionality is for inverse problems, wherein experimentally measured BC mixing state and absorption are input, and the morphology of BC is returned. Further, if absorption is measured at multiple wavelengths, the package facilitates the estimation of the imaginary refractive index of coating materials by combining the forward and inverse procedures. Our framework thus provides a computationally inexpensive source for calculation of absorption by BC and can be used to constrain light absorption throughout the atmospheric lifetime of BC.

54 ENVIRONMENTAL SCIENCES↗

Using quasar X-ray and UV flux measurements to constrain cosmological model parameters

ABSTRACT Risaliti and Lusso have compiled X-ray and UV flux measurements of 1598 quasars (QSOs) in the redshift range 0.036 ≤ z ≤ 5.1003, part of which, z ∼ 2.4 − 5.1, is largely cosmologically unprobed. In this paper we use these QSO measurements, alone and in conjunction with baryon acoustic oscillation (BAO) and Hubble parameter [H(z)] measurements, to constrain cosmological parameters in six different cosmological models, each with two different Hubble constant priors. In most of these models, given the larger uncertainties, the QSO cosmological parameter constraints are mostly consistent with those from the BAO + H(z) data. A somewhat significant exception is the non-relativistic matter density parameter Ωm0 where QSO data favour Ωm0 ∼ 0.5 − 0.6 in most models. As a result, in joint analyses of QSO data with H(z) + BAO data the 1D Ωm0 distributions shift slightly towards larger values. A joint analysis of the QSO + BAO + H(z) data is consistent with the current standard model, spatially-flat ΛCDM, but mildly favours closed spatial hypersurfaces and dynamical dark energy. Since the higher Ωm0 values favoured by QSO data appear to be associated with the z ∼ 2 − 5 part of these data, and conflict somewhat with strong indications for Ωm0 ∼ 0.3 from most z < 2.5 data as well as from the cosmic microwave background anisotropy data at z ∼ 1100, in most models, the larger QSO data Ωm0 is possibly more indicative of an issue with the z ∼ 2 − 5 QSO data than of an inadequacy of the standard flat ΛCDM model.

79 ASTRONOMY AND ASTROPHYSICS↗

Dust distribution around low-mass planets on converging orbits

Context : Super-Earths can form at large orbital radii and migrate inward due to tidal interactions with the circumstellar disk. In this scenario, convergent migration may occur and lead to the formation of resonant pairs of planets. Aim : We explore the conditions under which convergent migration and resonance capture take place, and what dynamical consequences can be expected on the dust distribution surrounding the resonant pair. Methods : We combine hydrodynamic planet–disk interaction models with dust evolution calculations to investigate the signatures produced in the dust distribution by a pair of planets in mean-motion resonances. Results : We find that convergent migration takes place when the outer planet is the more massive. However, convergent migration also depends on the local properties of the disk, and divergent migration may result as well. For similar disk parameters, the capture in low degree resonances (e.g., 2:1 or 3:2) is preferred close to the star where the resonance strength can more easily overcome the tidal torques exerted by the gaseous disk. Farther away from the star, convergent migration may result in capture in high degree resonances. The dust distribution shows potentially observable features typically when the planets are trapped in a 2:1 resonance. In other cases, with higher degree resonances (e.g., 5:4 or 6:5) dust features may not be sufficiently pronounced to be easily observable. Conclusions : The degree of resonance established by a pair of super-Earths may be indicative of the location in the disk where capture occurred. There can be significant differences in the dust distribution around a single super-Earth and a pair of super-Earths in resonance.

79 ASTRONOMY AND ASTROPHYSICS↗

System and method for modeling and forecasting cyclical demand systems with dynamic controls and dynamic incentives

Systems and methods for modeling and forecasting cyclical demand systems in the presence of dynamic control or dynamic incentives. A method for modeling a cyclical demand system comprises obtaining historical data on one or more demand measurements over a plurality of demand cycles, obtaining historical data on incentive signals over the plurality of demand cycles, constructing a model using the obtained historical data on the one or more demand measurements and the incentive signals, wherein constructing the model comprises specifying a state-space model, specifying variance parameters in the model, and estimating unknown variance parameters.

24 POWER TRANSMISSION AND DISTRIBUTION↗

torch-einshard v1.0

torch-einshard is a Python library for describing local and distributed PyTorch tensor computations with compact, einsum-like notation. Its expressions name logical axes, specify how they are sharded across a PyTorch DeviceMesh, and represent partial reductions. The library automatically performs contractions, permutations, reshaping, splitting, gathering, reduction, reduce-scatter, and repartitioning while preserving autograd. Additional features include sharding-aware FFTs, tensor rolls, halo exchange, sliding windows, 1D–3D convolutions, uneven-shard handling, parameter initialization and gradient management, and cost-based execution planning. It is designed for scientific machine learning and large-model workloads, including tensor-, sequence-, and spatial-parallel MLPs, attention, convolutions, and spectral operations. Compared with manually combining torch.einsum and distributed collectives, torch-einshard expresses both the mathematical operation and data placement in one readable formula. This reduces boilerplate and synchronization errors, keeps forward and backward communication consistent, and allows the library to select optimized collective strategies without changing model code.

Morozov, Dmitriy [Lawrence Berkeley National Labor↗

SNIa Cosmology Analysis Results from Simulated LSST Images: From Difference Imaging to Constraints on Dark Energy

Abstract The Vera Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to process ∼10 6 transient detections per night. For precision measurements of cosmological parameters and rates, it is critical to understand the detection efficiency, magnitude limits, artifact contamination levels, and biases in the selection and photometry. Here we rigorously test the LSST Difference Image Analysis (DIA) pipeline using simulated images from the Rubin Observatory LSST Dark Energy Science Collaboration Data Challenge (DC2) simulation for the Wide-Fast-Deep survey area. DC2 is the first large-scale (300 deg 2 ) image simulation of a transient survey that includes realistic cadence, variable observing conditions, and CCD image artifacts. We analyze ∼15 deg 2 of DC2 over a 5 yr time span in which artificial point sources from Type Ia supernova (SNIa) light curves have been overlaid onto the images. The magnitude limits per filter are u = 23.66 mag, g = 24.69 mag, r = 24.06 mag, i = 23.45 mag, z = 22.54 mag, and y = 21.62 mag. The artifact contamination levels are ∼90% of all detections, corresponding to ∼1000 artifacts deg –2 in g band, and falling to 300 deg –2 in y band. The photometry has biases <1% for magnitudes 19.5 < m < 23. Our DIA performance on simulated images is similar to that of the Dark Energy Survey difference-imaging pipeline on real images. We also characterize DC2 image properties to produce catalog-level simulations needed for distance bias corrections. We find good agreement between DC2 data and simulations for distributions of signal-to-noise ratio, redshift, and fitted light-curve properties. Applying a realistic SNIa cosmology analysis for redshifts z < 1, we recover the input cosmology parameters to within statistical uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Tomographic optical emission spectroscopy of an atmospheric pressure plasma jet and surface ionization waves on planar and structured surfaces

Here in this paper, an approach for 3D plasma structure diagnostics using tomographic optical emission spectroscopy (Tomo-OES) of a nanosecond pulsed atmospheric pressure plasma jet (APPJ) is presented. In contrast to the well-known Abel inversion, Tomo-OES does not require cylindrical symmetry to recover 3D distributions of plasma light emission. Instead, many 2D angular projections are measured with intensified cameras and the multiplicative algebraic reconstruction technique is used to recover the 3D distribution of light emission. This approach solves the line-of-sight integration problem inherent to optical diagnostics, allowing recovery of localized OES information within the plasma that can be used to better infer plasma parameters within complex plasma structures. Here, Tomo-OES was applied to investigate an APPJ operated with helium in ambient air and impinging on planar and structured dielectric surfaces. Surface charging caused the guided streamer from the APPJ to transition to a surface ionization wave (SIW) that propagated along the surface. The SIW experienced variable geometrical and electrical material properties as it propagated, leading to 3D configurations that were non-symmetric and spatially complex. Light emission from He, N$_2^+$, and N2 were imaged at ten angular projections and the respective time-resolved 3D emission distributions in the plasma were then reconstructed. The spatial resolution of each tomographic reconstruction was 7.4 µm and the temporal resolution was 5 ns, sufficient to observe the guided streamer and the effects of the structured surface on the SIW. Emission from He showed the core of the jet and emission from N$_2^+$ and N2 indicated effects of entrainment of ambient air. Penning ionization of N2 created a ring or outer layer of N$_2^+$ that spatially converged to form the 'plasma bullet' or spatially diverged across a surface as part of a SIW. The SIW entered trenches of size 150 µm, leading to decreases in plasma light emission in regions above the trenches. The plasma light emission was higher in some regions with trenches, possibly due to effects of field enhancement.

plasma interactions with complex surfaces↗

MPC-Based Local Voltage Control Strategy of DGs in Active Distribution Networks

The increasing penetration of distributed generators (DGs) exacerbates the voltage violations in active distribution networks (ADNs). The commonly used centralized voltage control may aggravate the burden of computation and communication. Local voltage control can realize fast response to frequent voltage fluctuations. However, effective local control for DGs still needs to be further investigated. In this paper, a model predictive control (MPC)-based local voltage control strategy of DGs is proposed to cope with voltage violation problems and enhance the adaptation to the DG volatility. DGs are regulated with local Q-V control curves to respond to voltage fluctuations rapidly. Based on the dynamic network partition, a unified control curve for DGs is set in the same area to realize efficient operation. To further improve the performance of voltage control, the parameter tuning of control curves is regularly conducted based on decentralized inter-area coordination with the rolling of the control horizon. Finally, case study is verified on the modified IEEE 33-node distribution system. The proposed local control strategy of DGs can effectively improve the voltage profile and reduce power losses of the whole network without high computation and communication burden.

active distribution network (ADN)↗

Generation of realistic input parameters for simulating atmospheric point-spread functions at astronomical observatories

High-fidelity simulated astronomical images are an important tool in developing and measuring the performance of image-processing algorithms, particularly for high precision measurements of cosmic shear – correlated distortions of images of distant galaxies due to weak gravitational lensing caused by the large-scale mass distribution in the Universe. For unbiased measurements of cosmic shear, all other sources of correlated image distortions must be modeled or removed. One such source is the correlated blurring of images due to optical turbulence in the atmosphere, which dominates the point-spread function (PSF) for ground-based instruments. In this work, we leverage data from weather forecasting models to produce wind speeds and directions, and turbulence parameters, that are realistically correlated with altitude. To study the resulting correlations in the size and shape of the PSF, we generate simulated images of the PSF across a ~10 square-degree field of view – the size of the camera focal plane for the Vera C. Rubin Observatory in Chile – using weather data and historical seeing for a geographic location near the Observatory. We make quantitative predictions for two-point correlation functions (2PCF) that are used in analyses of cosmic shear. We observe a strong anisotropy in the two-dimensional 2PCF, which is expected based on observations in real images, and study the dependence of the orientation of the anisotropy on dominant wind directions near the ground and at higher altitudes.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of Strain and Strain Rate on Dynamic Grain Growth and Subgrain Evolution During Plastic Deformation of an Interstitial-Free Steel at 850 ° C

Here, the effects of strain and strain rate on dynamic grain growth (DGG) and subgrain evolution are reported for an interstitial-free steel deformed at 850 ° C. Microstructures produced during tension tests at true-strain rates of 10 -4 and to 10 -3 s -1 true strains ranging from 0.02 to 0.2 were preserved following deformation. These were characterized using electron backscatter diffraction (EBSD), including the application of spherical harmonic transform indexing to produce high-angular-resolution EBSD (HR-EBSD) data. HR-EBSD data resolved the small misorientation angles of subgrain boundaries while imaging much larger data fields than possible with previously available techniques. The resulting data confirmed that steady-state flow stress is inversely proportional to the average subgrain size and that subgrain boundary misorientation angle increases with strain. The following new observations are reported. The rate of DGG increased with respect to time but decreased with respect to strain as strain rate increased. This behavior is rationalized through a simple model using separate rate parameters for the effects of time and strain. Subgrain size was not constant during steady-state deformation, but decreased slowly with increasing strain. Subgrain size distributions and subgrain boundary misorientation angle distributions were measured, and both remained approximately log-normal during steady-state deformation. Subgrain evolution demonstrated no dependence on parent grain size, crystallographic orientation, or Taylor factor. These new data suggest that steady-state flow stress is more likely controlled by the dislocation density internal to subgrains than by the spacing between subgrain boundaries.

dynamic grain growth↗

Investigation of 3D printed lightweight hybrid composites via theoretical modeling and machine learning

Hybrid composites combine two or more different fillers to achieve multifunctional or advanced material properties, such as lightweight and enhanced mechanical properties. The properties of the composites significantly depend on their microstructures, which can be tailored via advanced 3D printing processes. Understanding the process-structure-property relationships is critical to enable the design and engineering of novel hybrid composites for applications in aerospace, automotive, and protective coatings. Here, for this work, we develop 3D printable and lightweight hybrid composites and leverage the conventional design of experiments, a theoretical hybrid model, and an image-driven machine learning (ML) method to investigate their mechanical behaviors. The hybrid composites are formulated with elastomer matrix, microfillers, and thin-shell particles, enabling a significant degree of design freedom of microstructures with densities and mechanical properties varying up to 70% and 91%, respectively. Our statistical analysis indicates that the 3D printing path direction and the microfibers fraction are dominating process parameters with contribution percentages of 45.3% and 57.7% on the specific stiffness and strength, respectively. A hybrid mechanics model is developed based on a simple Weibull distribution function and classical single-filler models to effectively capture the variations in mechanical properties, however, it overestimates the values due to its statistical constraints and idealization of experimental uncertainty. The image-driven ML model leverages the microscale images directly without losing the structural details, shows more accurate predictions with experimental data, and has 48.6% lower root mean square error than the theoretical model.

3D printing↗

Thermal loading analysis of the ring injection dump for the Spallation Neutron Source facility

The ring injection dump (RID) is the largest beam dump in the Spallation Neutron Source facility at Oak Ridge National Laboratory and accepts a fraction of the beam from the ring that is not captured in the ring during injection. Thermo-fluid modeling of the RID components was performed during the initial design by using the software ANSYS CFX to characterize the power rating for the RID. From the calculations, the power rating was lowered from 200 kW to 150 kW, at 1.3 GeV beam energy due to concerns regarding the heating of the concrete structure. The Proton Power Upgrade project scope included reevaluating the power limit to the RID in anticipation that increasing it would provide more operational flexibility. The focus was on the shielding because it was the limiting factor, and data collected from thermocouples installed on the structure were used to benchmark the analysis. A three-stage process was adopted in the latest study. First, the already existing steady-state calculations were validated and extended with the latest Monte Carlo N-Particle source term calculations. Second, transient calculations were conducted to capture the dynamic state of the system in response to the energy dumped to the RID during the following operational period of the beam from 2005–2019 (14 years), and sensitivity analyses of crucial parameters were performed to benchmark the model with reasonable accuracy. Finally, an idealized transient cycle analysis was conducted with realistic duty factors to predict the temperature distribution at higher beam powers. Overall, the results confirmed the accuracy of the original steady state model but allowed the development of a new validated transient model with higher accuracy for future analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Computational fluid dynamics analysis of char conversion in Sandia’s pressurized entrained flow reactor

Design and analysis of practical reactors utilizing solid feedstocks rely on reaction rate parameters that are typically generated in lab-scale reactors. Evaluation of the reaction rate information often relies on assumptions of uniform temperature, velocity, and species distributions in the reactor, in lieu of detailed measurements that provide local information. This assumption might be a source of substantial error, since reactor designs can impose significant inhomogeneities, leading to data misinterpretation. Spatially resolved reactor simulations help understand the key processes within the reactor and support the identification of severe variations of temperature, velocity, and species distributions. In this work, Sandia’s pressurized entrained flow reactor is modeled to identify inhomogeneities in the reaction zone. Tracer particles are tracked through the reactor to estimate the residence times and burnout ratio of introduced coal char particles in gasifying environments. The results reveal a complex mixing environment for the cool gas and particles entering the reactor along the centerline and the main high-speed hot gas reactor flow. Furthermore, the computational fluid dynamics (CFD) results show that flow asymmetries are introduced through the use of a horizontal gas pre-heating section that connects to the vertical reactor tube. Computed particle temperatures and residence times in the reactor differ substantially from the idealized plug flow conditions typically evoked in interpreting experimental measurements. Furthermore, experimental measurements and CFD analysis of heat flow through porous refractory insulation suggest that for the investigated conditions (1350 °C, <20 atm), the thermal conductivity of the insulation does not increase substantially with increasing pressure.

47 OTHER INSTRUMENTATION↗

Room temperature buffer gas beam of metastable state titanium atoms

We produce beams of neutral titanium (Ti) atoms in their metastable 3d3(F4)4sa5F5 state by laser ablation into He, N2, and Ar buffer gases. The high temperatures associated with the ablation process populate the a5F5 level without the need for optical pumping. Remarkably, we observe that Ti atoms in the a5F5 state survive thousands of collisions with He and Ar buffer gas atoms without being quenched to lower-energy states. We study the yield of Ti atoms when ablated into buffer gases of varying species and pressure, quantify quenching rates and diffusion cross sections based on simple models, and provide insight into optimal design parameters for an ablation cell. Using a 3.3 cm ablation cell with interchangeable exit apertures, we produce metastable atom beams and quantify their brilliance and velocity distributions as functions of buffer gas pressure.

Schrott, Jackson↗