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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Application of a Distributed Element Roughness Model to Additively Manufactured Internal Cooling Channels

Abstract Design for cooling effectiveness in turbine blades relies on accurate models for dynamic losses and heat transfer of internal cooling passages. Metal additive manufacturing (AM) has expanded the design space for these configurations, but can give rise to large-scale roughness features. The range of roughness length scales in these systems makes morphology resolved computational fluid dynamics (CFD) impractical. However, volumetric roughness models can be leveraged, as they have computational costs orders of magnitude lower. In this work, a discrete element roughness model (DERM), based on the double-averaged Navier–Stokes equations, is presented and applied to additively manufactured rough channels, representative of gas turbine blade cooling passages. Unique to this formulation of DERM is a generalized sheltering-based treatment of drag, a two-layer model for spatially averaged Reynolds stresses, and explicit treatment of dispersion. Six different AM rough surface channel configurations are studied, with roughness trough to peak sizes ranging from 15% to 60% nominal channel passage half-width, and the roughness Reynolds number ranges from Rek = 60 to 300. DERM predictions for spatially and temporally averaged mean flow quantities are compared to previously reported direct numerical simulation results. Good agreement in the mean velocity profiles, stress balances, and drag partitions are observed. While DERM models are typically calibrated to specific deterministic roughness morphologies at comparatively small roughness Reynolds numbers, the present more generalized DERM formulation has wider applicability. Here, it is demonstrated that the model can accommodate random roughness of large scale, typical of AM.

Engineering↗
DOE OSTI2021

Simulations of biomass pyrolysis using glued-sphere CFD-DEM with 3-D intra-particle models

The shape of biomass particles influences intra-particle heat transfer, particle-particle collisions, interphase momentum, heat transfer, and pyrolysis yield. In this research, the effects of these multi-scale phenomena are captured using a glued-sphere computational fluid dynamic discrete element method (CFD-DEM) with 3-D intra-particle models. The glued-spheres resolved the shapes of particles and intra-particle temperature and species distributions. The implementation was validated using the pyrolysis data of a spherical particle and a cylindrical particle. Then, the influences of biomass shapes and sizes on pyrolysis were investigated. The results revealed the significant influence of surface areas and mass distributions on interphase heat transfer and chemical reactions. The intra-particle heat transfer is neglectable for a Biot number smaller than 0.41. Finally, the pyrolysis of irregularly shaped biomass particles in a fixed bed, fluidized bed, and spouted bed was simulated. Analysis of conversion distributions showed the largest deviations in the fixed bed and the slowest conversion rate in the spouted bed. A good balance between uniform and fast conversion is achieved in the fluidized bed. The multi-scale method proposed in this research provides a tool for the evaluation of large irregular-shaped biomass particle pyrolysis in different types of reactors.

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

Physics-Informed Neural Networks for Heat Transfer Problems

Abstract Physics-informed neural networks (PINNs) have gained popularity across different engineering fields due to their effectiveness in solving realistic problems with noisy data and often partially missing physics. In PINNs, automatic differentiation is leveraged to evaluate differential operators without discretization errors, and a multitask learning problem is defined in order to simultaneously fit observed data while respecting the underlying governing laws of physics. Here, we present applications of PINNs to various prototype heat transfer problems, targeting in particular realistic conditions not readily tackled with traditional computational methods. To this end, we first consider forced and mixed convection with unknown thermal boundary conditions on the heated surfaces and aim to obtain the temperature and velocity fields everywhere in the domain, including the boundaries, given some sparse temperature measurements. We also consider the prototype Stefan problem for two-phase flow, aiming to infer the moving interface, the velocity and temperature fields everywhere as well as the different conductivities of a solid and a liquid phase, given a few temperature measurements inside the domain. Finally, we present some realistic industrial applications related to power electronics to highlight the practicality of PINNs as well as the effective use of neural networks in solving general heat transfer problems of industrial complexity. Taken together, the results presented herein demonstrate that PINNs not only can solve ill-posed problems, which are beyond the reach of traditional computational methods, but they can also bridge the gap between computational and experimental heat transfer.

