Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Nonadiabatic Dynamics in a Continuous Circularly Polarized Laser Field with Floquet Phase-Space Surface Hopping

Non-adiabatic chemical reactions involving continuous circularly polarized light (cw CPL) have not attracted as much attention as dynamics in unpolarized/linearly polarized light. However, including circularly (in contrast to linearly) polarized light allows one to effectively introduce a complex-valued time-dependent Hamiltonian, which offers a new path for control or exploration through the introduction of Berry forces. Here, we investigate several inexpensive semiclassical approaches for modeling such nonadiabatic dynamics in the presence of a time-dependent complex-valued Hamiltonian, beginning with a straightforward instantaneous adiabatic fewest-switches surface hopping (IA-FSSH) approach (where the electronic states depend on position and time), continuing to a standard Floquet fewest switches surface hopping (F-FSSH) approach (where the electronic states depend on position and frequency), and ending with an exotic Floquet phase-space surface hopping (F-PSSH) approach (wherethe electronic states depend on position, frequency, and momentum). Using a set of model systems with time-dependent complex-valued Hamiltonians, we show that the Floquet phase- space adiabats are the optimal choice of basis as far as accounting for Berry phase effects and delivering accuracy. Furthermore, the F-PSSH algorithm sets the stage for modeling nonadiabatic dynamics under strong externally pumped circular polarization in the future.

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

Multirate partitioned Runge–Kutta methods for coupled Navier–Stokes equations

Earth system models are complex integrated models of atmosphere, ocean, sea ice, and land surface. Coupling the components can be a significant challenge due to the difference in physics, temporal, and spatial scales. Further, this study explores multirate partitioned Runge-Kutta methods for the fluid-fluid interaction problem and demonstrates its parallel performance by using the PETSc library. We consider compressible Navier-Stokes equations with gravity coupled through a rigid-lid interface. Our large-scale numerical experiments reveal that multirate partitioned Runge-Kutta coupling schemes (1) can conserve total mass; (2) have second-order accuracy in time; and (3) provide favorable strong- and weak-scaling performance on modern computing architectures. We also show that the speedup factors of multirate partitioned Runge-Kutta methods match theoretical expectations over their base (single-rate) method.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

nn-PINNs: Non-Newtonian physics-informed neural networks for complex fluid modeling

Time- and rate-dependent material functions in non-Newtonian fluids in response to different deformation fields pose a challenge in integrating different constitutive models into conventional computational fluid dynamic platforms. Considering their relevance in many industrial and natural settings alike, robust data-driven frameworks that enable accurate modeling of these complex fluids are of great interest. The main goal is to solve the coupled Partial Differential Equations (PDEs) consisting of the constitutive equations that relate the shear stress to the deformation and fully capture the behavior of the fluid under various flow protocols with different boundary conditions. Here, in this work, we present non-Newtonian physics-informed neural networks (nn-PINNs) for solving systems of coupled PDEs adopted for complex fluid flow modeling. The proposed nn-PINN method is employed to solve the constitutive models in conjunction with conservation of mass and momentum by benefiting from Automatic Differentiation (AD) in neural networks, hence avoiding the mesh generation step. nn-PINNs are tested for a number of different complex fluids with different constitutive models and for several flow protocols. These include a range of Generalized Newtonian Fluid (GNF) empirical constitutive models, as well as some phenomenological models with memory effects and thixotropic timescales. nn-PINNs are found to obtain the correct solution of complex fluids in spatiotemporal domains with good accuracy compared to the ground truth solution. We also present applications of nn-PINNs for complex fluid modeling problems with unknown boundary conditions on the surface, and show that our approach can successfully recover the velocity and stress fields across the domain, including the boundaries, given some sparse velocity measurements.

42 ENGINEERING↗
DOE OSTI2023

Grounding our Understanding of the Impacts of Boreal Forest Expansion on Shallow Cumulus Clouds with a Simple Modeling Framework

