Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “monotonicity”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Temperature-Dependent Configurational Entropy Calculations for Refractory High-Entropy Alloys

The cluster expansion formalism for alloys is used to construct surrogate models for three refractory high-entropy alloys (NbTiVZr, HfNbTaTiZr, and AlHfNbTaTiZr). These cluster expansion models are then used along with Monte Carlo methods and thermodynamic integration to calculate the configurational entropy of these refractory high-entropy alloys as a function of temperature. Many solid solution alloy design guidelines are based on the ideal entropy of mixing, which increases monotonically with N , the number of elements in the alloy. However, our results show that at low temperatures, the configurational entropy of these materials is largely independent of N , and the assumption described above only holds in the high-temperature limit. This suggests that alloy design guidelines based on the ideal entropy of mixing require further examination.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Data-Driven Constitutive Model for the Inelastic Response of Metals: Application to 316H Steel

Here, predictions of the mechanical response of structural elements are conditioned by the accuracy of constitutive models used at the engineering length-scale. In this regard, a prospect of mechanistic crystal-plasticity-based constitutive models is that they could be used for extrapolation beyond regimes in which they are calibrated. However, their use for assessing the performance of a component is computationally onerous. To address this limitation, a new approach is proposed whereby a surrogate constitutive model (SM) of the inelastic response of 316H steel is derived from a mechanistic crystal plasticity-based polycrystal model tracking the evolution of dislocation densities on all slip systems. The latter is used to generate a database of the expected plastic response and dislocation content evolution associated with several instances of creep loading. From the database, a SM is developed. It relies on the use of orthogonal polynomial regression to describe the evolution of the dislocation content. The SM is then validated against predictions of the dead load creep response given by the polycrystal model across a range of temperatures and stresses. When the SM is used to predict the response of 316H during complex non monotonic loading, extrapolating to new loading conditions, it is found that predictions compare particularly well against those from the physics-based polycrystal model.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Assessment of CMIP6 Performance and Projected Temperature and Precipitation Changes Over South America

We evaluate the performance of a large ensemble of Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6) over South America for a recent past reference period and examine their projections of twenty-first century precipitation and temperature changes. The future changes are computed for two time slices (2040–2059 and 2080–2099) relative to the reference period (1995–2014) under four Shared Socioeconomic Pathways (SSPs, SSP1–2.6, SSP2–4.5, SSP3–7.0 and SSP5–8.5). The CMIP6 GCMs successfully capture the main climate characteristics across South America. However, they exhibit varying skill in the spatiotemporal distribution of precipitation and temperature at the sub-regional scale, particularly over high latitudes and altitudes. Future precipitation exhibits a decrease over the east of the northern Andes in tropical South America and the southern Andes in Chile and Amazonia, and an increase over southeastern South America and the northern Andes—a result generally consistent with earlier CMIP (3 and 5) projections. However, most of these changes remain within the range of variability of the reference period. In contrast, temperature increases are robust in terms of magnitude even under the SSP1–2.6. Future changes mostly progress monotonically from the weakest to the strongest forcing scenario, and from the mid-century to late-century projection period. There is an increase in the seasonality of the intra-annual precipitation distribution, as the wetter part of the year contributes relatively more to the annual total. Furthermore, an increasingly heavy-tailed precipitation distribution and a rightward shifted temperature distribution provide strong indications of a more intense hydrological cycle as greenhouse gas emissions increase. The relative distance of an individual GCM from the ensemble mean does not substantially vary across different scenarios. We found no clear systematic linkage between model spread about the mean in the reference period and the magnitude of simulated sub-regional climate change in the future period. Overall, these results could be useful for regional climate change impact assessments across South America.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Elevated temperature ductility dip in an additively manufactured Al-Cu-Ce alloy