Engineering↗
DOE OSTI2025

A novel transformation of the ice sheet Stokes equations and some of its properties and applications

We introduce a novel transformation of the Stokes equations into a form closely resembling the shallow Blatter–Pattyn equations. The two forms differ by only a few additional terms, while their variational formulations differ only by a single term in each horizontal direction. Specifically, the variational formulation of the Blatter–Pattyn model drops the vertical velocity in the second invariant of the strain rate tensor. Here we make use of the new transformation in two ways. First, we consider incorporating the transformed equations into a code that can be very easily converted from a Stokes to a Blatter–Pattyn model, and vice versa, by switching these terms on or off. This may be generalized so that the Stokes model is switched on adaptively only where the Blatter–Pattyn model loses accuracy. Second, the key role played by the vertical velocity in the Blatter–Pattyn approximation motivates new approximations. Two examples are presented. These require a mesh that enables the discrete continuity equation to be invertible for the vertical velocity in terms of the horizontal velocity components. Examples of such meshes, such as the first-order P1–E0 mesh and the second-order P2–E1 mesh, are given in both 2D and 3D. However, the transformed Stokes model has the same type of gravity forcing as the Blatter–Pattyn model, determined by the ice surface slope, thereby forgoing some of the mesh generality of the traditional formulation of the Stokes model.

58 GEOSCIENCES↗
DOE OSTI2022

Sas20 is a highly flexible starch-binding protein in the Ruminococcus bromii cell-surface amylosome

Ruminococcus bromii is a keystone species in the human gut that has the rare ability to degrade dietary resistant starch (RS). This bacterium secretes a suite of starch-active proteins that work together within larger complexes called amylosomes that allow R. bromii to bind and degrade RS. Starch adherence system protein 20 (Sas20) is one of the more abundant proteins assembled within amylosomes, but little could be predicted about its molecular features based on amino acid sequence. Here, we performed a structure–function analysis of Sas20 and determined that it features two discrete starch-binding domains separated by a flexible linker. We show that Sas20 domain 1 contains an N-terminal β-sandwich followed by a cluster of α-helices, and the nonreducing end of maltooligosaccharides can be captured between these structural features. Furthermore, the crystal structure of a close homolog of Sas20 domain 2 revealed a unique bilobed starch-binding groove that targets the helical α1,4-linked glycan chains found in amorphous regions of amylopectin and crystalline regions of amylose. Affinity PAGE and isothermal titration calorimetry demonstrated that both domains bind maltoheptaose and soluble starch with relatively high affinity (K d ≤ 20 μM) but exhibit limited or no binding to cyclodextrins. Finally, small-angle X-ray scattering analysis of the individual and combined domains support that these structures are highly flexible, which may allow the protein to adopt conformations that enhance its starch-targeting efficiency. Taken together, we conclude that Sas20 binds distinct features within the starch granule, facilitating the ability of R. bromii to hydrolyze dietary RS.

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

Route Optimization for Energy Efficient Airport Shuttle Operations - A Case Study from Dallas Fort Worth International Airport

