Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetics Parameters”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Cellular Automaton Modeling of Dynamic Recrystallization of Nimonic 80A Superalloy Based on Inhomogeneous Distribution of Dislocations Inside Grains

A modified cellular automaton (CA) model in the time and space scale is developed to simulate the dynamic recrystallization (DRX) mechanism of Nimonic 80A nickel-base superalloy during hot deformation. Some necessary thermo-mechanical parameters used as input data in CA approach such as work hardening and softening coefficient are fitted according to the previous hot compression experiments. Based on the experimental results, nucleation kinetics is modified as a function of Zener-Holloman parameter. A dislocation density model at transient recrystallization boundary is proposed, in which the dislocation density is a function of temperature, strain rate and annealed dislocation density. The simulated results are in good agreement with experimental results, proving that the modified CA model is capable of predicting both flow behavior and microstructural evolution of Nimonic 80A during hot deformation. According to the experimental and simulated results, the DRX behavior is discussed in detail. It demonstrates that the flow stress and microstructural evolution of DRX have a relationship with the deformation temperature and strain rate.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Lutetium(III) aqua ion: On the dynamical structure of the heaviest lanthanoid hydration complex

The structure and dynamics of the lutetium(III) ion in aqueous solution have been investigated by means of a polarizable force field molecular dynamics (MD). An 8-fold square antiprism (SAP) geometry has been found to be the dominant configuration of the lutetium(III) aqua ion. Nevertheless, a low percentage of 9-fold complexes arranged in a tricapped trigonal prism (TTP) geometry has been also detected. Dynamic properties have been explored by carrying out six independent MD simulations for each of four different temperatures: 277 K, 298 K, 423 K, 632 K. The mean residence time of water molecules in the first hydration shell at room temperature has been found to increase as compared to the central elements of the lanthanoid series in agreement with previous experimental findings. Water exchange kinetic rate constants at each temperature and activation parameters of the process have been determined from the MD simulations. The obtained structural and dynamical results suggest that the water exchange process for the lutetium(III) aqua ion proceeds with an associative mechanism, in which the SAP hydration complex undergoes temporary structural changes passing through a 9-fold TTP intermediate. Such results are consistent with the water exchange mechanism proposed for heavy lanthanoid atoms.

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

Bulk Attenuation in the Earth’s Mantle Due to Phase Transitions

The calculated bulk attenuation due to phase transition from mineral physics data is reported here. With relaxation time less than 1 s, the calculated value for a pyrolite mantle is consistent with the inverted bulk attenuation of the upper mantle from seismic observations. The two important mechanisms of phase transitions, diffusion-controlled and kinetics-controlled, have different relaxation time as indicated by the models here. The diffusion controlled is more likely to contribute to the observed seismic bulk attenuation than the kinetic-controlled process based on the available diffusivity and kinetics data. The correlation between the bulk attenuation and relaxation time emphasizes the importance of a number of parameters in the mineral physics database such as Fe-Mg diffusivity and kinetics in olivine-wadsleyite-ringwoodite-perovskite, Mg-Ca-Al-Si diffusivity and kinetics in pyroxene-garnet-Ca perovskite, some of which are still unknown.

58 GEOSCIENCES↗
DOE OSTI2015

INTEGRAL FIELD SPECTROSCOPY OF AGN ABSORPTION OUTFLOWS: MRK 509 AND IRAS F04250–5718

