Search NASA⌕ Search

SEARCH · Search NASA

Results for “comparisons”

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

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

At least 253 records · Page 14

First implementation of gyrokinetic exact linearized Landau collision operator and comparison with models

Gyrokinetic simulations are fundamental to understanding and predicting turbulent transport in magnetically confined fusion plasmas. Previous simulations have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius (FLR) effects. We have implemented the linearized Fokker–Planck collision operator with exact field-particle terms and full FLR effects in a gyrokinetic code (GENE). The new operator, referred to as “exact” in this paper, allows the accuracy of model collision operators to be assessed. The conservative Landau form is implemented because its symmetry underlies the conservation laws and the H-theorem, and enables numerical methods to preserve this conservation, independent of resolution. The implementation utilizes the finite-volume method recently employed to discretize the Sugama collision model in GENE, allowing direct comparison between the two operators. Results show that the Sugama model appears accurate for the growth rates of trapped electron modes (TEMs) driven only by density gradients, but appreciably underestimates the growth rates as the collisionality and electron temperature gradient increase. The TEM turbulent fluxes near the nonlinear threshold using the exact operator are similar to the Sugama model for the n e = d ln T e /d ln n e = 0 case, but substantially larger than the Sugama model for the n e = 1 case. The FLR effects reduce the growth rates increasingly with wavenumber deepening a “valley” at the intermediate binormal wavenumber as the unstable mode extends from the TEM regime to the electron temperature gradient instability regime. Application to the Hinton–Rosenbluth problem shows that zonal flows decay faster as the radial wavenumber increases and the exact operator yields weaker decay rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cross-code comparison of the impact of the fill tube on high yield implosions on the National Ignition Facility

Fill tubes are used to inject deuterium and tritium fuel into inertial confinement fusion capsules fielded on the National Ignition Facility. These fill tubes have been shown to have a detrimental effect on capsule performance, primarily by introducing a low-density pathway into the central fuel region that enables the jetting of ablation material into the hot spot. Due to the complexity of the highly nonlinear flow associated with the fill tube and the challenge of diagnosing the evolution of the fill tube jet late in the implosion experiments, the uncertainty in how this perturbation source evolves is great. In this work, we report on the results of a detailed code comparison performed to understand uncertainties in computational modeling of the impact of fill tubes on implosion performance. The study employed two radiation-hydrodynamics codes, HYDRA and xRAGE, which employ very different meshing strategies and hydrodynamics solvers, as well as two radiation transport methodologies, discrete ordinates and multi-group diffusion. Our results demonstrate generally good agreement between codes through most of the implosion although they indicate sensitivity to opacity averaging methods. Late in the implosion, differences arise in the distribution and amount of contaminant although these differences have a remarkably small impact on the amount of yield reduction due to the fill tube. While these results demonstrate sensitivity in fill tube modeling to algorithmic choices, the observed differences between codes are small relative to known sensitivities due to expected variations in the fill tube geometry. Finally, we have developed a methodology for performing multi-group diffusion simulations that show good agreement with the more accurate discrete ordinates method

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extended magnetohydrodynamics simulations of thin-foil Z-pinch implosions with comparison to experiments

Cylindrical foil liners, with foil thicknesses on the order of 400 nm, are often used in university-scale Z-pinch experiments (~1 MA in 100 ns) to study physics relevant to inertial confinement fusion efforts on larger-scale facilities (e.g., the MagLIF effort on the 25-MA Z facility at Sandia National Laboratories). The use of ultrathin foil liners typically requires a central support rod to maintain the structural integrity of the liner target assembly prior to implosion. The radius of this support rod sets a limit on the maximum convergence ratio achievable for the implosion. In recent experiments with a support rod and a preimposed axial magnetic field, helical instability structures in the imploding foil plasma were found to persist as the foil plasma stagnated on the rod and subsequently expanded away from the rod. We have now used the 3D extended magnetohydrodynamics simulation code PERSEUS (which includes Hall physics) to study these experiments. The results suggest that it is the support rod which is responsible for the helical structures persisting beyond stagnation. Furthermore, we find that as the radius of the support rod decreases (i.e., as the convergence ratio increases), the integrity and persistence of the helical modes diminish. In the limit with no support rod, we find that the structure of the final stagnation column is governed by the structure of the central precursor plasma column. Furthermore, these simulation results and their comparisons to experiment are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of fixed charge and polarizable models for predicting the structural, thermodynamic, and transport properties of molten alkali chlorides

