Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simplified analysis”

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 127 records · Page 7

Calculation of mode lifetimes in weakly anharmonic solids using self-consistent ensemble eigenstates of the Liouvillian

Recent schemes for the efficient calculation of vibrational mode frequencies and lifetimes in weakly anharmonic systems such as solid crystal structures have shown promise, with simple relationships being observed between low order moments of the Liouvillian and the observed lifetime. While only shown empirically, two parameters seemed sufficient to determine the lifetime to high accuracy over a large range of temperatures for a few simple low degree of freedom models and a Lennard–Jones solid. In this paper we analytically reproduce the numerically observed relationship using a few generally applicable assumptions about the density of states. This advance is made by studying the autocorrelation based on an “ensemble eigenstate” of the Liouvillian. We propose the use of this ensemble eigenstate as the basis for autocorrelations to simplify calculations in richer systems. Finally, an analysis is also performed of the relationship between canonical and microcanonical autocorrelation functions and a long standing error in the long time behavior of such functions originally presented by Sen, Sinkovits, and Chakravarti is corrected.

36 MATERIALS SCIENCE↗

Unconventional nonlinear Hall effects in twisted multilayer 2D materials

We present the first investigation of unusual nonlinear Hall effects in twisted multilayer 2D materials. Contrary to expectations, our study shows that these nonlinear effects are not merely extensions of their monolayer counterparts. Instead, we find that stacking order and pairwise interactions between neighboring layers, mediated by Berry curvatures, play a pivotal role in shaping their collective nonlinear optical response. By combining large-scale Real-Time Time-Dependent Density Functional Theory (RT-TDDFT) simulations with model Hamiltonian analyses, we demonstrate a remarkable second-harmonic transverse response in hexagonal boron nitride four-layers, even in cases where the total Berry curvature cancels out. Furthermore, our symmetry analysis of the layered structures provides a simplified framework for predicting nonlinear responses in multilayer materials in general. Our investigation challenges the prevailing understanding of nonlinear optical responses in layered materials and opens new avenues for the design and development of advanced materials with tailored optical properties.

36 MATERIALS SCIENCE↗

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spring 2020 Dissertation Update [Slides]

An update is provided on the dissertation underway and what has been learned thus far. Work thus far: Developed analytic models for each region in the spent fuel cask – Used to identify and explain physical processes which create features in detailed casks; Developed simplified computational models to identify details not seen in analytic models; SC were calculated in the fuel region; The difference in SC’s between the analytic model and the simplified computational model were identified. Outstanding issues: Analysis of analytic models in stainless steel and carbon steel – These materials are thin and have few features (just the slope); Discrepancies between absorption SC’s in fuel; Create a test problem to show the effects of the high energy resonances in the fuel region. Future Work: Sensitivity analysis needs to be continued through the cask – The detailed model will be added to the remaining materials; Sensitivity analysis paper; and, Addressing outstanding issues.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗

Absorptive corrections to the forward elastic proton-proton analyzing power A N

In the RHIC Polarized Atomic Hydrogen Gas Jet Target (HJET) measurements of the protonproton analyzing power A N (t), absorptive corrections to the electromagnetic form factor were not considered in the data analysis. Recently, these corrections were evaluated theoretically. However, the results were presented in form of integral expressions which were calculated numerically. In this note, the expressions are simplified to make them more suitable for the HJET data analysis. Possible effect of the absorptive corrections in the A N (t) measurement is discussed.

43 PARTICLE ACCELERATORS↗

Hydrological analysis in watersheds with a variable-resolution global climate model (VR-CESM)

