Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics 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 1,423 records · Page 79

A high-fidelity multiphysics model for the new solid oxide iron-air redox battery part I: Bridging mass transport and charge transfer with redox cycle kinetics

A high-fidelity two-dimensional axial symmetrical multi-physics model is described in this paper as an effort to simulate the cycle performance of a recently discovered solid oxide metal-air redox battery (SOMARB). The model collectively considers mass transport, charge transfer and chemical redox cycle kinetics occurring across the components of the battery, and is validated by experimental data obtained from independent research. In particular, the redox kinetics at the energy storage unit is well represented by Johnson-Mehl-Avrami-Kolmogorov (JIVIAK) and Shrinking Core models. The results explicitly show that the reduction of Fe3O4 during the charging cycle limits the overall performance. Distributions of electrode potential, overpotential, Nernst potential, and H-2/H2O-concentration across various components of the battery are also systematically investigated. (C) 2015 Elsevier B.V. All rights reserved.

Jin, XF↗

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗

Hidden Sector Dark Matter Models for the Galactic Center Gamma-Ray Excess

The gamma-ray excess observed from the Galactic Center can be interpreted as dark matter particles annihilating into Standard Model fermions with a cross section near that expected for a thermal relic. Although many particle physics models have been shown to be able to account for this signal, the fact that this particle has not yet been observed in direct detection experiments somewhat restricts the nature of its interactions. One way to suppress the dark matter's elastic scattering cross section with nuclei is to consider models in which the dark matter is part of a hidden sector. In such models, the dark matter can annihilate into other hidden sector particles, which then decay into Standard Model fermions through a small degree of mixing with the photon, Z, or Higgs bosons. After discussing the gamma-ray signal from hidden sector dark matter in general terms, we consider two concrete realizations: a hidden photon model in which the dark matter annihilates into a pair of vector gauge bosons that decay through kinetic mixing with the photon, and a scenario within the generalized NMSSM in which the dark matter is a singlino-like neutralino that annihilates into a pair of singlet Higgs bosons, which decay through their mixing with the Higgs bosons of the MSSM.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Near-vacuum hohlraums for driving fusion implosions with high density carbon ablatorsa)

Recent experiments at the National Ignition Facility [M. J. Edwards et al., Phys. Plasmas 20, 070501 (2013)] have explored driving high-density carbon ablators with near-vacuum hohlraums, which use a minimal amount of helium gas fill. These hohlraums show improved efficiency relative to conventional gas-filled hohlraums in terms of minimal backscatter, minimal generation of suprathermal electrons, and increased hohlraum-capsule coupling. Given these advantages, near-vacuum hohlraums are a promising choice for pursuing high neutron yield implosions. Long pulse symmetry control, though, remains a challenge, as the hohlraum volume fills with material. Here, two mitigation methodologies have been explored, dynamic beam phasing and increased case-to-capsule ratio (larger hohlraum size relative to capsule). Unexpectedly, experiments have demonstrated that the inner laser beam propagation is better than predicted by nominal simulations, and an enhanced beam propagation model is required to match measured hot spot symmetry. Ongoing work is focused on developing a physical model which captures this enhanced propagation and on utilizing the enhanced propagation to drive longer laser pulses than originally predicted in order to reach alpha-heating dominated neutron yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An improved empirical model of electron and ion fluxes at geosynchronous orbit based on upstream solar wind conditions

In this study, a new empirical model of the electron fluxes and ion fluxes at geosynchronous orbit (GEO) is introduced, based on observations by Los Alamos National Laboratory (LANL) satellites. The model provides flux predictions in the energy range ~1 eV to ~40 keV, as a function of local time, energy, and the strength of the solar wind electric field (the negative product of the solar wind speed and the z component of the magnetic field). Given appropriate upstream solar wind measurements, the model provides a forecast of the fluxes at GEO with a ~1 h lead time. Model predictions are tested against in-sample observations from LANL satellites and also against out-of-sample observations from the Compact Environmental Anomaly Sensor II detector on the AMC-12 satellite. The model does not reproduce all structure seen in the observations. However, for the intervals studied here (quiet and storm times) the normalized root-mean-square deviation < ~0.3. It is intended that the model will improve forecasting of the spacecraft environment at GEO and also provide improved boundary/input conditions for physical models of the magnetosphere.

58 GEOSCIENCES↗

ANALYSIS OF VARIANCE-BASED MIXED MULTISCALE FINITE ELEMENT METHOD AND APPLICATIONS IN STOCHASTIC TWO-PHASE FLOWS

