Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “limitations”

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Populating the brown dwarf and stellar boundary: Five stars with transiting companions near the hydrogen-burning mass limit

We report the discovery of five transiting companions near the hydrogen-burning mass limit in close orbits around main sequence stars originally identified by the Transiting Exoplanet Survey Satellite (TESS) as TESS objects of interest (TOIs): TOI-148, TOI-587, TOI-681, TOI-746, and TOI-1213. Using TESS and ground-based photometry as well as radial velocities from the CORALIE, CHIRON, TRES, and FEROS spectrographs, we found the companions have orbital periods between 4.8 and 27.2 days, masses between 77 and 98 MJup , and radii between 0.81 and 1.66 RJup . These targets have masses near the uncertain lower limit of hydrogen core fusion (~73-96 MJup ), which separates brown dwarfs and low-mass stars. We constrained young ages for TOI-587 (0.2 ± 0.1 Gyr) and TOI-681 (0.17 ± 0.03 Gyr) and found them to have relatively larger radii compared to other transiting companions of a similar mass. Conversely we estimated older ages for TOI-148 and TOI-746 and found them to have relatively smaller companion radii. With an effective temperature of 9800 ± 200 K, TOI-587 is the hottest known main-sequence star to host a transiting brown dwarf or very low-mass star. We found evidence of spin-orbit synchronization for TOI-148 and TOI-746 as well as tidal circularization for TOI-148. These companions add to the population of brown dwarfs and very low-mass stars with well measured parameters ideal to test formation models of these rare objects, the origin of the brown dwarf desert, and the distinction between brown dwarfs and hydrogen-burning main sequence stars.

Nolan Grieves↗
NASA NTRS2020

Exploring the Limits of the Data-Model-Theory Synergy: “Hot” MW Transitions for Rovibrational IR Studies

In order to further improve the accuracy of rovibrational IR line lists generated from the “Best Theory +Reliable High-resolution Experiment” (BTRHE) strategy from 0.01-0.05 cm-1, or 300-1500 MHz, to ~10 MHz, we explore the current limits of the Data-Model-Theory synergy by examining the accuracy and consistency of existing data, then propose that “hot” bands in microwave (MW) spectra is the solution we need for future enhancements. The Ames SO2 J=0-20 rovibrational energy levels computed on the semi-empirically refined Ames-2 potential energy surface (PES) are fit to the Effective Hamiltonian (EH) model regularly used in the experimental infrared (IR) analysis for SO2 isotopologues. In the fitted EH(Ames) model, the rotational constants A/B/C and all 5 quartic centrifugal distortion constants display clear, systematic, and consistent patterns along the vibrational state energy or quanta. Such consistent patterns may facilitate the vibrational assignments for MW hot bands and extract more information from high temperature MW spectra. Some EH(Expt) analyses were carried out with the lowest order Coriolis Coupling term, C1. Their constants should not be directly compared with other EH(Expt) and EH(Ames) results. After excluding them, our  = EH(Ames)- EH(Expt) analyses for 5 isotopologues (626, 636, 646, 628 and 828) indicates some loss of accuracy and consistency starting from vibrational states as low as 22 or 1000 cm-1. Some EH parameters, e.g. K, may have relative deviations as large as 50-100% and totally lose any recognizable patterns. This simply means that current EH(Expt) models do not have the system-wide consistency we need to further refine the EH(Ames) and Ames rovibrational IR line lists. A large part of such defects are probably inherited from the limited precision of experimental line positions, i.e. 1E-3 ~ 1E-4 cm-1, or 3-30 MHz. This is confirmed in a series of truncation tests using the Ames data. Although the EH(Ames) consistency may help identify unreliable rovibrational EH(Expt) parameters, and make reliable predictions for minor isotopologues and unobserved vibrational bands, we believe only the highly accurate “hot” MW transitions can provide real enhancements for EH(Expt) accuracy and consistency. “Hot” MW spectra should play a more significant role in the future synergy of Data, Model, and Theory in the field of rovibrational IR studies.

Xinchuan Huang↗
NASA NTRS2022

Nutrient Limitation Predisposes A Cultivate of Burkholderia Contaminans From the ISS Water Processor Assembly to Biofilm Formation Under Simulated Microgravity.

