Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Semiparametric modeling: Correcting low-dimensional model error in parametric models

In this paper, a semiparametric modeling approach is introduced as a paradigm for addressing model error arising from unresolved physical phenomena. Our approach compensates for model error by learning an auxiliary dynamical model for the unknown parameters. Practically, the proposed approach consists of the following steps. Given a physics-based model and a noisy data set of historical observations, a Bayesian filtering algorithm is used to extract a time-series of the parameter values. Subsequently, the diffusion forecast algorithm is applied to the retrieved time-series in order to construct the auxiliary model for the time evolving parameters. The semiparametric forecasting algorithm consists of integrating the existing physics-based model with an ensemble of parameters sampled from the probability density function of the diffusion forecast. To specify initial conditions for the diffusion forecast, a Bayesian semiparametric filtering method that extends the Kalman-based filtering framework is introduced. In difficult test examples, which introduce chaotically and stochastically evolving hidden parameters into the Lorenz-96 model, we show that our approach can effectively compensate for model error, with forecasting skill comparable to that of the perfect model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Integrated optical and electrical modeling of plasmon-enhanced thin film photovoltaics: A case-study on organic devices

The nanoscale light control for absorption enhancement of organic photovoltaic (OPV) devices inevitably produces strongly non-uniform optical fields. These non-uniformities due to the localized optical modes are a primary route toward absorption enhancement in OPV devices. Therefore, a rigorous modeling tool taking into account the spatial distribution of optical field and carrier generation is necessary. Presented here is a comprehensive numerical model to describe the coupled optical and electrical behavior of plasmon-enhanced polymer:fullerene bulk heterojunction (BHJ) solar cells. In this model, a position-dependent electron-hole pair generation rate that could become highly non-uniform due to photonic nanostructures is directly calculated from the optical simulations. By considering the absorption and plasmonic properties of nanophotonic gratings included in two different popular device architectures, and applying the Poisson, current continuity, and drift/diffusion equations, the model predicts quantum efficiency, short-circuit current density, and desired carrier mobility ratios for bulk heterojunction devices incorporating nanostructures for light management. In particular, the model predicts a significant degradation of device performance when the carrier species with lower mobility are generated far from the collecting electrode. Consequently, an inverted device architecture is preferred for materials with low hole mobility. This is especially true for devices that include plasmonic nanostructures. Additionally, due to the incorporation of a plasmonic nanostructure, we use simulations to theoretically predict absorption band broadening of a BHJ into energies below the band gap, resulting in a 4.8% increase in generated photocurrent.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2020

GISS Model E2.2: A Climate Model Optimized for the Middle Atmosphere. Part 1: Model Structure, Climatology, Variability and Climate Sensitivity

We introduce a new climate model (GISS E2.2) that has been specially optimized for the Middle Atmosphere, and whose output is being contributed to the CMIP6 archive. The top of the model is at a geopotential altitude of 89 km, and parameterizations of moist convection and various forms of gravity wave drag based on tropospheric processes are chosen specifically for this optimization. We first evaluate the model in its configuration as a coupled atmosphere-chemistry model with respect to its simulation of the mean state of the Middle Atmosphere, from the mesosphere down through the Upper Troposphere/Lower Stratosphere. Then we assess its use as a coupled atmosphere-ocean climate model by exploring its mean ocean climatology. To evaluate its variability, we report on its simulation of the primary modes in the troposphere, stratosphere and ocean. Two climate change simulations are presented, the responses to instantaneous increases of 2xCO2 and 4xCO2, run with two different ocean models. Sensitivity studies are performed to illustrate the effect of parameterizations on the model results. We compare these results to the lower vertical resolution/top GISS Model E2.1, whose output has also been submitted to CMIP6. The different choices made for these models are explored. It is shown that important improvements in the circulation above and below the tropopause can be obtained when attention is paid to representation of Middle Atmosphere processes in climate model development.

D. Rind↗
DOE OSTI2017

Pairwise-interaction extended point-particle model for particle-laden flows

In this work we consider the pairwise interaction extended point-particle (PIEP) model for Euler–Lagrange simulations of particle-laden flows. By accounting for the precise location of neighbors the PIEP model goes beyond local particle volume fraction, and distinguishes the influence of upstream, downstream and laterally located neighbors. The two main ingredients of the PIEP model are (i) the undisturbed flow at any particle is evaluated as a superposition of the macroscale flow and a microscale flow that is approximated as a pairwise superposition of perturbation fields induced by each of the neighboring particles, and (ii) the forces and torque on the particle are then calculated from the undisturbed flow using the Faxén form of the force relation. The computational efficiency of the standard Euler–Lagrange approach is retained, since the microscale perturbation fields induced by a neighbor are pre-computed and stored as PIEP maps. Here we extend the PIEP force model of Akiki et al. with a corresponding torque model to systematically include the effect of perturbation fields induced by the neighbors in evaluating the net torque. Also, we use DNS results from a uniform flow over two stationary spheres to further improve the PIEP force and torque models. We then test the PIEP model in three different sedimentation problems and compare the results against corresponding DNS to assess the accuracy of the PIEP model and improvement over the standard point-particle approach. In the case of two sedimenting spheres in a quiescent ambient the PIEP model is shown to capture the drafting-kissing-tumbling process. In cases of 5 and 80 sedimenting spheres a good agreement is obtained between the PIEP simulation and the DNS. For all three simulations, the DEM-PIEP was able to recreate, to a good extent, the results from the DNS, while requiring only a negligible fraction of the numerical resources required by the fully-resolved DNS.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Testing a land model in ecosystem functional space via a comparison of observed and modeled ecosystem flux responses to precipitation regimes and associated stresses in a Central U.S. forest: Test Model in Ecosystem Functional Space