The stochastic partial differential systems have been widely used to model physical processes, where the inputs involve large uncertainties. Flows in random and heterogeneous porous media is one of the cases where the random inputs (e.g., permeability) are often modeled as a stochastic field with high-dimensional random parameters. To treat the high dimensionality and heterogeneity efficiently, model reduction is employed in both stochastic space and physical space. An analysis of variance (ANOVA)-based mixed multiscale finite element method MsFEM) is developed to decompose the high-dimensional stochastic problem into a set of lower-dimensional stochastic subproblems, which require much less computational complexity and significantly reduce the computational cost in stochastic space, and the mixed MsFEM can capture the heterogeneities on a coarse grid to greatly reduce the computational cost in the spatial domain. In addition, to enhance the efficiency of the traditional ANOVA method, an adaptive ANOVA method based on a new adaptive criterion is developed, where the most active dimensions can be selected to greatly reduce the computational cost before conducting ANOVA decomposition. This novel adaptive criterion is based on variance-decomposition method coupled with sparse-grid probabilistic collocation method or multilevel Monte Carlo method. The advantage of this adaptive criterion lies in its much lower computational overhead for identifying the active dimensions and interactions. A number of numerical examples in two-phase stochastic flows are presented and demonstrate the accuracy and performance of the adaptive ANOVA-based mixed Ms-FEM.

Wei, Jia↗

Machine Learning Tools Set for Natural Gas Fuel Cell System Design

This study is focusing on leveraging the system design tools set for the next-generation solid oxide fuel cell (SOFC) based natural gas fuel cell (NGFC) system. Conventionally, system design and optimization of NGFC systems rely heavily on traditional reduced order model (ROM) techniques and designers’ experience level. For overcoming the technical barriers of system design, multiple multi-physics models and machine learning (ML) tools have been utilized to automate the conceptual design process and enhance the reliability of solutions for the NGFC system. The proposed tools set includes a physics-informed ML tool for automated ROM construction that leverages advances in deep neural networks to significantly reduce ROM prediction error for the NGFC power island compared to traditional approaches. The constructed physics-informed ML ROM can be used in system design, and optimization tools set Institute for the Design of Advanced Energy Systems (IDAES) Process Systems Engineering (PSE) framework. The tools set also provides a user-friendly graphic user interface built within Jupyter Notebooks, and the complete tools set is open-source public available.

Wang, Dewei↗

Measurement and Modeling of Electron Temperature in Laboratory Photoionized Plasmas Relevant to Astrophysics

The electron temperature of photoionized plasmas characterizes the thermalization of photoelectrons, impacts the charge-state distribution, emissivity, and opacity through atomic recombination processes, and is needed to perform detailed comparisons with theory predictions. We discuss temperature measurements in laboratory photoionized plasmas and a comparison with model calculations done with several theory approximations and codes. These include a radiation-hydrodynamics simulation and two nonequilibrium heating and ionization models that tracked the evolution of the internal energy of the electrons. Furthermore, for the same physics model and X-ray flux time history, calculations were performed assuming steady-state or time-dependent conditions. The time history of steady-state results correlates with that of the X-ray flux, while that of the time-dependent cases does not, and it is qualitatively and quantitatively different from the steady-state case. Steady-state results significantly overestimated temperature measurements, while time-dependent results produced better approximations, which suggests the importance of transient effects in the experiment and also the need for time-resolved measurements.

Laboratory astrophysics↗

Measurement and modeling of proton-induced reactions on arsenic from 35 to 200 MeV

72 As is a promising positron emitter for diagnostic imaging that can be employed locally using a 72 Se generator. However, current reaction pathways to 72 Se have insufficient nuclear data for efficient production using regional 100–200 MeV high-intensity proton accelerators. In order to address this deficiency, stacked-target irradiations were performed at LBNL, LANL, and BNL to measure the production of the 72 Se/ 72 As positron emission tomography (PET) generator system via 75 As (p, x) between 35 and 200 MeV. This work provides the most well-characterized excitation function for 75 As ( p , 4 n ) 72 Se starting from threshold. Additional focus was given to report the first measurements of 75 As (p, x) 68 Ge and bolster an already robust production capability for the highly valuable 68 Ge/ 68 Ga PET generator. Thick target yield comparisons with prior established formation routes to both generators are made. In total, high-energy proton-induced cross sections are reported for 55 measured residual products from 75 As, nat Cu , and nat Ti targets, where the latter two materials were present as monitor foils. These results were compared with literature data as well as the default theoretical calculations of the nuclear model codes TALYS, COH, EMPIRE, and ALICE. Reaction modeling at these energies is typically unsatisfactory due to few prior published data and many interacting physics models. Therefore, a detailed assessment of the talys code was performed with simultaneous parameter adjustments applied according to a standardized procedure. Particular attention was paid to the formulation of the two-component exciton model in the transition between the compound and preequilibrium regions, with a linked investigation of level density models for nuclei off of stability and their impact on modeling predictive power. This paper merges experimental work and evaluation techniques for high-energy charged-particle isotope production in an extension to an earlier study of this kind.

