Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation 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 289 records · Page 16

Effects of Size and Shape on the Tolerances for Misalignment and Probabilities for Successful Oriented Attachment of Nanoparticles

Oriented attachment (OA) of nanoparticles is an important pathway of crystal growth, but tools for quantitatively modeling OA are lacking. Here we present several simple models that relate the probability of achieving OA to basic geometric parameters such as particle size, shape, and lattice periodicity. A Moiré-domain model is applied to understand twist-misorientations between parallel surfaces, and it predicts that the range of twist angles yielding perfect OA is inversely related to the width of the contact area. This is confirmed and further developed using a surface functional model, which predicts how crystallographic registration forces drive the emergence of complex orientational energy landscapes. The energy landscapes are predicted to possess local minima that can trap particles in imperfect alignments, and these local minima become deeper and more numerous as the contact area increases, making OA more challenging for large particles. Further, a second set of models is presented to understand the sequence of events by which two crystallographic faces become co-planer after collision. We use a ‘central force approximation’ to quantitatively predict the odds of attaining coalignment between various faces when particles collide with random misalignments, and we show that in the absence of biasing forces, the probability of attaining alignment on a given face is roughly proportional to its solid angle as viewed from the center of the particle. The model predicts that OA is most favorable between well-faceted particles and becomes exceedingly unlikely for large spherical particles that express many microfacets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Perspective—Mass Conservation in Models for Electrodeposition/Stripping in Lithium Metal Batteries

Electrochemical models at different scales and varying levels of complexity have been used in the literature to study the evolution of the anode surface in lithium metal batteries. This includes continuum, mesoscale (phase-field approaches), and multiscale models. In this paper, using a motivating example of a moving boundary model in one dimension, we show how battery models need proper formulation for mass conservation, especially when simulated over multiple charge and discharge cycles. The article concludes with some thoughts on mass conservation and proper formulation for multiscale models.

25 ENERGY STORAGE↗

BattPhase—A Convergent, Non-Oscillatory, Efficient Algorithm and Code for Predicting Shape Changes in Lithium Metal Batteries Using Phase-Field Models: Part I. Secondary Current Distribution

Electrochemical models at different scales and varying levels of complexity have been used in the literature to study the evolution of the anode surface in lithium metal batteries. This includes continuum, mesoscale (phase-field approaches), and multiscale models. Thermodynamics-based equations have been used to study phase changes in lithium batteries using phase-field approaches. However, grid convergence studies and the effect of additional parameters needed to simulate these models are not well-documented in the literature. In this paper, using a motivating example of a moving boundary model in one- and two-dimensions, we show how one can formulate phase-field models, implement algorithms for the same and analyze the results. An open-access code with no restrictions is provided as well. This article concludes with some thoughts on the computational efficiency of phase-field models for simulating dendritic growth.

25 ENERGY STORAGE↗

Reflectance Modeling

A classical description of the one dimensional radiative transfer treatment of vegetation canopies was completed and the results were tested against measured prairie (blue grama) and agricultural canopies (soybean). Phase functions are calculated in terms of directly measurable biophysical characteristics of the canopy medium. While the phase functions tend to exhibit backscattering anisotropy, their exact behavior is somewhat more complex and wavelength dependent. A Monte Carlo model was developed that treats soil surfaces with large periodic variations in three dimensions. A photon-ray tracing technology is used. Currently, the rough soil surface is described by analytic functions and appropriate geometric calculations performed. A bidirectional reflectance distribution function is calculated and, hence, available for other atmospheric or canopy reflectance models as a lower boundary condition. This technique is used together with an adding model to calculate several cases where Lambertian leaves possessing anisotropic leaf angle distributions yield non-Lambertian reflectance; similar behavior is exhibited for simulated soil surfaces.

Smith, J. A.↗

High-resolution heat-transfer-coefficient maps applicable to compound-curve surfaces using liquid crystals in a transient wind tunnel