The deformation and failure mechanisms of Al-9Cu-6Ce (wt%) based alloys fabricated with laser powder bed fusion were investigated from room temperature to 400°C. The yield and ultimate tensile strengths decreased monotonically with increase in temperature, but the tensile elongation dipped unexpectedly at elevated temperatures and exhibited a minimum at 300°C. The dip in tensile elongation occurred with a concomitant dip in strain-rate sensitivity (SRS) of deformation. The as-fabricated alloy microstructure was heterogeneous, and the heat affected zone (HAZ) underneath the melt pool boundary was prone to strain localization. At 300 °C, the reduced SRS promoted the progression of strain localization in the HAZ leading to failure initiation and the dip in tensile elongation. Furthermore, a higher SRS or strain-hardening rate at other temperatures improved the tensile elongation by slowing the progression of strain localization in the HAZ such that failure initiated by other mechanisms elsewhere in the microstructure. Notably, the tensile elongation was limited by the defect structure only in a narrow temperature range (150 - 200 °C) while at other temperatures it was limited by the inherent microstructural features. This investigation exemplifies unexpected deformation and failure mechanisms possible in heterogeneous microstructures that result from additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Development of grain-scale slip activity and lattice rotation fields in Inconel 718

Using a combination of in-situ high-resolution digital image correlation (HR-DIC), Heaviside-DIC method (H-DIC), and crystal plasticity finite element (CPFE), we investigate the evolution of intragranular lattice rotations and slip activity during monotonic and cyclic loading in a high performance, polycrystalline face centered cubic material. The CPFE employs a quasi-3D model microstructure, which is a highly resolved mirror representation of the experimental in-situ test sample. In agreement, the measurements and calculations reveal that most grains, regardless of their size and lattice orientation, develop intragranular lattice rotation gradients that span the grain. For a small cluster of grains on the deformed material, we perform HR-DIC analysis of slip lines to demonstrate agreement in the active slip systems and changes in this local slip activity across the individual grains. The combined analysis reveals that deforming grains are divided into sub-granular regions of uniform lattice rotation and these regions are most often associated with only one or two active slip systems. The gradient lines that divide them correspond to changes in the predominant slip system. The model is used to examine the evolution of intragranular lattice rotation in a single fully reversed tension-compression cycle. The calculations indicate that intragranular gradients intensify during the reverse loading path as nearest neighboring regions appear to shed lattice rotation, increasing the lattice rotation in some regions, while shutting down rotation in neighboring regions. So these findings provide insight into the irreversible changes that develop within deforming grains at the scale of the grain, particularly the heterogeneous development of intragranular lattice rotation in early stages of deformation, which could serve as precursors to localization.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Unveiling the origins of work-hardening enhancement and mechanical instability in laser shock peened titanium

Laser shock-peening (LSP) produces complex gradients in microstructures and residual stresses, which enables fabrication of metallic engineering components with superior mechanical properties. In this work, we use non-destructive and spatially resolved high-energy X-ray diffraction (HE-XRD) to investigate the influences of LSP on microstructure, surface topography, and residual stress for high purity titanium plates. LSP is found to produce large compressive in-plane residual stresses near the peened surface which monotonically decay to zero about 2.5 mm below the surface. These properties were also tracked during in-situ tensile loading, allowing stress partitioning and work-hardening rates to be measured as a function of depth. The surface region is found to have the highest work-hardening rate and remain mechanically stable up to sample failure. In addition, a crystal rigid rotation along a transversal axis of about 27 & nbsp; was observed at large applied strain near the surface region, which is attributed to the formation of denser dispersed shear bands. Furthermore, a sudden drop in tensile stress was observed at a depth of similar to 160 mu m from the LSP surface, indicating the formation of a localized mechanical instability zone (LMIZ). The origin of this LMIZ could be attributed to the increase in shear stress induced by a load transfer from LSP surface due to the formation of denser dispersed shear bands. These observations indicate a localized transition in deformation behavior has occurred from dispersed shear banding to LMIZ and homogeneous deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Assessment of metadynamic recrystallization in single copper particle impacts by focused ion beam tomography