43 PARTICLE ACCELERATORS↗

An Improved Microwave Semiempirical Model for the Dielectric Behavior of Moist Soils

Soil semiempirical dielectric models (SEMs) are powerful, and they are generally considered a useful hybrid of both empirical and physical models. In this work, the Wang-Schmugge dielectric model is improved to more accurately estimate the relative complex dielectric constants (CDCs) of moist soils. Instead of the Debye relaxation spectrum of liquid water located outside of the soil (i.e., free out-of-soil water) adopted in the Wang-Schmugge model, the Debye relaxation formula related to the free-water component inside the soil [i.e., free soil water (FSW)], which is correlated with the soil texture, is employed in the improved SEM. In addition, the effective conductivity loss term related to both soil texture and soil moisture is introduced to explain the ionic conductivity losses of FSW. Since the soil moisture influence is reduced at high frequencies, the effective conductivity loss term related to only the soil texture is also analyzed for 14-18 GHz. As in the Wang-Schmugge model, the relative CDC of bound soil water varies with the soil volumetric moisture content when the soil moisture is lower than the maximum bound water fraction in the new model, which takes a different approach than the Mironov mineralogy-based SEM. The proposed model obtains better fitting results than the three most widely employed SEMs. Lastly, the improved model exhibits a significantly improved accuracy with a higher correlation coefficient (R²), a closer 1:1 relationship, and a lower root-mean-square error, including in the L-band, and especially in the imaginary part of the L-band.

semiempirical dielectric model (SEM)↗

SYS-5620-Project - Achieving Robust Laser Performance utilizing Historical Shot Experiments - A Systems Engineering Proposal

The National Ignition Facility (NIF) at Lawrence Livermore National Laboratory precisely guides, amplifies, reflects, and focuses 192 powerful laser beams into a target about the size of a pencil eraser in a few billionths of a second, delivering more than 2 million joules of ultraviolet energy and 500 trillion watts of peak power. A crucial goal of the system is to trigger precise implosions of fuel capsules. This is achieved by delivering all 192 beams at user-specified times and locations on the target, minimizing any deviation from the requested performance. Power requirements can vary substantially on each experiment and the facility supports numerous amplifier pumping configurations and their attendant nonlinear effects. To achieve the tight performance required across such a broad array of configurations, constant comparison of measured and requested power delivery are tracked and long-term trends analyzed as a guide to understanding future performance. A very common question asked is given the current state of the laser today – how well would a similar experiment from the past perform today? Additionally, if one were to specify an alternate amplifier configuration using today’s model would and damage limits be exceeded and would there be an increase in performance. The physics model used to make these predictions and equipment protection checks is called the Virtual Beamline (VBL). VBL is used with an incoming desired pulse shape and energy to be delivered on target, and then does an iterative solve to predict the needed injected pulse in the front-end of the system to achieve this result. To effectively guide and predict future NIF experiment performance, laser scientists explore current amplifier configurations and compare them with historical data utilizing a tool called the Reverify Toolbox.

42 ENGINEERING↗

Rotating-fluid experiments with an atmospheric general circulation model

In order to determine features of rotating fluid flow that are dependent on the geometry, rotating annulus-type experiments are carried out with a numerical model in spherical coordinates. Rather than constructing and testing a model expressly for this purpose, it is found expedient to modify an existing general circulation model of the atmosphere by removing the model physics and replacing the lower boundary with a uniform surface. A regime diagram derived from these model experiments is presented; its major features are interpreted and contrasted with the major features of rotating annulus regime diagrams. Within the wave regime, a narrow region is found where one or two zonal wave numbers are dominant. The results reveal no upper symmetric regime; wave activity at low rotation rates is thought to be maintained by barotropic rather than baroclinic processes.

