Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phenomenological Models”

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

Progress Report on Model Development for the Transport of Aerosol through Microchannels

This report summarizes the current progress in the development of a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose is to introduce a generic, reliable numerical model for the prediction of aerosol transport, deposition, and plugging in leak paths while accounting for potential plugging formation, to a user community involving researchers, regulators, and industry. In that regard, a Graphical User Interface (GUI) was generated by integrating the individual MATLAB scripts that make up the model and adding additional features to aid the general understanding of the problem to a user. This report focuses on the model development and features included. Further, predictions from the GUI are compared with experimental data for validation. The strength of the GUI is the ability of a user to plug in basic parameters such as the initial pressure conditions and canister model specifications in order to obtain a first principal approximation of vital information such as the blowdown pressure differential, aerosol penetration, and deposition as a function of time. The user can obtain this without the know-how of the underlying MATLAB scripts and thus enables the model to be readily applied by regulators, industry, and shareholders to reduce the uncertainty in the off-site radiological consequences. The report also lists future tasks under this work scope. This includes ongoing efforts to include additional crack geometries (divergent and divergent-convergent slots) to the model; GUI development and improvement based on feedback from users, and integration of aerosol source term data from Sibling Pin tests into the model as we approach realistic canister stress corrosion cracking–induced aerosol release scenarios.

54 ENVIRONMENTAL SCIENCES↗

Baryon–antibaryon photoproduction cross sections off the proton

The GlueX experiment at Jefferson Lab has observed 𝑝⁢$\overline{p}$ and, for the first time, Λ$\overlineΛ$ and 𝑝$\overlineΛ$ photoproduction from a proton target at photon energies up to 11.6 GeV. The angular distributions are forward peaked for all produced pairs, consistent with Regge-like 𝑡-channel exchange. Asymmetric wide-angle antibaryon distributions show the presence of additional processes. In a phenomenological model, we find consistency with a double-𝑡-channel exchange process where antibaryons are created only at the middle vertex. The model matches all observed distributions with a small number of free parameters. In the hyperon channels, we observe a clear distinction between photoproduction of the Λ$\overlineΛ$ and 𝑝$\overlineΛ$ systems but general similarity to the 𝑝$\overline{p}$ system. We report both total cross sections and cross sections differential with respect to momentum transfer and the invariant masses of the created particle pairs. No narrow resonant structures were found in these reaction channels. The suppression of 𝑠$\overline{s}$ quark pairs relative to 𝑑$\overline{d}$ quark pairs is similar to what has been seen in other reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coupled Microbial-Conversion and Computational-Fluid-Dynamics (CFD) Models for Butanediol Production in Micro-Aerated Reactors

Microbial conversion of substrates to macromolecules has been widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries. These bioreactions are also being investigated in the production of low-value commodities such as biofuels [1]. Gas-liquid mass-transfer and transport-reaction coupling are important challenges when designing and scaling up these reactor systems. Experiments in wellmixed small-scale reactors have enabled characterization of microbial reactivity, while their coupling with macroscale transport remains relatively unexplored. In this work, we use a coupled metabolic-CFD model to study the action of a genetically engineered microbe Zymomonas mobilis [2] on sugars to produce 2,3-Butanediol (BDO). BDO is an important hydrocarbon intermediate that can be catalytically upgraded to several fuels and chemicals [3]. An important aspect to this particular microbial conversion is the need for micro-aerated environments as opposed to traditional aerobic fermentation. Slight variations in oxygen concentration can result in competing reaction pathways that disable BDO production. Hence, gas-liquid mass transfer and transport need to be optimized in large-scale reactors to maximize BDO production, for which CFD is a valuable tool. The aerobic-fermentation CFD model previously developed by the authors [5] for simulating bubble-column and airlift reactors at scale was used in this study. The Reynolds-averaged mass, momentum and energy transport equations for interpenetrating gas and liquid phase are solved in this model along with the transport and interphase mass transfer of oxygen. Our previous work used a phenomenological model for microbial oxygen uptake that neglected microbial growth and other reaction pathways. In this work, a detailed metabolic model enabled prediction of product formation and inhibition pathways. In order to manage computational cost, we used a subcycling technique [6] that takes advantage of the clear separation in transport (~ 200 sec) and reaction (~ 2-3 hours) timescales. The CFD model is first solved to steady state, after which the metabolic model is advanced at every cell in the computational domain using the local oxygen concentration. The CFD model is then run to achieve a new steady state that provides a new oxygen distribution for the metabolic model. This process, where reaction and fluid updates are interleaved together, is iterated until reactants are completely exhausted. This work will examine the performance of different reactor designs such as bubble column and airlift reactors at scale (250-500 m3). Oxygen mass-transfer coefficient and distribution are critically analyzed among reactors, and optimization studies pertaining to aeration is presented. Furthermore, it has been observed in experiments that high BDO production may be achieved by manipulating the aerobic environment over the course of reaction, such that oxygen concentration is high during the growth phase, and very low as sugar is depleted. This characteristic will be addressed by our simulations for which a timedependent scheduling strategy for aeration is presented that maximizes BDO production. [1] Humbird, D., Davis, R., and McMillan, J., Aeration costs in stirred-tank and bubble column bioreactors, Biochemical Engineering Journal, 127, 161—166, 2017 [2] Yang, S., Mohagheghi, A., Franden, M. A., Chou, Y.-C., Chen, X., Dowe, N., Himmel, M. E., and Zhang, M., Metabolic engineering of zymomonas mobilis for 2, 3-butanediol production from lignocellulosic biomass sugars. Biotechnology for biofuels, 9(1):189, 2016 [3] Kim, S. J., Sim, H. J., Kim, J. W., Lee, Y. G., Park, Y. C., and Seo, J. H., Enhanced production of 2,3-butanediol from xylose by combinatorial engineering of xylose metabolic pathway and cofactor regeneration in pyruvate decarboxylase-deficient Saccharomyces cerevisiae. Bioresource Technology, 245:1551–1557, 2017 [4] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

