Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum 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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Coupling of fully symmetric As phonon to magnetism in Ba(Fe 1-x Au x ) 2 As 2

In this work, we study the coupling of the fully symmetric vibration mode of arsenic atoms to magnetism in a Ba ( Fe 1 - x Au x ) 2 As 2 system by polarization-resolved Raman spectroscopy and neutron diffraction. In this system, there are two phase transitions: a tetragonal-to-orthorhombic structural phase transition at temperature T S and a magnetic phase transition into collinear spin-density wave (SDW) state at temperature T N ( ≤ T S ). T S and T N almost coincide in the pristine compound, whereas they differ by as much as 8 K for compounds with dilute gold substitution for iron. Raman coupling to the A g (As) phonon is forbidden for the X Y scattering geometry in the tetragonal phase above T S , whereas it becomes allowed in the orthorhombic phase below T S : The emerging mode's intensity indicates the lattice orthorhombicity. Furthermore, we find that upon cooling below T S , first, weak A g (As) phonon mode intensity appears in the X Y scattering geometry spectra; however, the mode's intensity is significantly enhanced in the magnetic phase below T N . The A g (As) phonon also shows an asymmetric line shape below T N and an anomalous linewidth broadening upon Au doping. We describe the anomalous behavior of the A g (As) mode in the X Y scattering geometry using a Fano model involving the A g (As) phonon interacting with the B 2 g ( D 4 h )-symmetry-like electron-hole continuum. We conclude that the temperature dependence of light coupling amplitude to the A g (As) phonon follows the evolution of the magnetic order parameter M ( T ) . We propose that the intensity enhancement of the A g (As) phonon in the X Y scattering geometry below T N is due to electronic anisotropy induced by the collinear SDW order parameter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Gauge-symmetrization method for energy-momentum tensors in high-order electromagnetic field theories

For electromagnetic field theories, canonical energy-momentum conservation laws can be derived from the underpinning spacetime translation symmetry according to the Noether procedure. However, the canonical energy-momentum tensors (EMTs) are neither symmetric nor gauge-symmetric (gauge invariant). The Belinfante-Rosenfeld (BR) method is a well-known procedure to symmetrize the EMTs, which also renders them gauge symmetric for first-order field theories. High-order electromagnetic field theories appear in the study of gyrokinetic systems for magnetized plasmas and the Podolsky system for the radiation reaction of classical charged particles. For these high-order field theories, gauge-symmetric EMTs are not necessarily symmetric and vice versa. In the present study, we develop a new gauge-symmetrization method for EMTs in high-order electromagnetic field theories. The Noether procedure is carried out using the Faraday tensor $F_{μν}$, instead of the 4-potential $A_{μ}$, to derive a canonical EMT $T^{μν}_{N}$. We show that the gauge-dependent part of $T^{μν}_{N}$ can be removed using the displacement-potential tensor $F^{σμν}$ ≡ $D^{σμ}A^{ν}/4π$, where $D^{σμ}$ is the antisymmetric electric displacement tensor. This method gauge-symmetrizes the EMT without necessarily making it symmetric, which is adequate for applications not involving general relativity. For first-order electromagnetic field theories, such as the standard Maxwell system, $F^{σμν}$ reduces to the familiar BR superpotential $S^{σμν}$, and the method developed can be used as a simpler procedure to calculate $S^{σμν}$ without employing the angular momentum tensor in 4D spacetime. When the electromagnetic system is coupled to classical charged particles, the gauge-symmetrization method for EMTs is shown to be effective as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The Lick AGN Monitoring Project 2016: Velocity-resolved Hβ Lags in Luminous Seyfert Galaxies

