Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic Modeling”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Uniform modelling of the stellar density of thirteen tidal streams within the Galactic halo

ABSTRACT We present the results of fitting a flexible stellar stream density model to a collection of thirteen streams around the Milky Way, using photometric data from DES, DECaLS, and Pan-STARRS. We construct density maps for each stream and characterize their tracks on the sky, width, and distance modulus curves along the length of each stream. We use these measurements to compute lengths and total luminosities of streams and identify substructures. Several streams show prominent substructures, such as stream broadening, gaps, large deviations of stream tracks, and sharp changes in stream densities. Examining the group of streams as a population, as expected we find that streams with globular cluster progenitors are typically narrower than those with dwarf galaxy progenitors, with streams around 100 pc wide showing overlap between the two populations. We also note the average luminosity of globular cluster streams is significantly lower than the typical luminosity of intact globular clusters. The likely explanation is that observed globular cluster streams preferentially come from lower luminosity and lower density clusters. The stream measurements done in a uniform manner presented here will be helpful for more detailed stream studies such as identifying candidate stream members for spectroscopic follow up and stellar stream dynamical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Molecular Structure of Adsorbed Water Phases in Silica Nanopores

The adsorption of water vapor in silica nanopores with different pore morphologies and surface hydrophilicities was studied to quantify the densities and thicknesses of the water sorption layers and deduce their molecular structures. Water adsorption to surface hydroxyls is described by a multilayer sorption model. At low pressure, the water adsorption isotherms are largely independent of pore size and the adsorbed amounts scale with the surface hydroxyl density. Adsorbed-phase densities corresponding to the adsorption of two water molecules per surface hydroxyl group are found in the first adsorbed water layer for a wide range of surface hydroxyl densities. The densities and layer thickness values found in narrow pores indicate that patchy adsorbed layers form if not enough water molecules exist for a full layer, which coexist with dry pore surface regions. This behavior indicates cooperative adsorption effects, i.e., a preference for the formation of hydrogen bonds between water molecules bound to surface hydroxyls. In narrow pores, pore condensation limits further growth of the sorption layer, and in larger pores and at the planar quartz surface, a second adsorbed water layer is formed, which can hold up to approximately 4 additional water molecules per surface hydroxyl group. The water molecules in these thicker adsorption layers arrange such that the sorption layer density is similar to the bulk water density. Pore confinement limitations on the sorption layers are observed in pores with radii of as large as 8 nm. Molecular dynamics modeling reveals two preferential orientations for water molecules adsorbing to surface hydroxyl groups and suggests an intralayer structuring in the adsorbed monolayer. Furthermore, adsorbed water molecules in the sorption layer are bonded to the surface hydroxyl group via the donation and acceptance of hydrogen bonds.

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

Supporting Bioproducts Industry Growth with a System Dynamics Decision-Support Tool

A bio-based economy requires chemical products as well as fuels to be produced from biomass. Although a variety of universities, government agencies, start-ups and established firms have engaged in bioproduct development, many projects have failed to reach the point of commercialization and commercialized bioproducts struggle to capture and maintain market share. To date, there has been no general research into the factors that contribute to bioproduct failure or success. This work presents the Bioproduct Transition Dynamics (BTD) system dynamics model, a decision-support tool that simulates the bioproduct development process from pre-piloting research through construction of the first commercial-scale plant, as well as the processes of obtaining funding from investors and government agencies. The core of the BTD is a feedback loop between bioproduct developers and funders, which relates development progress measured with indicators such as net present value to funders’ decisions to continue investing. External factors such as feedstock prices, market size and growth, and the existence of bioproduct consumers also influence a development project’s chances of receiving follow-on funding. Virtually any bioproduct can be represented with the BTD: direct replacements, performance-advantaged products, niche and commodity markets can all be modeled. The goal of the BTD project is to inform decisions made by developers in both established firms and start-ups, investors, government agencies, and other stakeholders interested in growing the nascent U.S. bioproducts industry. This talk will cover the general structure and functionality of the BTD, and present results from an analysis performed with the BTD to demonstrate its use as a decision-support tool and the insights it can provide. The goal of the analysis is to identify the most critical factors that lead to direct replacement and performance-advantaged bioproduct projects emerging successfully from the “Valley of Death”, and to determine if these factors differ between the two bioproduct types.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Self-Assembled Solute Networks in Crowded Electrolyte Solutions and Nanoconfinement of Charged Redoxmer Molecules