Geisler, J. E.↗

Emissions Prediction and Measurement for Liquid-Fueled TVC Combustor with and without Water Injection

An investigation is performed to evaluate the performance of a computational fluid dynamics (CFD) tool for the prediction of the reacting flow in a liquid-fueled combustor that uses water injection for control of pollutant emissions. The experiment consists of a multisector, liquid-fueled combustor rig operated at different inlet pressures and temperatures, and over a range of fuel/air and water/fuel ratios. Fuel can be injected directly into the main combustion airstream and into the cavities. Test rig performance is characterized by combustor exit quantities such as temperature and emissions measurements using rakes and overall pressure drop from upstream plenum to combustor exit. Visualization of the flame is performed using gray scale and color still photographs and high-frame-rate videos. CFD simulations are performed utilizing a methodology that includes computer-aided design (CAD) solid modeling of the geometry, parallel processing over networked computers, and graphical and quantitative post-processing. Physical models include liquid fuel droplet dynamics and evaporation, with combustion modeled using a hybrid finite-rate chemistry model developed for Jet-A fuel. CFD and experimental results are compared for cases with cavity-only fueling, while numerical studies of cavity and main fueling was also performed. Predicted and measured trends in combustor exit temperature, CO and NOx are in general agreement at the different water/fuel loading rates, although quantitative differences exist between the predictions and measurements.

Brankovic, A.↗

The FLUKA code for space applications: recent developments

The FLUKA Monte Carlo transport code is widely used for fundamental research, radioprotection and dosimetry, hybrid nuclear energy system and cosmic ray calculations. The validity of its physical models has been benchmarked against a variety of experimental data over a wide range of energies, ranging from accelerator data to cosmic ray showers in the earth atmosphere. The code is presently undergoing several developments in order to better fit the needs of space applications. The generation of particle spectra according to up-to-date cosmic ray data as well as the effect of the solar and geomagnetic modulation have been implemented and already successfully applied to a variety of problems. The implementation of suitable models for heavy ion nuclear interactions has reached an operational stage. At medium/high energy FLUKA is using the DPMJET model. The major task of incorporating heavy ion interactions from a few GeV/n down to the threshold for inelastic collisions is also progressing and promising results have been obtained using a modified version of the RQMD-2.4 code. This interim solution is now fully operational, while waiting for the development of new models based on the FLUKA hadron-nucleus interaction code, a newly developed QMD code, and the implementation of the Boltzmann master equation theory for low energy ion interactions. c2004 COSPAR. Published by Elsevier Ltd. All rights reserved.

Non-NASA Center↗

Object-Oriented/Data-Oriented Design of a Direct Simulation Monte Carlo Algorithm

Over the past decade, there has been much progress towards improved phenomenological modeling and algorithmic updates for the direct simulation Monte Carlo (DSMC) method, which provides a probabilistic physical simulation of gas Rows. These improvements have largely been based on the work of the originator of the DSMC method, Graeme Bird. Of primary importance are improved chemistry, internal energy, and physics modeling and a reduction in time to solution. These allow for an expanded range of possible solutions In altitude and velocity space. NASA's current production code, the DSMC Analysis Code (DAC), is well-established and based on Bird's 1994 algorithms written in Fortran 77 and has proven difficult to upgrade. A new DSMC code is being developed in the C++ programming language using object-oriented and data-oriented design paradigms to facilitate the inclusion of the recent improvements and future development activities. The development efforts on the new code, the Multiphysics Algorithm with Particles (MAP), are described, and performance comparisons are made with DAC.

Liechty, Derek S.↗

Gas-Particle Interaction Model Development in Plume Surface Interaction Erosion and Cratering