Tests were performed in a transient heat transfer tunnel in which the model under test was preheated prior to allowing room temperature air to be suddenly drawn over the model. The resulting movement of isothermal contours on the model is revealed using a surface coating of thermochromic liquid crystals that display distinctive colors at particular temperatures. A video record is obtained of a temperature and time data pair for all points on the model during a single test. Experiments on a duct model are reported in which the model was preheated using a hot air stream. A manner in which initial model temperature nonuniformities could be taken into account was investigated. The duct model was also tested with a steady-state measurement technique and results were compared with the transient measurements, but recognizing that differences existed between the upstream thermal boundary conditions. The steady-state and transient measurements were shown to be consistent with predicted values. The main advantage of this transient heat transfer technique using liquid crystals is that since the test model need not be actively heated, high-resolution measurements on surfaces with complex shapes may be obtained.

Jones, T. V.↗

High-resolution heat-transfer-coefficient maps applicable to compound-curve surfaces using liquid crystals in a transient wind tunnel

Tests were performed in a transient heat transfer tunnel in which the model under test was preheated prior to allowing room temperature air to be suddenly drawn over the model. The resulting movement of isothermal contours on the model is revealed using a surface coating of thermochromic liquid crystals that display distinctive colors at particular temperatures. A video record is obtained of a temperature and time data pair for all points on the model during a single test. Experiments on a duct model are reported in which the model was preheated using a hot air stream. A manner in which initial model temperature nonuniformities could be taken into account was investigated. The duct model was also tested with a steady-state measurement technique and results were compared with the transient measurements, but recognizing that differences existed between the upstream thermal boundary conditions. The steady-state and transient measurements were shown to be consistent with predicted values. The main advantage of this transient heat transfer technique using liquid crystals is that since the test model need not be actively heated, high-resolution measurements on surfaces with complex shapes may be obtained.

Jones, Terry V.↗

Properties of Backscattered Ultraviolet Reflectances in Cloudy Atmosphere

We use a Plane-Parallel Cloud (PPC) model to illustrate how Mie scattering from cloud particles interacts with Rayleigh scattering in the atmosphere and produces a complex wavelength dependence in the top-of-the-atmosphere (TOA) reflectances measured by satellite instruments that operate in the ultraviolet (UV) part of the spectrum. Comparisons of the PPC model-derived spectral dependence of reflectances with the Total Ozone Mapping Spectrometer (TOMS) measurements show surprisingly good agreement over a wide range of observational conditions. The PPC model results also are compared with the results of two other cloud models: Lambert Equivalent Reflectivity (LER) and Modified Lambert Equivalent Reflectivity (MLER) that have been used to analyze satellite data in the UV. These models assume that clouds are opaque Lambertian reflectors rather than Mie scattering particles. Although one of these models (MLER) agrees reasonably well with the data, the results from this model appear somewhat unphysical and may not be suitable for interpreting satellite data if one desires high accuracy. We also use the PPC model to illustrate how clouds can perturb tropospheric O3 absorption in complex ways that cannot be explained by models that treat them as reflecting surfaces rather than as volume scatterers.

Ahmad, Z.↗

Next-Generation Ecosystem Experiments–Arctic Flyer

An important challenge for Earth system models (ESMs) is to accurately represent land surface and subsurface processes and their complex interactions. This is true for all regions of the world, but is especially important for high-latitude Arctic ecosystems that are characterized by ice-rich landscapes where topography, hydrology, vegetation, and biogeochemistry are inexplicitly linked. To address this challenge, the Terrestrial Ecosystem Science (TES) program within the Department of Energy’s (DOE) Office of Biological and Environmental Research (BER) is supporting a next-generation ecosystem experiments project in the Arctic (NGEE–Arctic).

54 ENVIRONMENTAL SCIENCES↗

Modeling of Turbulent Flow in Electromagnetically Levitated Metal Droplets

