Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density distribution”

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,207 records · Page 67

Electron Cloud Simulations for the Electron-Ion Collider in Brookhaven National Laboratory

For high-intensity circular accelerators and storage rings, if the secondary electron emission yield (SEY) of the vacuum chamber surfaces is high, an EC (Electron Cloud) could build up with the passage of the circulating beam. The presence of EC can strongly affect the beam quality, such as transverse instabilities, transverse emittance growth, and beam loss. Furthermore, the heat load from EC can exceed the available cryogenic capability. For the RHIC superconducting (SC) arc magnets, to be used for the hadron storage ring of the Electron-Ion Collider (EIC), the dynamic heat load budget is 0.5 W/m to the 4.5 K stainless steel beam pipe. This can limit the maximum beam bunches or intensity of the EIC, hence diminishing the luminosity provided by the EIC. The EC has affected the beam instability and significantly contributed to cryogenic heat load in the Large Hadron Collider (LHC). Positron storage rings for which ECs have been a critical factor in the design and performance include KEKB in Japan and EC buildup remains one of the concerns for future high-intensity accelerators design. EC considerations have driven the SuperKEKB collider design and the positron damping ring for the proposed International Linear Collider (ILC). The LHC luminosity upgrade is contingent on reducing the bunch spacing to 25 ns; at this bunch spacing, severe EC buildup has been observed. The success of the upgrade is likely contingent on limiting EC buildup. To study the EC heat load, we did some EC simulations with PyECLOUD code for the dipole, quadrupole, sextupole magnets, and the warm (drift) section of the EIC hadron storage ring. PyECLOUD is an EC simulation code developed by CERN. The code has been validated and used to study EC in the LHC, SPS, and PS. To eliminate EC buildup, the sources of electrons must be minimized. First, we will reduce the production of electrons due to residual gas ionization by specifying the maximum gas density, which requires a vacuum chamber with low electron-stimulated desorption (ESD) yields and a sufficient (preferably distributed) pumping speed. Second, we should reduce the secondary electrons with a lower secondary electron yield (SEY) material.

43 PARTICLE ACCELERATORS↗

X-Ray Scattering Data Processing and Analysis

Semicrystalline polymers have a degree of crystallinity that depends on chain structure and crystallization conditions. Thermomechanical properties are influenced by the degree of crystallinity and crystalline morphology. A common method used to understand crystalline content, crystal and amorphous layer thickness, and layer thickness distributions, is X-ray scattering (XS). This report uses XS to first investigate crystalline orientation in high density polyethylene (HDPE) samples. Then, a single orientation is analyzed via XS in a temperature study to investigate crystalline behavior under changing temperatures. Lastly, block copolymer filaments extruded at varying speeds are analyzed. The X-ray data are used to gain a baseline understanding of crystalline behavior in these unique filaments.

36 MATERIALS SCIENCE↗

Thermal Conductivity of AM Fabricated Stainless Steel 316L Part

Thermal properties are an integral part of many diverse engineering applications, and additive manufacturing (AM), particularly laser powder bed fusion (LPBF) is shown to affect thermal properties due to the laser processing parameters. For 316L stainless steel, there is little prior research to determine the effects of the process on thermal properties. In this work, the temperature gradient is shown to create uniform, chess board-like distributions of grains in the build direction. These zones create dislocations which were visualized with EBSD techniques. Processing parameters cause hierarchal grain size distribution, with localized concentrations of small grains and large grains. Thermomechanical stresses in the rapid solidification increases dislocation density during grain formation. Previous research shows a higher density of dislocations decreases local thermal conductivity. Local and bulk thermal conductivity are shown in this work to have statistically lowered values to an average of 10-12 W/m-K compared to 14 W/m-K for conventional 316L.

36 MATERIALS SCIENCE↗

ATOMIC Simulations and Experimental Data for Basalt-like Compounds

