Biomass Particle Models with Realistic Morphology and Resolved Microstructure for Simulations of Intraparticle Transport Phenomena
Abstract not provided
09 BIOMASS FUELS↗
SEARCH · Engineering Papers
Results for “modelled”
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.
Model-to-model interface for multiscale materials modeling
A low-level model-to-model interface is presented that will enable independent models to be linked into an integrated system of models. The interface is based on a standard set of functions that contain appropriate export and import schemas that enable models to be linked with no changes to the models themselves. These ideas are presented in the context of a specific multiscale material problem that couples atomistic-based molecular dynamics calculations to continuum calculations of fluid ow. These simulations will be used to examine the influence of interactions of the fluid with an adjacent solid on the fluid ow. The interface will also be examined by adding it to an already existing modeling code, Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) and comparing it with our own molecular dynamics code.
74 ATOMIC AND MOLECULAR PHYSICS↗
Multiscale Modeling of Structurally-Graded Materials Using Discrete Dislocation Plasticity Models and Continuum Crystal Plasticity Models
A multiscale modeling methodology that combines the predictive capability of discrete dislocation plasticity and the computational efficiency of continuum crystal plasticity is developed. Single crystal configurations of different grain sizes modeled with periodic boundary conditions are analyzed using discrete dislocation plasticity (DD) to obtain grain size-dependent stress-strain predictions. These relationships are mapped into crystal plasticity parameters to develop a multiscale DD/CP model for continuum level simulations. A polycrystal model of a structurally-graded microstructure is developed, analyzed and used as a benchmark for comparison between the multiscale DD/CP model and the DD predictions. The multiscale DD/CP model follows the DD predictions closely up to an initial peak stress and then follows a strain hardening path that is parallel but somewhat offset from the DD predictions. The difference is believed to be from a combination of the strain rate in the DD simulation and the inability of the DD/CP model to represent non-monotonic material response.
Saether, Erik↗
A New BRDF Hotspot Model and its Use in Improving Radiative Transfer Model Efficiency and Application to VLIDORT-based PCRTM Model
The term “hotspot” refers to the sharp increase of reflectance occurring when incident (solar) and reflected (viewing) directions almost coincide in the backscatter direction. The accurate simulation of hotspot directional signatures is important for many remote sensing applications. The RossThick-LiSparse-Reciprocal (RTLSR) Bidirectional Reflectance Distribution Function (BRDF) model is widely used in radiative transfer simulations, and the hotspot model mostly used is from Maignan- Bréon but it typically requires large values of numerical quadrature and Fourier expansion terms in order to represent the hotspot accurately. To improve its use in atmospheric radiative transfer (RT) model simulations, we have developed a modified version based on the Maignan-Bréon’s hotspot BRDF model that converge much faster numerically, making it more practical for use in RT models that require Fourier expansion of BRDF to simulate the top-of-atmosphere (TOA) hotspot signatures. On basis of the vector linearized discrete ordinate radiative transfer model (VLIDORT), we have built a line-by-line based simulation system that can simulate the TOA radiance for atmospheres with 26 gases, different aerosols, dust and clouds, and surface with different BRDF. Using simulated results under diverse atmospheric conditions, we have developed VLIDORT-PCRTM model.
Xiaozhen Xiong↗
The DOE E3SM Model Version 2: Overview of the Physical Model and Initial Model Evaluation
Abstract This work documents version two of the Department of Energy's Energy Exascale Earth System Model (E3SM). E3SMv2 is a significant evolution from its predecessor E3SMv1, resulting in a model that is nearly twice as fast and with a simulated climate that is improved in many metrics. We describe the physical climate model in its lower horizontal resolution configuration consisting of 110 km atmosphere, 165 km land, 0.5° river routing model, and an ocean and sea ice with mesh spacing varying between 60 km in the mid‐latitudes and 30 km at the equator and poles. The model performance is evaluated with Coupled Model Intercomparison Project Phase 6 Diagnosis, Evaluation, and Characterization of Klima simulations augmented with historical simulations as well as simulations to evaluate impacts of different forcing agents. The simulated climate has many realistic features of the climate system, with notable improvements in clouds and precipitation compared to E3SMv1. E3SMv1 suffered from an excessively high equilibrium climate sensitivity (ECS) of 5.3 K. In E3SMv2, ECS is reduced to 4.0 K which is now within the plausible range based on a recent World Climate Research Program assessment. However, a number of important biases remain including a weak Atlantic Meridional Overturning Circulation, deficiencies in the characteristics and spectral distribution of tropical atmospheric variability, and a significant underestimation of the observed warming in the second half of the historical period. An analysis of single‐forcing simulations indicates that correcting the historical temperature bias would require a substantial reduction in the magnitude of the aerosol‐related forcing.
54 ENVIRONMENTAL SCIENCES↗
Solution to the sign problem in a frustrated quantum impurity model
In this work we solve the sign problem of a frustrated quantum impurity model consisting of three quantum spin-half chains interacting through an anti-ferromagnetic Heisenberg interaction at one end. We first map the model into a repulsive Hubbard model of spin-half fermions hopping on three independent one dimensional chains that interact through a triangular hopping at one end. We then convert the fermion model into an inhomogeneous one dimensional model and express the partition function as a weighted sum over fermion worldline configurations. By imposing a pairing of fermion worldlines in half the space we show that all negative weight configurations can be eliminated. This pairing naturally leads to the original frustrated quantum spin model at half filling and thus solves its sign problem.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Final Report: Modeling coupled ice sheet-ocean interactions in the Model for Prediction Across Scales (MPAS) and in DOE Earth System Models
Ice sheet-ocean interactions are involved in a potential climate tipping point in the climate system. Weakening of dynamical ocean barriers that currently hold warmer waters off the Antarctic continental shelf could lead to warming under Antarctica’s largest ice shelves, causing as much as a hundredfold increase in sub-ice-shelf melting. Increased ice-sheet melting could trigger dynamical instabilities in the ice sheet itself if sufficient buttressing from ice shelves is lost. This project developed and tested several modeling capabilities needed to investigate these phenomena. In collaboration with a number of international groups, notably the climate modeling group at Los Alamos National Laboratory (LANL), the following capabilities were added to the Energy Exascale Earth System Model (E3SM) and its ocean component (the Model for Prediction Across Scales Ocean, MPAS-Ocean): 1) a new vertical coordinate appropriate for ice-shelf cavities; 2) an initialization procedure appropriate for damping waves from imbalances in ocean pressure and that of the overlying ice shelf; 3) A suite of idealized test cases related to ice shelf-ocean interactions; 4) A set of diagnostics for tracking freshwater fluxes and heat budgets from ice-shelf melting; 5) infrastructure for computing ice sheet-ocean fluxes in the E3SM coupler; and 6) A set of analysis and visualization tools for E3SM ocean and sea-ice validation with a particular focus on Antarctic and Southern Ocean climate. Additionally, the project included participation and leadership in a number of international model intercomparison projects (MIPs). Most notably, the PI was a co-chair of the Marine Ice Sheet-Ocean Model Intercomparison Project (MISOMIP), designing the experiment, coordinating participation in the MIP and performing the bulk of the analysis. Although the project was cut short when the PI accepted a staff scientist position at LANL, the project played a role in seven publications and has segued smoothly into ongoing Department of Energy projects in the “ecosystem” surrounding E3SM.
54 ENVIRONMENTAL SCIENCES↗
Hierarchical calibration and validation for modeling bench-scale solvent-based carbon capture. Part 1: Non-reactive physical mass transfer across the wetted wall column
A hierarchical model calibration and validation are proposed here for quantifying the confidence level of mass transfer prediction using a computational fluid dynamics (CFD) model, where the solvent-based carbon dioxide (CO 2 ) capture is simulated and the predicted results are compared to the corresponding bench-scale experimental data. Two unit problems with an increasing level of complexity are proposed to break down the complex physical/chemical processes of solvent-based CO 2 capture into relatively simpler problems to separate the effects of physical transport and chemical reaction. This paper focuses on the calibration and validation of the first unit problem, i.e., the CO 2 mass transfer across a falling ethanolamine (MEA) film in the absence of chemical reaction. This problem is investigated both experimentally and numerically using non-reactive nitrous oxide (N 2 O) as a surrogate for CO 2 . To capture motion of the gas-liquid interface, a volume of fluid method is employed with a one-fluid formulation to compute the mass transfer between the two phases. Parallel bench-scale experiments are designed and conducted to validate and calibrate the CFD models using a general Bayesian calibration approach. Two important transport parameters, Henry's constant and gas diffusivity, are calibrated to produce the posterior distributions, which will be used as input for the second unit problem to address the chemical absorption of CO 2 across the MEA falling film, where both mass transfer and chemical reaction are involved. Finally, mass transfer coefficients predicted by CFD are also compared with those predicted by the traditional/empirical correlations under both steady and wavy falling film conditions.
54 ENVIRONMENTAL SCIENCES↗
Modelling of edge localised modes and edge localised mode control
Edge Localised Modes (ELMs) in ITER Q = 10 H-mode plasmas are likely to lead to large transient heat loads to the divertor. To avoid an ELM induced reduction of the divertor lifetime, the large ELM energy losses need to be controlled. In ITER, ELM control is foreseen using magnetic field perturbations created by in-vessel coils and the injection of small D2 pellets. ITER plasmas are characterised by low collisionality at a high density (high fraction of the Greenwald density limit). These parameters cannot simultaneously be achieved in current experiments. Therefore, the extrapolation of the ELM properties and the requirements for ELM control in ITER relies on the development of validated physics models and numerical simulations. In this paper, we describe the modelling of ELMs and ELM control methods in ITER. The aim of this paper is not a complete review on the subject of ELM and ELM control modelling but rather to describe the current status and discuss open issues.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
CFD Modeling and Image Analysis of Exhaled Aerosols due to a Growing Bronchial Tumor: towards Non-Invasive Diagnosis and Treatment of Respiratory Obstructive Diseases
Diagnosis and prognosis of tumorigenesis are generally performed with CT, PET, or biopsy. Such methods are accurate, but have the limitations of high cost and posing additional health risks to patients. In this study, we introduce an alternative computer aided diagnostic tool that can locate malignant sites caused by tumorigenesis in a non-invasive and low-cost way. Our hypothesis is that exhaled aerosol distribution is unique to lung structure and is sensitive to airway structure vari-ations. With appropriate approaches, it is possible to locate the disease site, determine the disease severity, and subsequently formulate a targeted drug delivery plan to treat the disease. This study numerically evaluated the feasibility of the proposed breath test in an image-based lung model with varying pathological stages of a bronchial squamous tumor. Large eddy simulations and a Lagran-gian tracking approach were used to model respiratory airflows and aerosol dynamics. Respira-tions of tracer aerosols of 1 µm at a flow rate of 20 L/min were simulated, with the distributions of exhaled aerosols recorded on a filter at the mouth exit. Aerosol patterns were quantified with multiple analytical techniques such as concentration disparity, spatial scanning and fractal analysis. We demonstrated that a growing bronchial tumor induced notable variations in both the airflow and exhaled aerosol distribution. These variations became more apparent with increasing tumor severity. The exhaled aerosols exhibited distinctive pattern parameters such as spatial probability, fractal dimension, and multifractal spectrum. Results of this study show that morphometric measures of the exhaled aerosol pattern can be used to detect and monitor the pathological states of respiratory diseases in the upper airway. The proposed breath test also has the potential to locate the site of the disease, which is critical in developing a personalized, site-specific drug de-livery protocol.
Aerosol breath test↗
Model-based estimation of wind fields over the ocean from wind scatterometer measurements. I - Development of the wind field model. II - Model parameter estimation
Techniques for the determination of near-surface mesoscale ocean wind fields on the basis of satellite scatterometer data are developed and demonstrated. The derivation of normal-boundary and parameterized-boundary-condition (PBC) wind-field models is outlined, and results from a simulation performed to estimate the model errors are presented in tables. It is shown that the PBC model provides accurate results while minimizing the number of unknowns. After a review of the principles of scatterometry and an analysis of scatterometer measurement noise, an objective function for the measurement parameters is developed and optimized on the basis of gradient search with initial values computed from pointwise wind estimates. The model is then applied to data from a simulation of the NASA Scatterometer (Li et al., 1984), and the results are presented in extensive graphs. The feasibility of model-based wind-field estimation and the appropriateness of the PBC model are demonstrated.
Long, David G.↗
Process-Based Climate Model Development Harnessing Machine Learning: II. Model Calibration From Single Column to Global
Explore the source record for details and available documents.
54 ENVIRONMENTAL SCIENCES↗
From minimum-viable-products to full models: a step-wise development of diagnostic forward models in support of design, analysis and modelling on the ST40 tokamak
Like most magnetic confined fusion experiments, the ST40 tokamak started off with a small subset of diagnostics and gradually increased the diagnostic set to include more complex and comprehensive systems. To make the most of each operational phase, forward models of various diagnostics are used and developed to aid design, provide consistency-checks during commissioning, test analysis methods, and build workflows to constrain high-level parameters to inform interpretation, theory and modelling. For new models and new analysis workflows, minimum-viable-products are released early, and their complexity is increased in a step-wise manner, facilitating the support of all programme phases on multiple parallel applications, while enabling learning opportunities and feedback loops. In this contribution we review the philosophy, scope and architecture of the framework under development. We discuss the details of some forward models, with examples on how they are used to aid diagnostic design, to investigate analysis methodologies through synthetic data, and how they are embedded in experimental analysis workflows. We compare previously published experimental results with new, more advanced analysis workflows employing more recent, detailed models and new diagnostic data, providing confirmation of the published material from the 2021–22 experimental campaign.
integrated data analysis↗
The interaction of reactants, intermediates and products with Cu ions in Cu-SSZ-13 NH3 SCR catalysts: An energetic and ab initio X-ray absorption modeling study
Copper-exchanged SSZ-13 (Cu-SSZ-13) provides excellent catalytic activity and hydrothermal stability in the selective catalytic reduction (SCR) of NOx by using NH3 as a reductant. In this contribution, we compare the most likely positions for Cu with a single charge compensating Al atom (ZCu), including the effect of the adsorption of reactants, intermediates and products of a NH3 SCR reaction by using first-principles calculations based on density functional theory. We find that the 6-membered ring (6MR) site is the most energetically favorable while the 8-membered ring (8MR) sites are less favorable with energy differences about 0.5 eV with respect to the 6MR site. Upon molecular adsorption, the energy differences between Cu in the 8MR and 6MR sites decreases and almost disappears, which means that the thermodynamic stability of ZCu with molecularly adsorbed species is nearly identical for Cu in different sites. Considering more complicated scenarios of NO (or CO) adsorption, we find that the co-adsorption of 2 NO (or 2 CO) molecules weakens their interaction with Cu, as well as their co-adsorption with OH and H2O. The X-ray absorption near edge structure (XANES) of Cu in Cu-SSZ-13 under different conditions was also modeled from first principles. We find that the pre-edge feature in the K-edge XANES of Cu in clean Cu-SSZ-13 with Cu in a 8MR site is induced by a 1s to 4s transition while there is no pre-edge peak when Cu is in a 6MR site. Molecular adsorption onto Cu in a 6MR or a 8MR site results in the same K-edge XANES. When NO is adsorbed onto Cu in SSZ-13, we find that the origin of the small pre-edge peak in the K-edge XANES lies in a 1s to 3d electron transition. On the other hand, when CO or N2 is adsorbed onto an isolated Cu atom in SSZ-13, a transition from the 1s state to the hybridized state of 3d and 4p results in a pre-edge peak in the K-edge of the XANES spectrum. It is also found that the edge position and intensity of the maximum peak of Cu K-edge XANES depend on the Cu locations and the adsorbed molecules. By studying the Cu K-edge XANES of ZCu with adsorption of species of different oxidizing power, we conclude that a higher oxidation state of Cu will result in a higher energy edge position in the Cu K-edge XANES. These results are expected to greatly aid our understanding of Cu-SSZ-13 SCR catalysts. This work was supported by institutional funds provided to JSM from the Voiland School of Chemical Engineering and Bioengineering. Financial support was also provided by the National Science Foundation GOALI program under contract Nos. CBET-1258717, CBET-1258715 and CBET-1258690. We thank Prof. Fabio Ribeiro, Mr. Atish Parekh, Prof. W. F. Schneider, Mr. Christopher Paolucci, Mr. Trunjoyo Anggara and Prof. Jeff Miller for stimulating discussions on the modeling of the XANES spectrum and Ms. Alyssa Hensley for her comments on the manuscript. J.Sz. acknowledges the financial support of this work by the U.S. Department of Energy (DOE), Office of Energy Efficiency and Renewable Energy, Vehicle Technologies Program. A portion of the computer time for the computational work was performed using EMSL, a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research and located at PNNL. PNNL is a multi-program national laboratory operated for the US DOE by Battelle.
Zhang, Renqin↗
Modeling and control of fuel distribution in a dual-fuel internal combustion engine leveraging late intake valve closings
In internal combustion engines, cycle-to-cycle and cylinder-to-cylinder variations of the combustion process have been shown to negatively impact the fuel efficiency of the engine and lead to higher exhaust emissions. The combustion variations are generally tied to differences in the composition and condition of the trapped mass throughout each cycle and across individual cylinders. Thus, advanced engines featuring exhaust gas recirculation, flexible valve actuation systems, advanced fueling strategies, and turbocharging systems are prone to exhibit higher variations in the combustion process. In this study, the cylinder-to-cylinder variations of the combustion process in a dual-fuel internal combustion engine leveraging late intake valve closing are investigated and a model to predict and address one of the root causes for these variations across cylinders is developed. The study is conducted on an inline six-cylinder heavy-duty dual-fuel engine equipped with exhaust gas recirculation, a variable geometry turbocharger, and a fully flexible variable intake valve actuation system. The engine is operated with late intake valve closure timings in a dual-fuel combustion mode in which a high reactivity fuel is directly injected into the cylinders and a low reactivity fuel is port injected into the cylinders. The cylinder-to-cylinder variations observed in the study have been associated with the maldistribution of the port-injected fuel, which is exacerbated at late intake valve timings. The resulting difference in indicated mean effective pressure between the cylinders ranges from 9% at an intake valve closing of 570° after top dead center to 38% at an intake valve closing of 620° after top dead center and indicates an increasingly uneven fuel distribution. The study leverages both experimental and simulation studies to investigate the distribution of the port-injected fuel and its impact on cylinder-to-cylinder variation. The effects of intake valve closing as well as the impact of intake runner length on fuel distribution were quantitatively analyzed, and a model was developed that can be used to accurately predict the fuel distribution of the port-injected fuel at different operating conditions with an average estimation error of 1.5% in cylinder-specific fuel flow. A model-based control strategy is implemented to adjust the fueling at each port and shown to significantly reduce the cylinder-to-cylinder variations in fuel distribution.
Kassa, Mateos↗
A numerical modeling study of a Montana thunderstorm. I - Model results versus observations involving nonelectrical aspects. II - Model results versus observations involving electrical aspects
Model results and the observed cloud behavior are examined in terms of nonelectrical and electrical aspects of a thunderstorm. The characteristics of the two-dimensional, time-dependent atmospheric electricity model used to simulate the cloud observations of July 19, 1981 in Miles City, Montana are described. The interactions of the dynamics and microphysics of the cloud with the charge separation mechanisms are analyzed. It is observed that the model accurately represents many of the observed characteristics of the cloud; however, the cloud base height, maximum liquid water content, and the time from first formation of precipitation until it reaches the ground are not accurately modeled. It is found that the model adequately represents the electrical field structure of the cloud and the electrical field strengths.
Helsdon, John H., Jr.↗
Wind-Tunnel Investigation of the Effects on the Aerodynamic Characteristics of Modifications to a Model of a Bomb Mounted on a Wing-Fuselage Model and to a Model of the Bomb Alone
An investigation was conducted in the Langley high-speed 7- by 10-foot tunnel to determine effects of modifications to a bomb model (particularly with regard to drag) when mounted on a wing-fuselage model and tested at Mach numbers from 0.70 to 1.10. In addition, the static longitudinal stability characteristics of several configurations of a larger scale model of the bomb alone were obtained over a Mach number range from 0.50 to 0.95. The results obtained for the wing-fuselage-bomb model indicate that large reductions in installation drag were obtained for the wing-fuselage-bomb model when the flat nose of the basic bomb was replaced by rounded or pointed noses of various calibers. Shortening the mounting pylon gave further decreases in the installation drag. The tests of the bomb alone indicated that only the flat-nose configurations were stable over the greater part of the Mach number range. Nose-shape modifications which improved the drag also caused the bombs to become unstable at low angles of attack. The stability of the low-drag bomb configurations could be improved by lengthening the cylindrical portion of the body behind the center of gravity.
King, Thomas J., Jr.↗
Electroweak baryogenesis in the exceptional supersymmetric standard model
We study electroweak baryogenesis in the E{sub 6} inspired exceptional supersymmetric standard model (E{sub 6}SSM). The relaxation coefficients driven by singlinos and the new gaugino as well as the transport equation of the Higgs supermultiplet number density in the E{sub 6}SSM are calculated. Our numerical simulation shows that both CP-violating source terms from singlinos and the new gaugino can solely give rise to a correct baryon asymmetry of the Universe via the electroweak baryogenesis mechanism.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