SEARCH · Engineering Papers

Results for “balance function”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,081 records · Page 60

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Determining the hydration energetics on carbon-supported Ru catalysts: An adsorption calorimetry and density functional theory study

Fundamental knowledge on the energetics at the interface between a water layer and a metal catalyst is essential so as to understand the roles that water can play in the synthesis, activation and regeneration of noble metal-based catalysts. Here, we report the direct measurement of the enthalpy of water adsorption (Δh ads ) on activated carbon (C) and activated C-supported Ru nanoparticles, which are promising catalyst as applied to the hydrogenation/hydrodeoxygenation (HDO) of oxygenates (phenolics, aldehydes, etc.). Specifically, the near-zero coverage enthalpy of water adsorption on a C-supported Ru catalyst is -75.3 ± 0.4 kJ/(mol water), suggesting favorable water–metal binding. This is much more exothermic than that on C, which has an enthalpy of adsorption of -50.3 ± 1.3 kJ/(mol water). Despite the favorable initial binding, the magnitudes of enthalpies of water condensation on C and Ru-C indicate that overall, their surfaces are both hydrophobic. Moreover, the experimentally-measured near-zero coverage water adsorption enthalpy at the Ru sites is in very good agreement with our density functional theory based calculations. At low coverages, we obtain a water binding energy of -61.7 kJ/(mol water), which increases to -78.1 kJ/(mol water) at saturation. Complementary results are also obtained from a thermal analysis, which employed a thermogravimetric analysis–differential scanning calorimetry–mass spectrometry (TG-DSC-MS), a spectroscopic investigation using ex situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) and a morphological evaluation with transmission electron microscopy (TEM). We point out that in carbon-supported metal catalysts, such as Ru-C, a strong hydration at near-zero coverage and relative weak water-surface interactions occurs upon saturation. Such heterogeneity is essential and crucial for catalytic hydrogenation/HDO reactions that involve balanced interactions among the water-rich reactant mixture and nonpolar organic products.

09 BIOMASS FUELS↗

DOE OSTI2022

Role Evaluation of Active Groups in Lignin on UV-Shielding Performance

The contribution of the active groups in lignin to ultraviolet (UV)-shielding performance still remains unclear. In this study, we employed 16 lignin model compounds and theoretical calculations to ascertain the contribution of different functional groups to UV-shielding performance. These results showed that the order of UV-shielding performance on functional groups was -C=C–CHO > -C=C–COOH > -C=C–CH 2 OH > --C=C > quinone > -COOH > -OCH 3 ≈ -OH. Moreover, the synergism between the para-substituted electron-donating group -OH and the electron-withdrawing group allyl on the aromatic ring played a major role in UV-shielding performance of lignin models, whereas the additive effect of meta/ortho-substituted -OCH 3 played a secondary role. Furthermore, we found that the hydrophile–lipophile balance (HLB) value of lignin models also greatly affected the sunscreen protection factor (SPF) of lignin model-based sunscreens. A closer HLB value of lignin to 8.89 corresponded with a higher SPF of the lignin model-based sunscreen. The commercial photoprotective ingredient ethylhexyl methoxycinnamate (EHMC) and lignin models had a synergistic effect on UV-shielding performance caused by the π–π stacking interaction between the parallel-displaced stacked aromatic rings of lignin models and EHMC. Additionally, the UV irradiation test showed that these sunscreens prepared using lignin models and EHMC exhibited good photostability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Continuous Counter‐Current Microfluidic Liquid–Liquid Extraction Achieved Using a Pair of Wettable Screen Meshes