Results from extensive molecular dynamics simulations of molten LiCl, NaCl, KCl, and RbCl over a wide range of temperatures are reported. Comparison is made between the “Polarizable Ion Model” (PIM) and the non-polarizable “Rigid Ion Model” (RIM). Densities, self-diffusivities, shear viscosities, ionic conductivities, and thermal conductivities are computed and compared with experimental data. In addition, radial distribution functions are computed from ab initio molecular dynamics simulations and compared with the two sets of classical simulations as well as experimental data. Here, the two classical models perform reasonably well at capturing structural and dynamic properties of the four molten alkali chlorides, both qualitatively and often quantitatively. With the singular exception of liquid density, for which the PIM is more accurate than the RIM, there are few clear trends to suggest that one model is more accurate than the other for the four alkali halide systems studied here.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Performance of microchannel plate based detectors for <25 keV x rays: Monte Carlo simulations and comparisons with experimental results

We present the results of Monte Carlo simulations of the microchannel plate (MCP) response to x rays in the 250 eV to 25 keV energy range as a function of both x-ray energy and impact angle and their comparisons with the experimental results from the X8A beamline at the National Synchrotron Light Source at Brookhaven National Laboratory. Incoming x rays interact with the lead glass of the microchannel plate, producing photoelectrons. Transport of the photoelectrons is neglected in this model, and it is assumed that photoelectrons deposit all their energy at the point they are created. This deposition leads to the generation of many secondary electrons, some fraction of which diffuse to the MCP pore surface where they can initiate secondary electron cascades in the pore under an external voltage bias. X-ray penetration through multiple MCP pore walls is increasingly important above 5 keV, and the effect of this penetration on MCP performance is studied. In agreement with past measurements, we find that the dependence of MCP sensitivity with angle relative to the pore bias changes from a cotangent dependence to angular independence and then proceeds to a secant dependence as the x-ray energy increases. We also find that with the increasing x-ray energy, the MCP gain sensitivity as a function of bias voltage decreases. The simulations also demonstrate that for x rays incident normal to the MCP surface, spatial resolution shows little dependence on the x-ray energy but degrades with the increasing x-ray energy as the angle of incidence relative to the surface normal increases. This agrees with experimental measurements. Simulation studies have also been completed for MCPs gated with a subnanosecond voltage pulse. In this work, we find that the optical gate profile width increases as the x-ray energy is increased above 5 keV, a consequence of increased x-ray penetration at energies >5 keV. Simulations of the pulsed dynamic range show that the dynamic range varies between ~100 and 1000 depending on x-ray energy and peak voltage.

47 OTHER INSTRUMENTATION↗

NaCl aggregation in water at elevated temperatures and pressures: Comparison of classical force fields

The properties of water vary dramatically with temperature and density. This can be exploited to control its effectiveness as a solvent. Thus, supercritical water is of keen interest as solvent in many extraction processes. The low solubility of salts in lower density supercritical water has even been suggested as a means of desalination. The high temperatures and pressures required to reach supercritical conditions can present experimental challenges during collection of required physical property and phase equilibria data, especially in salt-containing systems. Molecular simulations have the potential to be a valuable tool for examining the behavior of solvated ions at these high temperatures and pressures. However, the accuracy of classical force fields under these conditions is unclear. We have, therefore, undertaken a parametric study of NaCl in water, comparing several salt and water models at 200 bar–600 bar and 450 K–750 K for a range of salt concentrations. We report a comparison of structural properties including ion aggregation, hydrogen bonding, density, and static dielectric constants. All of the force fields qualitatively reproduce the trends in the liquid phase density. An increase in ion aggregation with decreasing density holds true for all of the force fields. The propensity to aggregate is primarily determined by the salt force field rather than the water force field. This coincides with a decrease in the water static dielectric constant and reduced charge screening. While a decrease in the static dielectric constant with increasing NaCl concentration is consistent across all model combinations, the salt force fields that exhibit more ionic aggregation yield a slightly smaller dielectric decrement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct comparison of wire, foil, and hybrid X-pinches on a 200 kA, 150 ns current driver