This data consists of simulations and experimental measurements of laser-induced breakdown spectroscopy (LIBS). The simulations are produced by ATOMIC, a general purpose plasma modeling and kinetics code that has been designed to compute emission (or absorption) spectra from plasmas [2]. The makeup of the plasma was considered to be divided into some proportion water, some proportion Martian atmosphere (CO2), and some proportion target (from the rock or object impacted by the laser), where these proportions add to 1. Based on expert knowledge, the proportion of water was kept in the range [0.0,0.5] and the proportion of atmosphere was kept in the range [0.02, 0.9]. Our overall suite of simulations contains six sets of simulations that differ in which elements were considered to make up the target. Within each set, we used uniformly drawn temperatures and log mass densities within pre-specified ranges. The temperature range was [0.5,1.5] eV and the log (base 10) mass density range was [-7,-4]. The proportion of water, atmosphere, and target were drawn from a symmetric Dirichlet distribution, but draws in which the propor- tion of water or atmosphere exceeded the pre-specified limits were rejected from the design. Up to eleven constituent elements (Si, Al, Fe, Mg, Ca, O, Ti, Mn, Na, K, P) were considered for the target, as they are the most common elements found in basalt compounds and were used in [1]. For each run, the proportions of the constituent elements making up the target were drawn from a symmetric Dirichlet distribution. We ran 1,350 simulations that included nonzero proportions of all eleven elements. We also ran simulations which excluded some of these elements. In particular, we ran 1,000 simulations that only included nonzero proportions for the six most common elements (Si, Al, Fe, Mg, Ca, O). We also ran five sets, each with 500 simulations, that only included nonzero proportions for five of the six most common elements (but where all sets included O). Thus, we generated a total of 4,850 spectra representing basalt-like compounds in which the target was comprised of oxygen and between four and ten other elements. The ATOMIC code produced spectra over a range of 240nm - 880nm that roughly mimics the range collected by the ChemCam instrument on the Mars rover Curiosity. Each spectra had 32,000 wavelengths split across three spectrometer ranges (to mimic ChemCam). The experimental data, described in [1], measures a prepared basalt sample. All files are kept in directories whose names indicate the set of elements considered for the target with file names numbered to indicate the line in the design files used to produce the simulation. The designs are provided as text files with names indicating their purpose. The experimental data is provided as a CSV file which contains a header with measurement information, followed by a collection of 50 shots across a collection of wavelengths, along with the computed median and mean across shots.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Environments of Green Pea Galaxies. I. The KISS Sample

We present environmental analyses for 13 KPNO International Spectroscopic Survey Green Pea (GP) galaxies. These galaxies were discovered via their strong [O ιιι] emission in the redshift range 0.29 < z < 0.42, and they are undergoing a major burst of star formation. A primary goal of this study is to understand what role the environment plays in driving the current star formation activity. By studying the environments around these extreme star-forming galaxies, we can learn more about what triggers their star formation processes and how they fit into the narrative of galaxy evolution. Using the Hydra multifiber spectrograph on the WIYN 3.5 m telescope, we mapped out the galaxy distribution around each of the GPs (out to ~15 Mpc at the redshifts of the targets). Using three density analysis methodologies chosen for their compatibility with the geometry of our redshift survey, we categorized the galaxian densities of the GPs into different density regimes. We find that the GPs in our sample tend to be located in low-density environments. We find no correlation between the density and the SFRs seen in the GPs. We conclude that the environments the GPs are found in are likely not the driving factor behind their extreme star formation activity.

79 ASTRONOMY AND ASTROPHYSICS↗

Verification of an energy-conserving semi-implicit electrostatic particle-in-cell scheme for modeling high-density plasma at scale

A verification study of a semi-implicit energy-conserving electrostatic particle-in-cell algorithm is presented. The algorithm relaxes the time-step and mesh-size constraints that require resolution of the plasma period and Debye length associated with traditional explicit momentum-conserving particle-in-cell algorithms. Physical implications and applicability of using the semi-implicit scheme for modeling high-density plasmas are discussed. Where possible, numerical results are compared against analytical solutions. The simulation results indicate that the algorithm is stable at time steps larger than twice the inverse plasma frequency and cell sizes larger than the Debye length. It is found that the algorithm gives adequate results, provided that the distribution function and the spatiotemporal scales dictating the physics of the problem are resolved. As such, the algorithm may provide a robust method for kinetic modeling of high-density plasmas at scale.