Testing complex land surface models has often proceeded by asking the question: does the model prediction agree with the observation? This approach has yet led to high-performance terrestrial models that meet the challenges of climate and ecological studies. Here we test the Community Land Model (CLM) by asking the question: does the model behave like an ecosystem? We pursue its answer by testing CLM in the ecosystem functional space (EFS) at the Missouri Ozark AmeriFlux (MOFLUX) forest site in the Central U.S., focusing on carbon and water flux responses to precipitation regimes and associated stresses. In the observed EFS, precipitation regimes and associated water and heat stresses controlled seasonal and interannual variations of net ecosystem exchange (NEE) of CO 2 and evapotranspiration in this deciduous forest ecosystem. Such controls were exerted more strongly by precipitation variability than by the total precipitation amount per se. A few simply constructed climate variability indices captured these controls, suggesting a high degree of potential predictability. While the interannual fluctuation in NEE was large, a net carbon sink was maintained even during an extreme drought year. Although CLM predicted seasonal and interanual variations in evapotranspiration reasonably well, its predictions of net carbon uptake were too small across the observed range of climate variability. Also, the model systematically underestimated the sensitivities of NEE and evapotranspiration to climate variability and overestimated the coupling strength between carbon and water fluxes. Its suspected that the modeled and observed trajectories of ecosystem fluxes did not overlap in the EFS and the model did not behave like the ecosystem it attempted to simulate. A definitive conclusion will require comprehensive parameter and structural sensitivity tests in a rigorous mathematical framework. We also suggest that future model improvements should focus on better representation and parameterization of process responses to environmental stresses and on more complete and robust representations of carbon-specific processes so that adequate responses to climate variability and a proper degree of coupling between carbon and water exchanges are captured.

58 GEOSCIENCES↗
DOE OSTI2014

Modeling time-coincident ultrafast electron transfer and solvation processes at molecule-semiconductor interfaces

Kinetic models based on Fermi's Golden Rule are commonly employed to understand photoinduced electron transfer dynamics at molecule-semiconductor interfaces. Implicit in such second-order perturbative descriptions is the assumption that nuclear relaxation of the photoexcited electron donor is fast compared to electron injection into the semiconductor. This approximation breaks down in systems where electron transfer transitions occur on 100-fs time scale. Here, we present a fourth-order perturbative model that captures the interplay between time-coincident electron transfer and nuclear relaxation processes initiated by light absorption. The model consists of a fairly small number of parameters, which can be derived from standard spectroscopic measurements (e.g., linear absorbance, fluorescence) and/or first-principles electronic structure calculations. Insights provided by the model are illustrated for a two-level donor molecule coupled to both (i) a single acceptor level and (ii) a density of states (DOS) calculated for TiO{sub 2} using a first-principles electronic structure theory. These numerical calculations show that second-order kinetic theories fail to capture basic physical effects when the DOS exhibits narrow maxima near the energy of the molecular excited state. Overall, we conclude that the present fourth-order rate formula constitutes a rigorous and intuitive framework for understanding photoinduced electron transfer dynamics that occur on the 100-fs time scale.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Test Results of the CERN HL-LHC low- Quadrupole Short Models MQXFS3c and MQXFS4

For the High Luminosity Upgrade of the CERN Large Hadron Collider, lower * quadrupole magnets based on advanced Nb3Sn conductors will be installed on each side of the ATLAS and CMS experiment insertion zones. As part of the technological developments needed to achieve the required field gradient of 132.6 T/m within a 150-mm aperture, short length model magnets, named MQXFS, are tested both at the CERN SM18 and FERMI Lab test facilities. The model magnets rely on two types of Nb3Sn conductors (RRP and PIT) and on an innovative bladders and keys design to provide mechanical support against the Lorentz forces. In 2016 and 2017, the powering tests of the first two models MQXFS3 (RRP) and MQXFS5 (PIT) proved that nominal performance (16.5 kA) could be reached with excellent memory of the quench current after thermal cycle. However both magnets showed a slow training behavior with clear observations of voltage disturbances before the quench. Besides MQXFS5 only could reach ultimate current (17.9 kA) whereas erratic behavior was observed on MQXFS3 due to conductor local degradation at the head of one of the coil. In 2018, this limiting coil was changed and the applied azimuthal pre-stress increased. If ultimate current could then be reached, no stable current could be maintained due to identified defect on the outer layer of the new coil. Finally the outcome of the test of the new model MQXFS4, featuring the final RRP conductors that will be used for the series production and variation on the inner layer quench heater designs are here reported in details.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Model reduction for hypersonic aerodynamics via conservative LSPG projection and hyper-reduction