Wire X-pinches (WXPs) have been studied comprehensively as fast (~1 ns pulse width), small (~1 μm) x-ray sources, created by twisting two or more fine wires into an “X” to produce a localized region of extreme magnetic pressure at the cross-point. Recently, two alternatives to the traditional WXP have arisen: the hybrid X-pinch (HXP), composed of two conical electrodes bridged by a thin wire or capillary, and the laser-cut foil X-pinch (LCXP), cut from a thin foil using a laser. We present a comparison of copper wire, hybrid, and laser-cut foil X-pinches on a single experimental platform: UC San Diego’s ~200 kA, 150 ns rise time GenASIS driver. All configurations produced 1–2 ns pulse width, ≤5 μm soft x-ray (Cu L-shell, ~1 keV) sources (resolutions diagnostically limited) with comparable fluxes. WXP results varied with linear mass and wire count, but consistently showed separate pinch and electron-beam-driven sources. LCXPs produced the brightest (~1 MW), smallest (≤5 μm) Cu K-shell sources, and spectroscopic data showed both H-like Cu K α lines indicative of source temperatures ≥2 keV, and cold K α (~8050 eV) characteristic of electron beam generated sources, which were not separately resolved on other diagnostics (within 1–2 ns and ≤200 μm). HXPs produced minimal K-shell emission and reliably single, bright, and small L-shell sources after modifications to shape the early current pulse through them. Benefits and drawbacks for each configuration are discussed to provide potential X-pinch users with the information required to choose the configuration best suited to their needs.

42 ENGINEERING↗

Comparison of the dielectric and magnetocaloric properties of bulk and film of GdFe 0.5 Cr 0.5 O 3

Here, we report a comparison of the magnetic, magnetocaloric, and dielectric properties of 50% iron substituted GdCrO 3 (GdFe 0.5 Cr 0.5 O 3 ) bulk pellet and 960 nm thick film of GdFe 0.5 Cr 0.5 O 3 (GFCO). The 960 nm film was synthesized on a platinized-silicon substrate by chem. soln. deposition and spin-coating methods. The X-ray diffraction scans of the bulk sample and the film as well as the morphology of the film as examined by the field-emission scanning electron microscope indicate phase-pure and polycrystalline nature of these samples. XPS was used to deermine the valence states of Gd, Fe, and Cr. The temp. dependence of the dielectric const. from 225 to 700 K shows peaks at T C = 525 K for the bulk and ≈450 K for the film due to ferroelectric to paraelectric transitions, since elec. polarization vs elec. field hysteresis loops are observed. at room temp. The dielectric studies in the bulk GFCO for T > T C indicate a relaxor-like behavior. The measurements of the magnetization (M) of the samples as a function of temp. (5-350 K) and magnetic field (H) up to 7 T (=70 kOe) depict hysteresis behavior at low temps. due to the canted antiferromagnetic order of Fe 3+ /Cr 3+ below the Néel temp. of ≈275 K. The M vs H isotherms at various temps. are used to det. and compare the magnetic entropy change (-ΔS) and relative cooling power (RCP) of the two samples, yielding (-ΔS) = 30.7 J/kg K (18.8 J/kg K) and RCP = 566.5 J/kg (375 J/kg) for the bulk (960 nm film) samples of GFCO at 7 K and 7 T, resp. The plot of RCP vs T shows that magnetic cooling for this system is most effective for T < 30 K. Comparatively smaller magnitudes of (-ΔS) and RCP for the film vis-a-vis the bulk sample of GFCO scale with its reduced magnetization. This suggests that further improvements in the quality of the films are needed to improve their magnetization and hence their magnetocaloric properties, possibly making them useful for on-chip cooling in miniaturized devices.

