Search NASA⌕ Search

SEARCH · Search NASA

Results for “discrete surface”

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 OSTI2023

Prediction of hemiwicking dynamics in micropillar arrays

Dynamic hemiwicking behavior is observable in both nature and a wide range of industrial applications ranging from biomedical devices to thermal management. We present a semi-analytical modeling framework (without empirical fitting coefficients) to predict transient capillary-driven hemiwicking behavior of a liquid through a nano/microstructured surface, specifically a micropillar array. In our model framework, the liquid domain is discretized into micropillar unit cells to enable the time marching of the hemiwicking front. A simplified linear pressure drop is assumed along the hemiwicking length such that the local meniscus curvature, contact angle, and effective liquid height are determined at each time step in our transient model. This semi-analytical model is validated with experimental data from our own experiments and from published literature for different fluids. Our model predicts hemiwicking dynamics with <20% error over a broad range of micropillar geometries with height-to-pitch ratio ranging between ≈0.34 and 6.7 and diameter-to-pitch ratio in the range of ≈0.25–0.7 and without any fitting parameters. For lower diameter-to-pitch ratio data points related to sparse micropillar array arrangements, we suggest modifications to the semi-analytical model. This work sheds light on complex and dynamic solid–liquid–vapor interfacial interactions which could serve as a guide for the design of textured surfaces for wicking enhancement in multi-phase thermal and mass transport technologies and applications.

Mechanics↗
DOE OSTI2020

Atmosphere-Biosphere Interaction Study (Field Campaign Report)

The objective of this campaign was to collect hyperspectral reflectance measurements and complementary optical and infrared data of land surface (i.e., vegetation and soil) throughout the growing season (i.e., from March to September) at high temporal frequency (e.g., < 30 minutes). To collect the measurements, the EcoSpec system, a tower-based sensor system housing seven optical and infrared sensors that was developed at Argonne National Laboratory, was deployed to the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility Southern Great Plains (SGP) atmospheric observatory (within the crop, south of the eddy correlation flux measurement system [ECOR] tower). The EcoSpec system consists of a spectroradiometer having 2,151 channels (350–2,500 nm), red-greenblue (RGB) camera, diffuse radiometer, thermal infrared (TIR) sensors, albedometers, Photochemical Reflectance Index (PRI) sensors, Normalized Difference Vegetation Index (NDVI) sensors, and wetness sensor. The system is designed to collect an array of land surface measurements, including hyperspectral reflectance of plants and soils, RGB photos capturing contextual information of land surface, radiant temperature of sky and land surface, incoming and reflected shortwave radiation, and incoming radiation components (direct and diffuse) as frequent as every minute from up to 12 discrete positions around the tower. This would allow us to capture inherent heterogeneity of the terrestrial ecosystem and to obtain representative spectral signatures of the plant and soil surfaces via averaging those measurements. The SGP site provides advanced atmospheric and other physical measurements that are essential for our investigation of ecosystem functions using hyperspectral reflectance measurements of land surfaces. Located in the mid-latitudes, the SGP site exhibits dynamic seasonality and a wide range of meteorological variability and flux properties that make it ideal for studying interactions between near-surface atmosphere and biosphere across a range of environmental conditions. The atmosphere, plants, and soil control terrestrial carbon and water cycles. To make more accurate climate forecasts, researchers need to understand ecosystem dynamics at the biosphere–atmosphere interface. Even though our ability to forecast ecosystem dynamics and climate at both the regional and global scales has advanced significantly, we still do not know how local phenomena affecting water and carbon fluxes—such as diurnal variations—interrelate with large-scale atmospheric/climate phenomena, and vice versa. The lack of a more thorough understanding of ecosystem dynamics and climate interactions is a key knowledge gap for the improvement of earth system models (ESMs).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Progress on the National Solar Radiation Data Base (NSRDB): A New DNI Computation

This study introduces a new technique to compute direct normal irradiance (DNI) for improving the National Solar Radiation Data Base (NSRDB). A finite-surface integration algorithm is developed to compute solar radiation in differential solid angles and efficiently infer its contribution to a surface perpendicular to the solar direction. A lookup table of cloud bi-directional transmittance distribution function (BTDF) is developed by use of the discrete ordinates radiative transfer (DISORT) model for possible solar and observing directions and various cloud optical and microphysical properties. In each solar incident direction, DNI is given by the cloud BTDFs from approximately 200 differential solid angles. The simulated DNI is calibrated and evaluated using surface observations at the National Renewable Energy Laboratory's (NREL's) Solar Energy Research Laboratory (SRRL) and the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) facility.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2024