Cyclotron resonance↗

Understanding the Effect of Electrochemical Properties and Microstructure on the Microgalvanic Corrosion of Mg Alloys via Phase-Field Simulations

In this work, we apply a phase-field model to gain insights into the corrosion behavior of Mg alloys. Specifically, we study the effect of electrochemical properties of the materials and their environment and the spatial distribution of second phases in a magnesium alloy. To this end, we perform sensitivity analyses in which we separately vary the exchange current density and corrosion potential of the second phase, as well as the electrolyte conductivity. We describe the effects of these parameters on the corrosion current and the resulting corrosion-front morphology. We also confirm that the ratio between the Wagner length, calculated for the anodic phase, and the initial width of the anodic domain is a good predictor of the morphological characteristic of the corroding surface. In addition, we examine the effect of the microstructure on the corrosion behavior in both 2D and 3D systems by varying the shape, size, and distribution of the second-phase regions. Our findings provide insights into the influence of the two-phase microstructure on the corrosion behavior in a magnesium alloy.

Electrochemistry↗

The Sizes of z ∼ 9−10 Galaxies Identified in the Brightest of Reionizing Galaxies (BoRG) Survey

Redshift z = 9–10 object selection is the effective limit of Hubble Space Telescope (HST) imaging capability, even when confirmed with Spitzer. If only a few photometry data points are available, it becomes attractive to add criteria based on their morphology in these J- and H-band images. One could do so through visual inspection, a size criterion, or alternate morphometrics. We explore a vetted sample of Brightest of Reionizing Galaxies (BoRG) z ∼ 9 and z ∼ 10 candidate galaxies and the object rejected by Morishita+ to explore the utility of a size criterion in z = 9–10 candidate selection. A stringent, point-spread function (PSF)-corrected effective radius criterion ([[Formula]]) would result in the rejection of 65%–70% of the interlopers visually rejected by Morishita et al. It may also remove up to ∼20% of bona fide brightest (L≫L{sup ∗}) z = 9 or 10 candidates from a BoRG selected sample based on the Mason et al. luminosity functions, assuming the Holwerda et al. z∼9 size–luminosity relation. We argue that including a size constraint in lieu of a visual inspection may serve in wide-field searches for these objects in, e.g., Euclid or HST archival imaging with the understanding that some brightest (L≫L{sup ∗}) candidates may be missed. The sizes of the candidates found by Morishita et al. follow the expected size distribution of z ∼ 9 for bright galaxies, consistent with the log normal in Shibuya et al. and single objects. Two candidates show high star formation surface density (Σ{sub SFR}>25M{sub ⊙} kpc{sup −2}) and all merit further investigation and follow-up observations.

79 ASTRONOMY AND ASTROPHYSICS↗

A Medium Voltage Three-Stage Power Converter Topology for Distribution Grid Scale Energy Storage Systems

This work proposes a three-stage converter topology with a medium frequency isolation transformer for direct integration of energy storage systems into medium voltage distribution grids. The distributed architecture of the topology, using standard AC/DC converters, has been developed with the aim of plug-and-play capability and voltage and power scalability. The medium frequency transformer facilitates improved power density, while its 20:1 voltage translation capability enables the use of limited voltage level (500 Vdc) of the commercially available batteries. This paper presents the topology structure and details how it can be scaled to reach the full range of the distribution grid voltages and the required power levels. A transformer design for 500 Hz rated frequency is discussed. The control scheme of the converter equipped with voltage support, frequency support, and black start capabilities is also presented. Control hardware-in-the-loop results using OPAL-RT, with converter switching models connected to a weak IEEE-4 node test feeder modified for the aforementioned grid services, are provided to show the topology feasibility.

Chinthavali, Madhu Sudhan↗

Parameter estimation with maximal updated densities