Continuous counter‐current microfluidic liquid–liquid extraction performs separations by flowing immiscible liquids in opposing directions within a single flow channel. In principle, this flow arrangement enables a large number of theoretical separation units in a small footprint, without using interstage valving, pumping, and phase separation. Despite its potential for excellent separation performance, this microfluidic scheme rarely appears in literature due to the requirement for capillary forces to be greater than hydrodynamic forces for stable flow. We present a novel microfluidic device and flow approaches that overcome this force‐balance challenge, enabling stable, long‐duration continuous counter‐current flow. Additionally, we cover a suite of methodologies for quantifying the performance of the microfluidic device, revealing the number of theoretical equilibrium stages achieved. The enabling technologies include a woven mesh screen‐based microfluidic device architecture that is easily fabricated outside of a clean room, surface functionalization strategies to promote conjugate (organic/aqueous) wettability, flow approaches to eliminate bubbles and carryover, and computer‐aided flow automation with optical measurement of extraction performance. The reported experiments lasted for over 36 h, terminated only at experiment conclusion, where the device still exhibited good performance. Automated Raman spectroscopy was used for solute quantitation of the ternary system tert‐butanol in a toluene/water matrix, a ternary system that was specifically chosen to analyze the device's performance with a small solute partition ratio and to enable in‐line Raman measurements of solute concentrations in both phases. The microfluidic device possessed a 55 mm contact length and a 38.5 µL internal volume. During counter‐current flow, we observed approximately 37 equilibrium stages (37 ± 13) based on a best‐fit of the solute fraction remaining in the aqueous phase using a Kremser Group Method analysis.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Structure−Property Relationships in a Series of Hybrid Halopyridinium Ruthenium(IV) Halides

Reported is the synthesis and characterization of a family of hybrid ruthenium halides consisting of haloruthenium octahedra (X = Cl, Br) charge balanced with halopyridinium (XPy; X= H, Cl, Br, I) organic cations which assemble via noncovalent interactions between ion pairs. Diffuse reflectance spectroscopy showed that compounds containing RuBr 6 2- octahedra displayed a lower bandgap (1.05 eV< x < 1.08 eV) compared to compounds with RuCl 6 2- octahedra (1.22 eV < x < 1.43 eV). Additionally, computational density functional theory (DFT) based natural bonding orbital (NBO) analysis and density of state (DOS) methods were used to characterize second-sphere non-covalent interaction strengths and elucidate molecular orbital perturbations to elucidate their influence on Ru-X orbital constructs and in turn rationalize band gap trends. Analyses showed ligand to metal charge transfer is responsible for the small bandgap energies and are unperturbed by second-sphere interactions. Here, this report serves as a platform for probing the relationship between the structural and photophysical properties within the haloruthenium family of low dimensional transition metal halide perovskites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Accuracy, transferability, and computational efficiency of interatomic potentials for simulations of carbon under extreme conditions

Large-scale atomistic molecular dynamics (MD) simulations provide an exceptional opportunity to advance the fundamental understanding of carbon under extreme conditions of high pressures and temperatures. However, the fidelity of these simulations depends heavily on the accuracy of classical interatomic potentials governing the dynamics of many-atom systems. Here, this study critically assesses several popular empirical potentials for carbon, as well as machine learning interatomic potentials (MLIPs), in their ability to simulate a range of physical properties at high pressures and temperatures, including the diamond equation of state, its melting line, shock Hugoniot, uniaxial compressions, and the structure of liquid carbon. Empirical potentials fail to accurately predict the behavior of carbon under high pressure–temperature conditions. In contrast, MLIPs demonstrate quantum accuracy, with Spectral Neighbor Analysis Potential (SNAP) and atomic cluster expansion (ACE) being the most accurate in reproducing the density functional theory results. ACE displays remarkable transferability despite not being specifically trained for extreme conditions. Furthermore, ACE and SNAP exhibit superior computational performance on graphics processing unit-based systems in billion atom MD simulations, with SNAP emerging as the fastest. In addition to offering practical guidance in selecting an interatomic potential with a fine balance of accuracy, transferability, and computational efficiency, this work also highlights transformative opportunities for groundbreaking scientific discoveries facilitated by quantum-accurate MD simulations with MLIPs on emerging exascale supercomputers.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Computational exploration of biomedical HfNbTaTiZr and Hf 0.5 Nb 0.5 Ta 0.5 Ti 1.5 Zr refractory high-entropy alloys