Sensitivity Analysis in the Presence of Intrinsic Stochasticity for Discrete Fracture Network Simulations

Abstract Large‐scale discrete fracture network (DFN) simulators are standard fare for studies involving the sub‐surface transport of particles since direct observation of real world underground fracture networks is generally infeasible. While these simulators have successfully been used in several engineering applications, estimates of output quantities of interest (QoI) — such as breakthrough time of particles reaching the edge of the system — suffer from two distinct types of uncertainty. A run of a DFN simulator requires several parameters to be set that dictate the placement and size of fractures, the density of fractures, and the overall permeability of the system; uncertainty on the proper parameters will lead to uncertainty in the QoI, called epistemic uncertainty. Furthermore, since these input settings to DFN simulators control the stochastic processes which place fractures and govern flow, understanding how this randomness affects the QoI requires several runs of the simulator at distinct random seeds. The uncertainty in the QoI attributed to different realizations (i.e., different seeds) of the same random process (i.e., identical input parameters) leads to a second type of uncertainty, called aleatoric uncertainty. In this paper, we perform a Sensitivity Analysis, which directly attributes the uncertainty observed in the QoI to the epistemic uncertainty from each input parameter and to the aleatoric uncertainty. Beyond the specific takeaways on which input variables influence uncertainty in the QoI the most, a major contribution of this paper is the introduction of a statistically rigorous workflow for characterizing the uncertainty in DFN flow simulations that exhibit heteroskedasticity.

58 GEOSCIENCES↗
DOE OSTI2020

A Bayesian Approach to Regional Decadal Predictability: Sparse Parameter Estimation in High-Dimensional Linear Inverse Models of High-Latitude Sea Surface Temperature Variability

Stochastic reduced models are an important tool in climate systems whose many spatial and temporal scales cannot be fully discretized or underlying physics may not be fully accounted for. One form of reduced model, the linear inverse model (LIM), has been widely used for regional climate predictability studies—typically focusing more on tropical or midlatitude studies. However, most LIM fitting techniques rely on point estimation techniques deriving from fluctuation–dissipation theory. In this methodological study we explore the use of Bayesian inference techniques for LIM parameter estimation of sea surface temperature (SST), to quantify the skillful decadal predictability of Bayesian LIM models at high latitudes. We show that Bayesian methods, when compared to traditional point estimation methods for LIM-type models, provide better calibrated probabilistic skill, while simultaneously providing better point estimates due to the regularization effect of the prior distribution in high-dimensional problems. We compare the effect of several priors, as well as maximum likelihood estimates, on 1) estimating parameter values on a perfect model experiment and 2) producing calibrated 1-yr SST anomaly forecast distributions using a preindustrial control run of the Community Earth System Model (CESM). Finally, we employ a host of probabilistic skill metrics to determine the extent to which an LIM can forecast SST anomalies at high latitudes. We find that the choice of prior distribution has an appreciable impact on estimation outcomes, and priors that emphasize physically relevant properties enhance the model’s ability to capture variability of SST anomalies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Meshfree methods on manifolds for hydrodynamic flows on curved surfaces: A Generalized Moving Least-Squares (GMLS) approach

We utilize generalized moving least squares (GMLS) to develop meshfree techniques for discretizing hydrodynamic flow problems on manifolds. We use exterior calculus to formulate incompressible hydrodynamic equations in the Stokesian regime and handle the divergence-free constraints via a generalized vector potential. This provides less coordinate-centric descriptions and enables the development of efficient numerical methods and splitting schemes for the fourth-order governing equations in terms of a system of second-order elliptic operators. Using a Hodge decomposition, we develop methods for manifolds having spherical topology. We show the methods exhibit high-order convergence rates for solving hydrodynamic flows on curved surfaces. The methods also provide general high-order approximations for the metric, curvature, and other geometric quantities of the manifold and associated exterior calculus operators. The approaches also can be utilized to develop high-order solvers for other scalar-valued and vector-valued problems on manifolds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Zonal-Based Emission Source Term Model for Predicting Particulate Emission Factors in Wildfire Simulations