Abstract The expansion of the boreal forest poleward is a potentially important driver of feedbacks between the land surface and Arctic climate. A growing body of work has highlighted the importance of differences in evaporative resistance between different possible future Arctic land covers, which in turn alters humidity and cloudiness in the boundary layer, for these feedbacks. While thus far this problem has been studied primarily with complex Earth system models, we turn to a locally focused, idealized model capable of diagnosing and testing the sensitivity of first-order processes connecting vegetation, the atmospheric boundary layer, and low clouds in this critical region. This allows us to benchmark the mechanisms and results at the center of predictions from larger-scale simulations. A surface dominated by broadleaf trees, characterized by higher albedo and lower surface evaporative resistance, drives cooling and moistening of the boundary layer relative to a surface of needleleaf trees, characterized by lower albedo and higher surface evaporative resistance. Differences in evaporative resistance between these hypothetical Arctic vegetation covers are of equal importance to changes in albedo for the initial response of the boundary layer to boreal expansion, even with our idealized approach. However, compensation between the elevation of the lifting condensation level (LCL) and more rapid growth of the mixed layer over higher evaporative resistance surfaces can minimize changes in the favorability of shallow clouds over different land cover types under some conditions. We then perform two tests on the sensitivity of this compensating effect, to changes in water availability, represented first by a reduction in boundary layer humidity and then by both a reduction in humidity and soil moisture available to our vegetation surface. Finally, given the importance of this potential LCL–mixed-layer height compensation in our idealized modeling results, we look to determine its relevance in observational data from a field campaign in boreal Finland. These observations do confirm that such a coupling plays an important role in cumulus-topped boundary layers over a needleleaf forest surface. While our results confirm some underlying mechanisms at the center of prior work with Earth system models, they also provide motivation for future work to constrain the impact of boreal forest expansion. This will include both large eddy simulations to examine the impact of processes and feedbacks not resolved by a mixed-layer model, as well as a more systematic evaluation and comparison of relevant observations at the site in Finland and sites from prior boreal field campaigns. Significance Statement Clouds and vegetation are both important components of the climate system that interact across a range of scales. These interactions are central to understanding how changes at the land surface feedback on climate. For example, if a forest expands or recedes, diagnosing how that will impact clouds will determine whether you predict warming or cooling temperatures from that shift in the forest area. These predictions are often made with complex Earth system models, but we look to a more idealized representation of the land–atmosphere system to diagnose how shallow clouds should respond to changes in surface properties with different scenarios of boreal forest expansion at a more foundational level. This both grounds our understanding of previous analysis and provides helpful direction for future studies of this relevant and impactful land cover change.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Surface Atmosphere Integrated Field Laboratory (SAIL) Science Plan

Mountains are the natural water towers of the world, effectively turning water vapor into readily available fresh water through precipitation, snowpack, and runoff. Unfortunately, Earth system models (ESMs) have persistently been unable to predict the timing and availability of water resources from mountains because the source(s) of model error are difficult to isolate in complex terrain with limited atmospheric or land-surface observations. Further complications arise from the gross scale mismatch between ESM grid box sizes and the relevant scales of mountainous hydrological processes. The mountain hydrometeorology community has repeatedly called for integrated atmospheric and land observations of water and energy budgets in complex terrain that span these scales to establish benchmarks against which scale-dependent models can be further developed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Soil Carbonyl Sulfide (OCS) Fluxes in Terrestrial Ecosystems: An Empirical Model