Refractory high entropy alloys (RHEAs) have been proven to be a potential candidate in the biomedical field due to their balanced mechanical properties and biocompatible composition. Recent experimental findings show that RHEAs like HfNbTaTiZr and Hf 0.5 Nb 0.5 Ta 0.5 Ti 1.5 Zr have good mechanical properties such as high polarization and wear resistance than others which establish them as potential materials for biomedical application. In this work, we performed first-principles density functional theory calculations on the mechanical and thermal properties of HfNbTaTiZr and Hf 0.5 Nb 0.5 Ta 0.5 Ti 1.5 Zr. The predicted lattice constant, density, Young's modulus, and Vickers hardness are consistent with the available experimental report, which verifies the accuracy of the applied model. The thermal coefficient of linear expansion of both RHEAs has been investigated by utilizing the Debye theory. The present methods could be applied to study other future RHEAs on exploration of their physical properties.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Center for Gas Separations (CGS)

The total energy consumption in the U.S. has been rising steadily for decades, and it currently amounts to ~98,000 TBtu/yr, with approximately 30% of this total attributable to the industrial sector. Reasonable estimates indicate that 45–55% of total industry energy consumption derives from chemical separations, and for example, over 120 TBtu/yr alone is used in carrying out olefin/paraffin separations via energy-intensive cryogenic distillation. Therefore, the pursuit of new, even radically different approaches to some of the most energy-intensive industrial separations processes is an imperative scientific pursuit for reducing energy consumption toward a more sustainable future. Adsorbent and membrane-based separations can require a fraction of the energy needed for distillation methods, and as such are considered promising solutions for balancing increasing energy demand in the U.S. with the need for a massive reduction in energy consumption. Although considerable research effort has been devoted to the design of materials capable of carrying out various gas separations, usually operating through size-selective, chemisorptive, or physisorptive mechanisms, it remains a great challenge to design materials that function adequately for real-world applications. Indeed, the chemical and physical differences between molecules in gas mixtures of interest are often small, and therefore it is necessary, through the use of nanoscience and synthetic chemistry, to engineer unprecedented molecular-level control in adsorbate–adsorbent interactions. The overarching mission of the Center for Gas Separations (CGS) was to discover fundamental innovations that have the potential to dramatically reduce the energy associated with critical gas separations. In particular, the CGS developed novel synthetic routes, guided by molecular chemistry principles, as well as advanced characterization and computational methods, that have enabled the discovery of new materials and membranes tailor-made to exhibit exceptional performance for a range of gas separations processes, as required in the clean use of fossil fuels and in reducing CO 2 emissions from industry. A challenge of this magnitude required the collaboration and synergy of a large team of researchers with expertise in materials synthesis, characterization, and computations. During the 11-year project period, the CGS created a range of new materials within the family of highly-tunable, porous solids known as metal–organic frameworks (MOFs). These new frameworks demonstrate novel mechanisms for key industrial gas separations, including revolutionary new cooperative adsorption processes that enable low-energy CO 2 and CO capture, and are capable of efficiently separating olefins from paraffins, O 2 from air, and the shape-selective separation of alkane isomers. In addition, the CGS developed new strategies for incorporating these materials into composite membranes toward highly efficient and selective membrane-based separations. As a testament to the success of the CGS, two start-up companies, Mosaic Materials,4 Inc. and Flux Technology, Inc., grew out of these research efforts, and these companies are seeking to commercialize MOF and composite membranes materials for key separations in industry, including large-scale CO 2 capture and hydrocarbon separations, respectively. Another company, framergy, Inc., licensed IP resulting from CGS research toward the commercialization of adsorbents for various energy-relevant applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Transition Metal Taggants in UO 2 from First Principles