This article details an effort to improve the understanding and prediction of turbulent flow inside a droplet of molten metal levitated in an electromagnetic field. It is shown that the flow field in a test case, a nickel droplet levitated under microgravity conditions, is in the transitional regime between laminar and turbulent flow. Past research efforts have used laminar, enhanced viscosity, and k-epsilon turbulence models to describe the flow. The method highlighted in our study is the renormalization group (RNG) algorithm. We show that an accurate description of the turbulent eddy viscosity is critical in order to obtain realistic velocity fields, and that the turbulent eddy viscosity cannot be uniform in levitated droplets. The RNG method does not impose isotropic length or time scales on the flow field, thus allowing such nonuniform features to be captured. A number of other materials processing applications exhibit similarly complex flow characteristics, such as highly recirculating, transitional, and free surface flows, for which this modeling approach may prove useful.

Berry, S.↗

Mechanisms of dissolution from gibbsite step edges elucidated by ab initio molecular dynamics with enhanced sampling

Surface reactions underpin our collective understanding of mineral dissolution. This impacts a range of predictive geochemical models (for weathering, the fate and transport of metals) as well as industrial utilization of minerals. Yet atomistic details are rarely known due to the complex mineral/fluid interfacial environment. There is significant need to understand the mechanistic details of dissolution reactions and how they depend on surface morphology and solution conditions. This work utilizes surface pit models in conjunction with changes to solution composition that mimic pH to explore surface detachment of aluminate from gibbsite, which is a primary source of Al in soils and within the industrial processing of aluminum. Ab initio molecular dynamics simulations with enhanced sampling has been used to explore the detailed process of the detachment, the results of which indicate two potential pathways that are differentiated based upon the extent of water hydration. Here, the heights of the energy barriers depend upon the local morphology which influence the number of bridges (quasi-)simultaneously broken (1 or 2) or the Al-O coordination of the neighboring aluminum atoms (5 or 6) at the armchair edge. pH effects are significant, with a nearly 50% reduction in barrier height under alkaline conditions that are relevant to geothermal fluids and Al extraction from minerals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessment of Buffet Forcing Function Development Process Using Unsteady Pressure Sensitive Paint

A wind tunnel test was conducted at the Ames Unitary Plan Wind Tunnel to characterize the transonic buffet environment of a generic launch vehicle forebody. The test examined a highly instrumented version of the Coe and Nute Model 11 test article first tested in the 1960s. One of the measurement techniques used during this test was unsteady pressure sensitive paint (uPSP) developed at the Arnold Engineering Development Complex. This optical measurement technique measured fluctuating pressures at over 300,000 locations on the surface of the model. The high spatial density of these measurements provided an opportunity to examine in depth the assumptions underpinning the development of buffet forcing functions (BFFs) used in the development of the Space Launch System vehicle. The comparison of discrete-measurement-based BFFs to BFFs developed by continuous surface pressure integration indicates that the current BFF development approach under predicts low frequency content of the BFFs while over predicting high frequency content. Coherence-based adjustments employed to reduce over prediction in the surface integration of discrete pressure measurements contribute to the inaccuracy of the BFFs and their implementation should be reevaluated.

Sekula, Martin K.↗

Arctic Amplification: Process Drivers and Sources of Uncertainty

Arctic amplification (AA) is a coupled atmosphere-sea ice-ocean process. This understanding has evolved from the early concept of AA, as a consequence of snow ice line progressions, through more than a century of research that has clarified the relevant processes and driving mechanisms of AA. The predictions made by early modeling studies, namely the fall/winter maximum, bottom-heavy structure, the prominence of surface albedo feedback, and the importance of stable stratification have withstood the scrutiny of multi-decadal observations and more complex models. Yet, the uncertainty in Arctic climate projections is larger than in any other region of the planet, making the assessment of high-impact, near-term regional changes difficult or impossible. Reducing this large spread in Arctic climate projections requires a quantitative process understanding. This presentation synthesizes current knowledge of AA and describes a set of recommendations to guide future research. It briefly reviews the history of AA science, summarizes observed Arctic changes, discusses modeling approaches and feedback diagnostics, and assesses the current understanding of the most relevant feedbacks to AA. These sections culminate in a conceptual model of the fundamental physical mechanisms causing AA and a collection of recommendations to accelerate progress towards reduced uncertainty in Arctic climate projections. Our conceptual model highlights the need to account for local feedback and remote process interactions within the context of the annual cycle to constrain projected AA.

