Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sub scale”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Improving initialization and evolution accuracy of cosmological neutrino simulations

Neutrino mass constraints are a primary focus of current and future large-scale structure (LSS) surveys. Non-linear LSS models rely heavily on cosmological simulations — the impact of massive neutrinos should therefore be included in these simulations in a realistic, computationally tractable, and controlled manner. Here, a recent proposal to reduce the related computational cost employs a symmetric neutrino momentum sampling strategy in the initial conditions. We implement a modified version of this strategy into the Hardware/Hybrid Accelerated Cosmology Code (HACC) and perform convergence tests on its internal parameters. We illustrate that this method can impart Script $\mathcal{O}$(1%) numerical artifacts on the total matter field on small scales, similar to previous findings, and present a method to remove these artifacts using Fourier-space filtering of the neutrino density field. Moreover, we show that the converged neutrino power spectrum does not follow linear theory predictions on relatively large scales at early times at the 15% level, prompting a more careful study of systematics in particle-based neutrino simulations. We also present an improved method for backscaling linear transfer functions for initial conditions in massive neutrino cosmologies that is based on achieving the same relative neutrino growth as computed with Boltzmann solvers. Our self-consistent backscaling method yields sub-percent accuracy in the total matter growth function. Comparisons for the non-linear power spectrum with the Mira-Titan emulator at a neutrino mass of m ν = 0.15 eV are in very good agreement with the expected level of errors in the emulator and in the direct N-body simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

UMass 2-Body WEC Techno-Economic Assessment

The University of Massachusetts (UMass) is developing a 2-body wave energy converter (WEC) device that is converting mechanical power into electricity using a mechanical motion rectifier that allows the system to couple to a flywheel. UMass has completed numerical modeling, wave tank testing, and PTO sub-system testing and needed assistance in developing a techno-economic model to enable optimization of their topology, comparison to a generic heaving point absorber topology, and guide the next steps in their development efforts. The core objective was to develop a techno-economic approach and modeling tool that allows benchmarking of the two topologies across a wide range of scales to evaluate their respective competitiveness in different application spaces. This data includes the final report as well as a supporting spreadsheet containing the data produced for this report.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Pore-Selective Fullerene Loading in a Porphyrin-Based Metal–Organic Framework Controls Photoinduced Charge-Transfer Dynamics

Building porous donor−acceptor networks based on host−guest interactions in metal organic frameworks (MOFs) provides unique opportunities for tuning charge separation in highly tailorable materials. Here we focus on installing electron-rich porphyrins and electron-deficient fullerene derivatives in the PCN-222 MOF using a solvent assisted ligand insertion (SALI) method. The fullerene is primarily bound in the large pore, where it is subject to distinct dielectric environments through dimethylformamide (DMF) and 1,4-dioxane solvent exchange. Following photoexcitation, sub-picosecond charge transfer involving initial exciplex population is observed, with different charge recombination pathways and lifetimes depending on solvent polarity through modulation of charge-transfer state energies. While the 1,4-dioxane environment yields charge recombination within 1 ns via local fullerene and porphyrin triplet state population, DMF results in charge recombination directly to the ground state on much longer time scales, including some lifetime components in the microsecond range. Fullerene loading influences these kinetics, and the potential for charge delocalization due to fullerene aggregation within the pores is evaluated by using molecular dynamics simulations.

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

Multiphysics simulation of microscale copper printing by confined electrodeposition using a nozzle array

3D printing of metals at the microscale and nanoscale is crucial to produce high-density interconnects and intricate structures in electronic devices. Conventional melting and sintering processes are not suitable for these scales due to a reliance on individual metal particles in the size range of tens of micrometers. Confined electrodeposition (CED) is an established alternative to conventional metal 3D printing processes in which an electrolyte is used to selectively induce deposition of the metal on the printing surface. However, commercialization and efficiency of this process have been limited due to a reliance on sub-micrometer nozzles to achieve desirable deposition rates and single nozzle to achieve uniformity of printed structures. Here, we address these challenges by computationally analyzing an array of microscale nozzles. We tailor the convection within the electrolyte to alter both deposition rate and geometric uniformity of the printed structures. The results show that for large nozzles the evaporation alone is not sufficient to obtain high deposition rates, yet an external pressure can be used to increase deposition and alter uniformity (thickness) of printed structures. Our results can be used to design and analyze new experiments toward parallel multi-nozzle deposition using CED toward high-throughput metal printing.