A Model to Account for the Effects of Load Ratio and Hydrogen Pressure on the Fatigue Crack Growth Behavior of Pressure Vessel Steels

A phenomenological model for estimating the effects of load ratio R and hydrogen pressure P H2 on the hydrogen-assisted fatigue crack growth rate (HA–FCGR) behavior in the transient and steady-state regimes of pressure vessel steels is described. The “transient regime” is identified with crack growth within a severely embrittled zone of intense plasticity at the crack tip. The “steady-state” behavior is associated with the crack growing into a region of comparatively lower hydrogen concentration located further away from the crack tip. The model treats the effects of R and P H2 as being functionally separable. In the transient regime, the effects of the hydrogen pressure on the HA–FCGR behavior were negligible but were significant in the steady-state regime. The hydrogen concentration in the steady-state region is modeled as being dependent on the kinetics of lattice diffusion, which is sensitive to pressure. Experimental HA–FCGR data from the literature were used to validate the model. The new model was shown to be valid over a wide range of conditions that ranged between –1 ≤ R ≤ 0.8 and 0.02 ≤ P H2 ≤ 103 MPa for pressure vessel steels.

36 MATERIALS SCIENCE↗

Neutron elastic scattering on calcium isotopes from chiral nuclear optical potentials

We formulate microscopic neutron-nucleus optical potentials from many-body perturbation theory based on chiral two- and three-body forces. The neutron self-energy is first calculated in homogeneous matter to second order in perturbation theory, which gives the central real and imaginary terms of the optical potential. Here, the real spin-orbit term is calculated separately from the density matrix expansion using the same chiral interaction as in the self-energy. Finally, the full neutron-nucleus optical potential is derived within the improved local density approximation utilizing mean-field models consistent with the chiral nuclear force employed. We compare the results of the microscopic calculations to phenomenological models and experimental data up to projectile energies of E = 200 MeV. Experimental elastic differential scattering cross sections and vector analyzing powers are generally well reproduced by the chiral optical potential, but we find that total cross sections are overestimated at high energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Adaptive four-level modeling of laser cooling of solids

Laser cooling of rare-earth doped solids has been demonstrated across a wide range of material platforms, inspiring the development of simple phenomenological models such as the four-level model to elucidate the universal properties of laser cooling under various operating conditions. However, these models usually require the input of full absorption spectra that must be provided experimentally or by additional complicated atomic modeling. In this Letter, we propose that a four-level model, when extended to admit effective energy levels adaptive to the pumping photon energy, can accurately predict the cooling efficiency as a function of temperature and pumping frequency using only few inputs such as the absorption coefficient measured at a single frequency and temperature. Our model exploits the quasi-equilibrium properties of the excitation of rare-earth ions for the determination of the effective four energy levels. The model is validated against published experimental results for a span of materials including ytterbium/thulium-doped glass and crystals. With the verified model, we derive explicit expressions for the optimal frequency and the operating bandwidth of pumping lasers. Our model significantly simplifies the modeling process of laser cooling and is expected to stimulate further development of optical refrigeration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spectral line-shape in collinear laser spectroscopy after atomic charge exchange