The International Space Station (ISS) Water Processor Assembly (WPA) experiences intermittent dormancy between water recycling events thus promoting biofilm formation within the system. In this work we aimed to gain a deeper understanding of the impact of nutrient limitation on bacterial growth and biofilm formation under microgravity in support of biofilm mitigation efforts in exploration water recovery systems. A representative species of bacteria that is commonly cultured from the ISS WPA was cultured in an WPA influent water ersatz formulation tailored for microbiology studies. Burkholderia contaminans was cultured under a simulated microgravity (SµG) treatment in a vertically rotating high-aspect rotating vessel (HARV), with a rotating control (R) in the horizontal plane at the determined optimal rpm of 15 along with a stationary (S) control. At different time points, the bacterial culture and ersatz were harvested for bacterial counts, transcriptomic and nutrient content analyses. Under the test conditions, the culture under SµG treatment consumed the essential nutrients faster than the R and S control cultures in the early stage of growth, thus approaching a nutrient limited growth condition earlier than the controls. The rapid uptake and subsequent depletion of essential nutrients was further illustrated in the transcriptomic response of the SµG culture when compared to the transcriptomic response of the R and S control conditions. The observed starvation response may serve as one element to explain a moderate enhancement of biofilm formation in the SµG treatment. One implication of this investigation is that biofilm mitigation in the ISS environment could be supported by ensuring a steady flow of water as a vehicle for essential nutrients within the WPA to avoid complete consumption which occurs in times of no flow lending to undesired biofilm formation.

Angie Diaz↗
NASA NTRS2023

Ch. 12. A Theoretical Approach to Management of Limited Attentional Resources to Support the m:N Operation in Advanced Air Mobility Ecosystem

Advanced air mobility (AAM) technologies incorporate increasingly autonomous systems that allow fully remote, independent, and intelligent operation of air vehicles to support the transportation of goods and passengers within and across urban and rural areas. With a myriad of automated technologies enabling the AAM ecosystem, the human operator’s role will likely be a passive supervisory monitor of the air vehicles, involving increasingly fewer humans (m) that manage many more autonomous systems (N), or m:N operations. Unfortunately, the human performance literature suggests that human operators will exhibit poor supervision of numerous autonomous agents due to the limits of attentional resources in the operators. In the general human information-processing model, a human operator exercises a limited pool of attentional resources to engage various information-processing stages including detecting, perceiving, comprehending, and predicting objects around them. Yamani and Horrey (2018) expanded the human information-processing model to characterize a tradeoff between information-processing demand and resource relief that automation brings in the context of automated driving. In their model, a driver interacting with an automated driving system is assumed to reallocate resources “freed” by automation to support other information-processing stages required for successful task performance. A future AAM ecosystem enabled by an orchestration of advanced automated systems, however, requires a single operator to interact with more than one air vehicle with varying levels and degrees of automated systems, making the traditional framework of human-automation interaction insufficient. To address this gap, we provide a review of the literature on situation assessment and trust, two constructs identified as critical for a fuller understanding of intimate and intricate interactions between a human operator and multiple air vehicles equipped with increasingly autonomous systems. Then, we propose an expansion of Yamani and Horrey’s (2018) model to motivate systematic research on the human operator’s role, identify factors that influence resource allocation and guide human-centered design of an interface supporting the m:N operation in the AAM environment.

Advanced Air Mobility↗
NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods Including Discontinuous Galerkin

As is well-known, popular high-order methods such as discontinuous Galerkin (DG) and flux reconstruction (FR) tend to generate oscillations near shocks and discontinuities, which can lead to negative pressure or density, ultimately causing code breakdown. For standard second-order methods, one approach to mitigating oscillations is to impose constraints on the calculated slopes, ensuring they do not become excessively steep. This limiting process is rooted in the idea of preserving monotonicity introduced by Van Leer (1974): when the data are monotone, limiting the slopes result in a monotone piecewise linear solution. However, a drawback of this approach is the loss of accuracy near extrema, where the non-monotone solution is somewhat flattened and achieves only first-order accuracy.