A recently developed measure-theoretic framework solves a stochastic inverse problem (SIP) for models where uncertainties in model output data are predominantly due to aleatoric (i.e., irreducible) uncertainties in model inputs (i.e., parameters). The subsequent inferential target is a distribution on parameters. Another type of inverse problem is to quantify uncertainties in estimates of “true” parameter values under the assumption that such uncertainties should be reduced as more data are incorporated into the problem, i.e., the uncertainty is considered epistemic. A major contribution of this work is the formulation and solution of such a parameter identification problem (PIP) within the measure-theoretic framework developed for the SIP. The approach is novel in that it utilizes a solution to a stochastic forward problem (SFP) to update an initial density only in the parameter directions informed by the model output data. In other words, this method performs “selective regularization” only in the parameter directions not informed by data. Further, the solution is defined by a maximal updated density (MUD) point where the updated density defines the measure-theoretic solution to the PIP. Another significant contribution of this work is the full theory of existence and uniqueness of MUD points for linear maps with Gaussian distributions. Data-constructed Quantity of Interest (QoI) maps are also presented and analyzed for solving the PIP within this measure-theoretic framework as a means of reducing uncertainties in the MUD estimate. We conclude with a demonstration of the general applicability of the method on two problems involving either spatial or temporal data for estimating uncertain model parameters. The first problem utilizes spatial data from a stationary partial differential equation to produce a MUD estimate of an uncertain boundary condition. The second problem utilizes temporal data obtained from the state-of-the-art ADvanced CIRCulation (ADCIRC) model to obtain a MUD estimate of uncertain wind drag coefficients for a simulated extreme weather event near the Shinnecock Inlet located in the Outer Barrier of Long Island, NY, USA.

42 ENGINEERING↗

The gravitational lensing imprints of DES Y3 superstructures on the CMB: a matched filtering approach

Low-density cosmic voids gravitationally lens the cosmic microwave background (CMB), leaving a negative imprint on the CMB convergence |$\kappa$|⁠. This effect provides insight into the distribution of matter within voids, and can also be used to study the growth of structure. We measure this lensing imprint by cross-correlating the Planck CMB lensing convergence map with voids identified in the Dark Energy Survey Year 3 (DES Y3) data set, covering approximately 4200 deg|$^2$| of the sky. We use two distinct void-finding algorithms: a 2D void-finder that operates on the projected galaxy density field in thin redshift shells, and a new code, Voxel, which operates on the full 3D map of galaxy positions. We employ an optimal matched filtering method for cross-correlation, using the Marenostrum Institut de Ciències de l’Espai N-body simulation both to establish the template for the matched filter and to calibrate detection significances. Using the DES Y3 photometric luminous red galaxy sample, we measure |$A_\kappa$|⁠, the amplitude of the observed lensing signal relative to the simulation template, obtaining |$A_\kappa = 1.03 \pm 0.22$| (⁠|$4.6\sigma$| significance) for Voxel and |$A_\kappa = 1.02 \pm 0.17$| (⁠|$5.9\sigma$| significance) for 2D voids, both consistent with Lambda cold dark matter expectations. We additionally invert the 2D void-finding process to identify superclusters in the projected density field, for which we measure |$A_\kappa = 0.87 \pm 0.15$| (⁠|$5.9\sigma$| significance). The leading source of noise in our measurements is Planck noise, implying that data from the Atacama Cosmology Telescope, South Pole Telescope and CMB-S4 will increase sensitivity and allow for more precise measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

Aggressively Dissipative Dark Dwarfs: The Effects of Atomic Dark Matter on the Inner Densities of Isolated Dwarf Galaxies

We present the first suite of cosmological hydrodynamical zoom-in simulations of isolated dwarf galaxies for a dark sector that consists of cold dark matter and a strongly dissipative subcomponent. The simulations are implemented in GIZMO and include standard baryons following the FIRE-2 galaxy formation physics model. The dissipative dark matter is modeled as atomic dark matter (aDM), which forms a dark hydrogen gas that cools in direct analogy to the Standard Model. Our suite includes seven different simulations of ∼10 10 M ⊙ systems that vary over the aDM microphysics and the dwarf’s evolutionary history. We identify a region of aDM parameter space where the cooling rate is aggressive and the resulting halo density profile is universal. In this regime, the aDM gas cools rapidly at high redshifts, and only a small fraction survives in the form of a central dark gas disk; the majority collapses centrally into collisionless dark “clumps,” which are clusters of subresolution dark compact objects. These dark clumps rapidly equilibrate in the inner galaxy, resulting in an approximately isothermal distribution that can be modeled with a simple fitting function. Even when only a small fraction (∼5%) of the total dark matter is strongly dissipative, the central densities of classical dwarf galaxies can be enhanced by over an order of magnitude, providing a sharp prediction for observations.

