Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MODEL”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Heterolytic Cleavage of Hydrogen by an Iron Hydrogenase Model: An Fe-H - - - H-N Dihydorgen Bond Characterized by Neutron Diffraction

Use of hydrogen as a fuel by [FeFe]-hydrogenase enzymes in nature requires heterolytic cleavage of the H-H bond into a proton (H+) and hydride (H-), a reaction that is also a critical step in homogeneous catalysts for hydrogenation of C=O and C=N bonds. An understanding of the catalytic oxidation of H2 by hydrogenases provides insights into the design of synthetic catalysts that are sought as cost-effective alternatives to the use of the precious metal platinum in fuel cells. Crystallographic studies on the [FeFe]-hydrogenase enzyme were critical to understanding of its reactivity, but the key H-H cleavage step is not readily observed experimentally in natural hydrogenases. Synthetic biomimics have provided evidence for H2 cleavage leading to hydride transfer to the metal and proton transfer to an amine. Limitations on the precise location of hydrogen atoms by x-ray diffraction can be overcome by use of neutron diffraction, though its use is severely limited by the difficulty of obtaining suitable crystals and by the scarcity of neutron sources. Here we show that an iron complex with a pendant amine in the diphosphine ligand cleaves hydrogen heterolytically under mild conditions, leading to [CpC5F4NFeH(PtBu2NtBu2H)]+BArF4-, [PtBu2NtBu2 = 1,5-di(tert-butyl)-3,7-di(tert-butyl)-1,5-diaza-3,7-diphosphacyclooctane; ArF = 3,5-bis(trifluoromethyl)phenyl]. The Fe-H- - - H-N moiety has a strong dihydrogen bond, with a remarkably short H • • • H distance of 1.489(10) Å between the protic N-Hδ+ and hydridic Fe-Hδ-. The structural data for [CpC5F4NFeH(PtBu2NtBu2H)]+ provide a glimpse of how the H-H bond is oxidized or generated in hydrogenase enzymes, with the pendant amine playing a key role as a proton relay. The iron complex [CpC5F4NFeH(PtBu2NtBu2H)]+BArF4- is an electrocatalyst for oxidation of H2 (1 atm) at 22 °C, so the structural data are obtained on a complex that is a functional model for catalysis by [FeFe]-hydrogenase enzymes. This research was supported as part of the Center for Molecular Electrocatalysis, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences. Pacific Northwest National Laboratory is operated by Battelle for the U.S. Department of Energy.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Simulation of Heterogeneous Atom Probe Tip Shape Evolution during Field Evaporation Using a Level Set Method and Different Evaporation Models

Accurate reconstruction of spatial positions of field evaporated ions from detector patterns depends on a correct understanding of the tip shape evolution and evaporation laws of component atoms. Reconstruction artifacts introduced in the heterogeneous tip can be contributed by the assumption on homogeneous evaporation behavior. A level set method based specimen shape evolution model is developed in this study to simulate the evaporation of synthetic layered-structured atom probe tip. The simulation results of the shape evolution by the present level set model agree well with the literature data using other numerical approaches, i.e. finite difference method. The asymmetric evolving shape predicted by the level set model demonstrates that the complex evaporation behavior of heterogeneous tips and interface curvature can potentially lead to the tip reconstruction distortion. Compared with other shape evaporation methods the new method provides considerably higher resolution and computational efficiency with the aid of the intrinsic subgrid scale accuracy.

atom probe tomography, field evaporation, tip shap↗
DOE OSTI2016

Addressing numerical challenges in introducing a reactive transport code into a land surface model: A biogeochemical modeling proof-of-concept with CLM-PFLOTRAN 1.0: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic discussion paper available at doi:10.5194/gmd-9-927-2016. The purpose is to document the simulations to allow verification, reproducibility, and follow-up studies. This dataset contains shell scripts to create the CLM-PFLOTRAN cases, specific input files for PFLOTRAN and CLM, outputs, and python scripts to make the figures using the outputs in the publication. Through these results, we demonstrate that CLM-PFLOTRAN can approximately reproduce CLM results in selected cases for the Arctic, temperate and tropic sites. In addition, the new framework facilitates mechanistic representations of soil biogeochemistry processes in the land surface model. This dataset includes four zipped files and one .pdf file.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Diagnosis of the Marine Low Cloud Simulation in the NCAR Community Earth System Model (CESM) and the NCEP Global Forecast System (GFS)-Modular Ocean Model v4 (MOM4) coupled model

We present a diagnostic analysis of the marine low cloud climatology simulated by two state-of-the-art coupled atmosphere-ocean models: the NCAR Community Earth System Model (CESM) and the NCEP Global Forecasting System (GFS). In both models, the shallow convection and boundary layer turbulence parameterizations have been recently updated: both models now use a mass-flux scheme for the parameterization of shallow convection, and a turbulence parameterization capable of handling Stratocumulus (Sc)-topped Planetary Boundary Layers (PBLs). For shallow convection, both models employ a convective trigger function based on the concept of convective inhibition and both include explicit convective overshooting/penetrative entrainment formulation. For Sc-topped PBL, both models treat explicitly turbulence mixing and cloud-top entrainment driven by cloud-top radiative cooling. Our focus is on the climatological transition from Sc to shallow Cumulus (Cu)-topped PBL in the subtropical eastern oceans. We show that in the CESM the coastal Sc-topped PBLs in the subtropical Eastern Pacific are well-simulated but the climatological transition from Sc to shallow Cu is too abrupt and happens too close to the coast. By contrast, in the GFS coupled simulation the coastal Sc amount and PBL depth are severely underestimated while the transition from Sc to shallow Cu is ³delayed² and offshore Sc cover is too extensive in the subtropical Eastern Pacific. We discuss the possible connections between such differences in the simulations and differences in the parameterizations of shallow convection and boundary layer turbulence in the two models.

Xiao, Heng↗
DOE OSTI2021

Drying model of a high salt content cementitious waste form: Effect of capillary forces and salt solution

Highlights: • Drying model for a high salt content cementitious waste form is developed. • Water vapor diffusion and capillary liquid flow are distinguished. • Capillary and salt solution effects are considered in vapor-liquid equilibrium. A water transport model coupling capillary liquid flow with vapor diffusion is developed to describe the drying process for a cementitious waste form with high salinity porewater. Vapor-liquid equilibrium is formulated as the driving force for vapor diffusion and the model accounts for pore capillary and high salinity effects on water thermodynamic activity. Pore filling and porewater surface tension as a function of pore size distribution and water saturation have been quantified for the material. Geochemical speciation modeling is used to simulate porewater activity as a function of composition over the range of saturation. The theoretical relationship between relative humidity and water saturation generally agrees with experimental measurement, and the developed model is capable of predicting drying rates under various external relative humidity conditions. The model was developed to be incorporated into reactive transport models considering the effects of drying such as salt redistribution and efflorescence.

36 MATERIALS SCIENCE↗
DOE OSTI2017

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↗
NASA NTRS2012

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↗
NASA NTRSDate not supplied

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↗
DOE OSTI2022

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↗
DOE OSTI2017

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↗
DOE OSTI2019

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↗
DOE OSTI2017

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↗
DOE OSTI2015

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↗
DOE OSTI2015

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↗
NASA NTRS1990

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.↗
DOE OSTI2024

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↗