numerical methods↗
DOE OSTI2020

Revisiting Theoretical Limits for One-Degree-of-Freedom Wave Energy Converters: Preprint

This work revisits the theoretical limits of one-degree of freedom wave energy converters (WECs). This paper considers the floating sphere used in the Ocean Energy Systems (OES) Task 10 Wave Energy Converter modeling and verification effort for analysis. Analytical equations are derived to determine bounds on the motion amplitude, time-averaged power, and power take-off (PTO) force. A unique result was found that shows the time-averaged power absorbed by a WEC can be defined solely by the inertial properties and the radiation hydrodynamic coefficients. In addition, a unique expression for the PTO force amplitude was derived which has provided upper and lower bounds when resistive control is used to maximize power generation. For complex conjugate control, this same expression can only provide a lower bound as there is theoretically no upper bound. These bounds are used to compare the performance of a floating sphere if it were to extract energy using surge or heave motion. The analysis shows that because of the differences in hydrodynamic coefficients of each oscillating mode there will be different frequency ranges that provide better power capture efficiency. The influence of a motion constraint on power absorption while also utilizing a nonideal power take-off is examined and found to reduce the losses associated with bidirectional energy flow. The expression to calculate the time-averaged power with a nonideal PTO is modified by the mechanical-to-electrical efficiency and the ratio of the PTO spring and damping coefficients. These were separated in the expression which allows for limits to be set on the possible values of PTO coefficients to ensure a net flow of power to the grid.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2021

Structure Formation Limits on Axion-Like Dark Matter

We derive structure formation limits on dark matter (DM) composed of keV-scale axion-like particles (ALPs), produced via freeze-in through the interactions with photons and Standard Model (SM) fermions. We employ Lyman-alpha (Ly-α) forest data sets as well as the observed number of Milky Way (MW) subhalos. We compare results obtained using Maxwell-Boltzmann and quantum statistics for describing the SM bath. It should be emphasized that the presence of logarithmic divergences complicates the calculation of the production rate, which can not be parameterized with a simple power law behaviour. The obtained results, in combination with X-ray bounds, exclude the possibility for a photophilic “frozen-in” ALP DM with mass below ~ 19. For the photophobic ALP scenario, in which DM couples primarily to SM fermions, the ALP DM distribution function is peaked at somewhat lower momentum and hence for such realization we find weaker limits on DM mass. Future facilities, such as the upcoming Vera C. Rubin observatory, will provide measurements with which the current bounds can be significantly improved to ~ 80.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Atomic Layer Deposition for Materials-Based H2 Storage: Opportunities and Limitations