cold dark matter↗

Critical Current Distributions of Recent Bi-2212 Round Wires

Bi 2 Sr 2 CaCu 2 O 8+x (Bi-2212) is the only high-field, high-temperature superconductor (HTS) capable of reaching a critical current density J c (16 T, 4.2 K) of 6500 A·mm –2 in the highly desirable round wire (RW) form. However, state-of-the-art Bi-2212 conductors still have a critical current density (J c ) to depairing current density (J d ) ratio around 20 to 30 times lower than that of state-of-the-art Nb–Ti or REBCO. Previously, we have shown that recent improvements in Bi-2212 RW J c are due to improved connectivity associated with optimization of the heat treatment process, and most recently due to a transition to a finer and more uniform powder manufactured by Engi-Mat. One quantitative measure of connectivity may be the critical current (I c ) distribution, since the local I c in a wire can vary along the length due to variable vortex-microstructure interactions and to factors such as filament shape variations, grain-to-grain connectivity variations and blocking secondary phase distributions. Modeling the experimental V-I transition measured on a low resistance shunt as a complex sum of voltage contributions of individual filament and wire sub-sections allows a numerical extraction of the I c distribution from the d 2 V/dI 2 treatment of the V-I curves. Here we compare ~0.1 m length I c distributions of Bi-2212 RWs with recent state-of-the-art very high-J c Engi-Mat powder and lower J c and older Nexans granulate powder. We do find that the I c spread for Bi-2212 wires is about twice the relative standard of high-J c Nb–Ti well below H irr . We do not yet see any obvious contribution of the Bi-2212 anisotropy to the I c distribution and are rather encouraged that these Bi-2212 round wires show relative I c distributions not too far from high-J c Nb–Ti wires.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

(U) A Linear Response Model Predicts Reactivity From a Density Profile

We tested the ability to predict the system reactivity, described by alpha, given a density profile using a simple linear system response. We generated a suite of 1-dimensional density profiles that consisted of nominal density, a discontinuity, and a decay. These profiles were prescribed a functional form and the mass was conserved in all cases. From these density profiles, we calculated the alpha value of the 3-dimensional system.We calculated a linear response function given a training set of the 1-dimensional density profiles, and the system reactivity described by alpha. We tested the robustness of the response function using the remaining test data. Our results showed very good agreement between the predicted and calculated test values, where the distribution of alpha differences was centered about zero and had a standard deviation of 0.005 gens/shake. The predicted and calculated alpha values did not significantly differ (t=-0.0009 p<0.99). We used Singular Value Decomposition (SVD) to reduce the matrix rank by retaining95% of the cumulative singular value contributions. This reduced the matrix rank by 91.7%. We generated the linear response matrix and calculated the difference between the predicted and calculated alpha values. Using the reduced order matrix, we showed good agreement between the predicted and calculated alpha values where the distribution of differences was centered near zero, the standard deviation was 0.006 gens/shake, and the statistical t-test showed good agreement (t=0.02, p<0.98). These results show a linear relationship between a series of 1-dimensional density profiles,where the mass was conserved, and the system reactivity. The next steps of this work will be to investigate the linear response using 2-dimensional density profiles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Experimental approach for development of a powder spreading metric in additive manufacturing

The powder spreading is a vital step of powder-based additive manufacturing (AM) processes. The quality of spread powder can considerably influence the properties of fabricated parts. Poorly packed powder beds with high surface roughness result in printed part layers with large porosity and low dimensional accuracy, leading to poor mechanical properties. Therefore, the powder spreadability and its dependence on process parameters and powder characteristics should be quantified to improve the efficiency of powder-based AM methods. This study proposes a novel dimensionless powder spreadability metric that can be commonly used in different powder-based AM processes. The quality of spread powder in terms of powder bed density and surface roughness was evaluated by adjusting the process parameters including recoating velocity and layer thickness, and powder characteristics including particle size distribution. In addition, the dynamic repose angle was proposed and examined as another powder spreadability metric. Furthermore, the results showed that these two proposed metrics were strongly correlated and lower recoating velocity and larger layer thickness led to higher spreadability and lower dynamic repose angle.