Ultraviolet (UV) absorption lines provide abundant spectroscopic information enabling the probe of the physical conditions in active galactic nucleus (AGN) outflows, but the outflow radii (and the energetics consequently) can only be determined indirectly. We present the first direct test of these determinations using integral field unit (IFU) spectroscopy. We have conducted Gemini IFU mapping of the ionized gas nebulae surrounding two AGNs, whose outflow radii have been constrained by UV absorption line analyses. In Mrk 509, we find a quasi-spherical outflow with a radius of 1.2 kpc and a velocity of ∼290 km s{sup −1}, while IRAS F04250–5718 is driving a biconical outflow extending out to 2.9 kpc, with a velocity of ∼580 km s{sup −1} and an opening angle of ∼70°. The derived mass flow rate ∼5 and >1 M{sub ⊙} yr{sup −1}, respectively, and the kinetic luminosity ≳1 × 10{sup 41} erg s{sup −1} for both. Adopting the outflow radii and geometric parameters measured from IFU, absorption line analyses would yield mass flow rates and kinetic luminosities in agreement with the above results within a factor of ∼2. We conclude that the spatial locations, kinematics, and energetics revealed by this IFU emission-line study are consistent with pre-existing UV absorption line analyses, providing a long-awaited direct confirmation of the latter as an effective approach for characterizing outflow properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Enabling Structure/Interface Regulation for High Performance Ni-Rich Cathodes

Here, further commercialization of Ni-rich layered cathodes is hindered by severe structure/interface degradation and kinetic hindrance that occur during electrochemical operation, which leads to safety risks and reduced range in electric vehicles (EVs). Herein, by selecting elements with different solubility properties, a multifunctional strategy that synchronously fabricates perovskite-type SrZrO 3 coating and Sr/Zr co-doping is employed to strengthen the structure/interface stability and the Li + transport mobility of LiNi 0.85 Co 0.10 Mn 0.05 O 2 (NCM). Perovskite-type SrZrO 3 protective layers formed on the particle surface can substantially mitigate the unexpected interfacial side reactions and surface phase transitions. In addition, a robust crystal framework is constructed by optimizing local O coordination through the introduction of strong Zr-O bonds. Notably, Li + diffusion kinetics is effectively improved due to expanded cell parameters and O-Li-O slab spacing with the incorporation of large-diameter Sr pillar ions, as revealed by X-ray diffraction. As a result, the Sr/Zr-modified NCM achieves a remarkable capacity retention of 99.4% after 200 cycles at 1 C, and a high rate capacity of 168.9 mAh g -1 at 10 C. This work opens new avenues to develop high-performance NCM cathodes with high energy and high power for EVs with long calendar life.

25 ENERGY STORAGE↗
DOE OSTI2016

Variations of characteristic time scales in rotating stratified turbulence using a large parametric numerical study

In this work, we study rotating stratified turbulence (RST) making use of numerical data stemming from a large parametric study varying the Reynolds, Froude and Rossby numbers, Re, Fr and Ro in a broad range of values. The computations are performed using periodic boundary conditions on grids of 10243 points, with no modeling of the small scales, no forcing and with large-scale random initial conditions for the velocity field only, and there are altogether 65 runs analyzed in this paper. The buoyancy Reynolds number defined as R B = ReFr 2 varies from negligible values to ≈ 105, approaching atmospheric or oceanic regimes. This preliminary analysis deals with the variation of characteristic time scales of RST with dimensionless parameters, focusing on the role played by the partition of energy between the kinetic and potential modes, as a key ingredient for modeling the dynamics of such flows. We find that neither rotation nor the ratio of the Brunt-Väisälä frequency to the inertial frequency seem to play a major role in the absence of forcing in the global dynamics of the small-scale kinetic and potential modes. Specifically, in these computations, mostly in regimes of wave turbulence, characteristic times based on the ratio of energy to dissipation of the velocity and temperature fluctuations, T V and T P , vary substantially with parameters. Their ratio γ=T V /T P follows roughly a bell-shaped curve in terms of Richardson number Ri. It reaches a plateau --on which time scales become comparable, γ≈0.6 -- when the turbulence has significantly strengthened, leading to numerous destabilization events together with a tendency towards an isotropization of the flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Magnetic resonance insights into the heterogeneous, fractal-like kinetics of chemically recyclable polymers