The incorporation of transition metals into nuclear fuel has gained attention both to improve fuel properties and as a possible nuclear forensics tool. Recent experimental studies by Ulrich et al. and Adorno Lopes et al. have investigated Ni and Fe as candidate transition metal dopants for potential nuclear forensics purposes and found that there is minimal alteration to key UO 2 fuel properties. In the present work, we performed density functional theory (DFT) investigations into possible defect structures for Ni and Fe incorporation. The dynamic stability of these defect structures was validated by calculating the phonon density of states. We found that transition metal incorporation likely occurs via substitution at U sites in the fluorite crystal structure with a nearby O vacancy for charge-balancing, which agrees with experimentally proposed structures. Additionally, this defect structure does not cause long-range alterations to the crystal parameters—an important consideration for use as a nuclear fuel taggant. Future needed work involves computational investigation of additional defect concentrations using larger supercells, additional transition metal charge states, and thermal effects. Such investigations can be carried forward into sintering models for a more complete understanding of the suitability of Ni and Fe as fuel taggants.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1983

Development and evaluation of a Fault-Tolerant Multiprocessor (FTMP) computer. Volume 2: FTMP software

The software developed for the Fault-Tolerant Multiprocessor (FTMP) is described. The FTMP executive is a timer-interrupt driven dispatcher that schedules iterative tasks which run at 3.125, 12.5, and 25 Hz. Major tasks which run under the executive include system configuration control, flight control, and display. The flight control task includes autopilot and autoland functions for a jet transport aircraft. System Displays include status displays of all hardware elements (processors, memories, I/O ports, buses), failure log displays showing transient and hard faults, and an autopilot display. All software is in a higher order language (AED, an ALGOL derivative). The executive is a fully distributed general purpose executive which automatically balances the load among available processor triads. Provisions for graceful performance degradation under processing overload are an integral part of the scheduling algorithms.

Lala, J. H.↗

NASA NTRS1985

The initial mass function and global rates of mass, momentum, and energy input to the interstellar medium via stellar winds

Published observational data are compiled and analyzed, using theoretical stellar-evolution models to determine the global rates of mass, momentum, and energy injected into the interstellar medium (ISM) by stellar winds. Expressions derived include psi = 0.00054 x (M to the -1.03) stars formed/sq kpc yr log M (where M is the initial mass function in solar mass units) and mass-loss = (2 x 10 to the -13th) x (L to the 1.25) solar mass/yr (with L in solar luminosity units). It is found that the wind/supernova injection of energy into the ISM and the mass loss from stars of 5 solar mass or more are approximately balanced by the dissipation of energy by cloud-cloud collisions and the formation of stars, respectively.

Van Buren, D.↗

NASA NTRS1987

Satellite sensing of aerosol absorption

A method is developed for remote sensing of aerosol absorption from satellite images of the earth's surface. The method is based on the measurement of the change in the upward radiances between a clear and a hazy day over a varying surface reflectance. For a zero change balance between brightening due to scattering and darkening due to absorption and scattering is reached. This balance is utilized via a radiative transfer model to derive the aerosol single-scattering albedo. A sensitivity study is performed, and the method is tested against laboratory measurements. It is suggested that for the case of haze introduced on top of an existing background aerosol and with a fair estimate of the scattering phase function, the error in the remotely sensed single-scattering albedo is in the range of 0.03-0.05. The main errors in the method arise from variations in the surface reflectance between the clear and the hazy days, uncertainty in the scattering phase function, and variation of the aerosol and gaseous absorption between these two days. If the satellite calibration varies with time, the measurements of single-scattering albedo can be substantially affected.

Kaufman, Yoram J.↗

NASA NTRS2005

Measurable characteristics of lysozyme crystal growth

The behavior of protein crystal growth is estimated from measurements performed at both the microscopic and molecular levels. In the absence of solutal flow, it was determined that a model that balances the macromolecular flux toward the crystal surface with the flux of the crystal surface well characterizes crystal growth observed using microscopic methods. Namely, it was determined that the model provides accurate estimates for the crystal-growth velocities upon evaluation of crystal-growth measurements obtained in time. Growth velocities thus determined as a function of solution supersaturation were further interpreted using established deterministic models. From analyses of crystal-growth velocities, it was found that the mode of crystal growth varies with respect to increasing solution supersaturation, possibly owing to kinetic roughening. To verify further the hypothesis of kinetic roughening, crystal growth at the molecular level was examined using atomic force microscopy (AFM). From the AFM measurements, it was found that the magnitude of surface-height fluctuations, h(x), increases with increasing solution supersaturation. In contrast, the estimated characteristic length, xi, decreases rapidly upon increasing solution supersaturation. It was conjectured that the magnitude of both h(x) and xi could possibly determine the mode of crystal growth. Although the data precede any exact theory, the non-critical divergence of h(x) and xi with respect to increasing solution supersaturation was nevertheless preliminarily established. Moreover, approximate models to account for behavior of both h(x) and xi are also presented.