Traditionally, watershed-scale hydrology is simulated by distributed hydrological models with offline meteorological forcing data, or by regional regional climate models that link atmospheric and land hydrology interactions. Global climate models (GCMs) are rarely used to study watershed-scale hydrology due to the relatively coarse grid resolution, computationally expensive downscaling, and simplified physical processes. Recently, however, watershed-scale hydrology analysis has become possible in GCMs due to the development of variable-resolution GCMs that dynamically couple the hydrological processes between atmosphere and land systems at fine resolutions in selected regions and coarse resolution elsewhere. In this study, we used the variable-resolution Community Earth System Model (VR-CESM) with refined-resolution (14 km) in the western U.S. and eastern China to evaluate watershed-scale hydrology. Furthermore, we compared the historical VR-CESM outputs with gauge measurements and other hydrological models (e.g., National Water Model in the U.S.) and calibrated the subsurface runoff capacities in four mountainous watersheds. An RCP8.5 projection from 2007 to 2050 is used to estimate the impact of changing precipitation and snow climatology on watershed hydrology. We also analyzed the long-term runoff variability and the possibility of extreme runoff events as simulated by the VR-CESM. Although calibration is not possible in larger-scale watersheds, VR-CESM simulates the long-term annual variability of watersheds and provides insights on climate change impacts on hydrology. We conclude that refined-resolution VR-CESM can be used for watershed-scale hydrology analysis to understand water resources and water balance, in addition to traditional watershed-scale hydrological models. It enables hydrological analysis at multiple watersheds in one simulation and can help understand the two-way dynamics between land surface hydrology and atmospheric processes, and is especially practical for projecting climate change impacts. However, it is challenging to apply VR-CESM for hydrologic analysis in regulated watersheds as human factors (e.g., pumping, irrigation, water diversion) have not been fully addressed in VR-CESM.

54 ENVIRONMENTAL SCIENCES↗

Nonlinear susceptibilities for weakly turbulent magnetized plasma: Electromagnetic formalism

This is a companion paper to the previous work [P. H. Yoon, Phys. Plasmas 31, 032309 (2024)] in which the nonlinear susceptibilities of weakly turbulent magnetized plasma are derived under a simplifying assumption of electrostatic interaction. The present paper extends the analysis to a general situation of electromagnetic interaction. The main novelty of the previous and present papers is that by employing the Bessel function addition theorem, the mathematical definitions for the susceptibilities are substantially simplified, a procedure that has not been discussed in the existing literature. In the present paper, a full set of Maxwell’s equations are considered in conjunction with the nonlinear Vlasov equation, which is solved by a perturbative method. The result is a fully general nonlinear susceptibility, given in tensorial form, which is applicable for weakly turbulent magnetized plasmas.

Materials Science↗

Catching the Killer: Dynamic Disorder Design Rules for Small–Molecule Organic Semiconductors

Small-molecule organic semiconductors are promising materials for applications ranging from solar cells to medical sensors. Their nearly infinite design space means that, in theory, it is possible to tailor the material to the exact specifications of a particular application. In reality; however, design rules to improve mobility (μ) remain elusive because of the complex calculations required to understand its limiters. Transient localization theory posits that charge carriers are slowed down by the collective phonon motions, called dynamic disorder, which localize charge carriers temporarily. Traditionally, a mode-by-mode analysis is performed to try and identify “killer” modes. Herein, the flaws are demonstrated with this analysis and present a streamlined simulation workflow that simplifies simulation of dynamic disorder and enables engineering of new molecular structures based on phonon dynamics. This workflow is combined with a novel visualization technique that enables per-atom insights into how μ is limited. Furthermore, this workflow is applied and analyzed to a series of -acenes and a series of BTBT-based molecules. Then design rules are identified in each series that identify possible ways to improve the μ beyond current experimental limits using phonon engineering.

36 MATERIALS SCIENCE↗

A quasi-meshfree method for constructing boundary-aware reproducing bases on geometrically complex domains using manifold geodesics

In many applications, physical domains are geometrically complex making it challenging to perform coarse-scale approximation. A defeaturing process is often used to simplify the domain in preparation for approximation and analysis at the coarse scale. In this paper, a methodology is presented for constructing a coarse-scale reproducing basis on geometrically complex domains given an initial fine-scale mesh of the fully featured domain. The initial fine-scale mesh can be of poor quality and extremely refined. The construction of the basis functions begins with a coarse-scale covering of the domain and generation of weighting functions with local support. Manifold geodesics are used to define distances within the local support for general applicability to non-convex domains. Conventional moving least squares is used to construct the coarse-scale reproducing basis. Applications in quasi-interpolation and linear elasticity are presented.

97 MATHEMATICS AND COMPUTING↗

Development of an accelerator-based neutron source to prototype Mo-99 production, part I: A liquid LBE windowless target