Collinear laser spectroscopy experiments on fast, neutral beams have been extensively used for studies on short-lived radioactive nuclei, taking advantage of its high sensitivity. The resulting resonance line-shape is known to show significant distortion, due to the energy exchange during the charge-exchange neutralization process, which can cause large systematic uncertainty in the determined centroid. A model for the line shape was constructed and simulated to be compared to measured Al, Si, and Ni hyperfine spectra. It is shown that the distortion is caused mainly by the transfer of electron into many different energy levels in the projectile atom and subsequent decays, rather than secondary inelastic collisions, which were often assumed in the line shape analysis before. Furthermore, the model can also be applied to other projectile–alkali pairs, providing a reliable line-shape with less fitting parameters than conventional phenomenological models.

Atom↗

Geometry-driven modeling of electron localization in InAs/GaAs double quantum dots

The coupled electronic states in two-dimensional (2D) and three-dimensional (3D) double quantum dot (DQD) systems are investigated using a phenomenological model applied to InAs/GaAs heterostructures. The single-band k • p effective potential approach previously proposed by our group is employed to numerically calculate the energy spectrum and spatial localization of a single electron, serving as an indicator of the coupling strength within the binary system. For identical quantum dots (QDs) in a DQD, the electronic states exhibit ideal coherence. We systematically vary the DQD geometry and the strength of the confinement potential (via an applied electric field) to examine the effects of symmetry breaking and the sensitivity of electron localization in both identical and nearly identical DQDs. Our results show that coherence in DQDs is highly sensitive to these subtle variations. This sensitivity can be harnessed to detect changes in the surrounding environment, such as fluctuations in chemical or electrical properties that affect the DQD system.

electron localization↗

Differential cross-sections for events with missing transverse momentum and jets measured with the ATLAS detector in 13 TeV proton-proton collisions

Measurements of inclusive, differential cross-sections for the production of events with missing transverse momentum in association with jets in proton-proton collisions at √s = 13 TeV are presented. The measurements are made with the ATLAS detector using an integrated luminosity of 140 fb –1 and include measurements of dijet distributions in a region in which vector-boson fusion processes are enhanced. They are unfolded to correct for detector resolution and efficiency within the fiducial acceptance, and are designed to allow robust comparisons with a wide range of theoretical predictions. A measurement of differential cross sections for the Z → νν process is made. The measurements are generally well-described by Standard Model predictions except for the dijet invariant mass distribution. Auxiliary measurements of the hadronic system recoiling against isolated leptons, and photons, are also made in the same phase space. Ratios between the measured distributions are then derived, to take advantage of cancellations in modelling effects and some of the major systematic uncertainties. These measurements are sensitive to new phenomena, and provide a mechanism to easily set constraints on phenomenological models. To illustrate the robustness of the approach, these ratios are compared with two common Dark Matter models, where the constraints derived from the measurement are comparable to those set by dedicated detector-level searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A transient site balance model for atomic layer etching

We present a transient site balance model of plasma-assisted atomic layer etching of silicon (Si) with alternating exposure to chlorine gas (Cl 2 ) and argon ions (Ar + ). Molecular dynamics (MD) simulation results are used to provide parameters for the model. The model couples the dynamics of a top monolayer surface region ('top layer') and a perfectly mixed subsurface region ('mixed layer'). The differential equations describing the rates of change of the Cl coverage in the two layers are transient mass balances. Model predictions include Cl coverages and rates of etching of various species from the surface as a function of Cl 2 or Ar + fluence. The simplified phenomenological model reproduces the MD simulation results well over a range of conditions. Comparing model predictions directly to experimental optical emission spectroscopy data, as reported in a previous paper, provides further evidence of the accuracy of the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spatiotemporal mapping of mesoscopic liquid dynamics