3D printing↗
DOE OSTI2023

Comparative Thermal Performance of Downdraft and Updraft Forced Convection in a Representative 3x3 PWR Fuel Assembly

A plethora of theoretical, numerical, and experimental investigations have relied on updraft forced convection of operating fuel assemblies for light water-cooled fuel assemblies. The inertial scales of turbulence in updraft flow through PWR fuel assemblies are primarily influenced by upstream mixing patterns within the lower plenum and consequent acceleration of flow through the lower core and support plates into the bottom nozzles. In secondary heat transfer applications such as steam generation, recuperation and economization, downdraft flow is utilized to retrieve maximal sensible heat from the primary coolant. Sub-channel volumes aligned vertically could also benefit from force of gravity in developing turbulent flows. A large reduction in pumping power, associated cost of operation and maintenance, and improved balance of plant is suggested to be possible through downdraft forced convection in PWR-type reactors. This article explores the impact of downdraft forced convection within the well-studied 3x3 sub-channel within a range of 20 hydraulic diameters across a representative intermediate flow mixing grid. A range of inflow Reynolds number from 10,000 to 60,000 is simulated over a range of equivalent heat loads of 20% to 100% spanning 20 hydraulic diameters. The fuel rods are represented in true geometric detail, and a constant-value heat profile is assumed for the investigation. At each inflow condition, the turbulence intensity, peak velocity, local heat transfer coefficient, peak wall heat flux, and peak wall temperature are collected across the sub-channel and compared with an updraft configuration evaluated at the same conditions.

Rao, Vivek↗
DOE OSTI2021

Adaptive Mesh Refinement for Parallel in Time Methods

The project applied the multigrid-reduction-in-time (MGRIT) algorithm to an existing sub-cycled adaptive mesh refinement (AMR) code to investigate the performance of flows dominated by inertial physics. Previous work demonstrated good performance from MGRIT+AMR applied to flows dominated by diffusive physics. Consistent with previous experience, inertial physics negatively affected convergence rates and performance. Efforts to circumvent this issue by appealing to the physics of turbulence were investigated. It has been demonstrated that scales can be effectively transferred between multigrid levels for a turbulent flow resulting in a) partial convergence observed and b) nearly identical results to sequential time-stepping. Performance improvements have not yet been demonstrated - attempts at coarsening the grid on coarser MG levels compromises the solution quality and leads to divergence. This report summarize the accomplishments for the time-frame from 10/5/2020 to 12/31/2020 with an informal no-cost-extension to 05/20/2021

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

IC Report for project “Developing nonlinear laser-plasma instability models for high-fidelity, multi-physics simulation capability for ICF/HED” [Slides]

We’ve achieved our technical goals: (1) Perform large-scale kinetic simulations of SRS and CBET a range of plasma and laser conditions; (2) Develop a new $\delta f$-Gaussian-mixture algorithm to represent non-Maxwellian distribution functions from particle trapping and time-dependent plasma response; (3) Develop physics-based preliminary nonlinear LPI models for coupling to LRT/ Mazinisin and AMP; and (4) Begin implementation of LPI sub-grid models into Mazinisin.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The EGS Collab Project: Status and Accomplishments