36 MATERIALS SCIENCE↗

Population balance modeling of polyurethane foam formation with pressure‐dependent growth kernel

Abstract Polyurethane foams are widely used materials often chosen for their useful characteristics such as low thermal conductivity, ease of application, and high strength‐to‐weight ratios. Computational models are needed to predict the dynamics of the flow and expansion, and the resulting material properties, to improve manufacturing processes. In this paper, a model for PMDI, a water‐blown polyurethane foam, is presented. By extending a kinetics‐based approach by adding bubble‐scale information via a population balance equation (PBE) using the quadrature method of moments, we can track bubble size distributions during foaming. We present results from a three‐dimensional computational fluid dynamics model using arbitrary Lagrangian–Eulerian interface tracking implemented in finite element software. The model compares favorably with experimental data, including dynamics, bubble distributions measured by both camera and diffusion wave spectroscopy, and post‐test bubble size from scanning electron microscopy and density measurements from x‐ray computed tomography.

Ortiz, Weston↗

Automatic and rapid calibration of urban building energy models by learning from energy performance database

Urban building energy modeling (UBEM) is attracting increasing attention in the energy modeling filed. Unlike modeling a single building using detailed building systems information, UBEM generally uses limited high-level building stock data to infer default assumptions about building characteristics and operations. Additionally, this practice inherently brings uncertainty to UBEM. This study introduced a novel method of automatic and rapid calibration of UBEM based on the annual electricity and natural gas energy use data by learning the correlations between crucial model input parameters and the building energy use from the reference building models. A case study was presented to calibrate 72 large office buildings built before 1978 in San Francisco. Seventeen model parameters were selected and Monte Carlo sampling was used to create 1000 samples that reasonably represent the parameter space. Then 1000 simulations were performed for the reference building model to create an energy performance database. The results showed that by learning from the energy performance database, it took less than four simulation runs on average to calibrate a building model. After the calibration, the distributions of each parameter were obtained to replace their single predefined default values. For example, the default lighting power density of 21.39 W/m 2 was calibrated to be 7.50 W/m 2 on average. The case study successfully demonstrated the effectiveness of the novel calibration method for UBEM in the mild climate. The method will be further tested in future for other climate zones and other building types.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Non-premixed counterflow methane flames in DC/AC/NS electric fields

The response of counterflow diffusion flames to sub-breakdown DC and AC electric fields, as well as their superposition with ns pulse discharge waveforms, is studied in the plane-to-plane electrode geometry. Sub-breakdown DC and low-frequency AC electric fields cause the flame displacement toward the grounded electrode, in the direction of the applied field, indicating that the body force on the positive ions exceeds that on the electrons and negative ions. As the AC frequency increases, the flame response becomes less pronounced, due to the reduction of the electrohydrodynamic (EHD) body force impulse over the AC half-period. The electric field in the electrode gap is determined by ps Electric Field Induced Second Harmonic (E-FISH) generation, with absolute calibration using sub-breakdown ns pulses overlapped with the measured electric field waveform. The results show that the electric field distribution across the flame in the current saturation regime follows the Laplacian field. This indicates that the space charge density in the gap is too low to distort the applied DC or AC field, consistent with the kinetic modeling predictions. Combining a nanosecond pulse discharge with a sub-breakdown DC field generates a diffuse plasma across the entire gap. Time-resolved and spatially resolved measurements of the electric field in the discharge indicate the ionization wave propagation between the electrodes. The present results do not exhibit a detectable flame displacement enhancement by ns discharge pulses combined with a sub-breakdown field, observed previously. Kinetic modeling calculations show that the absence of this effect in the plane-to-plane geometry is due to the rapid plasma self-shielding. This indicates that alternative electrode geometries limiting the self-shielding would be more effective for the plasma / electric field enhanced flameholding and flame stabilization applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