The transportation demands in our growing hydrogen (H2) economy requires robust storage systems. The current commercially implemented technology in fuel cell cars relies on the well-established technology of hydrogen gas compressed to 350-700 bar, depending on the application. The compressed gas tanks in use today are bulky and cost intensive. To address this challenge, material-based storage is one of the long-term alternatives considered and constitutes the focus of this talk. Material-based storage is broadly defined as hydrogen bound to solid materials, with its binding strength varying from physisorption to porous materials, such as zeolites and metal organic frameworks, to chemisorption in (complex) metal hydrides. The ultimate targets set by the U.S. Department of Energy for this technology include a system gravimetric capacity of 6.5 wt% and volumetric capacity of 40 g/L at 100 bar, operating temperatures ranging between -40 C and +40 C and adsorption/desorption timescales of < 5 min. Storage in the form of physisorbed or chemisorbed hydrogen has guided the materials research, which metal- organic framework and (complex) metal hydrides being the most promising materials classes. A variety of these materials have met one or more of the targets, but it has remained elusive for a single material system to meet all these stringent requirements. Nano-encapsulation and low-concentration chemical additives have previously been employed separately to overcome such challenges. Functionalization via atomic layer deposition (ALD), however, offers unique characteristics that make it suitable for both, nano-encapsulation and "doping" with low-concentration additives. This deposition technique has sub-monolayer thickness control, is highly conformal in high-surface area materials and is self-limiting, i.e., once the gas-phase precursor reacts with the available surface sites, the surface reactions stop. In this presentation, we will showcase examples where ALD, more generally vapor-phase functionalization, on (complex) metal hydrides and organic frameworks has improved the material properties for H2 storage. In our first study, Al2O3 was deposited on magnesium borohydride, Mg(BH4)2, at room temperature using trimethylaluminum (TMA) and water. From our findings, encouraging initial results were obtained: the H2 desorption temperature was lowered by 60-120 degrees C, the desorbed gravimetric H2 capacity at temperatures < 250 degrees C was doubled, and the desorption kinetics increased by a factor of ~6 compared to uncoated Mg(BH4)2. However, hydrolysis reactions caused by residual surface -OH groups from the ALD water-pulse degraded the sample substantially. Through this study, the use of TMA was observed to be highly reactive with the Mg(BH4)2 surface species, leading to the strategy of exposing the sample only to TMA. The relative mole fraction of the vapor-phase additive can be precisely controlled with the number of TMA pulses and pulse duration. By tuning these parameters, we show in our second study that 10 pulses of TMA at ambient conditions were able to decrease the H2 desorption temperature by ~100 degrees C while retaining <95 % of its H2 capacity. We applied this approach with other additives such as BBr3, TiCl4 and tetrahydrofuran, and demonstrate that this unique approach opens the door to a new class of molecular additives and catalysts which cannot easily be introduced with conventional mechano-chemical or solvent-based techniques for (complex) metal-hydrides. In the case of metal organic frameworks, ALD is a promising technique to functionalize the pores with metal atoms or functional groups able to tune the gas selectivity and binding energy, for which ~15 kJ/mol has been established as the optimal value for H2 storage in sorbent materials. Vapor-phase techniques such as ALD have numerous benefits over other functionalization tools for H2 storage materials opening innumerable opportunities to the field. To exploit these opportunities, the current limitations on room temperature and water-less ALD processes needs to be overcome which will greatly expand the possibilities of encapsulation and incorporation of additives for organic frameworks and (complex) metal hydrides.

atomic layer deposition↗
DOE OSTI2025

Power and Limitations of Linear Programming Decoder for Quantum LDPC Codes

Decoding quantum error-correcting codes is a key challenge in enabling fault-tolerant quantum computation. In the classical setting, linear programming (LP) decoders offer provable performance guarantees and can leverage fast practical optimization algorithms. Although LP decoders have been proposed for quantum codes, their performance and limitations remain relatively underexplored. In this work, we uncover a key limitation of LP decoding for quantum low-density parity-check (LDPC) codes: certain constant-weight error patterns lead to ambiguous fractional solutions that cannot be resolved through independent rounding. To address this issue, we incorporate a post-processing technique known as ordered statistics decoding (OSD), which significantly enhances LP decoding performance in practice. Our results show that LP decoding, when augmented with OSD, can outperform belief propagation with the same post-processing for intermediate code sizes of up to hundreds of qubits. These findings suggest that LP-based decoders, equipped with effective post-processing, offer a promising approach for decoding near-term quantum LDPC codes.

Gu, Shouzhen [Yale U.]↗
DOE OSTI2025

Power and Limitations of Linear Programming Decoder for Quantum LDPC Codes

Decoding quantum error-correcting codes is a key challenge in enabling fault-tolerant quantum computation. In the classical setting, linear programming (LP) decoders offer provable performance guarantees and can leverage fast practical optimization algorithms. Although LP decoders have been proposed for quantum codes, their performance and limitations remain relatively underexplored. In this work, we uncover a key limitation of LP decoding for quantum low-density parity-check (LDPC) codes: certain constant-weight error patterns lead to ambiguous fractional solutions that cannot be resolved through independent rounding. To address this issue, we incorporate a post-processing technique known as ordered statistics decoding (OSD), which significantly enhances LP decoding performance in practice. Our results show that LP decoding, when augmented with OSD, can outperform belief propagation with the same post-processing for intermediate code sizes of up to hundreds of qubits. These findings suggest that LP-based decoders, equipped with effective post-processing, offer a promising approach for decoding near-term quantum LDPC codes.

Gu, Shouzhen [Yale U.]↗
DOE OSTI2025

Understanding the limitations and potential of micro tensile testing of tungsten and needs for crystal plasticity modelling