Muramidase/chemistry↗

NASA NTRS2009

Laura Users Manual: 5.1-41601

This users manual provides in-depth information concerning installation and execution of LAURA, version 5. LAURA is a structured, multi-block, computational aerothermodynamic simulation code. Version 5 represents a major refactoring of the original Fortran 77 LAURA code toward a modular structure afforded by Fortran 95. The refactoring improved usability and maintainability by eliminating the requirement for problem-dependent re-compilations, providing more intuitive distribution of functionality, and simplifying interfaces required for multiphysics coupling. As a result, LAURA now shares gas-physics modules, MPI modules, and other low-level modules with the FUN3D unstructured-grid code. In addition to internal refactoring, several new features and capabilities have been added, e.g., a GNU-standard installation process, parallel load balancing, automatic trajectory point sequencing, free-energy minimization, and coupled ablation and flowfield radiation.

Mazaheri, Alireza↗

NASA NTRS2009

LAURA Users Manual: 5.2-43231

Mazaheri, Alireza↗

NASA NTRS2014

An Investigation of Low Earth Orbit Internal Charging

Low Earth orbit is usually considered a relatively benign environment for internal charging threats due to the low flux of penetrating electrons with energies of a few MeV that are encountered over an orbit. There are configurations, however, where insulators and ungrounded conductors used on the outside of a spacecraft hull may charge when exposed to much lower energy electrons of some 100's keV in a process that is better characterized as internal charging than surface charging. For example, the minimal radiation shielding afforded by thin thermal control materials such as metalized polymer sheets (e.g., aluminized Kapton or Mylar) and multilayer insulation may allow electrons of 100's of keV to charge underlying materials. Yet these same thermal control materials protect the underlying insulators and ungrounded conductors from surface charging currents due to electrons and ions at energies less than a few keV as well as suppress the photoemission, secondary electron, and backscattered electron processes associated with surface charging. We investigate the conditions required for this low Earth orbit "internal charging" to occur and evaluate the environments for which the process may be a threat to spacecraft. First, we describe a simple one-dimensional internal charging model that is used to compute the charge accumulation on materials under thin shielding. Only the electron flux that penetrates exposed surface shielding material is considered and we treat the charge balance in underlying insulation as a parallel plate capacitor accumulating charge from the penetrating electron flux and losing charge due to conduction to a ground plane. Charge dissipation due to conduction can be neglected to consider the effects of charging an ungrounded conductor. In both cases, the potential and electric field is computed as a function of time. An additional charge loss process is introduced due to an electrostatic discharge current when the electric field reaches a prescribed breakdown strength. For simplicity, the amount of charge lost in the discharge is treated as a random percentage of the total charge between a set maximum and minimum amount so a user can consider partial discharges of insulating materials (small loss of charge) or arcing from a conductor (large loss of charge). We apply the model to electron flux measurements from the NOAA-19 spacecraft to demonstrate that charging can reach levels where electrostatic discharges occur and estimate the magnitude of the discharge.

NeergaardParker, Linda↗

NASA NTRS2016

The Effects of Long Duration Bed Rest on Brain Functional Connectivity and Sensorimotor Functioning