The EGS Collab project, supported by the US Department of Energy, is addressing challenges in implementing enhanced geothermal systems (EGS). This includes improving understanding of the stimulation of crystalline rock to create appropriate flow pathways, and the ability to effectively simulate both the stimulation and the flow and transport processes in the resulting fracture network. The project is performing intensively monitored rock stimulation and flow tests at the 10-m scale in an underground research laboratory. Data and observations from the field test are compared to simulations to understand processes and to build confidence in numerical modeling of the processes. In Experiment 1, we examined hydraulic fracturing an underground test bed at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, at a depth of approximately 1.5 km. We drilled eight sub-horizontal boreholes in a well-characterized phyllite. Six of the boreholes were instrumented with many sensor types to allow careful monitoring of stimulation events and flow tests, and the other two boreholes were used for water injection and production. We performed a number of stimulations and flow tests in the testbed. Our monitoring systems allowed detailed observations and collection of numerous data sets of processes occurring during stimulation and during dynamic flow tests. Long-term ambient temperature and chilled water flow tests were performed in addition to many tracer tests to examine system behavior. Data were rapidly analyzed, allowing adaptive control of the tests. Numerical simulation was used to answer key experimental design questions, to forecast fracture propagation trajectories and extents, and to analyze and evaluate results. Many simulations were performed in near-real-time in conjunction with the field experiments, with more detailed process study simulations performed on a longer timeframe. Experiment 2 will examine hydraulic shearing in a test bed being built at the SURF at a depth of about 1.25 km in amphibolite under a different set of stress and fracture conditions than Experiment 1. Five sets of fracture orientations were considered in design, and three orientations seem to be consistently observed.

EGS Collab↗
DOE OSTI2026

Evolution of size-selected Pt cluster catalysts on prototypical oxide supports

The current quest for new pathways into sustainable, efficient and durable energy conversion technologies makes the used for a fundamental understanding of the atomic-scale phenomena underlying catalytic processes ever more pressing. In this context, characterizing catalyst particles in situ provides valuable information about the evolution of their composition, structure, oxidation state and charge state during an ongoing process. To disentangle the influence of individual parameters – temperature, pressure, gas composition, cluster size, as well as support acidity, redox state and defect density – it is crucial to control them precisely and separately in experiments. At the example of size-selected Pt n clusters – i.e. sub-nm particles defined to the exact number of atoms – on flat oxide supports, we follow their rich evolution phenomena via (synchrotron-based) X-ray photoelectron spectroscopy (XPS) and scanning tunnelling microscopy (STM) during temperature ramps and in various gas environments. Here, we present our experience with these highly defined, yet complicated-to-create samples in ultra-high vacuum (UHV) and at mbar pressures. We discuss their stability during transport to synchrotrons and under reaction conditions on three prototypical oxide supports and show the various phenomena that can be disentangled. Ptn clusters on the non-reducible silicon dioxide, SiO 2 , remain size-selected and show remarkable stability. Doping strongly influences the size-dependent binding energy shifts and changes the Pt response to oxidative and reaction conditions, which we attribute to different cluster geometries: p-type doping leads to wetting, enhanced sinter resistance and a diminished response to oxidative environments, compared to clusters on n-type doped samples. We compare the system with our previous findings of a similar change in dimensionality for Pt 20 clusters on the reducible ceria, CeO 2 (111), support, induced by modulation of its O vacancy density. The Pt n /CeO 2 system is particularly interesting for strategies to stabilize and redisperse Pt clusters dynamically. Finally, we study the evolution of Pt n clusters on another reducible oxide support, magnetite, Fe 3 O 4 (001), in 0.1 mbar alternating redox conditions at RT and elevated temperature. In analogy to findings previously reported for Pt/TiO 2 (110), the clusters either become encapsulated by a thin oxide film via strong metal–support interaction (SMSI) or deeply buried in the magnetite. Overall, our approach of following the evolution of size-selected clusters on oxide supports leads to fundamental atomic-scale insights on nano-scale catalyst materials, on our path to sustainable, dynamic and self-repairing catalysts.