Redoxmers are electrochemically active organic molecules storing charge and energy in electrolyte fluids circulating through redox flow batteries (RFBs). Such molecules typically have solvent-repelling cores and solvent-attracting pendant groups introduced to increase solubility in liquid electrolytes. These two features can facilitate nanoscale aggregation of the redoxmer molecules in crowded solutions. In some cases, this aggregation leads to the emergence of continuous networks of solute molecules in contact, and the solution becomes microscopically heterogeneous. In this study, we used small-angle X-ray scattering (SAXS) and molecular dynamics modeling to demonstrate formation of such networks and examine structural factors controlling this self-assembly. We also show that salt ions become excluded from these solute aggregates into small pockets of electrolytes, where these ions strongly associate. This confinement by exclusion is also likely to occur to charged redoxmer molecules in a “sea” of neutral precursors coexisting in the same solution. Here, we demonstrate that the decay lifetime of the confined charged molecules in such solutions can increase several fold compared to dilute solutions. We attribute this behavior to a “microreactor effect” on reverse reactions of the confined species during their decomposition.

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

INITIAL EXPLORATION OF A NOVEL TRANSIENT ARREST SYSTEM INVOLVING FUEL HEATING

A preliminary analysis on a novel accident response system to diminish the severity of supercritical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compare to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 \$ to 1.3 \$, it was found that the an optimal system response could reduce peak fuel temperatures during the transient by 3.5\% to 5\%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

INITIAL EXPLORATION OF A NOVEL TRANSIENT ARREST SYSTEM INVOLVING FUEL HEATING (Presentation)

A preliminary analysis on a novel accident response system to diminish the severity of supercritical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compare to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 \$ to 1.3 \$, it was found that the an optimal system response could reduce peak fuel temperatures during the transient by 3.5% to 5%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Effect of lot microstructure variations on detonation performance of the triaminotrinitrobenzene (TATB)-Based insensitive high explosive PBX 9502

PBX 9502 is an important insensitive high explosive due to its combination of safety properties and detonation performance. It is a polymer-bonded formulation consisting of 95 wt.% 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as the high explosive crystal, bound with Kel F-800 (FK-800), a co-polymer of chlorotrifluoroethylene and vinylidene-fluoride. Two different types are used, one known as virgin that uses only pristine manufactured TATB, and a second known as recycled that has 50 wt.% of its TATB reclaimed from machining scraps of previously pressed virgin PBX 9502. Recycled lots have a higher percentage of fine particles compared to virgin lots, due to the fracturing and damage sustained by TATB crystals during pressing. We examine the influence that this TATB microstructure difference between virgin and recycled PBX 9502 has on detonation performance properties. New rate-stick geometry diameter effect, detonation front shape and cylinder expansion test data are obtained for two previously uncharacterized virgin lots and one recycled lot of PBX 9502. This is combined with previously published data for one virgin and one recycled lot for evaluation purposes. Detonation shock dynamics model calibrations are conducted on each of the five lots to provide an assessment of the detonation timing characteristics of virgin versus recycled PBX 9502 lots. For two of the virgin lots, detonations propagate slower in the rate-stick geometry than those in the two recycled lots. However, the other virgin lot tested has propagation rates comparable to that of the recycled lots for larger diameter rate-sticks. The latter result, though, is shown to be geometry dependent and depends on the range of detonation curvatures accessed in different geometries. New copper-confined, cylinder expansion tests are conducted on each of the three virgin and two recycled lots to obtain detonation product Jones-Wilkins-Lee equations of state, enabling an assessment of the metal push capabilities for each of the five lots. We find that the metal push capabilities, characterized by the evolution of the heat of detonation with volume, are similar between the virgin and recycled lots. Thus, changes in microstructure between different PBX 9502 lots seemingly affect the rate of reaction, but not the overall energy content.

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

Metallization of colloidal crystals