We study single Cu-on-Cu impacts relevant to cold spray deposition and quantitatively analyze the metadynamic recrystallization (mDRX) that takes place after the impact by virtue of lingering impact adiabatic heat. Unlike prior studies, the current full 3D tomographic analysis of the mDRX volume shows that mDRX is extremely common in such impacts, although it is often missed when examining 2D sections. We also report an unexpected trend: there is a “sour spot” for mDRX at velocities about 20–40 % above the velocity for particle adhesion. This non-monotonic trend is contrary to the expectations based on increasing adiabatic heating with velocity. With a schematic model, we show that the trend can be explained on the basis of heat transfer: cooling of the heat-affected region is limited by transport through the bonded regions at the particle-substrate interface. Thus, bonding has a prominent role in the heat dissipation process and the best bonded particles most rapidly bulk quench, avoiding mDRX. Here, the developed semi-empirical model aligns with the experimental findings and may help inform microstructural evolution during cold spray and post-spray processing.

FIB-SEM tomography↗
DOE OSTI2025

Temperature and dose effects on dislocation loops in self-ion irradiated high-purity iron

Body-centered cubic (BCC) Fe-based alloys are promising candidate materials for advanced nuclear reactors. However, a detailed understanding of irradiation induced dislocation loop microstructure development remains unresolved. It is a widespread belief that 〈001〉 loops become increasingly favorable over ½〈111〉 loops as irradiation temperature rises above ∼300 C. Unfortunately, the temperature effects on 〈001〉 loop have been primarily examined in in-situ irradiation on TEM thin foils but poorly explored on bulk Fe due to exceedingly limited experimental studies on bulk specimens, raising concerns about the potential influence of TEM thin foil artifacts on observed results. Here, in this study, we conducted experiments on ultra-high purity BCC Fe specimens irradiated with 6.7–8 MeV Fe ions over a wide temperature range on bulk samples. We investigated the effects of temperature (T irr = 250–500 °C), dose rate (10⁻⁵ to 10⁻³ dpa/s), and dose (0.35 to 3.5 dpa) on the formation and evolution of 〈001〉 and ½〈111〉 loops as well as cavity (void) formation. Post-irradiation Burgers vector analysis via g•b method on dislocation segments and loops revealed that 〈001〉 loop fraction does not show a monotonic positive correlation with irradiation temperature. Combined with previous and current theoretical as well as experimental findings, we explore the temperature effects on all existing models of 〈001〉 loop formation. We conclude that the prevailing reports regarding the dominance of 〈001〉 loops in Fe at elevated temperatures are mainly attributable to the loss of glissile ½〈111〉 clusters in TEM thin foil experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Elucidating the effect of preheating temperature on melt pool morphology variation in Inconel 718 laser powder bed fusion via simulation and experiment

In laser powder bed fusion (L-PBF) additive manufacturing, the mechanical performance, microstructure and defects of fabricated parts are closely associated with the melt pool morphology, e.g., its dimension and shape through the building process. Past studies have largely focused on how the process parameters such as laser power and scan speed affect melt pool characteristics. In this study, the melt pool morphology variation as a function of preheating temperature in the conduction, transition, and keyhole regimes and the underlying mechanisms in each regime are investigated through ex-situ sample characterization and computation thermal fluid dynamics (CtFD) simulation. Single tracks with different combinations of laser power and scan speed are deposited on an Inconel 718 bare plate preheated to a temperature range of 100–500 °C in the experiment. Significant changes are observed in melt pool morphology as a function of preheating temperature from optical measurements of melt track cross sections. The depth of melt pool in the three regimes increases monotonically with preheating temperature, e.g., at 500 °C, the experimental melt pool depth is increased by 49% in conduction regime, 34% in transition regime and 33% in keyhole regime, respectively, while the variation of melt pool width in each regime does not all follow an increasing trend but depends on the melt pool regimes. Melt pool width variation in the conduction and transition regimes is found to depend on the enhanced heat conduction directly related to temperature dependent thermal properties. Through validated CtFD simulations, it is found that in the keyhole regime the evaporation mass, recoil pressure, and laser drilling effect is enhanced with higher preheating temperature, which gives rise to a deeper melt pool. The simulations also reveal that preheating temperature significantly elongates the melt track length due to the increased flow rate and strong recoil pressure that accelerates the backward flow.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Estimating source-sink distributions and fluxes of reactive nitrogen and sulfur within a mixed forest canopy

