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 991 records · Page 55

Transverse spin asymmetries and the electron Yukawa coupling at an FCC-ee

We show that measurements of single transverse-spin asymmetries at an 𝑒 + ⁢𝑒 − collider can enhance the sensitivity to the electron Yukawa coupling, possibly enabling observation of the Standard Model value for this quantity. We demonstrate that the significance in both the 𝑏⁢$\overline{𝑏}$ and semileptonic 𝑊⁢𝑊 final states can be enhanced by factors of up to 3 compared to inclusive cross-section determinations of this coupling for transversely polarized electrons. If the positrons can be simultaneously longitudinally polarized, even at the level of 30%, the significance can be enhanced by a factor of 5 or more. The method utilizes quantum interference between the Higgs signal and the continuum background and is also applicable in other 𝑊⁢𝑊 and 𝑍⁢𝑍 final states.

Phenomenology↗

An asymmetric electron-scattering photosphere around optical tidal disruption events

We report a star crossing the tidal radius of a supermassive black hole will be spectacularly ripped apart with an accompanying burst of radiation. A few tens of such tidal disruption events have now been identified in optical wavelengths, but the exact origin of the strong optical emission remains inconclusive. Here we report polarimetric observations of three tidal disruption events. The continuum polarization appears independent of wavelength, while emission lines are partially depolarized. These signatures are consistent with photons being scattered and polarized in an envelope of free electrons. An almost axisymmetric photosphere viewed from different angles is in broad agreement with the data, but there is also evidence for deviations from axial symmetry before the peak of the flare and significant time evolution at early times, compatible with the rapid formation of an accretion disk. By combining a super-Eddington accretion model with a radiative transfer code, we simulate the polarization degree as a function of disk mass and viewing angle and we show that the predicted levels are compatible with the observations for extended reprocessing envelopes of ~1,000 gravitational radii. Spectropolarimetry therefore constitutes a new observational test for tidal disruption event models, and opens an important new line of exploration in the study of tidal disruption events.

79 ASTRONOMY AND ASTROPHYSICS↗

Co-Based superalloy morphology evolution: A phase field study based on experimental thermodynamic and kinetic data

Cobalt-based superalloys with two phase $\gamma$/$\gamma '$ microstructures offer great promise as candidates for next-generation high-temperature alloys for applications, such as turbine blades. It is essential to understand the thermodynamic and kinetic factors that influence the microstructural evolution of these alloys in order to optimize the alloy compositions and processing steps with a goal to improve their coarsening, creep and rafting behavior. We are using a continuum phase field approach to study the diffusion process and to predict the equilibrium shapes of Co-Al-W $\gamma '$ precipitates. In order to obtain quantitatively predictive capabilities, we extract chemical free energies for the $\gamma$/$\gamma '$ phases based on CALculation of PHAse Diagrams (CALPHAD) thermodynamic data and diffusion mobilities for Co alloys based on CALPHAD kinetic data. We also use experimental or first-principles data for other quantities, such as misfit strain and interface information, for the parameterization of our model. A particular focus of our study is to understand how different energy balances, misfit strain and kinetics affect the coarsening and rafting behavior of $\gamma '$ precipitates, and the sensitivity of the final precipitate shape to materials parameters. Here, we find that the equilibrium shape of the precipitate results from a delicate competition between chemical, interfacial, and elastic energies, and it is very sensitive to changes in model parameters. Here, we examine how modeling input parameters affect the equilibrium shape of precipitates and relate these parameters to experimentally available values.

36 MATERIALS SCIENCE↗

Star formation rate and stellar mass calibrations based on infrared photometry and their dependence on stellar population age and extinction

The stellar mass (M $\star$ ) and the star formation rate (SFR) are among the most important features that characterize galaxies. Measuring these fundamental properties accurately is critical for understanding the present state of galaxies, their history, and future evolution. Infrared (IR) photometry is widely used to measure the M $\star$ and SFR of galaxies because the near-IR traces the continuum emission of the majority of their stellar populations (SPs), and the mid/far-IR traces the dust emission powered by star-forming activity. This work explores the dependence of the IR emission of galaxies on their extinction, and the age of their SPs. It aims to provide accurate and precise IR-photometry SFR and M $\star$ calibrations that account for SP age and extinction while providing quantification of their scatter. We used the CIGALE spectral energy distribution (SED) fitting code to create model SEDs of galaxies with a wide range of star formation histories, dust content, and interstellar medium properties. We fit the relations between M $\star$ and SFR with IR and optical photometry of the model-galaxy SEDs with the Markov chain Monte Carlo (MCMC) method. As an independent confirmation of the MCMC fitting method, we performed a machine-learning random forest (RF) analysis on the same data set. The RF model yields similar results to the MCMC fits, thus validating the latter. This work provides calibrations for the SFR using a combination of the WISE bands 1 and 3, or the JWST NIR-F200W and MIRI-F2100W. It also provides mass-to-light ratio calibrations based on the WISE band-1, the JWST NIR-F200W, and the optical u - r or g - r colors. These calibrations account for the biases attributed to the SP age, while they are given in the form of extinction-dependent and extinction-independent relations. The proposed calibrations show robust estimations while minimizing the scatter and biases throughout a wide range of SFRs and stellar masses. The SFR calibration offers better results, especially in dust-free or passive galaxies where the contributions of old SPs or biases from the lack of dust are significant. Similarly, the M $\star$ calibration yields significantly better results for dusty and high-SFR galaxies where dust emission can otherwise bias the estimations.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum Criticality of Antiferromagnetism and Superconductivity with Relativity