Measurements of carbonyl sulfide (OCS) enable independent estimates of regional stomatal conductance provided that non-stomatal OCS fluxes are well constrained. OCS is taken up through plant leaves, following the same pathway as CO 2 ; in contrast to CO 2 , OCS is irreversibly destroyed in plant leaves and plants do not typically exhibit OCS emissions. Ecosystem uptake of OCS can indicate changes in stomatal opening. Here we present an empirical model to assess the potential impact of soil OCS exchange, the non-Stomatal OCS exchange Empirical Model (SOCSEM, version 0). We created biome-specific response curves characterizing soil OCS exchange and restricted the model design to require only knowledge of soil moisture and surface temperature because remote sensing observations are available for these two features. Comparing the model to field-based chamber observations reveal deviations that can be attributed to missing complexity of the ground surface (having excluded litter and plants without regulated stomata), shortwave radiation, or the soil environment. For agricultural regions with known net emissions, we use remotely-sensed surface temperature data and demonstrate that data resolution can affect anticipated fluxes. We further investigate the influence of regions with unknown soil OCS responses, for example, Arctic tundra. We compare our model to a process-based and respiration-based soil OCS exchange model that has been implemented in a land surface model. Further field study of tropical and arctic ecosystems in conjunction with studies of non-stomatal surfaces in addition to soil (e.g., bryophytes) will increase confidence in applying OCS as a regional tracer for stomatal conductance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗
DOE OSTI2023

X-ray Absorption Spectroscopy Studies of a Molecular CO 2 -Reduction Catalyst Deposited on Graphitic Carbon Nitride

Metal-ligand complexes have been extensively explored as well-defined molecular catalysts in small molecule activation reactions such as carbon dioxide (CO 2 ) reduction. Many hybrid photocatalysts have been prepared by coupling such complexes with photoactive surfaces for use in solar CO 2 reduction. In this work, we employ X-ray absorption near edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) spectroscopies, density functional theory (DFT) and computational XANES modeling to interrogate the structure of a hybrid photocatalyst consisting of a macrocyclic cobalt complex deposited on graphitic carbon nitride (C 3 N 4 ). Results show that the cobalt complex binds on C 3 N 4 through surface OH or NH 2 groups. By refining the local geometry and binding sites of this well-defined molecular cobalt complex on C 3 N 4 , here we established an important benchmark for modeling a large class of molecular catalysts that can be adapted to in situ/operando studies and further enhanced by applying chemometrics-based approaches and machine learning methods of XANES data analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Perspectives on Machine Learning-assisted Plasma Medicine: Towards Automated Plasma Treatment

Cold atmospheric plasmas (CAPs) have shown great promise for medical applications through their synergistic chemical, electrical, and thermal effects, which can induce therapeutic outcomes. However, safe and reproducible plasma treatment of complex biological surfaces poses a major hurdle to the widespread adoption of CAPs for medical applications. Predictive modeling of the mutual interactions between the plasma and biological surfaces and, thus, systematic approaches to quantify and predict plasma treatment outcomes remain largely elusive due to the lack of mechanistic understanding of plasma-surface interactions that can span across vastly different length-and time-scales. In addition, real-time sensing capabilities in biomedical CAP devices are often limited, which can be detrimental to plasma treatment due to the intrinsic plasma and surface variability during the treatment, as well as sensitivity to external perturbations. All of these challenges can make reproducible and effective plasma treatment of biological surfaces difficult to realize, which is further compounded by errors due to human operation of hand-held CAP devices. Machine learning and data-driven approaches can be particularly useful in addressing these challenges in three major ways: (i) data-driven modeling of hard-to-model plasma-surface interactions and plasma treatment outcomes; (ii) learning data analytics for plasma and surface diagnostics in real-time; and (iii) developing predictive controllers that enable reliable and effective CAP treatments. Furthermore, this paper discusses the promise of machine learning to accelerate plasma medicine research in these areas, toward machine learning-assisted and automated CAP treatment of complex biological surfaces.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Characterizing the Small-Scale Dynamical, Ice Microphysical, and Residual Aerosol Properties of Mid-Latitude Cold Clouds: A Pilot Study (Final Report)