Air travel and requisite surface traffic supporting passenger arrival/departure constitutes a significant portion of travel and emissions in cities with large airports. An airport trip can segment into three parts namely: i) travel from a location in the city to the airport; ii) travel from a parking lot or rental car center to the terminal (i.e., within the airport premises), and iii) travel inside the terminal. Depending on the airport access mode all or a part of these legs comprise a traveler’s journey to the airport. The priority of airport ground transport management teams is to provide passengers with a seamless travel experience within the airport, so it is understandable that within airport shuttle routes might not be optimized for minimizing energy consumption. Solutions that meet the dual objective of reducing energy consumption from airport shuttle operations without compromising on passenger travel experience are key to improving system efficiency. There is currently a dearth of research and tools that can inform airports in making such decisions. Addressing this need, this research effort puts forth an optimization model that generates optimal shuttle routes for a given set of constraints, and a discrete-event simulator that evaluates the optimal solutions in a stochastic environment to understand the tradeoffs between passenger wait times, and within airport shuttle energy consumption. The proposed set of tools are tested in the context of optimizing airport shuttles routes within the Dallas Fort Worth International Airport (DFW). In addition to shuttle spatial positioning, and passenger demand information, high-fidelity vehicle data was collected using data loggers installed on DFW shuttles. Results show that 20% energy reduction in shuttle operations is possible with a modest two-minute increase in average passenger wait times. The tools developed in this research effort are designed to be generalizable and can help optimize shuttle operations planning at any major airport.

air travel↗
DOE OSTI2021

Assessment of a tomography-informed polyhedral discrete element modelling approach for complex-shaped granular woody biomass in stress consolidation

The design of handling equipment for granular biomass primarily requires experiments. Discrete element models (DEM) can provide designers with detailed insight into the behaviour of granular materials. However, granular biomass comprising complex-shaped particles is difficult to model with DEM. A tomography-informed DEM approach and an exhaustive assessment of the ability of the approach to predict the bulk behaviour of milled pines using experimental data is presented. Nano-CT scan was conducted to obtain 3D particle surface geometries as the basis for particle shape approximation by a polyhedral model and a sphero-polyhedral model. These models were applied in the simulation of a compressibility test. Our parametric study showed that particle Young's modulus and restitution coefficient are the two main properties influencing the simulated bulk behaviour of the DEM particles and that the level of approximating DEM particle shape for real particles is critical for accurately replicating the bulk behaviour of milled pines. The polyhedral model demonstrated better suitability than its sphero-polyhedral counterpart for modelling pine particles. The polyhedral model calibrated in a compressibility test was then applied in the simulations of a friction test without additional parameter tuning. Interlocking was dominant in the shear of bulk pine particles. Remarkably, the polyhedral model predicted the frictional behaviour of the pine particles when compared to the experimental results. The limitations of the model, as well as possible ways for enhancement, are discussed. This work provided novel insights into the suitability of complex-shaped DEM models for granular woody biomass.

Biomass flow↗
DOE OSTI2020

Efficient optimization method for finding minimum energy paths of magnetic transitions

Here, efficient algorithms for the calculation of minimum energy paths of magnetic transitions are implemented within the geodesic nudged elastic band (GNEB) approach. While an objective function is not available for GNEB and a traditional line search can, therefore, not be performed, the use of limited memory Broyden–Fletcher–Goldfarb–Shanno (LBFGS) and conjugate gradient algorithms in conjunction with orthogonal spin optimization (OSO) approach is shown to greatly outperform the previously used velocity projection and dissipative Landau–Lifschitz dynamics optimization methods. The implementation makes use of energy weighted springs for the distribution of the discretization points along the path and this is found to improve performance significantly. The various methods are applied to several test problems using a Heisenberg-type Hamiltonian, extended in some cases to include Dzyaloshinskii–Moriya and exchange interactions beyond nearest neighbours. Minimum energy paths are found for magnetization reversals in a nano-island, collapse of skyrmions in two-dimensional layers and annihilation of a chiral bobber near the surface of a three-dimensional magnet. The LBFGS-OSO method is found to outperform the dynamics based approaches by up to a factor of 8 in some cases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Mapping snow depth and volume at the alpine watershed scale from aerial imagery using Structure from Motion