Tungsten is the leading plasma facing material candidate due to its exceptional properties. Understanding the response to neutron irradiation is crucial for the lifetime evaluation of tungsten. The limited space in nuclear reactors and the high levels of radioactivity of the specimens after irradiation are significant barriers to accurately measuring the mechanical properties after neutron irradiation. Testing of micro tensile specimens is one approach to reduce the total amount of irradiated material needed for a set of mechanical testing experiments. However, micro samples are not necessarily measuring the bulk material properties, as size effects produce higher measured mechanical properties than is observed in engineering size counterparts. Determining the minimal sample size for reliable bulk property measurement, combined with modeling efforts, is essential. We conducted room temperature tensile tests on tungsten micro tensile specimens fabricated with focused ion beam, plasma focused ion beam and femto-second laser ablation system to dimension of 2×2×7?µm³, 7×7×18?µm³, and 80×100×233?µm³ (width×thickness×gauge-length), respectively. Elevated tensile testing was performed up to 475°C on 5×5×18?µm³ specimens fabricated with focused ion beam. The smallest specimens exhibited a high degree of ductility and strength, whereas the largest specimens demonstrated behavior akin to bulk tungsten. To investigate the effects introduced by the micro specimen fabrication processes and to obtain the necessary bulk dislocation density for the crystal plasticity model, we employed X-ray diffraction depth profiling. This measurement was performed using different X-ray sources with varying penetration depths. Finally, the potential and limitations of micro mechanical tests for tungsten will be discussed.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Power and Limitations of Linear Programming Decoder for Quantum LDPC Codes

Decoding quantum error-correcting codes is a key challenge in enabling fault-tolerant quantum computation. In the classical setting, linear programming (LP) decoders offer provable performance guarantees and can leverage fast practical optimization algorithms. Although LP decoders have been proposed for quantum codes, their performance and limitations remain relatively underexplored. In this work, we uncover a key limitation of LP decoding for quantum low-density parity-check (LDPC) codes: certain constant-weight error patterns lead to ambiguous fractional solutions that cannot be resolved through independent rounding. To address this issue, we incorporate a post-processing technique known as ordered statistics decoding (OSD), which significantly enhances LP decoding performance in practice. Our results show that LP decoding, when augmented with OSD, can outperform belief propagation with the same post-processing for intermediate code sizes of up to hundreds of qubits. These findings suggest that LP-based decoders, equipped with effective post-processing, offer a promising approach for decoding near-term quantum LDPC codes.

Gu, Shouzhen [Yale U.]↗
DOE OSTI2026

An improved direct limit on the muon electric dipole moment

A limit on the permanent electric dipole moment (EDM) of the positive muon is presented based on data from the Fermilab Muon g-2 Experiment taken between 2019 and 2020. The tracking detectors measure the average vertical decay angle of positrons from muon decays, enabling a search for an interaction between a possible muon EDM $d_μ$ and the lab-frame magnetic field. The result, $d_μ= (-0.35 \pm 0.19_{\mathrm{stat}} \pm 0.34_{\mathrm{sys}}) \times10^{-19}~e\cdot$cm, is consistent with zero and sets a new direct limit on the muon EDM of $|d_μ|<1.10\times10^{-19}~e\cdot$cm at the 95 percent confidence level.