The vertical source-sink distribution of air pollutants within and above forested canopies is necessary for describing the biological, physical, and chemical processes influencing the soil-vegetation-atmosphere exchange. Here, this study implemented inverse modeling methods to estimate the source-sink and flux profiles of reactive nitrogen (N) and sulfur (S) compounds from measurements of the mean concentration profiles of ammonia (NH 3 ), nitric acid (HNO 3 ), sulfur dioxide (SO 2 ), and particulate ammonium (NH 4 + ), nitrate (NO 3 − ), and sulfate (SO 4 2− ) at a forest site in the southern Appalachian Mountains. Three inverse approaches utilizing different approximations to scalar transport within the canopy were developed and evaluated against sensible heat flux measurements. The Eulerian model (EUL), which incorporates vertical velocity skewness, performed well in reproducing the turbulent heat fluxes and was subsequently used to calculate the chemical source-sink and flux profiles. Above-canopy fluxes of NH 3 were downward, indicating that the forest was a net sink of NH 3 . The soil/litter layer was both a source and a sink for NH 3 but the exchange rate at the forest floor was small. Fluxes of HNO 3 , SO 2 , NO 3 − , NH 4 + , and SO 4 2- were uni-directional (deposition only) between the air and the canopy/ground and increased monotonically from the forest floor to the canopy top. Crown foliage dominated the uptake of reactive N and S during the growing season, accounting for 80–90% of the total canopy-scale flux. Fluxes and canopy-ground partitioning estimated using the resistance-based Surface Tiled Aerosol and Gas Exchange (STAGE) model were generally comparable to EUL. The comparison highlights the need for improved parameterizations of litter exchange and NH 3 compensation points in resistance models for forest ecosystems. The findings here benefit the application of critical loads in forest ecosystems and guide further development of resistance-based exchange models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Impact of U-10Mo HALEU fuel element tolerances on the Massachusetts Institute of Technology reactor safety and operational performance – Thermal hydraulics

The U.S. is coordinating efforts for the conversion of six U.S. High Performance Research Reactors (USHPRR) including one critical facility from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts of Institute of Technology Reactor (MITR), and achieve similar performance, high-assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being evaluated. The impact of the fabrication specification and tolerances was assessed following the preliminary design of the MITR LEU fuel elements using the U-10Mo monolithic alloy. This research focuses on the analysis of fabrication specification impact on thermal hydraulics (TH) characteristic of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., coolant channel gap thickness, fuel plate thickness, etc.). The analyses are performed based on an all-fresh LEU fuel conversion plan identified in a preliminary safety analysis report submitted to the Nuclear Regulatory Commission. The reactor power margin to the onset of nucleate boiling (ONB) is assessed under the limiting safety system settings (LSSS), where a scram occurs, to ensure there is sufficient margin to the reactor safety limit, which is defined by the onset of flow instability that occurs after the ONB. The best estimate plus uncertainty approach is employed to analyze this TH characteristic, which yields realistic results while maintaining adequate conservatism, utilizing a statistical uncertainty propagation method with the STAT7 code. The TH characteristic is analyzed as a function of the variability of the specification parameters resulting from the fabrication process. The main findings of this study show that the MITR core can meet the TH safety and operational requirements at the all-LEU initial core startup (cycle 1), selected transition cycles (most reactive cycle and most limiting cycle: cycle 3 and 5, respectively) and equilibrium (cycle 14) cores under all limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core power margin to ONB on those specification parameters that have the most direct impact on TH performance is non-linear but monotonically decreasing within the specification tolerances. The third order polynomial fit curves are reported in detail for selected limiting cases and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion↗
DOE OSTI2022

“Self-trapping” in solar cell hybrid inorganic-organic perovskite absorbers

In the simplest picture, a “self-trapped” polaron forms when an excess electron or hole deforms a crystal lattice, creating a potential well with bound states. Properties of self-trapped polarons in methylammonium lead iodide perovskite (MAPbI 3 ), which is widely used as solar cell absorber, are of great interest, and are a subject of ongoing investigations and debates concerning the existence of large polarons with the co-presence of metastable self-trapping. Herein, we employ a self-interaction-free density functional theory method to investigate the stability of small polarons in tetragonal MAPbI 3 phase. The electron small polaron is found to be unstable, while the hole small polaron is found to be metastable at realistic operation temperatures of solar cells. Further, the hole polaron is found to have a hole band close to the conduction band, which in conjunction with its metastability suggests that small polarons will have an appreciable effect on charge-carrier recombinations in MAPbI 3 . Further, in this study, we posit that the existence of the metastable polarons in addition to the large polarons may explain the experimentally observed non-monotonic temperature dependence of bimolecular charge-carrier recombination rate in tetragonal MAPbI 3 phase.