Moving toward a circular plastics economy is a vital aspect of global resource management. Chemical recycling of plastics ensures that high-value monomers can be recovered from depolymerized plastic waste, thus enabling circular manufacturing. However, to increase chemical recycling throughput in materials recovery facilities, the present understanding of polymer transport, diffusion, swelling, and heterogeneous deconstruction kinetics must be systematized to allow industrial-scale process design, spanning molecular to macroscopic regimes. To develop a framework for designing depolymerization processes, we examined acidolysis of circular polydiketoenamine elastomers. We used magnetic resonance to monitor spatially resolved observables in situ and then evaluated these data with a fractal method that treats nonlinear depolymerization kinetics. This approach delineated the roles played by network architecture and reaction medium on depolymerization outcomes, yielding parameters that facilitate comparisons between bulk processes. These streamlined methods to investigate polymer hydrolysis kinetics portend a general strategy for implementing chemical recycling on an industrial scale.

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

An analysis of the uncertainty and bias in DCE-MRI measurements using the spoiled gradient-recalled echo pulse sequence

Purpose: The pharmacokinetic parameters derived from dynamic contrast-enhanced (DCE) MRI have been used in more than 100 phase I trials and investigator led studies. A comparison of the absolute values of these quantities requires an estimation of their respective probability distribution function (PDF). The statistical variation of the DCE-MRI measurement is analyzed by considering the fundamental sources of error in the MR signal intensity acquired with the spoiled gradient-echo (SPGR) pulse sequence. Methods: The variance in the SPGR signal intensity arises from quadrature detection and excitation flip angle inconsistency. The noise power was measured in 11 phantoms of contrast agent concentration in the range [0–1] mM (in steps of 0.1 mM) and in onein vivo acquisition of a tumor-bearing mouse. The distribution of the flip angle was determined in a uniform 10 mM CuSO{sub 4} phantom using the spin echo double angle method. The PDF of a wide range of T1 values measured with the varying flip angle (VFA) technique was estimated through numerical simulations of the SPGR equation. The resultant uncertainty in contrast agent concentration was incorporated in the most common model of tracer exchange kinetics and the PDF of the derived pharmacokinetic parameters was studied numerically. Results: The VFA method is an unbiased technique for measuringT1 only in the absence of bias in excitation flip angle. The time-dependent concentration of the contrast agent measured in vivo is within the theoretically predicted uncertainty. The uncertainty in measuring K{sup trans} with SPGR pulse sequences is of the same order, but always higher than, the uncertainty in measuring the pre-injection longitudinal relaxation time (T1{sub 0}). The lowest achievable bias/uncertainty in estimating this parameter is approximately 20%–70% higher than the bias/uncertainty in the measurement of the pre-injection T1 map. The fractional volume parameters derived from the extended Tofts model were found to be extremely sensitive to the variance in signal intensity. The SNR of the pre-injection T1 map indicates the limiting precision with which K{sup trans} can be calculated. Conclusions: Current small-animal imaging systems and pulse sequences robust to motion artifacts have the capacity for reproducible quantitative acquisitions with DCE-MRI. In these circumstances, it is feasible to achieve a level of precision limited only by physiologic variability.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Design of a multi beam klystron cavity from its single beam parameters

The klystron is a well-known microwave amplifier which uses kinetic energy of an electron beam for amplification of the RF signal. There are some limitations of conventional single beam klystron such as high operating voltage, low efficiency and bulky size at higher power levels, which are very effectively handled in Multi Beam Klystron (MBK) that uses multiple low purveyance electron beams for RF interaction. Each beam propagates along its individual transit path through a resonant cavity structure. Multi-Beam klystron cavity design is a critical task due to asymmetric cavity structure and can be simulated by 3D code only. The present paper shall discuss the design of multi beam RF cavities for klystrons operating at 2856 MHz (S-band) and 5 GHz (C-band) respectively. The design approach uses some scaling laws for finding the electron beam parameters of the multi beam device from their single beam counter parts. The scaled beam parameters are then used for finding the design parameters of the multi beam cavities. Design of the desired multi beam cavity can be optimized through iterative simulations in CST Microwave Studio.

42 ENGINEERING↗
DOE OSTI2014

AN OBSERVED CORRELATION BETWEEN THERMAL AND NON-THERMAL EMISSION IN GAMMA-RAY BURSTS