Colloidal crystals formed by size-asymmetric binary particles coassemble into a wide variety of colloidal compounds with lattices akin to ionic crystals. Recently, a transition from a compound phase with a sublattice of small particles to a metal-like phase in which the small particles are delocalized has been predicted computationally and observed experimentally. In this colloidal metallic phase, the small particles roam the crystal maintaining the integrity of the lattice of large particles, as electrons do in metals. A similar transition also occurs in superionic crystals, termed sublattice melting. Here, we use energetic principles and a generalized molecular dynamics model of a binary system of functionalized nanoparticles to analyze the transition to sublattice delocalization in different coassembled crystal phases as a function of temperature (T), number of grafted chains on the small particles, and number ratio between the small and large particles n s : n l . Here we find that n s : n l is the primary determinant of crystal type due to energetic interactions and interstitial site filling, while the number of grafted chains per small particle determines the stability of these crystals. We observe first-order sublattice delocalization transitions as T increases, in which the host lattice transforms from low- to high-symmetry crystal structures, including A20 → bct → bcc, A d → bct → bcc, and bcc → bcc/fcc → fcc transitions and lattices. Analogous sublattice transitions driven primarily by lattice vibrations have been seen in some atomic materials exhibiting an insulator-metal transition also referred to as metallization. We also find minima in the lattice vibrations and diffusion coefficient of small particles as a function of n s : n l , indicating enhanced stability of certain crystal structures for n s : n l values that form compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Control methods for mitigating flow oscillations in a supercritical CO 2 recompression closed Brayton cycle

A dynamic model of a 10 MWe supercritical CO 2 (sCO 2 ) recompression closed Brayton power generation cycle is used to investigate control methods for mitigating oscillations in flow conditions and power demand load, especially as the load ramp rate is increased and the load setpoint is more closely tracked. The focus is on control of the main compressor inlet temperature (MCIT) since oscillations are prompted by the strong nonlinearities in sCO 2 fluid properties near the critical point and by changes in system inventory for efficient load management. In this study, the oscillatory behavior during load ramping is first shown and then two potential control solutions for substantially reducing oscillations are given: 1) the use of the main compressor inlet guide vanes (IGV) and 2) bypassing a portion of the sCO 2 around the cooler to maintain the MCIT. These control methods reduce the effect of cycle feedback from interactions between the MCIT and load. When implementing either control method, a load setpoint ramp of 7.5%/min was maintainable while achieving closer load setpoint tracking with limited oscillations. Finally, the impact of sensor noise on the control was examined since derivative action was used in the IGV and cooler sCO 2 bypass control. Finally, while some effect was evident, signal noise was not problematic to the control.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

The effect of local chemical ordering on Frank-Read source activation in a refractory multi-principal element alloy

In this work, we investigate the operation of Frank-Read (FR) sources in a refractory multi-principal element alloy (MPEA). Simulations of discrete dislocation motion in MPEAs is enabled by the development of a phase field dislocation dynamics model that treats the atomic-scale fluctuations in lattice energies across the glide plane, present due to local ordering in the chemical composition within the nominally random MPEA atomic structure. We consider, through simulation, a range of length scales over which ordering occurs, varying from short-range lengths, a few times dislocation core width, to long-range lengths, an order of magnitude longer than the core. Characteristic of this body-centered cubic MPEA, the simulations also include screw/edge character dependence in glide resistance, as informed by atomic scale simulation. The critical stresses to activate the source for the same source size are found to be statistically distributed, as a direct consequence of the underlying variation in lattice energy. Furthermore, analysis of critical state for activating edge and screw FR sources in the MPEA reveals that FR source operation occurs via a two-step mechanism, involving athermal kink-pair formation, unlike the conventional FR source operation in a material with no composition fluctuations. This mechanism lowers the average critical stress required to activate the FR source and causes the statistical dispersion in critical stress to depend on the range of composition ordering. More importantly, it leads to a more severe dependence of source strength on FR source length than predicted by line tension alone.

36 MATERIALS SCIENCE↗
DOE OSTI2025

EMPDF : inferring the Milky Way mass with data-driven distribution function in phase space

We introduce the emPDF (empirical distribution function), a novel dynamical modelling method that infers the gravitational potential from kinematic tracers with optimal statistical efficiency under the minimal assumption of steady state. emPDF determines the best-fitting potential by maximizing the similarity between instantaneous kinematics and the time-averaged phase-space distribution function (DF), which is empirically constructed from observation upon the theoretical foundation of oPDF (Han et al. 2016). This approach eliminates the need for presumed functional forms of DFs or orbit libraries required by conventional DF- or orbit-based methods. emPDF stands out for its flexibility, efficiency, and capability in handling observational effects, making it preferable to the popular Jeans equation or other minimal assumption methods, especially for the Milky Way (MW) outer halo where tracers often have limited sample size and poor data quality. We apply emPDF to infer the MW mass profile using Gaia DR3 data of satellite galaxies and globular clusters, obtaining enclosed masses of M (,r) = 26±8, 46±8, 90±13⁠, and 149±40 x 10 10 M ⊙ at r = 30, 50, 100⁠, and 200 kpc, respectively. These are consistent with the updated constraints from simulation-informed DF fitting (Li et al. 2020). While the simulation-informed DF offers superior precision owing to the additional information extracted from simulations, emPDF is independent of such supplementary knowledge and applicable to general tracer populations. emPDF is currently implemented for tracers with complete 6D kinematics within spherical potentials, but it can potentially be extended to address more general problems.