Patrick C Taylor↗

Bridging Empirical and Physical Approaches for Landslide Monitoring and Early Warning

Rainfall-triggered landslides typically occur and are evaluated at local scales, using slope-stability models to calculate coincident changes in driving and resisting forces at the hillslope level in order to anticipate slope failures. Over larger areas, detailed high resolution landslide modeling is often infeasible due to difficulties in quantifying the complex interaction between rainfall infiltration and surface materials as well as the dearth of available in situ soil and rainfall estimates and accurate landslide validation data. This presentation will discuss how satellite precipitation and surface information can be applied within a landslide hazard assessment framework to improve landslide monitoring and early warning by considering two disparate approaches to landslide hazard assessment: an empirical landslide forecasting algorithm and a physical slope-stability model. The goal of this research is to advance near real-time landslide hazard assessment and early warning at larger spatial scales. This is done by employing high resolution surface and precipitation information within a probabilistic framework to provide more physically-based grounding to empirical landslide triggering thresholds. The empirical landslide forecasting tool, running in near real-time at http://trmm.nasa.gov, considers potential landslide activity at the global scale and relies on Tropical Rainfall Measuring Mission (TRMM) precipitation data and surface products to provide a near real-time picture of where landslides may be triggered. The physical approach considers how rainfall infiltration on a hillslope affects the in situ hydro-mechanical processes that may lead to slope failure. Evaluation of these empirical and physical approaches are performed within the Land Information System (LIS), a high performance land surface model processing and data assimilation system developed within the Hydrological Sciences Branch at NASA's Goddard Space Flight Center. LIS provides the capabilities to quantify uncertainty from model inputs and calculate probabilistic estimates for slope failures. Results indicate that remote sensing data can provide many of the spatiotemporal requirements for accurate landslide monitoring and early warning; however, higher resolution precipitation inputs will help to better identify small-scale precipitation forcings that contribute to significant landslide triggering. Future missions, such as the Global Precipitation Measurement (GPM) mission will provide more frequent and extensive estimates of precipitation at the global scale, which will serve as key inputs to significantly advance the accuracy of landslide hazard assessment, particularly over larger spatial scales.

Kirschbaum, Dalia↗

Satellite Derived Land Surface Temperature for Model Assimilation

Studies have shown that land surface temperature (LST) tendencies are sensitive to the surface moisture availability which is a function of soil moisture and vegetation. The assimilation of satellite derived LST tendencies into the surface energy budget of mesoscale models has shown promise in improving the representation of the complex effects of both soil moisture and vegetation within the models for short term simulations. LST derived from geostationary satellites has the potential of providing the temporal and spatial resolution needed for an LST assimilation process. This paper presents an analysis comparing the LST derived from GOES-8 infrared measurements with LST calculated by the MM5 numerical model. The satellite derived LSTs are calculated using a physical split window approach using channels 4 and 5 of GOES-8. The differences in the LST data sets, especially the tendencies, are presented and examined. Quantifying the differences between the data sets provide insight of possible weaknesses in the model parameterizations affecting the surface energy budget calculations and an indication of the potential effectiveness o f assimilating LST into the models.

Suggs, Ronnie J.↗

Dynamic PRA-Based Estimation of PWR Coping Time Using a Surrogate Model for Accident Tolerant Fuel