Time series mapping of snow volume in the mountains at global scales and at resolutions needed for water resource management is an unsolved challenge to date. Snow depth mapping by differencing surface elevations from airborne lidar is a mature measurement approach filling the observation gap operationally in a few regions, primarily in mountain headwaters in the Western United States. The same concept for snow depth retrieval from stereo- or multi-view photogrammetry has been demonstrated, but these previous studies had limited ability to determine the uncertainties of photogrammetric snow depth at the basin scale. For example, assessments used non-coincident or discrete points for reference, masked out vegetation, or compared a subset of the fully snow-covered study domain. Here, using a unique data set with simultaneously collected airborne data, we compare snow depth mapped from multi-view Structure from Motion photogrammetry to that mapped by lidar at multiple resolutions over an entire mountain basin (300 km 2 ). After excluding reconstruction errors (negative depths), SfM had lower snow-covered area (~27%) and snow volume (~16%) compared to lidar. The reconstruction errors were primarily in areas with vegetation, shallow snow (< 1 m), and steep slopes (> 60°C). Across the overlapping snow extent, snow depths compared well to lidar with similar mean values (< 0.03 m difference) and snow volume (± 5%) for output resolutions of 3 m and 50 m, and with a normalized median absolute deviation of 0.19 m. Our results indicate that photogrammetry from aerial images can be applied in the mountains but would perform best for deeper snowpacks above tree line.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Continuous snow depth, ground interface temperature and shallow soil temperature measurements from 2021-10-1 to 2022-6-14, Seward Peninsula, Alaska

The dataset contains co-located snow depth, ground interface temperature, and shallow soil temperature measured at 98 discrete locations in a watershed located along the Nome-Teller road at mile marker 27 (referred to as T27), and at 53 discrete locations on a hillslope located along the Kougarok road at mile marker 64 (referred to as K64), in Seward Peninsula, Alaska. The dataset aims to understand the local heterogeneity of snow depth, snow temperature, and soil temperature dynamics and their interactions in a discontinuous permafrost region. The dataset is also valuable to train and evaluate machine learning and physical models to predict snow depth or the impact of snow depth on ground surface temperature. At each location, temperatures above and below the ground surface were measured by a pair of vertically deployed distributed temperature profiling probes designed based on Dafflon et al (2022). The probes have high precision temperature sensors spaced at 5 or 10 cm. The mean daily snow depth was estimated by identifying the pair of consecutive sensors with maximum drop of daily temperature high-frequency fluctuations. The ground interface temperature was measured by the sensor located 1-5 cm above the ground surface at 15-minute intervals. The shallow soil temperature was measured by the sensor located 1-5 cm below the ground surface at 15-minute intervals. The dataset includes a description of the probe locations in the "Probe_locations_*.csv" file and the 3 data files (Snow_depths_*.csv, Ground_interface_temperatures_*.csv, Shallow_soil_temperatures_*.csv). * is either T27 or K64, which are the two study sites. In each data file, each column corresponds to a measurement location. Metadata files include data descriptions (_dd.csv) for tabular data. All included files are listed and described in xxxx_flmd.csv. A more detailed description of data processing, along with an updated dataset incorporating multiple seasons and improved snow depth estimation is available at https://doi.org/10.15485/2480365 (Wang et al., 2025a, Wang et al., 2025b).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 OSTI2020

FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials

This paper presents the application of a fast Fourier transform (FFT) based method to solve two phase field models designed to simulate crack growth of strongly anisotropic materials in the brittle regime. By leveraging the ability of the FFT-based solver to generate solutions with higher-order and global continuities, we design two simple algorithms to capture the complex fracture patterns (e.g. sawtooth, and curved crack growth) common in materials with strongly anisotropic surface energy via the multi-phase-field and high-order phase-field frameworks. A staggered operator-split solver is used where both the balance of linear momentum and the phase field governing equations are formulated in the periodic domain. The unit phase field of the initial failure region is prescribed by the penalty method to alleviate the sharp material contrast between the initial failure region and the base material. The discrete frequency vectors are generalized to estimate the second and fourth order gradients such that the Gibbs effect near shape interfaces or jump conditions can be suppressed. Furthermore, a preconditioner is adopted to improve the convergence rate of the iterative linear solver. Three numerical experiments are used to systematically compare the performance of the FFT-based method in the multi-phase-field and high-order phase-field frameworks.