Recent observations by the Fermi Gamma-ray Space Telescope have confirmed the existence of thermal and non-thermal components in the prompt photon spectra of some gamma-ray bursts (GRBs). Through an analysis of six bright Fermi GRBs, we have discovered a correlation between the observed photospheric and non-thermal γ-ray emission components of several GRBs using a physical model that has previously been shown to be a good fit to the Fermi data. From the spectral parameters of these fits we find that the characteristic energies, E {sub p} and kT, of these two components are correlated via the relation E {sub p}∝T {sup α} which varies from GRB to GRB. We present an interpretation in which the value of the index α indicates whether the jet is dominated by kinetic or magnetic energy. To date, this jet composition parameter has been assumed in the modeling of GRB outflows rather than derived from the data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Numerical analysis of experimental studies of methane hydrate formation in a sandy porous medium

Here we analyse numerically an earlier experimental study that involved the formation of methane hydrates by the excess water method in a small reactor filled with a sandy porous medium, and seek to address questions about the type of the hydration reaction and the phase heterogeneity in the resulting hydrate-bearing sand. Using a fine discretization describing the reactor assembly, the experimental process is faithfully replicated numerically. The multi-stage process of hydrate formation is subdivided in 7 steps. The experimental data from the continuously-monitored pressure and temperature during each step are used for comparison against the numerical predictions, the identification of the dominant processes and the determination of the associated parameters through a history-matching process that minimizes deviations between observations and simulation results. The results of this first-ever study on this subject demonstrate unequivocally that the hydration reaction is a kinetic (as opposed to an equilibrium) process, and that the spatial distributions of the various phases (aqueous, gas and hydrate) at the end of the formation process are strongly heterogeneous. This has serious implications in simulation studies of hydrate dissociation that assume uniform initial phase saturation distributions. Finally, the history-matching process indicates that (a) the system behaviour is sensitive to some flow parameters (porosity and irreducible water saturation) only during the first water injection, (b) it is insensitive to the sand intrinsic permeability during all steps of the study, and (c) thermal processes dominate after the first water injection, yielding estimates of the thermal properties of the sand and of time-variable key parameters of the kinetic reaction.

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

Energy injection in gamma-ray burst afterglows

We present multi-wavelength observations and modeling of gamma-ray bursts (GRBs) that exhibit a simultaneous re-brightening in their X-ray and optical light curves, and are also detected at radio wavelengths. We show that the re-brightening episodes can be modeled by injection of energy into the blastwave and that in all cases the energy injection rate falls within the theoretical bounds expected for a distribution of energy with ejecta Lorentz factor. Our measured values of the circumburst density, jet opening angle, and beaming-corrected kinetic energy are consistent with the distribution of these parameters for long-duration GRBs at both z∼1 and z≳6, suggesting that the jet launching mechanism and environment of these events are similar to that of GRBs that do not have bumps in their light curves. However, events exhibiting re-brightening episodes have lower radiative efficiencies than average, suggesting that a majority of the kinetic energy of the outflow is carried by slow-moving ejecta, which is further supported by steep measured distributions of the ejecta energy as a function of Lorentz factor. We do not find evidence for reverse shocks over the energy injection period, implying that the onset of energy injection is a gentle process. We further show that GRBs exhibiting simultaneous X-ray and optical re-brightenings are likely the tail of a distribution of events with varying rates of energy injection, forming the most extreme events in their class. Future X-ray observations of GRB afterglows with Swift and its successors will thus likely discover several more such events, while radio follow-up and multi-wavelength modeling of similar events will unveil the role of energy injection in GRB afterglows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The effects of pre-plasma scale length on the relativistic electron beam directionality