We carried out spectroscopic monitoring of 21 low-redshift Seyfert 1 galaxies using the Kast double spectrograph on the 3 m Shane telescope at Lick Observatory from 2016 April to 2017 May. Targeting active galactic nuclei (AGNs) with luminosities of λL$_{λ}$(5100 Å) ≈ 10$^{44}$ erg s$^{-1}$ and predicted Hβ lags of ~20–30 days or black hole masses of 10$^{7}$–10$^{8.5}$M$_{⊙}$, our campaign probes luminosity-dependent trends in broad-line region (BLR) structure and dynamics as well as to improve calibrations for single-epoch estimates of quasar black hole masses. Here we present the first results from the campaign, including Hβ emission-line light curves, integrated Hβ lag times (8–30 days) measured against V-band continuum light curves, velocity-resolved reverberation lags, line widths of the broad Hβ components, and virial black hole mass estimates (10$^{7.1}$–10$^{8.1}$M$_{⊙}$). Our results add significantly to the number of existing velocity-resolved lag measurements and reveal a diversity of BLR gas kinematics at moderately high AGN luminosities. AGN continuum luminosity appears not to be correlated with the type of kinematics that its BLR gas may exhibit. Follow-up direct modeling of this data set will elucidate the detailed kinematics and provide robust dynamical black hole masses for several objects in this sample.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Multiscale modeling of microwave-heated multiphase systems

Modeling of microwave heating of multiphase systems requires resolving the electromagnetic field from the single particle to the entire cavity. We introduce a multiscale methodology for computationally affordable simulations of microwave heating of multiphase systems consisting of one phase dispersed in a continuum phase. The methodology homogenizes the original multiphase system, by taking advantage of the large separation of length scales, to calculate the effective permittivity, effective thermal conductivity, and volumetric power absorbed in each phase. The methodology is rigorously assessed against particle-resolved detailed numerical simulations of model systems. We demonstrate a significant reduction in computational cost while retaining the accuracy of detailed simulations. Lastly, the present methodology enables high throughput exploration of novel multiphase system designs utilizing microwave heating.

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

Topology of the $O(3)$ non-linear sigma model under the gradient flow

The $O(3)$ non-linear sigma model (NLSM) is a prototypical field theory for QCD and ferromagnetism, featuring topological qualities. Though the topological susceptibility should vanish in physical theories, lattice simulations of the NLSM find that the topological susceptibility diverges in the continuum limit. We introduce a $\theta$-term and analyze the topological charge as a function of $\theta$ under the gradient flow. Our results show that divergence persists in the presence of the flow, even at non-zero $\theta$.

Thomas, Stuart↗
DOE OSTI2022

A Micromorphic Length-Scale Coupling Framework for the Determination of Higher-Order Constitutive Models and the Multi-Scale Simulation of Heterogeneous Materials [Thesis]

Heterogeneous materials and materials with complex microstructures pose a unique challenge in the development of accurate models of their response to external stimuli. These difficulties principally arise due to the difficulty in characterizing and modeling the constituents and their interactions. It is usually possible, though non-trivial, to construct an explicit representation of the microstructure (a direct numeric simulation or DNS) but the method by which the complex modes of deformation and other processes can be homogenized to a reduced-order approximation is frequently unclear. Many homogenization approaches are ad-hoc and lack a strong justification beyond ease of computation. Furthermore, the homogenization approach can, in some cases, not utilize the full breadth of information available in the computation of the macro-scale stresses and deformation measures. It is also noteworthy that homogenization, by its very nature, will tend to obfuscate details of processes occurring at the lower length-scale. It is therefore of interest to include as much information as possible in the construction of the reduced-order model so as to be predictive in a variety of loading environments. We here present a length-scale bridging technique based upon the micromorphic continuum mechanics of Eringen which incorporates volume and surface area averages as a part of its construction. This approach enforces the balance equations at the micro-scale and then studies the effect of the spatially varying nature on the macro-scale. This leads, naturally, to further balance equations which are solved at the macro-scale. This work details a homogenization framework which arises naturally from the micromorphic construction of Eringen attempting to introduce no definitions beyond which are justifiable from micro-structural considerations. One of the results of this effort is the so-called “micromorphic filter” which has been developed and applied to several DNS to demonstrate its effectively. In order to determine the macroscopic degrees of freedom we utilize the special case of an overlap coupling technique where the macro-scale is fully constrained to the micro-scale. This enables us to study the resulting material properties as expected but also allows us to further study the boundary conditions on the additional degrees of freedom at the macro-scale.

42 ENGINEERING↗
DOE OSTI2021

UCGretina geant4 simulation of the GRETINA Gamma-Ray Energy Tracking Array