42 ENGINEERING↗
DOE OSTI2021

Review of heat transfer enhancement techniques in two-phase flows for highly efficient and sustainable cooling

We show that two-phase internal flow is ubiquitous to many systems due to its ability to transfer large amounts of heat effectively. Recent advances motivated by sustainability have pushed towards augmenting the heat transfer between the tube and the two-phase working fluid. Augmenting heat transfer in two-phase flows enables process intensification, and compactness leading to reduced cost and material consumption. This review investigates the effect of heat transfer enhancement techniques on different fluids including R134a, R407C, R22, etc … The techniques are categorized into passive and active and focus on durable enhancement methods which do not include surface structuring or coating. Passive techniques utilize geometrical and surface modifications to induce better mixing. Several passive techniques are analyzed including fins, inserts, and more. In analyzing these techniques, the effect of changing geometric parameters is studied, and the limitations are highlighted. Further, active techniques that enable heat transfer enhancement through applying power are also investigated. Electrohydrodynamic approaches are analyzed and compared to passive techniques. Additionally, the implications of heat exchange enhancement in the context of global energy are discussed. An energy analysis discretized by the use sector is conducted, with an emphasis on predicting the future potential for renewable and sustainable development. Moreover, although active approaches provide the highest heat transfer enhancement, due to implementation difficulties, passive techniques are more frequently adopted. The review ends with discussion of next-generation heat exchangers which will rely more on additive manufacturing due to its flexibility and reduced cost per part as the technology matures.

Nusselt number↗
DOE OSTI2022

Development of lean, efficient, and fast physics-framed deep-learning-based proxy models for subsurface carbon storage

In this work, we present deep-learning-based surrogate models for CCUS developed with four different algorithms and a physics-framed two-phase flow problem involving displacement of water by CO 2 . The deep-learning models were trained using 3D datasets describing the pressure plume, CO 2 saturation plume, and water extraction rate generated by numerical simulation. The hyperparameters defining the architecture of the neural networks were optimized to determine the slimmest network size and training parameters that give the most efficient performance at the least training cost. To develop a robust model that closely mimics the governing physical laws, the discretized form of the two-phase fluid transport equation was used to formulate the supervised deep-learning task. The algorithms investigated in this study predicted the data to above 95% accuracy, with the multi-layer perceptron model demonstrating the best performance by balancing training speed, prediction time, and prediction accuracy with lean network capacity. Furthermore, the surrogate models simultaneously predict reservoir pressure and CO 2 saturation in every grid block, including the surface well extraction rate and bottomhole pressure, at all simulation times for a given static model realization in just a few seconds on a standard desktop computer. A key outcome of this study is that limits can be placed on network design parameters to avoid over designing neural networks, with associated efficiencies in training and prediction times. This is very useful because large volumes of data may be generated in CCUS projects and over-design of neural network architectures imposes penalties that are antithetical to the goal of near-real time forecasting.

58 GEOSCIENCES↗
DOE OSTI2020

Performance of Closed Loop Venturi Aspirated Aeration System: Experimental Study and Numerical Analysis with Discrete Bubble Model

In wastewater treatment plants, aeration plays a significant role as it encourages aerobic respiration of microbes, which are necessary to break down carbonaceous matter in the waste stream. This process can account for the majority of energy use in wastewater treatment plants. The aeration process is also necessary in odor control in lagoons and in the aquaculture industry. Generally, the aeration process is accomplished with compressors or blowers which can be of low efficiency due to ideal gas laws. This study introduces the idea of increasing aeration efficiency by looping water from a reservoir through a piping network which includes a venturi aspirator at its inlet. For this purpose, an experimental study has been conducted in a laboratory setup with a pump which pulls water from a tank, passes it through a Venturi aspirator and a helical piping network intended to increase bubble residence time, before depositing it back into the bulk fluid tank. This same process is modeled computationally, using a discrete bubble method (DBM), with good agreement with experiments. The overall purpose here is to determine the optimal configurations for standard aeration efficiency (SAE) and the standard oxygen transfer rate (SOTR). A parametric study has been implemented using the DBM based on different hydraulic and flow parameters. The model is also used to predict the SAE of a hypothetical aeration system. Results indicate that it is possible to achieve SAE values in the range of surface aerators or submersed jet aerators using the proposed aeration system with less complex components, thereby decreasing overall costs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