36 MATERIALS SCIENCE↗

Code-to-code comparison and validation of the radiation-hydrodynamics capabilities of the FLASH code using a laboratory astrophysical jet

The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and behavior of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities in the publicly available FLASH radiation-hydrodynamics code, we present detailed code-to-code comparisons between FLASH and the HYDRA code developed at Lawrence Livermore National Laboratory using previously published HYDRA simulations from Grava et al. [Phys. Rev. E 78, 016403 (2008)]. That study describes a laser experiment that produced a jet-like feature that the authors compare to astrophysical jets. Importantly, the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment included detailed x-ray interferometric measurements of electron number densities and a time-integrated measurement of the soft x-ray spectrum. Despite markedly different methods for treating the computational mesh, and different equations of state and opacity models, the FLASH results resemble the results from HYDRA and, most importantly, the experimental measurements of electron density. Having validated the FLASH code in this way, we use the code to further investigate and understand the formation of the jet seen in the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment and discuss its relation to the Wan et al. [Phys. Rev. E 55, 6293 (1997)] experiment at the NOVA laser.

Orban, Chris (ORCID:0000000312004538)↗

System-Level Comparison of Sodium and Salt Systems in Support of the Gen3 Liquids Pathway

This paper presents the comparison of the annual performance and the Levelised Cost of Energy (LCOE) of two high- temperature solar power tower configurations using a 565MWth chloride salt receiver and a 540 MWth liquid sodium receiver. Each power tower consists of a sodium or salt receiver, a 175 m high concrete tower, and a two-tank chloride salt system for thermal energy storage (TES). The power plants were simulated using SolarTherm to estimate their energy per year, capacity factor, and LCOE. The simulation was verified against the System Advisor Model (SAM), showing that the LCOE calculation differences are within ±1.0%. Results showed that the sodium-receiver system yields 69.60 USD/MWh real LCOE, versus 78.55 USD/MWh for the chloride salt system. Results also showed that sodium wins here due to a smaller and more efficient receiver, enabled by higher flux limits and reduced field size, leading to less mirror reflective area and site improvement cost.

capacity factor↗

Large-database cross-verification and validation of tokamak transport models using baselines for comparison

State-of-the-art 1D transport solvers ASTRA and TRANSP are verified, then validated across a large database of semi-randomly selected, time-dependent DIII-D discharges. Various empirical models are provided as baselines to contextualize the validation figures of merit using statistical hypothesis tests. For predicting plasma temperature profiles, no statistically significant advantage is found for the ASTRA and TRANSP simulators over a baseline empirical (two-parameter) model. For predicting stored energy, a significant advantage is found for the simulators over a baseline empirical model based on confinement time scaling. Uncertainty in the results due to diagnostic and profile fitting uncertainties is approximated and determined to be insignificant due in part to the large quantity of discharges employed in the study. Advantages are discussed for validation methodologies like this one that employ (1) large databases and (2) baselines for comparison that are specific to the intended use-case of the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of tungsten versus molybdenum for double shell capsules using machine learning design optimization