In this article, Molybdenum-99 (Mo-99)’s decay product, technetium-99 (Tc-99 m), is one of the most critical isotopes for medical diagnostics. To provide U.S. domestic supply of Mo-99 without using high-enriched uranium (HEU), a subcritical uranium target assembly (UTA) is irradiated by an accelerator-based neutron source to create Mo-99 through fission. This study discusses the development of the accelerator-based neutron source. The high-energy electrons from the accelerator irradiate a liquid lead-bismuth eutectic (LBE) target to produce neutrons. Part I of this work focuses on numerical and experimental analysis towards the development of a liquid LBE windowless target. Unlike the existing windowless targets in literature, the current design creates a vertical free surface for a beam to irradiate. First, a hydrodynamic analysis of the LBE windowless target is performed. Simplified analytical calculations are assisted by 2D computational fluid dynamics (CFD) simulations to design the target, with the focus on eliminating recirculation zones and avoiding cavitation. With the optimized geometry, the experimental study is performed to investigate the flow hydrodynamics using liquid LBE. The experiments (1) compare pressure drop in the system to correlation predictions; (2) visualize the free surface liquid LBE flow from the beam view; (3) validate the LBE flow profile using temperature sensitive paint from the side view; and (4) validate the liquid LBE film thickness using gamma densitometer measurements. Second, the power handling capability of the designed windowless target is investigated. The divider plate in the current design is susceptible to overheating due to the thin LBE film in front. As LBE erosion and corrosion is likely to occur at an LBE velocity of 2.0 m/s and temperature above 500 °C, a power limit of 10 kW of beam power was established to prevent this corrosion from occurring, which is calculated by a Nusselt number correlation. The divider plate surface temperature at 10 kW agrees well with the 3D CFD simulation results. Part I demonstrates the fundamental physics in liquid LBE windowless target design and associated testing. A companion paper, Part II will demonstrate how to couple this windowless target into the Mo-99 production system, including an accelerator system operating under an ultra-high vacuum and the UTA cooled by water at room temperature.

43 PARTICLE ACCELERATORS↗

Ergodic Lagrangian dynamics in a superhero universe

We present a fictional scenario that, while undeniably whimsical, provides the foundation for a unique exercise in extended problem solving, physics analysis, and quantitative model development. Starting with the foundational premise of the Wild Cards shared-world superhero universe, we demonstrate how a variety of concepts appropriate to the advanced undergraduate level—ergodicity, functional analysis, Lagrangian mechanics, and the ever-important simplifying approximation—can be combined into a rich, coherent mathematical model. The goal of this case study is to develop a useful pedagogical exercise in exploring an open-ended research question that presents, at first glance, no clear path forward. Being both eclectic and lengthy, this exercise offers a unique way for students to apply their core physics and mathematics education. It is perhaps best used within a senior honors seminar or within a brief (e.g., January term) elective class.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Charged Microdroplets in Air: Characterization and Use in Surface Functionalization and Nanomaterials Preparation (Final Technical Report)

The major theme of this project is the exploration of ions and charged droplets in the open air. The specific goals are (1) focusing ions in the open air, (2) separating and measuring the sizes of charged droplets, (3) examining chemical reactions occurring within droplets, (4) examining chemical reactions occurring between charged droplets and chemical vapors, and (5) characterizing surfaces post droplet deposition. The analysis of ions in the open air greatly simplifies instrumentation as it removes the need for large, costly vacuum pumps. Major strides in this project include the sizing of electrospray droplets using structured illumination microscopy, 3D printed drift tube ion mobility (both standalone and coupled with mass spectrometry), and accelerated Suzuki reactions in Leidenfrost droplets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Integrators

Many electromagnetic measurements detect signals that are the time derivative of the actual quantity of interest. Examples include so-called “B-dot” and “D-dot” detectors that are used for measuring transient magnetic and electric fields. Integration of these signals to recover the quantity of interest is performed either by hardware integrators in the signal line, or by software coding in the analysis programming. Hardware integration is usually done using a classic resistive-capacitive (RC) circuit, or internally in the detector by induction (L/R) or stray capacitance (RC). These hardware integration methods are only approximate, and must be corrected for distortion. This is usually called “droop correction.” In this note, we will examine this correction in some detail. Transient analysis via Laplace transforms will be used to simplify the mathematics. A rudimentary understanding of this method is assumed.

42 ENGINEERING↗

Complete Demonstration of a Prototype Version of FORCE User Interface and Conduct Analyst Survey Collecting Feedback on Interface Features and Usability