A physics/chemistry-based numerical model for predicting the emission of fine particles from wildfires is proposed. This model implements the fundamental mechanisms of soot formation in a combustion environment: soot nucleation, surface growth, agglomeration, oxidation, and particle fragmentation. These mechanisms occur on a scale too fine for the discretization of most wildfire models, which need to simulate landscape-scale dynamics. As a result this model implements a zonal approach, where the computed soot particle distribution is partitioned into process zones within a single resolved grid cell. These process zones include: an inception zone (for nucleation), a heating zone (for coagulation, surface growth, and fragmentation), a reaction zone (for oxidation), and a quenched zone (for atmospheric processes). Governing mechanisms are applied to the appropriate zones to predict total particle growth and emission. Additionally, the proposed model is implemented into HIGRAD/FIRETEC, a physics-based wildfire simulation code which couples interactions between fire, fuels, atmosphere, and topography on a landscape scale. Fire simulations among grasslands and conifer forests are performed and compared against experimental data for emission factors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗
DOE OSTI2020

An adaptive space-time phase field formulation for dynamic fracture of brittle shells based on LR NURBS

Here, we present an adaptive space-time phase field formulation for dynamic fracture of brittle shells. Their deformation is characterized by the Kirchhoff–Love thin shell theory using a curvilinear surface description. All kinematical objects are defined on the shell’s mid-plane. The evolution equation for the phase field is determined by the minimization of an energy functional based on Griffith’s theory of brittle fracture. Membrane and bending contributions to the fracture process are modeled separately and a thickness integration is established for the latter. The coupled system consists of two nonlinear fourth-order PDEs and all quantities are defined on an evolving two-dimensional manifold. Since the weak form requires C 1 -continuity, isogeometric shape functions are used. The mesh is adaptively refined based on the phase field using Locally Refinable (LR) NURBS. Time is discretized based on a generalized-α method using adaptive time-stepping, and the discretized coupled system is solved with a monolithic Newton–Raphson scheme. Finally, the interaction between surface deformation and crack evolution is demonstrated by several numerical examples showing dynamic crack propagation and branching.

42 ENGINEERING↗
DOE OSTI2022

Reaction–drift–diffusion models from master equations: application to material defects

We present a general method to produce well-conditioned continuum reaction–drift–diffusion equations directly from master equations on a discrete, periodic state space. We assume the underlying data to be kinetic Monte Carlo models (i.e. continuous-time Markov chains) produced from atomic sampling of point defects in locally periodic environments, such as perfect lattices, ordered surface structures or dislocation cores, possibly under the influence of a slowly varying external field. Our approach also applies to any discrete, periodic Markov chain. Here, the analysis identifies a previously omitted non-equilibrium drift term, present even in the absence of external forces, which can compete in magnitude with the reaction rates, thus being essential to correctly capture the kinetics. To remove fast modes which hinder time integration, we use a generalized Bloch relation to efficiently calculate the eigenspectrum of the master equation. A well conditioned continuum equation then emerges by searching for spectral gaps in the long wavelength limit, using an established kinetic clustering algorithm to define a proper reduced, Markovian state space.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Highly dispersed Co deposited on Al 2 O 3 particles via CoCp 2 + H 2 ALD

Highly dispersed cobalt atoms were deposited on porous alumina particles using atomic layer deposition (ALD) with a CoCp 2 /H 2 chemistry at approximately 7 wt%. H 2 did not completely reduce the cyclopentadienyl organic ligands bound to deposited Co atoms at ALD reaction conditions. A sharp decline in Co deposited per cycle for two or more ALD cycles indicates that much of the Al 2 O 3 surface is sterically blocked from further CoCp 2 deposition after the first CoCp 2 exposure. Temperature programmed reduction confirmed that the adsorbed precursor organic ligands persist after H 2 exposures during ALD and temperatures as high as 500 °C are required to fully reduce the organic ligands to CH 4 . High resolution, element sensitive imaging showed that Co atoms were dispersed on the Al 2 O 3 surface and could deposit in previously unobserved multiple growth morphologies, specifically layers that were continuous over several angstroms or discrete nanoparticles. Density functional theory calculations were used to examine Co atom adsorption, show the altered haptic binding of cracked Cp ligands, and to calculate the thermodynamics of Cp ligand decomposition. The lateral steric hindrance between organic ligands bound to deposited Co atoms, Cp ligand decomposition mechanism, and local Al 2 O 3 surface termination all likely determine the observed Co growth morphology during initial ALD cycles.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Kinetic pathways of fast lithium transport in solid electrolyte interphases with discrete inorganic components