UCGretina, a geant4 simulation of the GRETINA gamma-ray tracking array of highly-segmented high-purity germanium detectors is described. We have developed a model of the array, in particular of the Quad Module and the capsules, that gives good agreement between simulated and measured photopeak efficiencies over a broad range of gamma-ray energies and reproduces the shape of the measured Compton continuum. Both of these features are needed in order to accurately extract gamma-ray yields from spectra collected in in-beam gamma-ray spectroscopy measurements with beams traveling at v / c ≳0.3 at the National Superconducting Cyclotron Laboratory and the Facility for Rare Isotope Beams. In the process of developing the model, we determined that millimeter-scale layers of passive germanium surrounding the active volumes of the simulated crystals must be included in order to reproduce measured photopeak efficiencies. In conclusion, we adopted a simple model of effective passive layers and developed heuristic methods of determining passive-layer thicknesses by comparison of simulations and measurements for a single crystal and for the full array. Prospects for future development of the model are discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

An iterative method to deblend AGN-Host contributions for Integral Field spectroscopic observations

ABSTRACT We present a new iterative deblending method to separate the host galaxy (HG) and their Active Galactic Nuclei (AGNs) emission with the use of Integral Field spectroscopic (IFS) data. The method decomposes the resolved HG emission from the unresolved AGN emission by modelling the two-dimensional surface brightness (SB) profile of the point-spread function (PSF) and the two-dimensional SB HG continuum simultaneously per each monochromatic slide. Our method does not require any prior information about the observed SB profile or a detailed fitting of the PSF, making it ideal for the automatic analysis of large galaxy samples. In this work, we test the quality of our method, its advantages, and its disadvantages. We test our method by using a set of IFS mock data cubes to quantify the reliability of our deblending process and further compare our method with the qdblend3d analysis tool. Furthermore, we applied our method to three data cubes selected from the MaNGA survey according to the dominance of either its HG or its AGN. We show that our deblending method is capable of disengaging the bright, non-resolved AGN emission from the HG continuum and its narrow emission lines. However, the decoupling depends on how well the IFS spatially resolves the PSF, and on the relative flux intensity of the HG-AGN. Therefore, the method is ideal for disentangling the bright-flux contribution from AGN-dominated spectra.

Ibarra-Medel, H. (ORCID:0000000297906313)↗
DOE OSTI2022

New narrow resonances observed in the unbound nucleus F 15

The structure of the unbound 15 F nucleus is investigated using the inverse kinematics resonant scattering of a radioactive 14 O beam impinging on a CH 2 target. The analysis of 1 H( 14 O,p) 14 O and 1 H( 14 O,2p) 13 N reactions allowed the confirmation of the previously observed narrow 1/2 - resonance, near the two-proton decay threshold, and the identification of two new narrow 5/2 - and 3/2 - resonances. The newly observed levels decay by 1p emission to the ground of 14 O, and by sequential 2p emission to the ground state (g.s.) of 13 N via the 1 - resonance of 14 O. Gamow shell model (GSM) analysis of the experimental data suggests that the wave functions of the 5/2 - and 3/2 - resonances may be collectivized by the continuum coupling to nearby 2p- and 1p- decay channels. Finally, the observed excitation function 1 H( 14 O, p) 14 O and resonance spectrum in 15 F are well reproduced in the unified framework of the GSM.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Semi-Dirac Fermions in a Topological Metal

Topological semimetals with massless Dirac and Weyl fermions represent the forefront of quantum materials research. In two dimensions, a peculiar class of fermions that are massless in one direction and massive in the perpendicular direction was predicted 16 years ago. These highly exotic quasiparticles—the semi-Dirac fermions—ignited intense theoretical and experimental interest but remain undetected. Using magneto-optical spectroscopy, we demonstrate the defining feature of semi-Dirac fermions— B 2 / 3 scaling of Landau levels—in a prototypical nodal-line metal ZrSiS. In topological metals, including ZrSiS, nodal lines extend the band degeneracies from isolated points to lines, loops, or even chains in the momentum space. With calculations and theoretical modeling, we pinpoint the observed semi-Dirac spectrum to the crossing points of nodal lines in ZrSiS. Crossing nodal lines exhibit a continuum absorption spectrum but with singularities that scale as B 2 / 3 at the crossing. Our work sheds light on the hidden quasiparticles emerging from the intricate topology of crossing nodal lines and highlights the potential to explore quantum geometry with linear optical responses. Published by the American Physical Society 2024

Shao, Yinming (ORCID:0000000228910028)↗
DOE OSTI2023