The effects of pre-plasmas on the electron beam directionality was experimentally and numerically investigated. Single material and layered targets made of Ti and/or CH were used to simultaneously measure high-energy (≥3 MeV) electrons along two directions, pre-pulse energy and pre-plasma density. The electron directionality is quantified by using a new parameter, the electron energy ratio of the total kinetic energies along the two directions. Measurements and radiation–hydrodynamic (RH) simulations show that a large (≥3.5 μm) plasma scale length at the critical surface enhances electrons along the laser axis, and such pre-plasma conditions could only be achieved with the CH targets. Particle-in-cell simulations were performed on the RH generated pre-plasmas from Ti and CH targets, and the results show that the CH target provided conditions for higher forward momentum gains by electrons. First, the CH target allowed longer distances for electrons to interact with laser. Second, the intense laser pulse modified the critical surface, but the resulting surface differed. The CH target resulted in a smooth surface where a retro-reflection was observed while the Ti target resulted in a rippled surface that scattered the reflected light. As results, the CH electrons gained higher forward momentum via a direct-laser-acceleration in the counter propagating laser fields. So the results presented in this article show a way of controlling the high-energy electron directionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Mean force kinetic theory: A convergent kinetic theory for weakly and strongly coupled plasmas

A new closure of the BBGKY hierarchy is developed, which results in a convergent kinetic equation that provides a rigorous extension of plasma kinetic theory into the regime of strong Coulomb coupling. The approach is based on a single expansion parameter which enforces that the exact equilibrium limit is maintained at all orders. Because the expansion parameter does not explicitly depend on the range or the strength of the interaction potential, the resulting kinetic theory does not suffer from the typical divergences at short and long length scales encountered when applying the standard kinetic equations to Coulomb interactions. Here, the approach demonstrates that particles effectively interact via the potential of mean force and that the range of this force determines the size of the collision volume. When applied to a plasma, the collision operator is shown to be related to the effective potential theory [S. D. Baalrud and J. Daligault, Phys. Rev. Lett. 110, 235001 (2013)]. In addition to the collision operator, this systematic derivation reveals a second term that is related to the excess (nonideal) components of the pressure and internal energy in the hydrodynamic limit. The relationship between this and previous kinetic theories is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A quantitative phase-field model for gas bubble evolution in UO2

Due to the large formation energy of vacancies and noble gas atoms in the form of interstitials or substitutional atoms in nuclear fuel (UO2), the thermodynamic equilibrium concentrations of these species are very low in the nuclear fuel matrix even at very high temperatures, which imposes difficulties upon the quantitative study of bubble evolution via the phase-field method. In this study, a quantitative phase-field model is proposed to deal with this problem. The system’s free energy density is derived according to the principles of thermodynamics, with consideration of the elastic effect and with the use of real material parameters from experiments. The model is useful for the study of the kinetics of gas bubble growth with very dilute concentrations of vacancy and gas atoms in the matrix. This model is applied to study single bubble growth and multiple bubble growth under various concentrations of vacancy and gas atoms and at various temperatures. The elastic effect and the effects of the generation rate of vacancies and gas atoms on bubble growth are analyzed.

low defect concentration, nuclear fuel, gas bubble↗
DOE OSTI2015

Mathematical modeling of positron emission tomography (PET) data to assess radiofluoride transport in living plants following petiolar administration

Background: Ion transport is a fundamental physiological process that can be studied non-invasively in living plants with radiotracer imaging methods. Fluoride is a known phytotoxic pollutant and understanding its transport in plants after leaf absorption is of interest to those in agricultural areas near industrial sources of airborne fluoride. Here we report the novel use of a commercial, high-resolution, animal positron emission tomography (PET) scanner to trace a bolus of [¹⁸F]fluoride administered via bisected petioles of Brassica oleracea , an established model species, to simulate whole plant uptake of atmospheric fluoride. This methodology allows for the first time mathematical compartmental modeling of fluoride transport in the living plant. Radiotracer kinetics in the stem were described with a single-parameter free- and trapped-compartment model and mean arrival times at different stem positions were calculated from the free-compartment time-activity curves. Results: After initiation of administration at the bisected leaf stalk, [¹⁸F] radioactivity climbed for approximately 10 minutes followed by rapid washout from the stem and equilibration within leaves. Kinetic modeling of transport in the stem yielded a trapping rate of 1.5 +/- 0.3%/min (mean +/- s.d., n = 3), velocity of 2.2 +/- 1.1 cm/min, and trapping fraction of 0.8 +/- 0.5%/cm. Conclusion: Quantitative assessment of physiologically meaningful transport parameters of fluoride in living plants is possible using standard positron emission tomography in combination with petiolar radiotracer administration. Movement of free fluoride was observed to be consistent with bulk flow in xylem, namely a rapid and linear change in position with respect to time. Trapping, likely in the apoplast, was observed. Future applications of the methods described here include studies of transport of other ions and molecules of interest in plant physiology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Oxidation and pyrolysis of methyl propyl ether