Double shell targets are an alternative ignition platform for inertial confinement fusion. One design consideration for double shell targets is the choice of inner shell material to help trap radiation emitted by the hot fuel to aid ignition. Materials such as molybdenum and tungsten are of interest for the inner shell layer of the targets. While molybdenum has a lower density that could inhibit instability growth and allow for radiography and code benchmarking, tungsten has a higher density that could provide more compression and confinement. These tradeoffs have been explored using optimized designs for each material. Our previous work [Vazirani et al., “Coupling 1D xRAGE simulations with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 28, 122709 (2021); Vazirani et al., “Coupling multi-fidelity xRAGE with machine learning for graded inner shell design optimization in double shell capsules,” Phys. Plasmas 30, 062704 (2023); and Vazirani et al., “Bayesian batch optimization for molybdenum versus tungsten inertial confinement fusion double shell target design,” Stat. Anal. Data Min. 17, e11698 (2024)] resulted in a multi-fidelity Bayesian optimization framework to find yield-optimized double shell target geometries. By leveraging simulations of varying fidelities (one-dimensional and two-dimensional) to inform one another, the multi-fidelity optimization was able to optimize a design in the highest fidelity with significantly fewer simulations than would be used in a systematic parameter scan. In this work, we apply the multi-fidelity Bayesian optimization to explore the optimized designs of double shell targets with molybdenum and tungsten inner shells as well as the physics producing the high performing implosions. A physics exploration of all the simulations used in this study shows trends in designs that contribute to high yields, ion temperatures, and fuel areal densities. Comparison of molybdenum and tungsten simulations shows that they can produce similar implosion conditions with different geometries, which would be important to study in experiments. Graded density layers produce varying performances with the two materials but continue to be of interest for future studies along with studies of doped inner shell materials and applied surface roughness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct comparison of gyrokinetic and fluid scrape-off layer simulations

Typically, fluid simulations are used for tokamak divertor design. However, fluid models are only valid if the scrape-off layer (SOL) is highly collisional. This assumption is valid in many present-day experiments but is questionable in the upstream SOL of some high-power scenarios envisioned for burning plasmas and fusion pilot plants. This paper reports on comparisons between fluid and kinetic simulations of the SOL for upstream parameters and geometry representative of the Spherical Tokamak for Energy Production fusion pilot plant. The SOLPS-ITER (fluid) and Gkeyll (gyrokinetic) codes are operated in a two-dimensional axisymmetric mode, which replaces turbulence with ad-hoc diffusivities. In kinetic simulations, we observe that the ions in the upstream SOL experience significant mirror trapping. This substantially increases the upstream temperature and has important implications for impurity dynamics. We show that the mirror force, which is excluded in SOLPS’s fluid equations, enhances the electrostatic potential drop along the field line in the SOL. We also show that the assumption of equal main ion and impurity temperatures, which is made in commonly used fluid codes, is invalid for the regimes explored here. The combination of these effects results in superior confinement of impurities to the divertor region in kinetic simulations, consistent with our earlier predictions [Kotschenreuther et al., in 29th IAEA 29 Fusion Energy Conference (IAEA, London, UK, 2023)]. This effect can be dramatic, reducing the midplane impurity density by orders of magnitude. These results indicate that in lower collisionality SOL’s the tolerable downstream impurity densities may be higher than would be predicted by fluid simulations, allowing for higher radiated power while avoiding unacceptable core contamination. Our results highlight the importance of kinetic simulations for divertor design and optimization for fusion pilot plants.

Computational fluid dynamics↗

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A comparison of the dielectric function for p-point water models

We compute the longitudinal dielectric function ɛl(k) for p-point water models, p = 3 (SPC/E, TIP3P), p = 4 (TIP4P-EW, OPC), and p = 5 (TIP5P-E) from the charge–charge fluctuation function SZZ(k), obtained from two methods: average of the charge in k-space (method:sum) and from the pair distribution function in Fourier space (method:integral). The latter requires a continuation to small k-values, which we thoroughly discuss. We conclude with a detailed comparison of the longitudinal dielectric function for the different p-point models and its main characteristics: small k and large k limits, position of the poles, and its (complex) zeros.