ProvLight: Efficient Workflow Provenance Capture on the Edge-to-Cloud Continuum

Modern scientific workflows require hybrid infrastructures combining numerous decentralized resources on the IoT/Edge interconnected to Cloud/HPC systems (aka the Computing Continuum) to enable their optimized execution. Understanding and optimizing the performance of such complex Edge-to-Cloud workflows is challenging. Capturing the provenance of key performance indicators, with their related data and processes, may assist in understanding and optimizing workflow executions. However, the capture overhead can be prohibitive, particularly in resource-constrained devices, such as the ones on the IoT/Edge.To address this challenge, based on a performance analysis of existing systems, we propose ProvLight, a tool to enable efficient provenance capture on the IoT/Edge. We leverage simplified data models, data compression and grouping, and lightweight transmission protocols to reduce overheads. We further integrate ProvLight into the E2Clab framework to enable workflow provenance capture across the Edge-to-Cloud Continuum. This integration makes E2Clab a promising platform for the performance optimization of applications through reproducible experiments.We validate ProvLight at a large scale with synthetic workloads on 64 real-life IoT/Edge devices in the FIT IoT LAB testbed. Evaluations show that ProvLight outperforms state-of-the-art systems like ProvLake and DfAnalyzer in resource-constrained devices. ProvLight is 26—37x faster to capture and transmit provenance data; uses 5—7x less CPU; 2x less memory; transmits 2x less data; and consumes 2—2.5x less energy. ProvLight [1] and E2Clab [2] are available as open-source tools.

Rosendo, Daniel↗
DOE OSTI2022

An Optical Spectrum of the Diffuse Galactic Light from BOSS and IRIS

We present a spectrum of the diffuse Galactic light (DGL) between 3700 and 10,000 Å, obtained by correlating optical sky intensity with far-infrared dust emission. We use nearly 250,000 blank-sky spectra from BOSS/SDSS-III together with IRIS-reprocessed maps from the IRAS satellite. The larger sample size compared to SDSS-II results in a factor-of-2 increase in signal to noise. We combine these data sets with a model for the optical/far-infrared correlation that accounts for self-absorption by dust. The spectral features of the DGL agree remarkably well with the features present in stellar spectra. There is evidence for a difference in the DGL continuum between the regions covered by BOSS in the northern and southern Galactic hemispheres. We interpret the difference at red wavelengths as the result of a difference in stellar populations, with mainly old stars in both regions, but a higher fraction of young stars in the south. There is also a broad excess in the southern DGL spectrum over the prediction of a simple radiative transfer model, without a clear counterpart in the north. We interpret this excess, centered at ~6500 Å, as evidence for luminescence in the form of extended red emission. The observed strength of the 4000 Å break indicates that at most ~7% of the dust-correlated light at 4000 Å can be due to blue luminescence. Our DGL spectrum provides constraints on dust scattering and luminescence, independent of measurements of extinction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Circuit Quantum Electrodynamics in Hyperbolic Space: From Photon Bound States to Frustrated Spin Models

Circuit quantum electrodynamics is one of the most promising platforms for efficient quantum simulation and computation. In recent groundbreaking experiments, the immense flexibility of superconducting microwave resonators was utilized to realize hyperbolic lattices that emulate quantum physics in negatively curved space. Here we investigate experimentally feasible settings in which a few superconducting qubits are coupled to a bath of photons evolving on the hyperbolic lattice. We compare our numerical results for finite lattices with analytical results for continuous hyperbolic space on the Poincaré disk. We find good agreement between the two descriptions in the long-wavelength regime. We show that photon-qubit bound states have a curvature-limited size. We propose to use a qubit as a local probe of the hyperbolic bath, for example, by measuring the relaxation dynamics of the qubit. We find that, although the boundary effects strongly impact the photonic density of states, the spectral density is well described by the continuum theory. We show that interactions between qubits are mediated by photons propagating along geodesics. We demonstrate that the photonic bath can give rise to geometrically frustrated hyperbolic quantum spin models with finite-range or exponentially decaying interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Description of the Proton-Decaying 0 2 + Resonance of the α Particle

