A multiple-grid-time-step lattice Boltzmann method for transport phenomena with dissimilar time scales: Application in dendritic solidification
Explore the source record for details and available documents.
Elaheh Dorari↗
SEARCH · Search NASA
Results for “scaling methods”
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.
Intercomparison of methods of coupling between convection and large‐scale circulation: 1. Comparison over uniform surface conditions
Explore the source record for details and available documents.
C. L. Daleu↗
Scaling Building Energy Audits through Machine Learning Methods on Novel Drone Image Data
Building energy audits are time-consuming and labor-intensive. This paper describes a new method using machine learning (ML) techniques on novel data sources (drone images) to improve the identification of building characteristics and retrofit opportunities, and thereby reduce the effort for audits. The new ML method includes: (1) Building footprint extraction using line extraction, polygonization, and polygon-merging, (2) Building envelope extraction using PIX4d modeling software to reconstruct a building 3D model, (3) Visualization tool for viewing images from the 3D model, (4) Window-to-wall ratio (WWR) using state-of-art deep neural network semantic segmentation, (5) Envelope thermal anomaly detection using an unsupervised machine learning clustering algorithm, and (6) Rooftop energy equipment detection based on an object detection algorithm. The testing of this method involved a comparison of additional ML-generated information overlaid on current ‘state-of-practice’ audit and remote assessment baselines using evaluation metrics: labor time and associated cost, marginal benefits of using ML-generated information in workflows for audits and remote assessments, integration potential with existing processes and tools, and replicability/scalability of the method. In two test buildings in California that had comprehensive drawings and meter data available, the ML method effectively generated a building footprint, envelope, rooftop equipment, WWR, and locations of envelope thermal anomalies. Projected target segments of the ML method are sites with minimal drawings and energy data, and underserved sectors such as multistoried housing, disadvantaged communities, and schools for which the ML method can enable identification of building asset characteristics and prioritization of envelope retrofits and decentralized energy equipment retrofits.
Singh, Reshma↗
Application of aerial photography to the study of small scale upper ocean phenomena
The industrial waste dumped 180 n. miles south of Galveston was monitored in July 1977 by water sampling, hydrographic measurements, acoustic tracking on board two vessels, and by aerial photography. The plume of the waste diffused vertically and horizontally. Photodensitometry of aerial photos of the plume showed lateral dispersion of the plume in agreement with two other methods: acoustic tracking of the waste suspensoid and transmissometer sampling. In addition, the method showed small scale features like the lateral and longitudinal variations in the photodensity, indicating the waste concentration. This waste concentration showed periodic changes in its axial distance, with the spectral peak at about 160 m wave length. It shows a sharp increase at the windward edge of the plume as do the acoustic records. This phenomenon is explained in terms of the shearing current near the surface together with vertical diffusion. The periodic change along the axis is explained in terms of the Langmuir circulation and in terms of internal ship waves.
Ichiye, T.↗
Project Methods to Enhance Wellbore Cement Integrity with Microbially-Induced Calcite Precipitation (MICP) (Final Scientific/Technical Report)
The goal of this project was to develop improved methods for sealing compromised wellbore cement in leaking oil and gas wells, thereby reducing the risk of unwanted upward fluid migration. Novel methods for improving wellbore integrity, such as microbially induced calcite precipitation (MICP), can reduce leakage potential, improve the safety of fossil fuel extraction, improve the public perception of hydraulic fracturing, and promote environmentally-prudent unconventional oil and gas development. Microbes, with the urease enzyme, can catalyze the chemical reaction of urea hydrolysis to induce the precipitation of calcium carbonate which can be used as a cementitious material to seal leakage pathways. In this project, methods to promote robust bio-composite cementitious materials were designed and tested in the laboratory. Scale-up of those methods were tested in meso-scale reactor systems and in field applications. In this report, in Section One, we describe laboratory efforts to develop injection strategies to promote precipitation in wellbore analogs and determine the strength of the bio-composite cements as compared to fine cement. In Section Two, we describe the efforts to scale up the work and study the use of materials that can be used in field application, for example exploring the use of calcium chloride ice melt or urea fertilizer as source chemicals. In Section Three, the three field trials (methods and results) performed as part of the project are described and summarized. At the end of the report is a comprehensive summary and conclusion section which highlights the key findings of the project. The work performed during this project significantly advanced the technology readiness level (TRL) of the MICP wellbore sealing strategy.
02 PETROLEUM↗
Methods to Enhance Wellbore Cement Integrity with Microbially-Induced Calcite Precipitation (MICP) (Final Scientific/ Technical Report)
The goal of this project was to develop improved methods for sealing compromised wellbore cement in leaking oil and gas wells, thereby reducing the risk of unwanted upward fluid migration. Novel methods for improving wellbore integrity, such as microbially induced calcite precipitation (MICP), can reduce leakage potential, improve the safety of fossil fuel extraction, improve the public perception of hydraulic fracturing, and promote environmentally-prudent unconventional oil and gas development. Microbes, with the urease enzyme, can catalyze the chemical reaction of urea hydrolysis to induce the precipitation of calcium carbonate which can be used as a cementitious material to seal leakage pathways. In this project, methods to promote robust bio-composite cementitious materials were designed and tested in the laboratory. Scale-up of those methods were tested in meso-scale reactor systems and in field applications. In this report, in Section One, we describe laboratory efforts to develop injection strategies to promote precipitation in wellbore analogs and determine the strength of the bio-composite cements as compared to fine cement. In Section Two, we describe the efforts to scale up the work and study the use of materials that can be used in field application, for example exploring the use of calcium chloride ice melt or urea fertilizer as source chemicals. In Section Three, the three field trials (methods and results) performed as part of the project are described and summarized. At the end of the report is a comprehensive summary and conclusion section which highlights the key findings of the project. The work performed during this project significantly advanced the technology readiness level (TRL) of the MICP wellbore sealing strategy.
03 NATURAL GAS↗
Data‐driven variational method for discrepancy modeling: Dynamics with small‐strain nonlinear elasticity and viscoelasticity
Abstract The effective inclusion of a priori knowledge when embedding known data in physics‐based models of dynamical systems can ensure that the reconstructed model respects physical principles, while simultaneously improving the accuracy of the solution in the previously unseen regions of state space. This paper presents a physics‐constrained data‐driven discrepancy modeling method that variationally embeds known data in the modeling framework. The hierarchical structure of the method yields fine scale variational equations that facilitate the derivation of residuals which are comprised of the first‐principles theory and sensor‐based data from the dynamical system. The embedding of the sensor data via residual terms leads to discrepancy‐informed closure models that yield a method which is driven not only by boundary and initial conditions, but also by measurements that are taken at only a few observation points in the target system. Specifically, the data‐embedding term serves as residual‐based least‐squares loss function, thus retaining variational consistency. Another important relation arises from the interpretation of the stabilization tensor as a kernel function, thereby incorporating a priori knowledge of the problem and adding computational intelligence to the modeling framework. Numerical test cases show that when known data is taken into account, the data driven variational (DDV) method can correctly predict the system response in the presence of several types of discrepancies. Specifically, the damped solution and correct energy time histories are recovered by including known data in the undamped situation. Morlet wavelet analyses reveal that the surrogate problem with embedded data recovers the fundamental frequency band of the target system. The enhanced stability and accuracy of the DDV method is manifested via reconstructed displacement and velocity fields that yield time histories of strain and kinetic energies which match the target systems. The proposed DDV method also serves as a procedure for restoring eigenvalues and eigenvectors of a deficient dynamical system when known data is taken into account, as shown in the numerical test cases presented here.
Masud, Arif↗
Large-Scale High-Resolution Coastal Mangrove Forests Mapping Across West Africa With Machine Learning Ensemble and Satellite Big Data
Coastal mangrove forests provide important ecosystem goods and services, including carbon sequestration, biodiversity conservation, and hazard mitigation. However, they are being destroyed at an alarming rate by human activities. To characterize mangrove forest changes, evaluate their impacts, and support relevant protection and restoration decision making, accurate and up-to-date mangrove extent mapping at large spatial scales is essential. Available large-scale mangrove extent data products use a single machine learning method commonly with 30 m Landsat imagery, and significant inconsistencies remain among these data products. With huge amounts of satellite data involved and the heterogeneity of land surface characteristics across large geographic areas, finding the most suitable method for large-scale high-resolution mangrove mapping is a challenge. The objective of this study is to evaluate the performance of a machine learning ensemble for mangrove forest mapping at 20 m spatial resolution across West Africa using Sentinel-2 (optical) and Sentinel-1 (radar) imagery. The machine learning ensemble integrates three commonly used machine learning methods in land cover and land use mapping, including Random Forest (RF), Gradient Boosting Machine (GBM), and Neural Network (NN). The cloud-based big geospatial data processing platform Google Earth Engine (GEE) was used for pre-processing Sentinel-2 and Sentinel-1 data. Extensive validation has demonstrated that the machine learning ensemble can generate mangrove extent maps at high accuracies for all study regions in West Africa (92%–99% Producer’s Accuracy, 98%–100% User’s Accuracy, 95%–99% Overall Accuracy). This is the first-time that mangrove extent has been mapped at a 20 m spatial resolution across West Africa. The machine learning ensemble has the potential to be applied to other regions of the world and is therefore capable of producing high-resolution mangrove extent maps at global scales periodically.
coastal environment↗
Guest-Host Interactions in Clathrate Hydrates: Benchmark MP2 and CCSD(T)/CBS Binding Energies of CH4, CO2 and H2S in (H2O)20 Cages
We present benchmark binding energies of naturally occurring gas molecules CH4, CO2, and H2S in the small cage, namely the pentagonal dodecahedron (512) (H2O)20, which is one of the constituent cages of the 3 major lattices (structures I, II and H) of clathrate hydrates. These weak interactions require higher levels of electron correlation and converge slowly with increasing basis set to the Complete Basis Set (CBS) limit, necessitating the use of large basis sets up to the augcc- pV5Z and subsequent correction for Basis Set Superposition Error (BSSE). For the host hollow (H2O)20 cages we have identified a most stable isomer with binding energy of -200.8 ± 2.1 kcal/mol at the CCSD(T)/CBS limit (-199.2 ± 0.5 kcal/mol at the MP2/CBS limit). Additionally, we report converged second order Moller-Plesset (MP2) CBS binding energies for the encapsulation of guests in the (H2O)20 cage of -4.3 ± 0.1 for CH4@(H2O)20, -6.6 ± 0.1 for CO2@(H2O)20 and -8.5 ± 0.1 kcal/mol for H2S@(H2O)20, respectively. For CH4@(H2O)20, exhibiting the weakest encapsulation affinity among the three, we report CCSD(T)/aug-cc-pVTZ binding energies and, based on them, a CCSD(T)/CBS estimate of -4.75 ± 0.1 kcal/mol. To the best of our knowledge, the CCSD(T)/aug-cc-pVTZ calculation for CH4@(H2O)20 is the largest one reported to date (168 valence electrons, 1978 basis functions and the correlation of 84 doubly occupied and 1873 virtual orbitals) and required a scalable implementation of the (T) module on 6144 nodes (350208 cores) of the “Cori” supercomputer at the National Energy Research Supercomputing Center (NERSC) for a total execution time of 195 minutes (for the (T) part). These efficient scalable implementations of highly correlated methods offer the capability to obtain long-lasting benchmarks of intermolecular interactions in complex systems. They also provide a path towards parametrizing classical potentials needed to study the dynamical and transport properties in these complex systems as well as assess the accuracy of lower scaling electronic structure methods such as Density Functional Theory (DFT) and MP2 including its spin-biased variants.
Heindel, Joseph↗
Study of self-interaction-errors in barrier heights using locally scaled and Perdew-Zunger self-interaction methods
We study the effect of self-interaction errors on the barrier heights of chemical reactions. For this purpose we use the well-known Perdew-Zunger [J. P. Perdew and A. Zunger, Phys. Rev. B,23, 5048 (1981)] self-interaction-correction(PZSIC), as well two variations of the recently developed, locally scaled self-interaction correction (LSIC) [Zope et al.,J. Chem. Phys.151, 214108 (2019)] to study the barrier heights of the BH76 benchmark dataset. Our results show that both PZSIC and especially the LSIC methods improve the barrier heights relative to the local density approximation(LDA). The version of LSIC that uses the iso-orbital indicator z as a scaling factor gives a more consistent improvement than an alternative version that uses an orbital-dependent factor based on the ratio of orbital densities to the total electron density. We show that LDA energies evaluated using the self-consistent and self-interaction-free PZSIC densities can be used to assess density-driven errors. Furthermore, the LDA reaction barrier errors for the BH76 set are found to contain significant density-driven errors for all types of reactions contained in the set, but the corrections due to adding SIC to the functional are much larger than those stemming from the density for the hydrogen transfer reactions and of roughy equal size for the non-hydrogen transfer reactions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Measuring Gas Transport and Sorption in Large Intact Geologic Specimens via the Piezometric Method
Uptake of noble gases into heterogeneous geologic core samples was measured using a piezometric methodology. In addition to measuring accessible porosity—as with gas pycnometry—by monitoring the rate of pressure decay, this method can also be used to estimate the gas effective diffusivity in the sample. In contrast to previous applications of this method, where milligram quantities of fractured grains are characterized, here approximately kilogram core samples were left intact when tested. In doing so, a more representative sample of the heterogeneous field geology is provided. Additionally, alteration of the pore structure and connectivity during sample preparation is avoided. Furthermore, to scale the piezometric method from milligrams to kilograms, the system was designed to operate at medium vacuum (1 to 100 Pa) to restrict transport in pores less than approximately 60 µm to large Knudsen numbers. To test the system performance, two samples of interest were selected: a rhyolitic welded tuff from Blue Canyon Dome at the Energetic Materials Research and Testing Center and a zeolitized non-welded rhyolitic tuff from the Nevada National Security Site. Three noble gases were utilized in this test series; Argon and xenon as they are of direct interest to nuclear monitoring efforts and helium as it is a weakly adsorbing reference standard. Additionally, mercury intrusion porosimetry measurements were made on subsamples of the core to compare the observed porosity by the two methods and to discuss gas transport rates in the context of the measured pore distribution.
58 GEOSCIENCES↗
Area-Preserving Mapping of 3D Carotid Ultrasound Images Using Density-Equalizing Reference Map
Carotid atherosclerosis is a focal disease at the bifurcations of the carotid artery. To quantitatively monitor the local changes in the vessel-wall-plus-plaque thickness (VWT) and compare the VWT distributions for different patients or for the same patients at different ultrasound scanning sessions, a mapping technique is required to adjust for the geometric variability of different carotid artery models. In this work, we propose a novel method called density-equalizing reference map (DERM) for mapping 3D carotid surfaces to a standardized 2D carotid template, with an emphasis on preserving the local geometry of the carotid surface by minimizing the local area distortion. The initial map was generated by a previously described arc-length scaling (ALS) mapping method, which projects a 3D carotid surface onto a 2D non-convex L-shaped domain. A smooth and area-preserving flattened map was subsequently constructed by deforming the ALS map using the proposed algorithm that combines the density-equalizing map and the reference map techniques. This combination allows, for the first time, one-to-one mapping from a 3D surface to a standardized non-convex planar domain in an area-preserving manner. Evaluations using 20 carotid surface models show that the proposed method reduced the area distortion of the flattening maps by over 80% as compared to the ALS mapping method.
60 APPLIED LIFE SCIENCES↗
In-Process Monitoring and Structural Health Monitoring of Large-Scale Additive Manufacturing Using Acoustic Emission Technique
ORNL collaborated with MISTRAS Group, Inc. to investigate acoustic emission (AE) as a structural health monitoring (SHM) method for large-scale additive manufacturing (AM). Large-scale AM is being adapted as method of producing large structures in a short lead time and cost-effective way. With the growing advancement in AM techniques and application, machine monitoring and part qualification is highly needed. There has been leading research focused on the manufacturing, feedstock material but minimum research on the SHM, defect detection, and nondestructive evaluation (NDE) for AM. Scanning large structure using conventional nondestructive testing (NDT) techniques, such as ultrasound or X-ray, and searching for potential defects can be very time consuming, challenging and cost prohibitive. AE is a passive technique that can be used to monitor and locate defect progression in large structure by distributing group of sensors around the part. This project utilized AE technique and system manufactured/designed by MISTRAS Group to monitor large-scale AM equipment (i.e. Big Area Additive Manufacturing (BAAM) system located at the Oak Ridge National Laboratory – Manufacturing Demonstration Facility (ORNL-MDF) and the printed parts it produces. The AE system provided valuable insight on defect development/progression during and post-printing process.
36 MATERIALS SCIENCE↗
Progressive Damage Response of 3D Woven Composites via the Multiscale Recursive Micromechanics Solution with Tailored Fidelity
Progressive failure simulations have been performed for orthogonal 3D woven composites consisting of RTM6 resin matrix and AS4 carbon fibers. The Multiscale Recursive Micromechanics approach has been used, which, while being computationally efficient, captures the primary effects of the microstructure at each considered length scale. This approach also enables use of any micromechanics theory at any length scale, and herein, the fidelity of the chosen theories across the scales has been tailored to strike a balance with computational efficiency. The Mori-Tanaka method is employed at the lowest length scale, the Generalized Method of Cells is used at intermediate scales, and the High-Fidelity Generalized Method of Cells is used at the highest woven composite repeating unit cell scale. Furthermore, two different damage models, also with different levels of fidelity and efficiency, have been used for the resin material at the lowest length scale. Results for the mechanical behavior in response to loading in various directions are compared for the two damage models and with available test data.
3D woven↗
Optimizing the optimizer for physics-informed neural networks and Kolmogorov-Arnold networks
Physics-Informed Neural Networks (PINNs) have revolutionized the computation of PDE solutions by integrating partial differential equations (PDEs) into the neural network’s training process as soft constraints, becoming an important component of the scientific machine learning (SciML) ecosystem. More recently, physics-informed Kolmogorv-Arnold networks (PIKANs) have also shown to be effective and comparable in accuracy with PINNs. In their current implementation, both PINNs and PIKANs are mainly optimized using first-order methods like Adam, as well as quasi-Newton methods such as BFGS and its low-memory variant, L-BFGS. However, these optimizers often struggle with highly nonlinear and non-convex loss landscapes, leading to challenges such as slow convergence, local minima entrapment, and (non)degenerate saddle points. In this study, we investigate the performance of Self- Scaled BFGS (SSBFGS), Self-Scaled Broyden (SSBroyden) methods and other advanced quasi-Newton schemes, including BFGS and L-BFGS with different line search strategies. These methods dynamically rescale updates based on historical gradient information, thus enhancing training efficiency and accuracy. We systematically compare these optimizers – using both PINNs and PIKANs – on key challenging PDEs, including the Burgers, Allen-Cahn, Kuramoto-Sivashinsky, Ginzburg-Landau, and Stokes equations. Additionally, we evaluate the performance of SSBFGS and SSBroyden for Deep Operator Network (DeepONet) architectures, demonstrating their effectiveness for data-driven operator learning. Our findings provide state-of-the-art results with orders-of-magnitude accuracy improvements without the use of adaptive weights or any other enhancements typically employed in PINNs. More broadly, our work reveal insights into the effectiveness of quasi-Newton optimization strategies in significantly improving the convergence and accurate generalization of PINNs and PIKANs.
97 MATHEMATICS AND COMPUTING↗
Spectral Methods for Computational Fluid Dynamics
As a tool for large-scale computations in fluid dynamics, spectral methods were prophesized in 1944, born in 1954, virtually buried in the mid-1960's, resurrected in 1969, evangalized in the 1970's, and catholicized in the 1980's. The use of spectral methods for meteorological problems was proposed by Blinova in 1944 and the first numerical computations were conducted by Silberman (1954). By the early 1960's computers had achieved sufficient power to permit calculations with hundreds of degrees of freedom. For problems of this size the traditional way of computing the nonlinear terms in spectral methods was expensive compared with finite-difference methods. Consequently, spectral methods fell out of favor. The expense of computing nonlinear terms remained a severe drawback until Orszag (1969) and Eliasen, Machenauer, and Rasmussen (1970) developed the transform methods that still form the backbone of many large-scale spectral computations. The original proselytes of spectral methods were meteorologists involved in global weather modeling and fluid dynamicists investigating isotropic turbulence. The converts who were inspired by the successes of these pioneers remained, for the most part, confined to these and closely related fields throughout the 1970's. During that decade spectral methods appeared to be well-suited only for problems governed by ordinary diSerential eqllations or by partial differential equations with periodic boundary conditions. And, of course, the solution itself needed to be smooth. Some of the obstacles to wider application of spectral methods were: (1) poor resolution of discontinuous solutions; (2) inefficient implementation of implicit methods; and (3) drastic geometric constraints. All of these barriers have undergone some erosion during the 1980's, particularly the latter two. As a result, the applicability and appeal of spectral methods for computational fluid dynamics has broadened considerably. The motivation for the use of spectral methods in numerical calculations stems from the attractive approximation properties of orthogonal polynomial expansions.
Zang, T. A.↗
Irrigation Characterization Improved by the Direct Use of SMAP Soil Moisture Anomalies Within a Data Assimilation System
Prior soil moisture data assimilation (DA) efforts to incorporate human management features such as agricultural irrigation has only shown limited success. This is partly due to the fact that observational rescaling approaches for bias correction used in soil moisture DA systems are less effective when unmodeled processes such as irrigation are the dominant source of systematic biases. In this article, we demonstrate an alternative approach, i.e. anomaly correction for overcoming this limitation. Unlike the rescaling approaches, the proposed method does not scale remote sensing soil moisture retrievals to the model climatology, but it extracts the temporal variability information from the retrievals. The study demonstrates this approach through the assimilation of soil moisture retrievals from the Soil Moisture Active Passive mission into the Noah land surface model. The results demonstrate that DA using the anomaly correction method can better capture the effect of irrigation on soil moisture in agricultural areas while providing comparable performance to the DA integrations using rescaling approaches in non-irrigated areas. These findings emphasize the need to reduce inconsistencies between remote sensing and the models so that assimilation methods can employ information from remote sensing more directly to develop representations of unmodeled processes such as irrigation.
soil moisture↗
Drop clustering and drop size correlations from holographic imagery suggest cloud droplet spectral broadening via entrainment-mixing
The question of how droplets rapidly grow large enough to initiate collision-coalescence has persisted for decades. Many theories explaining the production of sufficiently large drops (i.e., those in the “bottleneck” size range; ∼ 25–50 µm diameters) involve drop clustering on millimeter scales. A novel method is introduced to evaluate drop clustering trends particle-by-particle (i.e., the number/proximity of neighboring drops for given droplets; defined as drop clustering fields) which are diagnosed relative to drops within their shared drop environments – in contrast to previous studies which diagnose drop clustering of defined sample volumes, or in terms of absolute length scales. Specifically, this study evaluates the statistical likelihood that drops of a given size are associated with either a significant number of neighboring drops, or are significantly isolated from neighboring drops. Observations are acquired from the HOLODEC during the Cloud System Evolution in the Trades campaign, which sampled subtropical marine clouds. The HOLODEC measures drop size distributions and the 3D spatial coordinates of droplets. Results show drops within the bottleneck size range (diameters of ∼ 25–50 µm) are most likely to be significantly isolated from neighboring drops. This “isolated large drop trend” is primarily observed at subsaturated conditions, suggesting entrainment is the contributing factor. Holograms associated with this trend are more likely to have broader drop size distributions, larger maximum drop sizes and overly regions where precipitation reaches the lowest altitudes from the sampled cloud, suggesting entrainment-mixing drop size distribution broadening is a relevant precipitation-initiation mechanism.
D'Alessandro, John J. [Univ. of Washington, Seattl↗