Aguillard, D. P. [Michigan U.] (ORCID:000000025556↗
DOE OSTI2020

Kinetic Limitations in Single-Crystal High-Nickel Cathodes

High-nickel cathodes attract immense interest for use in lithium-ion batteries, to boost Li-storage capacity while reducing cost. For overcoming the intergranular-cracking issue in polycrystals, single-crystals are considered an appealing alternative, but aggravating concerns on compromising the electrochemical kinetics of electrodes. We report here a quantitative assessment of electrochemical reaction in single-crystal LiNi 0.8 Mn 0.1 Co 0.8 O 2 using operando hard X-ray microscopy/spectroscopy, revealing a strong dependence of redox kinetics on the state of charge (SOC) - being sluggish at low SOC while increasing rapidly as SOC increases both in the bulk electrode and individual particles. The observation is corroborated by transport measurements and finite-element simulation, indicating that the sluggish kinetics in single-crystals is governed by ionic transport at low SOC and may be alleviated through synergistic interaction with polycrystals integrated into a same electrode. Finally, with the elucidated mechanistic origin and alleviation solution of kinetic limitations in single-crystal high-Ni cathodes, this work constitutes a step forward in enabling their practical deployment.

25 ENERGY STORAGE↗
DOE OSTI2020

Electronic couplings for singlet fission: Orbital choice and extrapolation to the complete basis set limit

Abstract For the search for promising singlet fission candidates, the calculation of the effective electronic coupling, which is required to estimate the singlet fission rate between the initially excited state (S 0 S 1 ) and the multiexcitonic state ( 1 TT, two triplets on neighboring molecules, coupled into a singlet), should be sufficiently reliable and fast enough to explore the configuration space. We propose here to modify the calculation of the effective electronic coupling using a nonorthogonal configuration interaction approach by: (a) using only one set of orbitals, optimized for the triplet state of the molecules, to describe all molecular electronic states, and (b) only taking the leading configurations into consideration. Furthermore, we also studied the basis set convergence of the electronic coupling, and we found, by comparison to the complete basis set limit obtained using the cc‐pV n Z series of basis sets, that both the aug‐cc‐pVDZ and 6–311++G** basis sets are a good compromise between accuracy and computational feasibility. The proposed approach enables future work on larger clusters of molecules than dimers.

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

Quantitative XPS of plutonium: Evaluation of the Pu4f peak shape, relative sensitivity factors and estimated detection limits

We report that high-resolution X-ray photoelectron spectra have been acquired from sputter cleaned and in situ oxidized samples of α-plutonium and a face-centred cubic δ-plutonium–gallium alloy. The differences in the Pu4f peak shape in these two materials have been investigated, and the poorly screened satellite peaks have been quantified. Curve fitting parameters for the Pu4p 3/2 , Pu4d 5/2 , Pu4f and Pu6p 3/2 photoelectron peaks are reported, and relative sensitivity factors have been determined. The Pu4f curve fit model has been applied to data acquired using different spectrometers and alloys. Examples of quantification of the plutonium spectra are provided and minimum detection limits are calculated for common impurities in plutonium metal.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating the impact of the ionomer equivalent weight on a platinum group metal-free PEMFC cathode via oxygen limiting current

Leveraging the interactions between ionomer and catalyst can increase the performance of proton exchange membrane fuel cells. The impacts of the equivalent weight (EW) of perfluorosulfonic acid–based ionomers on the platinum group metal-free electrode structure and fuel cell performance have not been fully explored. Four membrane electrode assemblies (MEAs) were prepared by using a commercial Fe–N–C catalyst, two perfluorosulfonic acid ionomers with different EWs, that is, Aquivion 720 (A720) and Nafion 1100 (N1100), and two ionomer-to-catalyst (I/C) ratios. The four MEAs were characterized to understand the impact of the ionomer EW and content on the capacitance, proton conductivity, and mass transport on the cathode. The mass transport resistance was measured for the first time using a new oxygen reduction reaction limiting current method enabling to couple the effects of oxygen diffusion with liquid water generation. Low EW ionomer combined with a moderate I/C results in improved performance due to its enhanced proton conductivity. However, when used at high I/C, it can cause severe water flooding at high current density due to the enhanced liquid water uptake, especially at high relative humidity, resulting in lower catalyst utilization and higher mass transport resistance.

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

Well-Balanced Second-Order Convex Limiting Technique for Solving the Serre–Green–Naghdi Equations

In this article, we introduce a numerical method for approximating the dispersive Serre–Green–Naghdi equations with topography using continuous finite elements. The method is an extension of the hyperbolic relaxation technique introduced in Guermond et al. (J Comput Phys 450:110809, 2022). It is explicit, second-order accurate in space, third-order accurate in time, and is invariant-domain preserving. It is also well balanced and parameter free. Special attention is given to the convex limiting technique when physical source terms are added in the equations. The method is verified with academic benchmarks and validated by comparison with laboratory experimental data.

97 MATHEMATICS AND COMPUTING↗