The development of a predictive simulation capability for extra-terrestrial Plume Surface Interaction (PSI) environments is undertaken by the Fluid Dynamics Branch at NASA/MSFC under the Game Changing Development (GCD) Program funded by NASA’s Space Technology Mission Directorate (STMD). Predictive simulation capability for propulsive Martian and Lunar landing PSI requires accurate modeling of the complex high-speed plume impingement and resulting gas/particle-cloud and ejecta stream formation. Physics modeling gap analysis during formulation of the PSI project scope identified two particular phenomena of first order importance in gas-particle interactions that lacked existing modeling formulations. The first arises from the lack of models capturing the unsteady drag forces imparted on particles by the rapidly expanding supersonic compressible flow. The second addresses the fluctuating forces and dispersions imparted on both the particle and the fluid resulting from the interference and wake turbulence generated by close proximity particle clouds, dubbed the Particle Turbulent Kinetic Energy(PTKE). Their first order significance has been identified in experiments, but simulation models currently do not exist for either effect. The development of models and the eventual implementation into the Eulerian Gas-Granular Flow Solver (GGFS) simulation tools was constructed as one of four major tasks of the PSI project. In this process, data on particle kinetics and gas-particle interactions are collected from carefully designed experiments of particles embedded in jets. The effects observed in the experiments are then replicated in high-fidelity particle-resolved CFD simulations to inform the formulation of improvements to particle phase drag models for implementations in the more efficient Eulerian-Lagrangian CFD simulations. The resulting models are ultimately ported to the Eulerian-Eulerian models applied for most efficient simulations in PSI production application tools. This paper describes the activities and accomplishments of the past year in the gas-particle interaction modeling task of the PSI project.

Jeff West↗

The SM and NLO Multileg Working Group: Summary Report

After years of waiting, and after six Les Houches workshops, the era of LHC running is finally upon us, albeit at a lower initial center-of-mass energy than originally planned. Thus, there has been a great sense of anticipation from both the experimental and theoretical communities. The last two years, in particular, have seen great productivity in the area of multi-parton calculations at leading order (LO), next-to-leading order (NLO) and next-to-next-to-leading order (NNLO), and this productivity is reflected in the proceedings of the NLM group. Both religions, Feynmanians and Unitarians, as well as agnostic experimenters, were well-represented in both the discussions at Les Houches, and in the contributions to the write-up. Next-to-leading order (NLO) is the first order at which the normalization, and in some cases the shape, of perturbative cross sections can be considered reliable. This can be especially true when probing extreme kinematic regions, as for example with boosted Higgs searches considered in several of the contributions to this writeup. A full understanding for both standard model and beyond the standard model physics at the LHC requires the development of fast, reliable programs for the calculation of multi-parton final states at NLO. There have been many advances in the development of NLO techniques, standardization and automation for such processes and this is reflected in the contributions to the first section of this writeup. Many calculations have previously been performed with the aid of semi-numerical techniques. Such techniques, although retaining the desired accuracy, lead to codes which are slow to run. Advances in the calculation of compact analytic expressions for Higgs + 2 jets have resulted in the development of much faster codes, which extend the phenomenology that can be conducted, as well as making the code available to the public for the first time. A prioritized list of NLO cross sections was assembled at Les Houches in 2005 and added to in 2007. This list includes cross sections which are experimentally important, and which are theoretically feasible (if difficult) to calculate. Basically all 2-3 cross sections of interest have been calculated, with the frontier now extending to 2 {yields} 4 calculations. Often these calculations exist only as private codes. Since 2007, two additional calculations have been completed: t{bar t}b{bar b} and W+3 jets, reflecting the advance of the NLO technology to 2 {yields} 4 processes. In addition, the cross section for b{bar b}b{bar b} has been calculated for the q{bar q} initial state with the gg initial state calculation in progress. Final states of such complexity usually lead to multi-scale problems, and the correct choice of scales to use can be problematic not only at LO, but also at NLO. The size of the higher order corrections and of the residual scale dependence at NLOcan depend strongly on whether the considered cross section is inclusive, or whether a jet veto cut has been applied. Depending on the process, dramatically different behavior can be observed upon the application of a jet veto. There is a trade-off between suppressing the NLO cross section and increasing the perturbative uncertainty, with application of a jet veto sometimes destroying the cancellation between infra-red logs of real and virtual origin, and sometimes just suppressing large (and very scale-sensitive) tree-level contributions. So far, there is no general rule predicting the type of behavior to be expected, but this is an important matter for further investigation. From the experimental side, an addition to the above wish-list that will be crucial is the determination of the accuracy to which each of the calculations needs to be known. This is clearly related to the experimental accuracy at which the cross sections can be measured at the LHC, and can determine, for example, for what processes it may be necessary to calculate electo-weak corrections, in addition to the higher order QCD corrections. On the theoretical side, it would also be interesting to categorize the impact of a jet veto on the size and stability of each of the NLO cross sections. The technology does exist to carry out a calculation for W/Z production at NNLO (QCD) and at NLO (EW). This process was placed on the wish-list in 2007 and it is unfortunate that the combined calculation has not yet been carried out, as this precision benchmark will be very useful and important at the LHC.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