Shevy, Noah [Ames National Laboratory, Department ↗

Comparison of Irradiated TRISO Fuel Radioactivity from Multiple Advanced Reactor Designs

The SCALE code system was used to model, deplete, and compare several different TRISO-fueled reactor designs: a helium-cooled prismatic reactor, a helium-cooled pebble-bed reactor (PBR), and a Fluoride-Lithium-Beryllium (FLIBE) molten-salt-cooled PBR. The purpose of this comparison was to understand how differences in the reactor designs affect the radioactivity of the fuel after discharge and whether those differences are significant. First, the various reactor designs were build and depleted in the TRITON module for each design and fuel enrichment. Then, the TRITON outputs were used to create burn-up dependent reactor libraries. These libraries were then used by ORIGEN to determine the activities of discharged fuel for each reactor, which were compared to generic Westinghouse 17x17 fuel. Overall, the results showed that short term decays are dominated by reactors with higher operating powers, and the reactor type, initial fuel enrichment, and maximum burn-up are of only secondary importance. Although this analysis only focuses on activities in becquerels, these dependencies are consistent with the expected behavior of decay heat. However, analysis of long- term time periods post-irradiation shows that the reactor type and maximum burn-up have strong impacts on the activities; initial fuel enrichment has a secondary impact while operating power is inconsequential. These results would be useful for analyses, such as dose assessment and modeling in post- release scenarios; normal fuel handling operations; and spent fuel transport, storage and disposal. Of particular interest, the results in this report show that analyses that focus on spent nuclear fuel of advanced reactors need to consider each parameter carefully. Unsurprisingly, if the correct operating power is not used in short term analyses, the results will not be correct. Perhaps unexpectedly, however, if the correct reactor type is not used, then the long term results will also be incorrect, especially for areas such as permanent disposal. Even though this report focuses on the total activity of nuclear fuel, it provides initial results on the effects of various input parameters and also provides a framework to extend the work into other analyses of spent fuel from advanced reactors, especially those employing TRISO fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Design and Analysis of the Open-Surface Slow Li Flow Divertor and Comparison to the Fast Li Flow Divertor

In the ongoing U.S. project, “Liquid Metal Plasma Facing Components,” sponsored by the U.S. Department of Energy, efforts have been taken to develop two open-surface divertor designs for the Fusion Nuclear Science Facility using liquid lithium (Li) as a heat and particle flux removal media. The main focus of this study is the design and analysis of a slow (~1 mm/s) and thin (<1 mm) open-surface Li flow divertor with a Li-cooled substrate, which is then compared with an earlier design of a fast (up to 10 m/s) and thick (~0.5 cm) Li flow divertor with the substrate cooled with helium. The slow Li flow divertor design is based on the original LiWall concept developed at the Princeton Plasma Physics Laboratory. Such a thin and slow Li layer can remove the particle flux by reducing the recycling flux, while the heat flux is removed mainly through the heat sink located beneath. In the present study, the heat sink is provided through a Li cooling flow inside the substrate of reduced activation ferritic/martensitic steel. By performing a multiphysics analysis with COMSOL that included liquid-metal magnetohydrodynamics (MHD), heat transfer, and structural mechanics, the impact of various factors on the divertor heat removal capability, such as Li flow velocity, MHD effects, and inlet velocity boundary condition, were examined. Based on comparisons of the two divertor designs, it was shown that the fast-flow divertor significantly outperformed the slow-flow design, whose heat removal capability was limited to ~1 to 2 MW/m2.

Jiang, Yuchen [ORNL] (ORCID:0000000283280284)↗

A comparison of Eulerian and semi-Lagrangian approaches for modeling stream water quality

This paper describes and compares some of the advantages and limitations of Eulerian and Lagrangian approaches to water quality modeling and introduces a mixed Eulerian-Lagrangian (or semi-Lagrangian) methodology that captures the strengths of both approaches. The semi-Lagrangian modeling approach is applied to advection-dominated rivers, and flexibly ensures unconditional stability for all time step durations and grid segmentations. The semi-Lagrangian modeling approach is demonstrated by applying it to estimate the dissolved oxygen concentrations in the Sava River in Slovenia, focusing on aspects of the methodology and findings that would be of broad interest to managers of water quality in fluvial water bodies. Results of comparisons of the semi-Lagrangian model with the Eulerian-based QUAL2K model in steady and non-steady scenarios demonstrate that while both models are fully capable of producing satisfactory results when optimally configured, the semi-Lagrangian approach offers accuracy and stability without sensitivity to the interaction of time step size and computational grid segmentation scheme.

Sava River↗