We study a quantum phase transition from a massless to massive Dirac fermion phase in a new two-dimensional bipartite lattice model of electrons that is amenable to sign-free quantum Monte Carlo simulations. Importantly, interactions in our model are not only invariant under SU(2) symmetries of spin and charge like the Hubbard model, but they also preserve an Ising-like electron spin-charge flip symmetry. From unbiased fermion bag Monte Carlo simulations with up to 2304 sites, we show that the massive fermion phase spontaneously breaks this Ising symmetry, picking either antiferromagnetism or superconductivity, and that the transition at which both orders are simultaneously quantum critical belongs to a new “chiral spin-charge symmetric” universality class. We explain our observations using effective potential and renormalization group calculations within the framework of a continuum field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Monte Carlo technique to model performance of streak camera-based time-resolving x-ray spectrometers

A Monte Carlo technique has been developed to simulate the expected signal and the statistical noise of x-ray spectrometers that use streak cameras to achieve the time resolution required for ultrafast diagnostics of laser-generated plasmas. The technique accounts for statistics from both the photons incident on the streak camera’s photocathode and the electrons emitted by the photocathode travelling through the camera’s electron optics to the sensor. We use the technique to optimize the design of a spectrometer, which deduces the temporal history of electron temperature of the hotspot in an inertial confinement fusion implosion from its hard x-ray continuum emission spectra. The technique is general enough to be applied to any instrument using an x-ray streak camera.

Stoupin, S. (ORCID:0000000226225270)↗

The Coronal Line Spectrum of the Iron Nova V3666 Ophiuchi

Optical and infrared spectroscopy are presented for the classical nova V3666 Ophiuchi (Nova Ophiuchi 2018 No.3) from 11 months after its outburst. At that epoch, V3666 Oph displayed an [Fe x] λ0.6375 emission line stronger than Hα and nearly a factor of 3 larger, relative to Hα, than in any nova observed previously. The [Fe x] feature, with an equivalent width exceeding 1300 Å, was enormously strong relative to the continuum emission as well. V3666 Oph also reached a state of high excitation for an Fe ii type nova, displaying coronal lines of [Fe xi], [S ix], and [Ar x], and permitted features that can be identified with the C vi (8–7), (9–8), and (10–9) transitions. In addition to [Fe x] and [Fe xi], V3666 Oph also displayed bright lines of [Fe vii], indicating that the prominent iron features were not confined to a single or narrow range of ionization states. Modeling the emission line spectrum of V3666 Oph with the photoionization code Cloudy indicated that the bulk of the emission line gas was matter-bounded, severely reducing the amount of low ionization gas and resulting in coronal lines that were very strong in comparison to lines such as H i and He i, and the familiar low excitation forbidden lines of nitrogen and oxygen. Such matter-bounded regions are observationally rare in novae since they are both faint and short-lived, offering a plausible explanation for why [Fe x] λ0.6375 had not be seen at this strength previously.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Field-tuned quantum renormalization of spin dynamics in the honeycomb lattice Heisenberg antiferromagnet YbCl 3

The basis for our understanding of quantum magnetism has been the study of elegantly simple model systems. However, even for the antiferromagnetic honeycomb lattice with isotropic spin interactions–one of the simplest model systems–a detailed understanding of quantum effects is still lacking. Here, using inelastic neutron scattering measurements of the honeycomb lattice material YbCl 3 , we elucidate how quantum effects renormalize the single-magnon and multimagnon excitations and how this renormalization can be tuned and ultimately driven to the classical limit by applying a magnetic field. Additionally, our work reveals that the quantum effects tuned by the magnetic field not only renormalize the magnetic excitations but also induce a distinctive sharp feature inside the multimagnon continuum. From a more general perspective, this result demonstrates that structures within magnetic continua can occur over a wide experimental parameter space and can be used as a reliable means of identifying quantum phenomena.