The study of liquid dynamics at mesoscopic scales is still strewn with difficulty due to limitations in theory and experiment. Historically, significant attention has been given to the analysis of space-time correlation functions and their frequency-Fourier transforms at a few discrete wave numbers. The massive computing power afforded by modern high performance computing clusters and the advent of a wide-angle neutron spin-echo spectrometer, however, have unlocked a more intuitive and fruitful approach to this problem. Using molecular dynamics simulations, here we demonstrate the benefits of spatiotemporally mapping intermediate scattering functions on a dense grid of correlation times and wave numbers. Four model systems are investigated: a Lennard-Jones liquid, a coarse-grained bead-spring polymer, a molten sodium chloride, and a poly(ethylene oxide) melt. We show that the spatiotemporal mapping approach is particularly useful for elucidating the mesoscopic dynamics in these liquids, where several underlying mechanisms, such as molecular relaxations, hydrodynamic modes, and nonhydrodynamic excitations, are potentially at play. Therefore, compared to the traditional method, direct visualization of density space-time correlation functions on two-dimensional color maps permits appraisals of complicated dynamical behavior at mesoscales in a global manner. For example, the scaling relations between space and time for different types of molecular motions can be straightforwardly identified on these plots, without any model-dependent analysis. Additionally, we show how theoretical ideas regarding collective mesoscopic dynamics, such as the classical hydrodynamic theory, the convolution approximation, and a recently proposed phenomenological model, can be discussed in terms of the global features of spatiotemporal maps of intermediate scattering functions. The new perspective offered by the spatiotemporal mapping method should prove useful for the study of liquid dynamics in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Data-driven modeling of dislocation mobility from atomistics using physics-informed machine learning

Dislocation mobility, which dictates the response of dislocations to an applied stress, is a fundamental property of crystalline materials that governs the evolution of plastic deformation. Traditional approaches for deriving mobility laws rely on phenomenological models of the underlying physics, whose free parameters are in turn fitted to a small number of intuition-driven atomic scale simulations under varying conditions of temperature and stress. This tedious and time-consuming approach becomes particularly cumbersome for materials with complex dependencies on stress, temperature, and local environment, such as body-centered cubic crystals (BCC) metals and alloys. In this paper, we present a novel, uncertainty quantification-driven active learning paradigm for learning dislocation mobility laws from automated high-throughput large-scale molecular dynamics simulations, using Graph Neural Networks (GNN) with a physics-informed architecture. We demonstrate that this Physics-informed Graph Neural Network (PI-GNN) framework captures the underlying physics more accurately compared to existing phenomenological mobility laws in BCC metals.

36 MATERIALS SCIENCE↗

Full-Scale Ab Initio Simulation of Magic-Angle-Spinning Dynamic Nuclear Polarization

Theoretical models aimed at describing magic-angle-spinning (MAS) dynamic nuclear polarization (DNP) NMR have great potential in facilitating the in silico design of DNP polarizing agents and formulations. These models must typically face a trade-off between the accuracy of a strict quantum mechanical description and the need for using realistically large spin systems, for instance, using phenomenological models. In this paper, we show that the use of aggressive state-space restrictions and an optimization strategy allows full-scale ab initio MAS-DNP simulations of spin systems containing thousands of nuclei. Our simulations are shown to reproduce experimental DNP enhancements quantitatively, including their MAS rate dependence, for both frozen solutions and solid materials. They also reveal the importance of a previously unrecognized structural feature found in some polarizing agents that helps minimize the sensitivity losses imposed by the spin diffusion barrier.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Excluded-volume model for quarkyonic matter: Three-flavor baryon-quark mixture

The single-flavor excluded-volume model based on the effective size of baryons reproduces the hard-soft density evolution of the equation of state (EoS) required by the recent studies of GW170817. This phenomenological model basically realizes the concept of quarkyonic matter which is introduced from large-$N_c$ gauge theory for dense matter. Further, enhanced nucleon interactions and dynamically generated quark degrees of freedom can reproduce the hard-soft evolution of the EoS. In this paper, we extend the excluded-volume model to a three-flavor system by considering electromagnetic charge and possible weak equilibrium in order to obtain a proper description for the hard-soft behavior of the EoS inferred from the gravitational waves observations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling iterative reconstruction and displacement field in the large scale structure