Our understanding of cirrus clouds has been advancing rapidly due to improvements in measurements made inside clouds, computer cloud models, and laboratory experimentation. This advancement in knowledge has made it clear that cirrus are complicated systems where the intricate interactions among dynamic air motions, solar and terrestrial radiation, ice crystal growth, and water and air chemistry can substantially influence cirrus structure and evolution. Measurements made inside clouds have shown that cirrus crystals are often complex with shapes that are often unclassifiable, especially at small sizes. Laboratory measurements have revealed that the ice crystal surface is complex, even at small sizes. Recent laboratory and satellite studies show that this surface complexity modifies the radiative properties of cirrus clouds in important ways. Moreover, observations and computer cloud modeling work show that the larger and smaller scale properties of cirrus are sensitive to the interplay between small-scale air motions and the aerosol articles upon which cirrus ice crystals form. Taken together, these studies indicate that accurately accounting for the radiative properties of cirrus clouds requires knowing the air motions and aerosol properties that influence the number concentration and crystal size, along with an understanding of the surface properties of vapor grown ice. Although this prior work has been illuminating, there have been few studies that provide detailed shape and surface complexity measurements for the small ice crystals that often compose cirrus clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Two-phase flow numerical analysis of electrode geometry for alkaline water electrolyzers

Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.

Alkaline↗
DOE OSTI2025

A neural master equation framework for multiscale modeling of molecular processes: application to atomic-scale plasma processes

Plasma-surface interactions (PSI) play a crucial role in microelectronics fabrication; however, their multiscale nature and array of complex, often unknown interactions make computational modeling of PSIs extremely difficult. To this end, we propose a general neural master equation (NME) framework that uses master equations to describe the dynamics of a molecular process, wherein neural networks learned from atomistic simulations represent unknown transitions between different system states. By leveraging the physics-based structure of master equations and data-driven state transitions, the NME framework promotes generalizability and physics interpretability, and can bridge disparate length and time scales. The framework is demonstrated for multiscale modeling of Si atomic layer etching and reactive ion etching, where the learned NME-based surface kinetic models exhibit good predictive and extrapolative capabilities for predicting experimentally relevant observables as a function of process parameters. The NME-based surface kinetic models obey physical constraints, which are violated in models based on neural ordinary differential equations. The proposed NME framework for multiscale modeling of molecular processes can pave the way for the discovery of new chemistries and materials in atomic-scale plasma processes.

Chemical engineering↗
DOE OSTI2020

Potential-pH diagrams considering complex oxide solution phases for understanding aqueous corrosion of multi-principal element alloys

The potential-pH diagram, a graphical representation of the thermodynamically predominant reaction products in aqueous corrosion, is originally proposed for the corrosion of pure metals. The original approach only leads to stoichiometric oxides and hydroxides as the oxidation products. However, numerous experiments show that non-stoichiometric oxide scales are prevalent in the aqueous corrosion of alloys. In the present study, a room temperature potential-pH diagram considering oxide solid solutions, as a generalization of the traditional potential-pH diagram with stoichiometric oxides, is constructed for an FCC single-phase multi-principal element alloy (MPEA) based on the CALculation of PHAse Diagram method. The predominant reaction products, the ions in aqueous solution, and the cation distribution in oxides are predicted. The oxide solid solution is stabilized by the mixing free energy (or mixing entropy) and the stabilizing effect becomes more significant as the temperature increases. Consequently, solid solution oxides are stable in large regions of the potential-pH diagram and the mixing free energy mostly affects the equilibrium composition of the stable oxides, while the shape of stable regions for oxides is mostly determined by the structure of the stable oxides. Agreements are found for Ni 2+ , Fe 2+ , and Mn 2+ between the atomic emission spectroelectrochemistry measurements and thermodynamic calculations, while deviations exist for Cr 3+ and Co 2+ possibly due to surface complexation with species such as Cl - and the oxide dissolution. By incorporating the solution models of oxides, the current work presents a general and more accurate way to analyze the reaction products during aqueous corrosion of MPEAs.

Materials Science↗
DOE OSTI2022

Hardware-Accelerated Ray Tracing of CAD-Based Geometry for Monte Carlo Radiation Transport