The transport of lithium ions in the solid electrolyte interphase (SEI) has been previously accepted to proceed in two steps: a fast pore diffusion through the outer, porous organic layer followed by a slow knock-off or vacancy diffusion in the inner, dense inorganic layer. The second step is believed to be the rate-limiting step during fast-charging. In this study, we have intentionally constructed a thicker SEI (SEI-rich) structure on the surface of monoclinic Nb 2 O 5 (H-Nb 2 O 5 ) by adding LiNO 3 into a conventional ethylene carbonate based electrolyte. The electrochemical performance of two electrodes, one SEI-rich and one with few SEI (SEI-lean), was found to be almost the same, including their fast-charging capability and cycling stability, despite the significant difference in their SEI structure. Importantly, analysis using cryogenic scanning/transmission electron microscopy showed the discrete decoration of individual inorganic particles (e.g., Li 2 O) and amorphous species (LiN x O y /organic components) over the surface of H-Nb 2 O 5 . These discrete inorganic particles are in contradiction to the formation of dense inner inorganic layer, which has been commonly postulated. Based on these findings, we propose a new mechanism for Li ion transport through the SEI: one-step pore diffusion, without the second step slow diffusion. This one-step pore diffusion process provides an extremely fast Li ion transport, and effectively removes the kinetic limitation of Li ion transport in the SEI for fast charging. Finally, these results strongly suggest that the influence of SEI structure on the transport kinetics of lithium ions is much less significant than previously accepted. These results offer a new understanding of possible lithium ion transport pathway within SEI and may have implications for the future designs of fast-charging battery materials.

25 ENERGY STORAGE↗
DOE OSTI2020

Effect of Particle Spreading Dynamics on Powder Bed Quality in Metal Additive Manufacturing

Powder spreading precedes creation of every new layer in powder bed additive manufacturing (AM). The powder spreading process can lead to powder layer defects such as porosity, poor surface roughness and particle segregation. Therefore, the creation of homogeneous layers is the first task for optimal part printing. Discrete element methods (DEM) powder spreading simulations are typically limited to a single layer and/or small number of particles. Therefore, results from such model configurations may not be generalized to multiple layer processes. In this study, a computationally efficient multi-layer powder spreading DEM simulation model is proposed. The model is calibrated experimentally using static Angle of Repose measurements. The adhesion model parameter, cohesive energy density is related to adhesive surface energy and strain energy release rate parameters. The model results show that interaction between particle and the powder spreading rake leads to noticeable variation in packing density, surface roughness, dynamic angle of repose (AOR), particle size distribution, and particle segregation. Finally, the powder model is experimentally validated using a recoater spreading rig to measure the dynamic AOR at spreading speeds consistent with recoating speeds and layer heights used in AM processes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Flowability of Crumbler Rotary Shear Size-Reduced Granular Biomass: An Experiment-Informed Modeling Study on the Angle of Repose

Biomass has potential as a carbon-neutral alternative to petroleum for chemical and energy products. However, complete replacement of fossil fuel is contingent upon efficient processes to eliminate undesirable characteristics of biomass, e.g., low bulk density, variability, and storage-induced quality problems. Mechanical size reduction via comminution is a processing operation to engineer favorable biomass flowability in handling. Crumbler rotary shear mill has been empirically demonstrated to produce more uniformly shaped particles with higher flowability than hammermilled biomass. This study combines modeling and experimentation to unveil fundamental understandings of the relation between granular particle characteristics and biomass flow behavior, which elucidate underlying mechanisms and guide selection of critical processing parameters. For this purpose, the impact of critical material attributes, including particle size (2–6 mm), particle shape (briquette, chip, clumped-sphere, cube, etc.), and surface roughness, on the angle of repose (AOR) of milled pine chips were investigated using discrete element method (DEM) simulations. Forest Concepts Crumbler rotary shear system is used to produce milled pine particles within the same size range considered in DEM simulations. AOR of different sets of these particles were measured experimentally to benchmark DEM results against experimental data. Specific energy consumption for the comminution of biomass with different particle size and moisture content are measured for technoeconomic analysis. Our results show that the smaller size (2 mm) of pine particle achieves better followability (i.e., smaller AOR) while the energy cost of comminution is significantly higher and bulk density is almost the same as the 6-mm pine particles. For the 2-mm particle size, Crumbles from veneer have better flow properties than Crumbles from chips. Contrarily, no significant difference was observed between the AOR of the two materials for the 6-mm particle size. Furthermore, from DEM simulations, mechanical interlocking between particles was found as a dominant factor in determining AOR of complex-shaped particles such as milled pine, which cannot be accurately captured by using simple particle shapes (e.g., mono-sphere) with a rolling resistance model. Conversely, clumped-sphere model alleviates this limitation without increasing computational cost significantly and can be used for accurate representation of biomass granular particles when simulating free-flow behavior.

09 BIOMASS FUELS↗
DOE OSTI2020

A comparison of lodgepole and spruce needle chemistry impacts on terrestrial biogeochemical processes during isolated decomposition