In this study, we propose an interpolation-based response surface surrogate methodology to manage a large number of scenarios in dynamic probabilistic risk assessment. It adopts the shape Dynamic Time Warping algorithm to cluster the interpolation neighborhood from time series sample data. The interpolation method was adapted from Taylor Kriging to allow a reduced-order model of the Taylor series. In order to demonstrate its applicability to complex issues in risk assessment for nuclear engineering, an example risk response surface to estimate emergency core cooling system (ECCS) criteria for triplex silicon carbide (SiC) accident-tolerant fuel was constructed. The response surface was exploited to estimate the cumulative failure probability of the fuel cladding structure due to the uncertainties in operator actions and safety systems. The functional failures were assessed based on a combination of individual layer failures computed by coupling Risk Analysis Virtual Environment software with a pressurized water reactor 1000-MW(electric) RELAP5 model and the in-house fuel performance assessment module. Results showed that SiC cladding failure probability spiked less than 1 min after a large-break loss-of- coolant accident whenever the current ECCS criteria for Zircaloy-4 (Zr-4) cladding was used. However, it still provides an increased safety margin of three orders of magnitude compared to Zr-4. This positive margin could be utilized to relax active ECCS requirements by allowing deviations of up to 450 s in its actuation time. The proposed surrogate methodology generated a response surface of SiC cladding failure probability reasonably well, with a significant savings of computation time. This methodology is expected to be useful in the analysis of system response with complex uncertainty sources.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Unsteady transonic potential flow over a flexible fuselage

A flexible fuselage capability has been developed and implemented within version 1.2 of the CAP-TSD code. The capability required adding time dependent terms to the fuselage surface boundary conditions and the fuselage surface pressure coefficient. The new capability will allow modeling the effect of a flexible fuselage on the aeroelastic stability of complex configurations. To assess the flexible fuselage capability several steady and unsteady calculations have been performed for slender fuselages with circular cross-sections. Steady surface pressures are compared with experiment at transonic flight conditions. Unsteady cross-sectional lift is compared with other analytical results at a low subsonic speed and a transonic case has been computed. The comparisons demonstrate the accuracy of the flexible fuselage modifications.

Gibbons, Michael D.↗

Evaluation of Planetary Boundary Layer Structure from NASA Global Modeling and Assimilation Office’s Next Retrospective Analysis Product MERRA-21C

The Planetary Boundary Layer (PBL) is a complex interface between Earth’s surface and the atmosphere with high spatiotemporal variability in its characteristics, and accurate simulation and observation of the PBL has proven to be a challenge. In this study, we evaluate PBL structures from NASA Global Modeling and Assimilation Office (GMAO)’s next retrospective analysis product: the Goddard Earth Observing System Retrospective Analysis of the early 21st Century (GEOS-R21C), with 25-km horizontal resolution. The GEOS-R21C incorporates a wide range of observing systems and various improvements over previous GMAO reanalysis. We compare the PBL thermodynamic structure as well as PBL Height (PBLH) from GEOS-R21C with various PBL observations. For PBLH, definitions are different depending on each observing system. PBLHs derived from radiosonde observations and Global Navigation Satellite System Radio Occultation (GNSS RO) are based on bulk Richardson number and refractivity gradient, respectively. The results of the evaluation of PBL structure from GMAO’s next retrospective analysis product will be presented. In addition, comparison between near-real-time (NRT) commercial and NASA commercial GNSS RO PBLH will also be discussed.

Eun-Gyeong Yang↗

A general analysis of propagation along multiple-layer superconducting stripline and microstrip transmission lines

A rigorous spectral-domain formulation for a superconducting stripline or microstrip transmission line with a multiple-layer dielectric substrate is presented. The formulation models the strip conductor as a surface current with an equivalent surface impedance, where the surface impedance is approximated in closed form when the strip is either much thinner or much thicker than a penetration depth. In either case the surface impedance is related to the complex conductivity of the material, which is calculated from a two-fluid model. Results are presented to show the slow-wave propagation and attenuation along both microstrip and stripline packages in a realistic multiple-layer configuration, which accounts for the field penetration into the superconducting ground planes.

Nghiem, David↗