A soft departure from jamming: the compaction of deformable granular matter under high pressures

Here, the high-pressure compaction of three-dimensional granular packings is simulated using a bonded particle model (BPM) to capture linear elastic deformation. In the model, grains are represented by a collection of point particles connected by bonds. A simple multibody interaction is introduced to control Poisson's ratio and the arrangement of particles on the surface of a grain is varied to model both high- and low-frictional grains. At low pressures, the growth in packing fraction and coordination number follows the expected behavior near jamming and exhibit friction dependence. As the pressure increases, deviations from the low-pressure power-law scaling emerge after the packing fraction grows by approximately 0.1 and results from simulations with different friction coefficients converge. These results are compared to predictions from traditional discrete element method simulations which, depending on the definition of packing fraction and coordination number, may only differ by a factor of two. As grains deform under compaction, the average volumetric strain and asphericity, a measure of the change in the shape of grains, are found to grow as power laws and depend heavily on the Poisson's ratio of the constituent solid. Larger Poisson's ratios are associated with less volumetric strain and more asphericity and the apparent power-law exponent of the asphericity may vary. The elastic properties of the packed grains are also calculated as a function of packing fraction. In particular, we find the Poisson's ratio near jamming is 1/2 but decreases to around 1/4 before rising again as systems densify.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Room temperature corrosion behaviour of plastically deformed AISI 304 stainless steel by friction stir welding in neutral and acidified chloride solutions

Severe plastic deformation of 304L stainless steel plate was carried out by friction stirring under an isothermal tool temperature at 825°C. The friction stirred zone showed extensive grain refinement, high fraction of low angle grain boundaries due to dislocation rearrangement, Σ3 and Σ9 special grain boundaries and discrete Fe–Cr–Mn-rich particles. Electrochemical polarisation measurements were done on the base metal (BM) and friction stir welded (FSW) specimens in the neutral and acidified 3.5% NaCl solutions. The BM showed marginally better pitting resistance in the neutral chloride solution than the FSW because of a larger fraction of special grain boundaries, and larger grain size. Whereas in acidified 3.5% NaCl solution, the FSW showed better corrosion resistance than the BM. Furthermore, the surface film of the FSW had a lower impedance and a lower defect concentration than the BM indicating higher diffusivity of point defects in the FSW.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Status of the Mu2e calorimeter readout electronics

The Mu2e experiment [1] at Fermilab will search for the neutrino-less coherent conversion of a muon into an electron in the field of a nucleus. Mu2e detectors comprise a straw tracker, an electromagnetic calorimeter and a veto for cosmic rays. The calorimeter employs 1348 Cesium Iodide crystals readout by silicon photomultipliers and fast front-end and digitization electronics. The front-end electronics consists of two discrete readout circuits (AMP-HV) for each crystal. These provide the amplification, shaping stage and linear regulation of the SiPM bias voltage and monitoring. The SiPM and front-end control electronics is implemented in a battery of mezzanine boards each equipped with an ARM processor that controls a group of 20 Amp-HV circuits distributing the low voltage and the high-voltage. The electronic is hosted in crates located on the external surface of calorimeter disks. The crates also host the waveform digitizer board (DIRAC) that performs digitization of the front end signals and transmit the digitized data to the Mu2e DAQ. Calorimeter electronic is hosted inside the cryostat and must sustain very high radiation and magnetic field so it was necessary to fully qualify it. The system design and quality assurance procedures will be reviewed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