The next generation of galaxy surveys like the Dark Energy Spectroscopic Instrument and Euclid will provide datasets orders of magnitude larger than anything available to date. Our ability to model nonlinear effects in late time matter perturbations will be a key to unlock the full potential of these datasets, and the area of initial condition reconstruction is attracting growing attention. Iterative reconstruction developed by Marcel Schmittfull et al. is a technique designed to reconstruct the displacement field from the observed galaxy distribution. The nonlinear displacement field and initial linear density field are highly correlated. Therefore, reconstructing the nonlinear displacement field enables us to extract the primordial cosmological information better than from the late time density field at the level of the two-point statistics. This paper will test to what extent the iterative reconstruction can recover the true displacement field and construct a perturbation theory model for the postreconstructed field. We model the iterative reconstruction process with Lagrangian perturbation theory (LPT) up to third order for dark matter in real space and compare it with N-body simulations. Here we find that the simulated iterative reconstruction does not converge to the nonlinear displacement field, and the discrepancy mainly appears in the shift term, i.e., the term correlated directly with the linear density field. On the contrary, our 3LPT model predicts that the iterative reconstruction should converge to the nonlinear displacement field. We discuss the sources of discrepancy, including numerical noise/artifacts on small scales, and present an ad hoc phenomenological model that improves the agreement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tunable plasmon-enhanced second-order optical nonlinearity in transition-metal dichalcogenide nanotriangles (Final Report for SC0008712)

The control of light-matter interaction at the nanoscale is a grand challenge that cuts across modern photonics, plasmonics, and optoelectronics. The development of simulation tools for predictive modeling of nanostructured light emitters and absorbers is critical to successfully tackling this grand challenge. Existing state-of-the-art tools are not suitable for the following reasons: a) Phenomenological models of materials, commonly employed in the full-wave solutions to Maxwell’s equations, are fit to experimental data for bulk materials and fall short when attempting to analyze nanostructures. b) Nonlinearities in optical response can arise solely from nanostructuring, and require time-dependent quantum electronic transport techniques to capture. c) Emergent materials with exotic and tunable physical properties are increasingly used, but their properties are very sensitive to the substrate, impurities, and edge termination, and require an accurate description of interband and intraband processes, as well as dissipation. d) Microscopic modeling of electronic transport usually employs simplified electrodynamics, and realistic models that couple full electronic transport with full-wave electrodynamics are needed for accurate photonics simulations. The objective of this project was to computational tools that will overcome the limitations faced by current modeling tools and provide an unprecedented level of accuracy for predictive modeling of light-matter interaction at the nanoscale.

36 MATERIALS SCIENCE↗

Updated parameters of the LArQL model

The need for a microscopic description of scintillation light generation in liquid argon becomes increasingly desirable with the upcoming operation of large scale LArTPCs in the next decade. While a detailed mathematical account of the process is still to be achieved, a phenomenological model for simultaneously treating ionization and scintillation, LArQL, has been successfully employed to describe the range of electric fields from 0 to 0.75 kV/cm and dE/dx from 2 to 40 MeV/cm providing the anti-correlation between the free ionization charge and scintillation light. A reanalysis of the original model parameter values has been performed within a global fit procedure and is presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Far-forward neutrinos at the Large Hadron Collider

We present a new calculation of the energy distribution of high-energy neutrinos from the decay of charm and bottom hadrons produced at the Large Hadron Collider (LHC). In the kinematical region of very forward rapidities, heavy-flavor production and decay is a source of tau neutrinos that leads to thousands of charged-current tau neutrino events in a 1 m long, 1 m radius lead neutrino detector at a distance of 480 m from the interaction region. In our computation, next-to-leading order QCD radiative corrections are accounted for in the production cross-sections. Non-perturbative intrinsic-kT effects are approximated by a simple phenomenological model introducing a Gaussian kT-smearing of the parton distribution functions, which might also mimic perturbative effects due to multiple initial-state soft-gluon emissions. The transition from partonic to hadronic states is described by phenomenological fragmentation functions. To study the effect of various input parameters, theoretical predictions for D± s production are compared with LHCb data on double-differential cross-sections in transverse momentum and rapidity. The uncertainties related to the choice of the input parameter values, ultimately affecting the predictions of the tau neutrino event distributions, are discussed. We consider a 3+1 neutrino mixing scenario to illustrate the potential for a neutrino experiment to constrain the 3+1 parameter space using tau neutrinos and antineutrinos. We find large theoretical uncertainties in the predictions of the neutrino fluxes in the far-forward region. Untangling the effects of tau neutrino oscillations into sterile neutrinos and distinguishing a 3+1 scenario from the standard scenario with three active neutrino flavours, will be challenging due to the large theoretical uncertainties from QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