Monte Carlo radiation transport (MCRT) methods have been used to simulate radiation environments for many decades by tracking individual particles through a model to accumulate statistical information. MCRT geometry is historically formed using the constructive solid geometry (CSG). Recently, significant work has been performed to support simulations using computer-aided design (CAD)-based tessellated surfaces to support highly complex geometries. Ray tracing acceleration data structures from the rendering and visualization community are applied to accelerate particle tracking in CAD-based models. Despite these efforts, CSG representations provide the superior performance in surface intersection operations during particle flight. Concurrently, pseudo Monte Carlo methods have become prevalent in rendering applications to support more realistic models for scattering media, motivating innovations that are advantageous for MCRT simulations. Finally, the authors’ work extends these innovations by employing Intel’s Embree ray tracing kernel within a geometry toolkit for Monte Carlo to improve the simulation performance using CAD-based models by factors of 1.5 to 2.

42 ENGINEERING↗
DOE OSTI2025

A comparison of surrogate constitutive models for viscoplastic creep simulation of HT-9 steel

Mechanistic microstructure-informed constitutive models for the mechanical response of polycrystals are a cornerstone of computational materials science. However, as these models become increasingly more complex – often involving coupled differential equations describing the effect of specific deformation modes – their associated computational costs can become prohibitive, particularly in optimization or uncertainty quantification tasks that require numerous model evaluations. To address this challenge, surrogate constitutive models that balance accuracy and computational efficiency are highly desirable. Data-driven surrogate models, that learn the constitutive relation directly from data, have emerged as a promising solution. In this work, we develop two local surrogate models for the viscoplastic response of a steel: a piecewise response surface method and a mixture of experts model. These surrogates are designed to adapt to complex material behavior, which may vary with material parameters or operating conditions. The surrogate constitutive models are applied to creep simulations of HT-9 steel, an alloy of considerable interest to the nuclear energy sector due to its high tolerance to radiation damage, using training data generated from viscoplastic self-consistent (VPSC) simulations. In conclusion, we define a set of test metrics to numerically assess the accuracy of our surrogate models for predicting viscoplastic material behavior, and show that the mixture of experts model outperforms the piecewise response surface method in terms of accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Hector V3.2.0: functionality and performance of a reduced-complexity climate model

Abstract. Hector is an open-source reduced-complexity climate–carbon cycle model that models critical Earth system processes on a global and annual basis. Here, we present an updated version of the model, Hector V3.2.0 (hereafter Hector V3), and document its new features, implementation of new science, and performance. Significant new features include permafrost thaw, a reworked energy balance submodel, and updated parameterizations throughout. Hector V3 results are in good general agreement with historical observations of atmospheric CO2 concentrations and global mean surface temperature, and the future temperature projections from Hector V3 are consistent with more complex Earth system model output data from the sixth phase of the Coupled Model Intercomparison Project. We show that Hector V3 is a flexible, performant, robust, and fully open-source simulator of global climate changes. We also note its limitations and discuss future areas for improvement and research with respect to the model's scientific, stakeholder, and educational priorities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Accelerating laser ray tracing in high fidelity physics simulations of laser melting using squeeze U-net

Laser melting is a core component of the ongoing industrial revolution, dubbed Industry 4.0, as lasers facilitate fast and precise melting and fusion in advanced manufacturing. There is a strong need to optimize the laser process using simulations. However, this has proven challenging as high fidelity simulations are needed for predictive modeling and this is currently prohibitively expensive even when run on hundreds of processors on high performance computers. The challenge is capturing complex physics of laser material interaction, fluid dynamics, thermal physics and material phase transformations at various length and time scales. To close this technological gap, we modified a squeeze U-net to accelerate the laser ray tracing component of such high fidelity models by ~4x–40x while preserving the core physics principle of conservation of energy with 97% accuracy. This approach enables the accurate modeling of global laser energy absorption as a function of local surface temperatures and complex surface topologies, which govern the reflection directions and energy losses of laser rays upon interacting with the material surface.

Computer science↗