14 SOLAR ENERGY↗
DOE OSTI2021

Effect of substrate tilt on sputter-deposited AuTa films

Gold-tantalum alloy films are of interest for biomedical and magnetically-assisted inertial confinement fusion applications. However, growth mechanisms of such Au-Ta alloy films are not well understood. Here, we systematically study the effect of substrate tilt on properties of films deposited by direct current magnetron sputtering from a AuTa alloy target. Experimental data is correlated with distributions of landing energies and incident angles of sputtered and backscattered species calculated by Monte Carlo simulations of ballistic collisions and gas phase atomic transport. Results reveal that the deposition rate and film density monotonically decrease with increasing substrate tilt. Properties are similar for films with a tilt of 40° and below: films are amorphous, with an average compressive stress of 0.5 GPa, a density of 16 g/cm 3 , and electrical resistivity of 200 μΩ cm. The critical substrate tilt angle for the onset of porosity, a transition to a reduced stress state, and an increase in electrical resistivity lies between 40 and 60°. Based on Monte Carlo simulations, these changes in film properties are correlated with an increase in the average impact angle of the ballistic component of atomic flux. Implications of these results to deposition onto non-planar substrates are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mixed Tin-Titanium Oxides by Atomic Layer Deposition on Planar Substrates: Physical and Electronic Structure

Ternary atomic layer deposition of Sn-doped TiO x (Sn:TiO x ) and Ti-doped SnO x (Ti:SnO x ) was performed using tetrakisdimethylamidotitanium (IV) and tin(IV) with water at 150 °C on (1 0 0) Si and quartz substrates. The physical structure of the films was investigated using scanning electron, atomic force, and Kelvin probe force microscopies as well as spectroscopic ellipsometry, Raman spectroscopy, and X-ray diffraction. All as-deposited films are amorphous and possessed flat and conformal topographies with low surface roughness (<1 nm). Annealed Sn:TiO x and Ti:SnO x crystallize into the tetragonal anatase TiO 2 and cassiterite SnO 2 structures, respectively, without large variations to the film topography. The chemical composition of the as-deposited and annealed films was investigated through X-ray photoelectron spectroscopy. While dopants are homogenously mixed, the ternary ALD to form Sn:TiO x and Ti:SnO x surprisingly does not follow the rule of mixtures, despite similarities in both the growth rates of binary SnO x and TiO x and the precursors used. This finding is also supported by spectroscopic ellipsometry. Here, we postulate the deviation from the rule of mixtures is due to non-ideal surface chemistry of the binary ALD of SnO x . The optical bandgap of all films was quantified by the Tauc analysis and the valence band position was investigated by valence band X-ray photoelectron spectroscopy. Intra-bandgap states above the valence band maximum were observed for as-deposited Ti:SnO x and Sn:TiO x and ascribed to 5 s states from Sn 2+ ions. The presence of these intra-bandgap states induces enhanced visible absorption and bathochromic shifts to the apparent bandgap for all as-deposited films. Annealing the films removes the intra-bandgap states and increases the optical bandgap. Simplified band diagrams created using the data available demonstrated that the conduction and valence band positions did not vary between the positions of pure TiO 2 or SnO 2 with monotonic increases to the dopant incorporation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Scanning tunneling microscope tip-induced formation of Bi bilayers on Bi 2 Te 3