This study investigates the isolated decomposition of spruce and lodgepole conifer needles to enhance our understanding of how needle litter impacts near-surface terrestrial biogeochemical processes. Harvested needles were exported to a subalpine meadow to enable a discrete analysis of the decomposition processes over 2 years. Initial chemistry revealed the lodgepole needles to be less recalcitrant with a lower carbon to nitrogen (C:N) ratio. Total C and N fundamentally shifted within needle species over time with decreased C:N ratios for spruce and increased ratios for lodgepole. Differences in chemistry correlated with CO 2 production and soil microbial communities. The most pronounced trends were associated with lodgepole needles in comparison to the spruce and needle-free controls. Increased organic carbon and nitrogen concentrations associated with needle presence in soil extractions further corroborate the results with clear biogeochemical signatures in association with needle chemistry. Interestingly, no clear differentiation was observed as a function of bark beetle impacted spruce needles vs those derived from healthy spruce trees despite initial differences in needle chemistry. These results reveal that the inherent chemistry associated with tree species has a greater impact on soil biogeochemical signatures during isolated needle decomposition. By extension, biogeochemical shifts associated with bark beetle infestation are likely driven more by changes such as the cessation of rhizospheric processes than by needle litter decomposition.

Bark Beetle Disturbance↗
DOE OSTI2022

Vibrational levels of a generalized Morse potential

A Generalized Morse Potential (GMP) is an extension of the Morse Potential (MP) with an additional exponential term and an additional parameter that compensate for MP’s erroneous behavior in the long range part of the interaction potential. Because of the additional term and parameter, the vibrational levels of the GMP cannot be solved analytically, unlike the case for the MP. We present several numerical approaches for solving the vibrational problem of the GMP based on Galerkin methods, namely, the Laguerre Polynomial Method (LPM), the Symmetrized LPM, and the Polynomial Expansion Method (PEM), and apply them to the vibrational levels of the homonuclear diatomic molecules B 2 , O 2 , and F 2 , for which high level theoretical near full configuration interaction (CI) electronic ground state potential energy surfaces and experimentally measured vibrational levels have been reported. Overall, the LPM produces vibrational states for the GMP that are converged to within spectroscopic accuracy of 0.01 cm –1 in between 1 and 2 orders of magnitude faster and with much fewer basis functions/grid points than the Colbert–Miller Discrete Variable Representation (CN-DVR) method for the three homonuclear diatomic molecules examined in this study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Discrete Orientations of Interfacial Waters Direct Crystallization of Mica-Binding Proteins

Understanding the basis of templated molecular assembly on a solid surface requires a fundamental comprehension of both short- and long-range aqueous response to the surface under a variety of solution conditions. Herein we provide a detailed picture of how the molecular-scale response to different mica surfaces yields distinct solvent orientations that produce quasi-static directional potentials onto which macromolecules can adsorb. We connect this directionality to observed (a)symmetric epitaxial alignment of designed proteins onto these surfaces, corroborate our findings with 3D atomic force microscopy experiments, and identify slight differences in surface structure as the origin of this effect. Furthermore, our work provides a detailed picture of the intrinsic electrolyte response in the vicinity of mineral interfaces, with clear predictions for experiment, and highlights the role of solvent on the predictive assembly of hierarchical materials on mineral surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Implementation of extrinsic cohesive zone model (ECZM) in 2D finite-discrete element method (FDEM) using node binding scheme

The combined finite-discrete element method (FDEM) has been widely used for rock fracturing simulations. Conventionally, FDEM is realized using the intrinsic cohesive zone model (ICZM); however, it has the drawback of artificial compliance and high computational expense. As a complement, the extrinsic cohesive zone model (ECZM) is seen to be realized in FDEM recently, whereas the node splitting scheme utilized is cumbersome. Here, within the framework of ICZM-based FDEM, we propose a node binding scheme to efficiently bind the pre-discretized finite elements and thus guarantee the continuum behavior of materials in the elastic stage. The yield surfaces, controlled by ECZM, are dynamically embedded by invoking the pre-inserted cohesive elements. The effectiveness and efficiency of the proposed approach are validated and tested by performing a suite of numerical experiments. Compared with ICZM-based FDEM, the proposed approach can correctly capture material deformation and reduce the computation cost. In contrast to the existing ECZM-based FDEM, the proposed approach can overcome the frequent and complex element topology updating. Finally, this work provides a novel perspective that fully inherits the advantages of both ICZM and ECZM, but circumvents their shortcomings, which guarantees a more efficient and effective simulation of brittle material evolution from continuum to discontinuum.

58 GEOSCIENCES↗