Long duration spaceflight has been associated with detrimental alterations in human sensorimotor functioning. Prolonged exposure to a head-down tilt (HDT) position during long duration bed rest can resemble several effects of the microgravity environment such as reduced sensory inputs, body unloading and increased cephalic fluid distribution. The question of whether microgravity affects other central nervous system functions such as brain functional connectivity and its relationship with behavior is largely unknown, but of importance to the health and performance of astronauts both during and post-flight. In the present study, we investigate the effects of prolonged exposure to HDT bed rest on resting state brain functional connectivity and its association with behavioral changes in 17 male participants. To validate that our findings were not due to confounding factors such as time or task practice, we also acquired resting state functional magnetic resonance imaging (rs-fMRI) and behavioral measurements from 14 normative control participants at four time points. Bed rest participants remained in bed with their heads tilted down six degrees below their feet for 70 consecutive days. Rs-fMRI and behavioral data were obtained at seven time points averaging around: 12 and 8 days prior to bed rest; 7, 50, and 70 days during bed rest; and 8 and 12 days after bed rest. 70 days of HDT bed rest resulted in significant increases in functional connectivity during bed rest followed by a reversal of changes in the post bed rest recovery period between motor cortical and somatosensory areas of the brain. In contrast, decreases in connectivity were observed between temporoparietal regions. Furthermore, post-hoc correlation analyses revealed a significant relationship between motor-somatosensory network connectivity and standing balance performance changes; participants that exhibited the greatest increases in connectivity strength showed the least deterioration in postural equilibrium with HDT bed rest. This suggests that neuroplastic processes may facilitate adaptation to the HDT bed rest environment. The findings from this study provide novel insights into the neurobiology and future risk assessments of long-duration spaceflight.

Cassady, K.↗

NASA NTRS2017

Monitoring Land Surface Albedo and Vegetation Dynamics Using High Spatial and Temporal Resolution Synthetic Time Series from Landsat and the MODIS BRDF/NBAR/Albedo Product

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warmingcooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500-meter Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF (Bidirectional Reflectance Distribution Function) / NBAR (Nadir BRDF-Adjusted Reflectance) / albedo products and 30-meter Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDFAlbedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30-meter Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30-meter albedos for the intervening daily time steps in this study. These enhanced daily 30-meter spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of 0.006. These synthetic time series provide much greater spatial detail than the 500 meter gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 kilometers by 14 kilometers) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF-Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30-meter resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces.

Vegetation Index↗

DOE OSTI2025

Sharp front tracking with geometric interface reconstruction

Here, this paper presents a novel sharp front-tracking method designed to address limitations in classical front-tracking approaches, specifically their reliance on smooth interpolation kernels and extended stencils for coupling the front and fluid mesh. In contrast, the proposed method employs exclusively sharp, localized interpolation and spreading kernels, restricting the coupling to the interfacial fluid cells–those containing the interface/front. This localized coupling is achieved by integrating a divergence-preserving velocity interpolation method with a piecewise parabolic interface calculation (PPIC) and a polyhedron intersection algorithm to compute the indicator function and local interface curvature. Surface tension is computed using the Continuum Surface Force (CSF) method, maintaining consistency with the sharp representation. Additionally, we propose an efficient local roughness smoothing implementation to account for surface mesh undulations, which is easily applicable to any triangulated surface mesh. Building on our previous work, the primary innovation of this study lies in the localization of the coupling for both the indicator function and surface tension calculations. By reducing the interface thickness on the fluid mesh to a single cell, as opposed to the 4–5 cell spans typical in classical methods, the proposed sharp front-tracking method achieves a highly localized and accurate representation of the interface. This sharper representation mitigates parasitic currents and improves force balancing, making it particularly suitable for scenarios where the interface plays a critical role, such as microfluidics, fluid-fluid interactions, and fluid-structure interactions. The proposed method is comprehensively validated and tested on canonical interfacial flow problems, including stationary and translating Laplace equilibria, oscillating droplets, and rising bubbles. The presented results demonstrate that the sharp front-tracking method significantly outperforms the classical approach in terms of accuracy, stability, and computational efficiency. Notably, parasitic currents are reduced by approximately two orders of magnitude and stable results are obtained for parameter ranges where classical front tracking fails to converge.

42 ENGINEERING↗