Astrophysics of Galaxies (astro-ph.GA)↗
DOE OSTI2023

Enhancing Organic Semiconductor Molecular Packing Using Perovskite Interfaces to Improve Singlet Fission

Singlet fission, a process by which one singlet exciton is converted into two lower energy triplet excitons, is sensitive to the degree of electronic coupling within a molecular packing structure. Variations in molecular packing can be detrimental to triplet formation and triplet–triplet separation, ultimately affecting the harvesting of triplets for electricity in organic photovoltaic devices. Here, six phase-pure molecular packing structures of 6,13-bis(triisopropylsilylethynyl)pentacene (TIPS-pentacene) with varying optoelectronic properties are isolated using 2D lead halide perovskites as tunable, crystalline surfaces for crystallization. Transient absorption spectroscopy reveals that while triplet formation is fast (<100 fs) regardless of template structure, the increased ordering in perovskite-templated samples speeds up triplet–triplet separation and recombination, providing evidence that the benefits of phase-purity offset minor variations in molecular packing. Molecular dynamics modeling of the interface reveals that perovskite-templating allows for closer packing of TIPS-pentacene molecules for all perovskite templates. With an extensive number of organic molecule-perovskite pairings, this work provides a methodology to use ordered, periodic surfaces to elucidate structure–property relationships of small organic molecules in order to adjust structural or optoelectronic responses, such as molecular packing and singlet fission.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase segregation mechanisms of small molecule-polymer blends unraveled by varying polymer chain architecture

As phase separation between the small-molecule semiconductor and the polymer binder is the key enabler of blend-based organic field-effect transistors (OFETs) fabricated by low-cost solution processing, it is crucial to understand the underlying phase separation mechanisms that determine the phase morphology, which significantly impacts device performance. Beyond the parameter space investigated in previous work, here we investigate the formation of blends by varying the branch architecture of the polymer binder and by shortening the solvent dry time using ultrasonic spray casting. The phase morphologies of the resulting blend films have been thoroughly characterized with a variety of techniques in three dimensions over multiple length scales, including AFM, energy-filtered transmission electron microscope, and neutron reflectivity, and have been correlated with electrical transport performance. From the results, we have inferred that the phase morphology is kinetically determined, limited by the inherent slow movement of polymer macromolecules. The kinetic picture, supported by molecular dynamics modeling, not only consistently explains our observations but also resolves inconsistencies in previous works. The achieved mechanistic understanding will guide further optimization of blend-based organic electronics, such as OFETs and organic photovoltaics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Topological or rotational non-Abelian gauge fields from Einstein-Skyrme holography

We report analytically known states at non-zero temperature which may serve as a powerful tool to reveal common topological and thermodynamic properties of systems ranging from the QCD phase diagram to topological phase transitions in condensed matter materials. In the holographically dual gravity theory, these are analytic solutions to a five-dimensional non-linear-sigma (Skyrme) model dynamically coupled to Einstein gravity. This theory is shown to be holographically dual to \( \mathcal{N} \) = 4 Super-Yang-Mills theory coupled to an SU(2)-current. All solutions are fully backreacted asymptotically Anti-de Sitter (AdS) black branes or holes. One family of global AdS black hole solutions contains non-Abelian gauge field configurations with positive integer Chern numbers and finite energy density. Larger Chern numbers increase the Hawking-Page transition temperature. In the holographically dual field theory this indicates a significant effect on the deconfinement phase transition. Black holes with one Hawking temperature can have distinct Chern numbers, potentially enabling topological transitions. A second family of analytic solutions, rotating black branes, is also provided. These rotating solutions induce states with propagating charge density waves in the dual field theory. We compute the Hawking temperature, entropy density, angular velocity and free energy for these black holes/branes. These correspond to thermodynamic data in the dual field theory. For these states the energy-momentum tensor, (non-)conserved current, and topological charge are interpreted.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