High-speed aerospace engineering applications rely heavily on computational fluid dynamics (CFD) models for design and analysis due to the expense and difficulty of flight tests and experiments. This reliance on CFD models necessitates performing accurate and reliable uncertainty quantification (UQ) of the CFD models. However, it is very computationally expensive to run CFD for hypersonic flows due to the fine grid resolution required to capture the strong shocks and large gradients that are typically present. Furthermore, UQ approaches are “many-query” problems requiring many runs with a wide range of input parameters. One way to enable computationally expensive models to be used in such many-query problems is to employ projection-based reduced-order models (ROMs) in lieu of the (high-fidelity) full-order model. In particular, the least-squares Petrov–Galerkin (LSPG) ROM (equipped with hyper-reduction) has demonstrated the ability to significantly reduce simulation costs while retaining high levels of accuracy on a range of problems including subsonic CFD applications. This allows computationally inexpensive LSPG ROM simulations to replace the full-order model simulations in UQ studies, which makes this many-query task tractable, even for large-scale CFD models. This work presents the first application of LSPG to a hypersonic CFD application. In particular, we present results for LSPG ROMs of the HIFiRE-1 in a three-dimensional, turbulent Mach 7.1 flow, showcasing the ability of the ROM to significantly reduce computational costs while maintaining high levels of accuracy in computed quantities of interest.

Fluid Dynamics↗
NASA NTRS1949

Longitudinal Stability and Control Characteristics of a Semispan Wind-Tunnel Model of the XF7U-1 Airplane and a Comparison with Complete-Model Wind-Tunnel Tests and Semispan-Model Wing-Flow Tests

An investigation was conducted on an 0.08-scale semispan model of the Chance Vought XF7U-1 airplane in the Langley high-speed 7- by 10-foot tunnel in the Mach number range from 0.40 to 0.97. The results are compared with those obtained with an 0.08-scale sting-mounted complete model tested in the same tunnel and with an 0.026-scale semispan model tested by the wing-flow method. The lift-curve slopes obtained for the 0.08-scale semispan model and the 0.026-scale wing-flow model were in good agreement but both were generally lower than the values obtained for the sting model. The results of an unpublished investigation have shown that tunnel-wall boundary-layer and strut-leakage effects can came the difference noted between the lift-curve slopes of the sting and the semispan data. Fair agreement was obtained among the data of the three models as regard the variation of pitching-moment coefficients with lift coefficient. The agreement between the complete and the semispan models was more favorable with the vertical fine on, because the wall-boundary-layer and strut leakage effects were less severe. In the Mach number range between 0.94 and 0.97, ailavator-control reversal was indicated in the wing-flow data near zero lift; Whereas, these same trends were indicated in the larger scale semispan data at somewhat higher lift coefficients.

Goodson, Kenneth W.↗
NASA NTRS2024

Modeling Atmospheric Brown Carbon in the GISS ModelE Earth System Model

Brown carbon (BrC) is an absorbing organic aerosol (OA), primarily emitted through biomass burning (BB), which exhibits light absorption unique to both black carbon (BC) and other organic aerosols. Despite many field and laboratory studies seeking to constrain BrC properties, the radiative forcing (RF) of BrC is still highly uncertain. To better understand its climate impact, we introduced BrC to the One-Moment Aerosol (OMA) module of the GISS ModelE Earth system model (ESM). We assessed ModelE sensitivity to primary BrC processed through a novel chemical aging scheme and to secondary BrC formed from biogenic volatile organic compounds (BVOCs). Initial results show that BrC typically contributes a top-of-the-atmosphere (TOA) radiative effect of 0.04 W m −2 . Sensitivity tests indicate that explicitly simulating BrC (separating it from other OAs), including secondary BrC, and simulating chemical bleaching of BrC contribute distinguishable radiative effects and should be accounted for in BrC schemes. This addition of prognostic BrC to ModelE allows greater physical and chemical complexity in OA representation with no apparent trade-off in model performance, as the evaluation of ModelE aerosol optical depth against Aerosol Robotic Network (AERONET) and Moderate Resolution Imaging Spectroradiometer (MODIS) retrieval data, with and without the BrC scheme, reveals similar skill in both cases. Thus, BrC should be explicitly simulated to allow more physically based chemical composition, which is crucial for more detailed OA studies like comparisons to in situ measurement campaigns. We include a summary of best practices for BrC representation within ModelE at the end of this paper.

brown carbon aerosols↗
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↗