The recent precise experimental determination of the monopole transition form factor from the ground state of 4 He to its $0$$^{+}_{2}$ resonance via electron scattering has reinvigorated discussions about the nature of this first excited state of the α particle. The $0$$^{+}_{2}$ state has been traditionally interpreted in the literature as the isoscalar monopole resonance (breathing mode) or, alternatively, as a particle-hole shell-model excitation. To better understand the nature of this state, which lies only ~410 keV above the proton emission threshold, we employ the coupled-channel representation of the no-core Gamow shell model. By considering the [ 3 H + p], [ 3 He + n], and [ 2 H + 2 H] reaction channels, we explain the excitation energy and monopole form factor of the $0$$^{+}_{2}$ state. Furthermore, we argue that the continuum coupling strongly impacts the nature of this state, which carries characteristics of the proton decay threshold.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Application of the core shell model for strengthening polymer filament interfaces

While surface segregation has been intensively studied in thin films, relatively little has been investigated in filament geometries. In contrast to thin films, the filament surface forms a continuum, with interactions both at the free air and internal filament interface, which makes it difficult to achieve an equilibrium configuration. This asymmetry has become particularly apparent with the increasing use of additive manufacturing where entire structures are constructed via layered deposition of filamentous layers from a moving thermal nozzle. In addition, this process is highly non-equilibrated with poor thermal retention and is responsible for insufficient filament–filament interfacial fusion and internal pore formation along the printing direction. These features all contribute to poor mechanical and structural properties, which are exacerbated in semi-crystalline polymers, such as Polylactic Acid (PLA), where the melting point, T m , is much higher than the glass transition, T g . Even though system is liquid the molecules become ordered within the crystal structure and this ordering interferes with interdiffusion, which in turn prevents proper interfacial fusion. Capitalizing on the differential in surface tension, we show that addition of less than 1%, or trace amounts, of longer chain polymers, an oriented core shell structure, will be formed with a self-limiting thickness of three time the radius of gyration, R g , of the tracer chains. The lower T g of the tracers and their orientation makes them poised for interdiffusion when subsequent layers are deposited producing structures with significant improvement of the mechanical properties and internal structural integrity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A description of structured waves in shock compressed particulate composites

Dynamic compression of composite materials is of scientific interest because the mechanical mismatch between internal phases challenges continuum theories. Typical assumptions about steady wave propagation and quasi-instantaneous state changes require reexamination along with the need for time-dependent models. To that end, data and models are presented here for the shock compression of an idealized particulate composite. To serve as a generic representative of this material class, a polymer matrix was filled with tungsten particles, ranging from 1 to 50 vol. %. This creates a simple microstructure containing randomly scattered particles with an extreme impedance mismatch to the binding matrix. These materials were parallel plate impact loaded by Al flyers traveling at 1.8–5.0 km/s. Velocimetry provided records of the equilibrium state and the compression wave structure for each case with trends quantified by an empirical fit. The same quantities were also studied as a function of the wave's propagation distance. A homogenized viscoelastic model then made it possible to progress from cataloging trends to making predictions. Starting from a Mie–Gruneisen equation of state, additional time varying terms were added to capture the transient response. After calibration, accurate predictions of the steady wave structure were possible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

The QPO in RE J1034+396 Originates in the Hot Corona

RE J1034+396 is one of the few active galactic nuclei (AGNs) with a significant quasiperiodic oscillation (QPO). The QPO has been observed in over 1 Ms of XMM-Newton observations spanning over a decade. We investigate the power spectral density function (PSD) of seven long (∼90 ks) XMM-Newton observations of the AGN RE J1034+396 in two energy bands. The soft (0.3–0.5 keV) band targets emission from the disk, while the hard (2–7 keV) band isolates the primary X-ray continuum emission from the corona. The QPO is significantly detected in the hard band of five of the seven observations. The best-fitting models indicate that the QPO detection in both bands is entirely attributable to the coronal emission with no additional contribution from the disk. This explains the strong coherence between the hard and soft bands at the QPO frequency. The covariance spectrum is consistent with this picture as the variability at QPO frequencies is attributed solely to fluctuations in the hot corona. The time lag as a function of energy is well described by a ∼2000 s intrinsic soft lag, resulting from the disk responding to emission from the corona, that undergoes phase wrapping at approximately the QPO frequency. By demonstrating that in this system the QPO arises in the corona, we provide new insights into the mechanisms generating QPOs.

Active galactic nuclei↗