In 2024 the US Department of Energy (DOE) Office of Nuclear Energy (NE) Integrated Energy System (IES) program continued to develop the Framework for Optimization of Resources and Economics (FORCE) analysis ecosystem into a more traditional toolset with simplified software installation, automated workflows, and interactive results visualization. The DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) Workbench continued to be leveraged for user input, application workflow and runtime environment, and interactive results visualization capabilities. This report documents the demonstration of a FORCE User Interface (UI) prototype and the results of a survey of analysts’ using the Holistic Energy Resource Optimization Network (HERON) tool in FORCE with the prototype UI.

97 MATHEMATICS AND COMPUTING↗

Simplifying Geospatial Workflows with GeoGridFusion

Growing demands to understand PV deployment and reliability across an expanding range of climates, along with increasing computational power, are driving the need for simplified computing tools that support large-scale analysis and geospatial workflows. While existing open-source libraries and PV system modeling tools offer extensive collections of empirical and physical models, they often lack the ability to scale to large geospatial datasets. The tool presented here, GeoGridFusion, enables PV modelers to store and utilize gridded geospatial data from sources such as NSRDB, PVGIS, and others. GeoGridFusion harmonizes diverse datasets, taxonomies, and nomenclatures, and includes utilities that support intuitive geospatial area selection.

14 SOLAR ENERGY↗

Chemical timescale effects on detonation convergence

Numerical simulations of detonation-containing flows have emerged as crucial tools for designing next- generation power and propulsion devices. As these tools mature, it is important for the combustion community to properly understand and isolate grid resolution effects when simulating detonations. To this end, the objective of this work is to provide a comprehensive analysis of the numerical convergence of unsteady detonation simulations, with focus on isolating the impacts of chemical timescale modifications on convergence characteristics in the context of operator splitting. With the aid of an AMReX-based adaptive mesh refinement flow solver-which enables resolutions up to ($\mathcal{O}$ (1000) cells-per-induction length-the convergence analysis is conducted using two kinetics configurations: (1) the simplified three-step Arrhenius-based model mechanism of Short and Quirk (1997), where chemical timescales in the detonation are modified by adjusting activation energies in the initiation and branching reactions, and (2) a detailed hydrogen- air mechanism, where the chemical timescales are adjusted by varying the ambient pressure. The convergence of unsteady self-sustained detonations in one-dimensional channels is then analyzed with reference to steady-state theoretical baseline solutions using these mechanisms. The goal of the analysis is to provide a detailed comparison of the effects of grid resolution on both macroscopic (peak pressures and wave speeds) and microscopic (wave structure) quantities of interest, drawing connections between the deviations from steady-state baselines and minimum chemical timescales. In particular, chemical timescale reductions were found to have minimal impact on the convergence of macroscopic properties. Furthermore, analyses of microscopic convergence trends, particularly in the reaction front location, revealed a key insight: maintaining the induction time while eliminating prohibitive chemical timescales through mechanism simplifications and combustion modeling can significantly enhance detonation convergence properties. Ultimately, this work uncovers resolution-dependent unsteady detonation convergence regimes and highlights the important role played by not only the chemical timescales, but also the ratio between the chemical timescale and induction time on the numerical convergence of the detonation wave structure.

Adaptive mesh refinement↗

Topological structure of complex predictions

Abstract Current complex prediction models are the result of fitting deep neural networks, graph convolutional networks or transducers to a set of training data. A key challenge with these models is that they are highly parameterized, which makes describing and interpreting the prediction strategies difficult. We use topological data analysis to transform these complex prediction models into a simplified topological view of the prediction landscape. The result is a map of the predictions that enables inspection of the model results with more specificity than dimensionality-reduction methods such as tSNE and UMAP. The methods scale up to large datasets across different domains. We present a case study of a transformer-based model previously designed to predict expression levels of a piece of DNA in thousands of genomic tracks. When the model is used to study mutations in the BRCA1 gene, our topological analysis shows that it is sensitive to the location of a mutation and the exon structure of BRCA1 in ways that cannot be found with tools based on dimensionality reduction. Moreover, the topological framework offers multiple ways to inspect results, including an error estimate that is more accurate than model uncertainty. Further studies show how these ideas produce useful results in graph-based learning and image classification.

Computer Science↗