We report the formation of Bi(111) bilayer islands and crater structures on Bi 2 Te 3 (111) surfaces induced by voltage pulses from a scanning tunneling microscope tip. Pulses above a threshold voltage (+3 V) produce craters ∼0.5μm in diameter, similar to the size of the tip. Redeposited material self-assembles into a network of atomically ordered islands with a lattice constant identical to the underlying Bi2Te3 surface. The island size monotonically decreases over several micrometers from the pulse site, until the pristine Bi 2 Te 3 surface is recovered. We assign these islands to Bi bilayer based on atomic resolution images, analysis of step heights, and tunneling spectroscopy. Here, the dependence of bilayer formation on bias polarity and the evidence for defect diffusion together suggest a mechanism driven by the interplay of field evaporation and tunneling-current-induced Joule heating.

Bi bilayer↗
DOE OSTI2023

Characteristics of upward-connecting-leader current leading to attachment in downward negative cloud-to-ground lightning strokes

For this work, we analyzed currents associated with upward connecting leaders (UCLs) initiated from the Kennedy Space Center Industrial Area Tower in two negative cloud-to-ground strokes that struck the tower. One stroke was also recorded using a high-speed video camera located 760 m from the KSC IAT. The NLDN-reported peak currents for the two strokes were – 31.7 and – 98.5 kA. During the UCL development phase the current waveforms exhibited a monotonically (quasi-exponentially) increasing “background” current overlaid with 10-μs scale pulses with a median amplitude of 51.1 A. The UCL current durations for the two strokes were 1039 and 449 μs, respectively. During the pre-attachment processes (UCL and slow front) the total negative charge effectively transferred to ground were 70.2 and 55 mC, respectively. For the stroke captured on high-speed video, the average line-charge-density for the 109-m long UCL was found to be 0.5 mC/m. The average UCL 2-D speed was 2.4 × 10 5 m/s, and it was observed to accelerate toward the downward leader prior to attachment. We observed that UCL-pulse amplitudes are larger, background currents are higher, and interpulse intervals are shorter at later times during UCLs, which can be attributed to the intensification of the local electric field due to the approaching downward negative leader. The median positive charge injected into the UCL by a pulse was 297 μC. The UCL associated with the higher peak-current stroke produced the highest injected pulse-charge values about three times sooner during its development, likely due to the 2–2.7 times faster average downward leader vertical speeds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Nematic cell alignment directs calcium waves in an epithelial monolayer

Tissues rely on supracellular signals to coordinate their cells over a long range. Two such tissuescale cues are calcium waves and patterns of cell-cell alignment or nematic order. During wound healing, for example, calcium waves propagate across a tissue to guide directed cell migration and reepithelialization. Defects in long-range cell-cell alignment, or nematic orientation, can act to localize morphogenetic events in a tissue. Although these two cues have been considered in isolation, we demonstrate a relationship in epithelial tissue between long-range calcium signaling and the cell’s nematic order: The speed of a wound-induced calcium wave depends monotonically on the angle between the wave vector and cell axis, with maximal wave speed occurring perpendicular to the tissue’s orientation. Including anisotropic di↵usive coupling between cells in a canonical reactiondi↵usion model recapitulates our measured calcium wave dynamics. Our model demonstrates how orientation defects can desynchronize information propagation across a tissue. A calcium wave front is bent around nematic defects, therefore cells the same distance from a wound can receive the calcium signal at di↵erent times. Our work elucidates how spatial patterns in global cell alignment can control collective communication via calcium signaling during development, wound healing, and disease.

Winterstrain, Annemarie C. [Brandeis University, W↗
DOE OSTI2024

Optimization-based, property-preserving algorithm for passive tracer transport

Here we present a new optimization-based property-preserving algorithm for passive tracer transport. The algorithm utilizes a semi-Lagrangian approach based on incremental remapping of the mass and the total tracer. However, unlike traditional semi-Lagrangian schemes, which remap the density and the tracer mixing ratio through monotone reconstruction or flux correction, we utilize an optimization-based remapping that enforces conservation and local bounds as optimization constraints. In so doing we separate accuracy considerations from preservation of physical properties to obtain a conservative, second-order accurate transport scheme that also has a notion of optimality. Moreover, we prove that the optimization-based algorithm preserves linear relationships between tracer mixing ratios. We illustrate the properties of the new algorithm using a series of standard tracer transport test problems in a plane and on a sphere.

97 MATHEMATICS AND COMPUTING↗