{phi} Meson Photoproduction on the Nucleon and 4He Targets

We study the reaction mechanism of phi-meson photoproduction on the nucleon and He-4 targets by using a dynamical model based on a Hamiltonian. In addition to the dominant contribution of the Pomeron exchange, various meson exchanges are considered in the t channel to describe the CLAS data in the low energy region root s =(1.97-2.84) GeV. The direct phi radiations are taken into account in the s- and u-channels. The backward structures at root s approximate to 2.1 and 2.3 GeV are well reproduced by the inclusion of the s-channel N(2000,5/2(+)) and N(2300,1/2(+))resonances, respectively. We also consider the final phi N interactions by the gluon exchange, the direct phi N interactions, and the box diagrams arising from the couplings with the pi N, rho N, K Lambda, and K Sigma channels. The effects of the final state interactions are found to be very weak. Then the resulting Hamiltonian is used to study the coherent gamma He-4 -> phi He-4 reaction within the distorted-weave impulse approximation. The calculated differential cross sections account for the LEPS data quite well.

Kim, Sang-Ho↗
DOE OSTI2021

Tasting nuclear pasta made with classical molecular dynamics simulations

Nuclear clusters or voids in the inner crust of neutron stars were predicted to have various shapes collectively nicknamed nuclear pasta. Furthermore, the recent review by López, Dorso and Frank summarized their systematic investigations into properties especially the morphological and thermodynamical phase transitions of the nuclear pasta within a Classical Molecular Dynamics model, providing further stimuli to find more observational evidences of the predicted nuclear pasta in neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

DownScaleBench for developing and applying a deep learning based urban climate downscaling- first results for high-resolution urban precipitation climatology over Austin, Texas

Abstract Cities need climate information to develop resilient infrastructure and for adaptation decisions. The information desired is at the order of magnitudes finer scales relative to what is typically available from climate analysis and future projections. Urban downscaling refers to developing such climate information at the city (order of 1 – 10 km) and neighborhood (order of 0.1 – 1 km) resolutions from coarser climate products. Developing these higher resolution (finer grid spacing) data needed for assessments typically covering multiyear climatology of past data and future projections is complex and computationally expensive for traditional physics-based dynamical models. In this study, we develop and adopt a novel approach for urban downscaling by generating a general-purpose operator using deep learning. This ‘DownScaleBench’ tool can aid the process of downscaling to any location. The DownScaleBench has been generalized for both in situ (ground- based) and satellite or reanalysis gridded data. The algorithm employs an iterative super-resolution convolutional neural network (Iterative SRCNN) over the city. We apply this for the development of a high-resolution gridded precipitation product (300 m) from a relatively coarse (10 km) satellite-based product (JAXA GsMAP). The high-resolution gridded precipitation datasets is compared against insitu observations for past heavy rain events over Austin, Texas, and shows marked improvement relative to the coarser datasets relative to cubic interpolation as a baseline. The creation of this Downscaling Bench has implications for generating high-resolution gridded urban meteorological datasets and aiding the planning process for climate-ready cities.

Singh, Manmeet (ORCID:0000000233747149)↗
DOE OSTI2022

Transfer learning for smart buildings: A critical review of algorithms, applications, and future perspectives

Smart buildings play a crucial role toward decarbonizing society, as globally buildings emit about one-third of greenhouse gases. In the last few years, machine learning has achieved a notable momentum that, if properly harnessed, may unleash its potential for advanced analytics and control of smart buildings, enabling the technique to scale up for supporting the decarbonization of the building sector. In this perspective, transfer learning aims to improve the performance of a target learner exploiting knowledge in related environments. The present work provides a comprehensive overview of transfer learning applications in smart buildings, classifying and analyzing 77 papers according to their applications, algorithms, and adopted metrics. The study identified four main application areas of transfer learning: (1) building load prediction, (2) occupancy detection and activity recognition, (3) building dynamics modeling, and (4) energy systems control. Furthermore, the review highlighted the role of deep learning in transfer learning applications that has been used in more than half of the analyzed studies. The paper also discusses how to integrate transfer learning in a smart building's ecosystem, identifying, for each application area, the research gaps and guidelines for future research directions.

Pinto, G↗