The ignition, oxidation, and pyrolysis chemistry of methyl propyl ether (MPE) was probed experimentally at several different conditions, and a comprehensive chemical kinetic model was constructed to help understand the observations, with many of the key parameters computed using quantum chemistry and transition state theory. Experiments were carried out in a shock tube measuring time variation of CO concentrations, in a flow tube measuring product concentrations, and in a rapid compression machine (RCM) measuring ignition delay times. The detailed reaction mechanism was constructed using the Reaction Mechanism Generator software. Sensitivity and flux analyses were used to identify key rate and thermochemical parameters, which were then computed using quantum chemistry to improve the mechanism. Validation of the final model against the 1-20 bar 600-1500 K experimental data is presented with a discussion of the kinetics. The model is in excellent agreement with most of the shock tube and RCM data. Strong non-monotonic variation in conversion and product distribution is observed in the flow-tube experiments as the temperature is increased, and unusually strong pressure dependence and significant heat release during the compression stroke is observed in the RCM experiments. These observations are largely explained by a close competition between radical decomposition and addition to O2 at different sites in MPE; this causes small shifts in conditions to lead to big shifts in the dominant reaction pathways. The validated mechanism was used to study the chemistry occurring during ignition in a diesel engine, simulated using Ignition Quality Test (IQT) conditions. At the IQT conditions, where the MPE concentration is higher, bimolecular reactions of peroxy radicals are much more important than in the RCM.

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

RANS Simulation of Variable Density Turbulent Round Jets with Coflow using xRAGE Hydrodynamic Code and BHR Turbulence Models

This work for the fiscal year 2023 (FY23) is a continuation of previous efforts to evaluate the BHR turbulence models for their ability to accurately simulate variable density turbulent round jets with coflow. As before, RANS simulations are carried out using the xRAGE hydrodynamic code. The following are some of the previous findings. Israel showed that i) the symmetry boundary conditions for the BHR models in axisymmetric simulations were in error, ii) three grids of different resolutions did not lead to converging solutions, and iii) BHR 3.1 simulation exhibited instabilities and did not reach a steady state. Saenz and Rauenzahn derived and implemented into xRAGE the correct BHR boundary conditions at the symmetry axis. Cline conducted sensitivity studies with various parameters including the BHR models (versions 2, 3.1 and 4), gravity, material pressure, specific heat, initial turbulent kinetic energy and initial turbulent length scale, and found that the largest factor impacting on simulation results was the BHR model version, followed by the initial turbulent length scale. In addition, freeze boundary conditions at the exit and the side wall of the computational domain were used to remove anomalous flow behavior. Cline adjusted the inlet jet velocity, initial turbulent kinetic energy and initial turbulent length scale to obtain the best reasonable match with the experimental data of Charonko and Prestridge. It was found that the BHR 2 and 3.1 models performed in a similar manner, but the BHR 2 model produced much lower levels of density-specific-volume covariance and turbulent kinetic energy. The main focus for the FY23 is to study the effects of computational parameters related to the boundary conditions, mesh, domain size and timestep size. The reasoning behind this is that, unless simulation results can be shown to be reasonably independent from the aforementioned computational parameters, it would be difficult to attribute any discrepancies between simulation and experimental results to turbulence models. This important aspect has largely been overlooked in the previous years, and therefore will be studied comprehensively here. Additionally, effects of varying the initial turbulent length scale will be examined because it was previously identified as a major factor affecting the flow fields.

42 ENGINEERING↗