Falling, Lorenz J. [Technical Univ. of Munich (Ger↗
DOE OSTI2022

Do chaotic field lines cause fast reconnection in coronal loops?

Over the past decade, Boozer has argued that three-dimensional (3D) magnetic reconnection fundamentally differs from two-dimensional reconnection due to the fact that the separation between any pair of neighboring field lines almost always increases exponentially over distance in a 3D magnetic field. According to Boozer, this feature makes 3D field-line mapping chaotic and exponentially sensitive to small non-ideal effects; consequently, 3D reconnection can occur without intense current sheets. We test Boozer's theory via ideal and resistive reduced magnetohydrodynamic simulations of the Boozer–Elder coronal loop model driven by sub-Alfvénic footpoint motions [A. H. Boozer and T. Elder, Phys. Plasmas 28, 062303 (2021)]. In this work, our simulation results significantly differ from their predictions. The ideal simulation shows that Boozer and Elder under-predict the intensity of current density due to missing terms in their reduced model equations. Furthermore, resistive simulations of varying Lundquist numbers show that the maximal current density scales linearly rather than logarithmically with the Lundquist number.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Toward UV models of kinetic mixing and portal matter. VII. A light dark photon in the 3 𝑐 ⁢3 𝐿 ⁢1 𝐴 ⁢1 𝐵 model

The kinetic mixing (KM) portal, by which the Standard Model (SM) photon mixes with a light dark photon arising from a new 𝑈⁢(1) 𝐷 gauge group, allows for the possibility of viable scenarios of sub-GeV thermal dark matter with appropriately suppressed couplings to the SM. This KM can only occur if particles having both SM and dark quantum numbers, here termed portal matter, also exist. The presence of such types of states and the strong suggestion of a need to embed 𝑈⁢(1) 𝐷 into a non-Abelian gauge structure not too far above the TeV scale based on the renormalization-group equation running of the 𝑈⁢(1)𝐷 gauge coupling is potentially indicative of an enlarged group linking together the visible and dark sectors. The gauge group 𝐺 = 𝑆⁢𝑈⁢(3) 𝑐 × 𝑆⁢𝑈⁢(3) 𝐿 × 𝑈⁢(1) 𝐴 × 𝑈⁢(1) 𝐵 = 3 𝑐 ⁢3 𝐿 ⁢1 𝐴 ⁢1 𝐵 is perhaps the simplest setup wherein the SM and dark interactions are partially unified in a non-Abelian fashion that is not a simple product group of the form 𝐺 = 𝐺 SM × 𝐺 𝐷 encountered frequently in earlier work. The present paper describes the implications and phenomenology of this type of setup.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

First Constraints on the Epoch of Reionization Using the Non-Gaussianity of the Kinematic Sunyaev-Zel’dovich Effect from the South Pole Telescope and Herschel-SPIRE Observations

We report results from an analysis aimed at detecting the trispectrum of the kinematic Sunyaev-Zel’dovich (kSZ) effect by combining data from the South Pole Telescope (SPT) and Herschel-SPIRE experiments over a 100 deg 2 field. The SPT observations combine data from the previous and current surveys, namely SPTpol and SPT-3G, to achieve depths of 4.5, 3, and 16 μ K - arcmin in bands centered at 95, 150, and 220 GHz. For SPIRE, we include data from the 600 and 857 GHz bands. We reconstruct the velocity-induced large-scale correlation of the small-scale kSZ signal with a quadratic estimator that uses two cosmic microwave background (CMB) temperature maps, constructed by optimally combining data from all the frequency bands. We reject the null hypothesis of a zero trispectrum at 10.3 σ level. However, the measured trispectrum contains contributions from both the kSZ and other undesired components, such as CMB lensing and astrophysical foregrounds, with kSZ being sub-dominant. We use the agora simulations to estimate the expected signal from CMB lensing and astrophysical foregrounds. After accounting for the contributions from CMB lensing and foreground signals, we do not detect an excess kSZ-only trispectrum and use this nondetection to set constraints on reionization. By applying a prior based on observations of the Gunn-Peterson trough, we obtain an upper limit on the duration of reionization of Δ z re , 50 < 4.5 (95% confidence level). We find these constraints are fairly robust to foregrounds assumptions. This trispectrum measurement is independent of, but consistent with, Planck’s optical depth measurement. This result is the first constraint on the epoch of reionization using the non-Gaussian nature of the kSZ signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

On the Orbital Elements of the Two-body Problem with Slowly Decreasing Mass: The Gyldén–Mestchersky Cases

We study the dynamics of two-body problems with gravitational parameters μ = μ(t) decreasing with the time. In particular we focus our attention on the behavior of the orbital elements in the first two Gyldén–Mestchersky cases: in the first case (GM-I) μ(t)=μ{sub 0}(1+αt){sup −1}, where μ{sub 0} is the constant initial value of the parameter and α a small positive constant. In the second case (GM-II) μ(t)=μ{sub 0}(1+2αt){sup −1/2}. In the GM-I problem, using a suitable transformation of the radius vector and the physical time, which reduces the original problem to a pure Kepler problem, we are able to obtain explicit analytical expressions of the orbital elements that hold true for all t⩾0 and therefore derive their asymptotic behavior. Moreover using the alternative scaled time these expressions show clearly their periodic and secular terms. For the GM-II problem a similar transformation to the GM-I problem leads to a perturbed Kepler problem with a radial direction perturbation. Although exact solutions in terms of elliptic functions may be obtained, the use of some integrals leads to analytical expressions of orbital elements that allow us to study their behavior. Finally, some comparisons between the time evolution of elements in both problems are also included.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Modeling and Optimizing Pumped Storage in a Multi-stage Large Scale Electricity Market under Portfolio Evolution

To leverage the fast-ramping capability of resources to provide great value to the grid, electricity system operators such as the Midcontinent Independent System Operator (MISO) continue to evolve their approaches for integrating energy storage resources, including pumpedstorage hydro (PSH), into the electricity markets. However, new challenges arise in modeling and optimizing these energy-limited resources across multiple market clearing processes and planning studies with uncertainties and imperfect information. For instance, current market practices of PSH owners specifying pumping/generating hours can result in sub-optimal generation dispatch. Letting grid operators optimize PSH with the consideration of multiple operating modes and energy limitation constraints can potentially bring economic benefits to both the system and the PSH owners. However, in multi-stage clearing process of electricity markets, utilizing the PSH flexibility to deal with realized uncertainties can cause deviation in the multi-stage scheduling processes. The resulting financial risks from the schedule deviation may not be acceptable to PSH owners. In addition, to effectively utilize this energy limited resource, the state of charge (SOC) constraints of PSH needs to be continuously optimized and the marginal cost of deviation need to reflect the expected cost to purchase or sell energy at future times to compensate for deviations. This project aims to develop a prototype enhanced PSH model and improved price signals in the multi-stage market clearing process with proper consideration of the unique characteristics of PSH, in order to better align underlying PSH capabilities with evolving grid needs, particularly including the needs for more frequent and larger cycling to manage variability and uncertainty from renewables.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

A Goldstone boson equivalence for inflation

The effective field theory of single-field inflation characterizes the inflationary epoch in terms of a pattern of symmetry breaking. An operator acquires a time-dependent vacuum expectation value, defining a preferred spatial slicing. In the absence of dynamical gravity, the fluctuations around the time-dependent background are described by the Goldstone boson associated with this symmetry breaking process. With gravity, the Goldstone is eaten by the metric, becoming the scalar metric fluctuation. In this paper, we will show that in general single-field inflation, the statistics of scalar metric fluctuations are given by the statistics of this Goldstone boson decoupled from gravity up to corrections that are controlled as an expansion in slow-roll parameters. This even holds in the presence of additional parameters, like the speed of sound, that naively enhance the impact of the gravitational terms. In the process, we derive expressions for leading and sub-leading gravitational corrections to all-orders in the Goldstone boson.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Energy-peak based method to measure top quark mass via B -hadron decay lengths

We develop a method for the determination of the top quark mass using the distribution of the decay length of the B -hadrons originating from its decay. This technique is based on our earlier observation regarding the location of the peak of the b quark energy distribution. Such “energy-peak” methods enjoy a greater degree of model-independence with respect to the kinematics of top quark production compared to earlier proposals. The CMS experiment has implemented the energy-peak method using associated b -jet energy as an approximation for b quark energy. The new method uses B -hadron decay lengths, which are related to b quark energies by convolution. The advantage of the new decay length method is that it can be applied in a way that evades jet-energy scale (JES) uncertainties. Indeed, CMS has measured the top quark mass using B -hadron decay lengths, but they did not incorporate the energy-peak method. Therefore, mismodeling of top quark transverse momentum remains a large uncertainty in their result. We demonstrate that, using energy-peak methods, this systematic uncertainty can become negligible. We show that with the current LHC data set, a sub-GeV statistical uncertainty on the top quark mass can be attained with this method. To achieve a comparable systematic uncertainty as is true for many methods based on exclusive or semi-inclusive observables using hadrons, we find that the quark-hadron transition needs to be described significantly better than is the case with current fragmentation functions and hadronization models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Simulation of dynamic crystal plasticity with a Lagrangian discontinuous Galerkin hydrodynamic method

Here we present a new Lagrangian modal discontinuous Galerkin (DG) hydrodynamic method that supports a dynamic dislocation based crystal plasticity model for simulating the mechanical behavior of crystallographic materials, both single crystal and polycrystalline, under dynamic conditions. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated by Taylor series polynomials of varying degree. These polynomials describe macro-scale hydrodynamic behavior while their evolution is determined by evaluating the dynamic crystal plasticity model at material points within the element. The dynamic crystal plasticity model is sensitive to the time increment size, with too large of time increments leading to instability in the model. To mitigate this, the temporal evolution of the dynamic crystal plasticity model is achieved with the combination of a sub-incrementing scheme with Heun’s third-order time integration scheme, which is also used to temporally evolve the governing equations. The implementation of the dynamic crystal plasticity model within the DG framework is tested using a 2D approximation of the Taylor impact test with a single crystal material, using quadratic elements that have faces that can bend. In addition to the standard continuous material modeling, we propose a new simulation method that would represent the heterogeneous behavior of polycrystalline microstructures within an element by varying the position and material properties of the material points within the element. This method is demonstrated using random orientation distributions on materials points that are arranged in both structured and random configurations.

42 ENGINEERING↗
DOE OSTI2022

Using Ensemble Data Assimilation to Estimate Transient Hydrologic Exchange Flow Under Highly Dynamic Flow Conditions

Abstract Quantifying dynamic hydrologic exchange flows (HEFs) within river corridors that experience high‐frequency flow variations caused by dam regulations is important for understanding the biogeochemical processes at the river water and groundwater interfaces. Heat has been widely used as a tracer to infer steady‐state flow velocities through analytical solutions of heat transport defined by the diurnal temperature signals. Under sub‐daily dynamic flow conditions, however, such analytical solutions are not applicable due to the violation of their fundamental assumptions. In this study, we developed a data assimilation‐based approach to estimate the sub‐daily flux under highly dynamic flow conditions using multi‐depth temperature observations at a 5‐min resolution. If the hydraulic gradient is measured, Darcy's law was used to calculate the flux with permeability estimated from temperature responses below the riverbed. Otherwise, flux was estimated directly by assimilating multi‐depth temperature data at 1‐ or 2‐hr time intervals assuming one‐dimensional flow and heat transport governing equation. By comparing estimated fluxes with model‐generated synthetic truth, we demonstrated that both schemes have robust performance in estimating fluxes under highly dynamic flow conditions. This data assimilation‐based flux estimation method was able to capture the vertical sub‐daily fluxes using multi‐depth high‐resolution temperature data alone, even in the presence of multi‐dimensional flow. This approach has been successfully applied to real field temperature data collected at the Hanford site, which experiences highly dynamic HEFs. Our study shows the promise of adopting distributed 1‐D temperature monitoring to capture spatial and temporal exchange dynamics in river corridors at a watershed scale or beyond.

54 ENVIRONMENTAL SCIENCES↗