36 MATERIALS SCIENCE↗

GLUE Code: A framework handling communication and interfaces between scales

Many scientific applications are inherently multiscale in nature. Such complex physical phenomena often require simultaneous execution and coordination of simulations spanning multiple time and length scales. This is possible by combining expensive small-scale simulations (such as molecular dynamics simulations) with larger scale simulations (such continuum limit/hydro solvers) to allow for considerably larger systems using task and data parallelism. However, the granularity of the tasks can be very large and often leads to load imbalance. Traditionally, we use approximations to streamline the computation of the more costly interactions and this introduces trade-offs between simulation cost and accuracy. In recent years, the available computational power and the advances in machine learning have made computing these scale-bridging interactions and multiscale simulations more feasible. One driving application has been plasma modeling in inertial confinement fusion (ICF), which is fundamentally multiscale in nature. This requires deep understanding of how to extrapolate microscopic information into macroscopically relevant scales. For example, in ICF one needs an accurate understanding of the connection between experimental observables and the underlying microphysics. The properties of the larger scales are often affected by the microscale behavior incorporated usually into the equations of state and ionic and electronic transport coefficients (Liboff, 1959; Rinderknecht et al., 2014; Rosenberg et al., 2015; Ross et al., 2017). Instead of incorporating this information using reliable molecular dynamics (MD) simulations, one often needs to use theoretical models, due to the inability of MD to reach engineering scales (Glosli et al., 2007; Marinak et al., 1998). One approach to resolve this issue is by coupling two MD simulations of different scales via force interpolation, e.g., the AdResS method (Krekeler et al., 2018; Nagarajan et al., 2013). Another approach, which we will pursue in the scope of this work, is by enabling scale bridging between MD simulations and meso/macro-scale models through the development and support of application programming interfaces that these different applications can interact through.

54 ENVIRONMENTAL SCIENCES↗

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗

An FFT-based micromechanical model for gradient enhanced brittle fracture

Damage models incorporated within FFT-based micromechanical methods have received much attention recently because of the need to better understand and predict brittle and ductile fracture. An important aspect of a damage model is non-local regularization, which removes the mesh dependence of the predictions that otherwise become physically unacceptable upon grid refinement. In this work, the Helmholtz-type equation for non-local gradient regularization of a damage model on a distorted grid is solved using an FFT-based approach. Further, the resulting system of equations is solved using the Jacobi iterative method. The model is applied to simulate brittle fracture of an intermetallic. The influence of the time and space discretization, the length-scale parameter, and intermetallic crystallographic orientation on crack evolution is studied.

36 MATERIALS SCIENCE↗

Continuum of magnetic excitations in the Kitaev honeycomb iridate D 3 LiIr 2 O 6

Inelastic neutron scattering (INS) measurements of powder D 3 ( 7 Li)( 193 Ir) 2 O 6 reveal low energy magnetic excitations with a scattering cross-section that is broad in ∣Q∣ and energy transfer. The magnetic nature of the excitation spectrum is demonstrated by longitudinally polarized neutron scattering. The total magnetic moment of 1.8(4)μ B /Ir inferred from the observed magnetic scattering cross-section is consistent with the effective moment inferred from magnetic susceptibility data and expectations for the J eff = 1/2 single ion state. The rise in the dynamic correlation function $\mathcal{S}(Q,w)$ for ℏω < 5 meV can be described by a simple model assuming nearest-neighbor anisotropic spin exchange, such as that found in the Kitaev model. Exchange disorder associated with the D site likely plays an important role in stabilizing the low T quantum fluctuating state.

Halloran, Thomas [Johns Hopkins Univ., Baltimore, ↗

Gluon moment and parton distribution function of the pion from 𝑁 𝑓 = 2 +1 +1 lattice QCD

We present the first calculation of the pion gluon moment from lattice QCD in the continuum-physical limit. The calculation is done using clover fermions for the valence action with three pion masses, 220, 310 and 690 MeV, and three lattice spacings, 0.09, 0.12, and 0.15 fm, using ensembles generated by MILC Collaboration with 2+1+1 flavors of highly improved staggered quarks (HISQ). On the lattice, we nonperturbatively renormalize the gluon operator in RI/MOM scheme using the cluster-decomposition error reduction (CDER) technique to enhance the signal-to-noise ratio of the renormalization constant. We extrapolate the pion gluon moment to the continuum-physical limit and obtain ⟨𝑥⟩ 𝑔 = 0.394⁢(58) stat+NPR ⁢(39) mixing in the $\overline{MS}$ scheme at 2 GeV, with first error being the statistical error and uncertainties in nonperturbative renormalization, and the second being a systematic uncertainty estimating the effect of ignoring quark mixing. Our pion gluon momentum fraction has a central value lower than two recent single-ensemble lattice-QCD results near physical pion mass but is consistent with the recent global fits by JAM and xFitter and with most QCD-model estimates.

Astronomy & Astrophysics↗

Dependence of Narrow-line Region Sizes on [O III ] Luminosity in Low-redshift Active Galactic Nuclei with Double-peaked Broad Balmer Emission Lines

In this paper, simple but interesting results are reported on the upper limits of narrow-line region (NLR) sizes of a small sample of 38 low-redshift (z < 0.1) active galactic nuclei (AGN) with double-peaked broad emission lines (double-peaked BLAGN), in order to check whether the NLR sizes in type-1 AGN (broad line) and type-2 AGN (narrow line) obey a similar empirical dependence on [O III ] luminosity. In order to correct the inclination effects on projected NLR sizes of type-1 AGN, the accretion disk origin is commonly applied to describe the double-peaked broad Hα line, leading to the determined inclination angles of central disk-like broad-line regions of 38 double-peaked BLAGN. Then, considering the fixed Sloan Digital Sky Survey (SDSS) fiber radius, the upper limits of the NLR sizes of the 38 double-peaked BLAGN can be estimated. Meanwhile, a strong linear correlation between continuum luminosity and [O III ] luminosity is applied to confirm that the [O III ] emissions of the 38 double-peaked BLAGN are totally covered in the SDSS fibers. Considering the reddening-corrected measured [O III ] luminosity, the upper limits of the NLR sizes of the 38 double-peaked BLAGN are within a 99.9999% confidence interval of the expected results from the empirical relation between NLR size and [O III ] luminosity in type-2 AGN. In the current understanding, there are no challenges to the unified model of AGN through the space properties of NLRs.

79 ASTRONOMY AND ASTROPHYSICS↗

Structure Factors of Neutron Matter at Finite Temperature

Here, we compute continuum and infinite volume limit extrapolations of the structure factors of neutron matter at finite temperature and density. Using a lattice formulation of leading-order pionless effective field theory, we compute the momentum dependence of the structure factors at finite temperature and at densities beyond the reach of the virial expansion. The Tan contact parameter is computed and the result agrees with the high momentum tail of the vector structure factor. All errors, statistical and systematic, are controlled for. This calculation is a first step towards a model-independent understanding of the linear response of neutron matter at finite temperature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A field–particle correlation analysis of a perpendicular magnetized collisionless shock

Using the field–particle correlation technique, we examine the particle energization in a three-dimensional (one spatial dimension and two velocity dimensions; 1D-2V) continuum Vlasov–Maxwell simulation of a perpendicular magnetized collisionless shock. The combination of the field–particle correlation technique with the high-fidelity representation of the particle distribution function provided by a direct discretization of the Vlasov equation allows us to ascertain the details of the exchange of energy between the electromagnetic fields and the particles in phase space. We identify the velocity-space signatures of shock-drift acceleration of the ions and adiabatic heating of the electrons arising from the perpendicular collisionless shock by constructing a simplified model with the minimum ingredients necessary to produce the observed energization signatures in the self-consistent Vlasov–Maxwell simulation. We are thus able to completely characterize the energy transfer in the perpendicular collisionless shock considered here and provide predictions for the application of the field–particle correlation technique to spacecraft measurements of collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Positive correlations between soil C stabilization and loss are regulated by roots, minerals, and microbes

Modern conceptual models of soil organic carbon (SOC) cycling place a strong focus on the microbe-mineral interactions that regulate C stabilization. However, the formation of ‘stable’ soil organic matter, which consists mainly of microbial residues bonded to mineral surfaces, is inextricably linked to C loss through microbial respiration. Therefore, what is the net impact of microbial metabolism on the total quantity of C held in the soil? To address this question, we constructed artificial root-soil systems to identify controls on C cycling across the plant-microbe-mineral continuum, simultaneously quantifying the formation of persistent mineral-associated C and C losses to respiration. Roots and minerals interacted to regulate these processes: while roots stimulated respiratory C losses and depleted mineral-associated C pools in low-activity clays, roots triggered formation of stable C in high-activity clays. Moreover, we observed a strong positive correlation between the formation of mineral-associated C and respiration. This suggests that the growth of persistent C pools which can resist subsequent disturbance comes at the expense of C loss from the system.

58 GEOSCIENCES↗