Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “control volume methods”

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 325 records · Page 18

Liquid rocket booster study. Volume 2, book 4, appendices 6-8: Reports of Rocketdyne, Pratt and Whitney, and TRW

For the pressure fed engines, detailed trade studies were conducted defining engine features such as thrust vector control methods, thrust chamber construction, etc. This was followed by engine design layouts and booster propulsion configuration layouts. For the pump fed engines parametric performance and weight data was generated for both O2/H2 and O2/RP-1 engines. Subsequent studies resulted in the selection of both LOX/RP-1 and O2/H2 propellants for the pump fed engines. More detailed analysis of the selected LOX/RP-1 and O2/H2 engines was conducted during the final phase of the study.

Source record↗

SU-E-I-90: Characterizing Small Animal Lung Properties Using Speckle Observed with An In-Line X-Ray Phase Contrast Benchtop System

Purpose: We demonstrate a novel X-ray phase-contrast (XPC) method for lung imaging representing a paradigm shift in the way small animal functional imaging is performed. In our method, information regarding airway microstructure that is encoded within speckle texture of a single XPC radiograph is decoded to produce 2D parametric images that will spatially resolve changes in lung properties such as microstructure sizes and air volumes. Such information cannot be derived from conventional lung radiography or any other 2D imaging modality. By computing these images at different points within a breathing cycle, dynamic functional imaging will be readily achieved without the need for tomography. Methods: XPC mouse lung radiographs acquired in situ with an in-line X-ray phase contrast benchtop system. The lung air volume is varied and controlled with a small animal ventilator. XPC radiographs will be acquired for various lung air volume levels representing different phases of the respiratory cycle. Similar data will be acquired of microsphere-based lung phantoms containing hollow glass spheres with known distributions of diameters. Image texture analysis is applied to the data to investigate relationships between texture characteristics and airspace/microsphere physical properties. Results: Correlations between Fourier-based texture descriptors (FBTDs) and regional lung air volume indicate that the texture features in 2D radiographs reveal information on 3D properties of the lungs. For example, we find for a 350 × 350 πm2 lung ROI a linear relationship between injected air volume and FBTD value with slope and intercept of 8.9×10{sup 5} and 7.5, respectively. Conclusion: We demonstrate specific image texture measures related to lung speckle features are correlated with physical characteristics of refracting elements (i.e. lung air spaces). Furthermore, we present results indicating the feasibility of implementing the technique with a simple imaging system design, short exposures, and low dose which provides potential for widespread use in laboratory settings for in vivo studies. This research was supported in part by NSF Award CBET1263988.

60 APPLIED LIFE SCIENCES↗

SU-E-T-436: Fluence-Based Trajectory Optimization for Non-Coplanar VMAT

Purpose: To investigate a fluence-based trajectory optimization technique for non-coplanar VMAT for brain cancer. Methods: Single-arc non-coplanar VMAT trajectories were determined using a heuristic technique for five patients. Organ at risk (OAR) volume intersected during raytracing was minimized for two cases: absolute volume and the sum of relative volumes weighted by OAR importance. These trajectories and coplanar VMAT formed starting points for the fluence-based optimization method. Iterative least squares optimization was performed on control points 24° apart in gantry rotation. Optimization minimized the root-mean-square (RMS) deviation of PTV dose from the prescription (relative importance 100), maximum dose to the brainstem (10), optic chiasm (5), globes (5) and optic nerves (5), plus mean dose to the lenses (5), hippocampi (3), temporal lobes (2), cochleae (1) and brain excluding other regions of interest (1). Control point couch rotations were varied in steps of up to 10° and accepted if the cost function improved. Final treatment plans were optimized with the same objectives in an in-house planning system and evaluated using a composite metric - the sum of optimization metrics weighted by importance. Results: The composite metric decreased with fluence-based optimization in 14 of the 15 plans. In the remaining case its overall value, and the PTV and OAR components, were unchanged but the balance of OAR sparing differed. PTV RMS deviation was improved in 13 cases and unchanged in two. The OAR component was reduced in 13 plans. In one case the OAR component increased but the composite metric decreased - a 4 Gy increase in OAR metrics was balanced by a reduction in PTV RMS deviation from 2.8% to 2.6%. Conclusion: Fluence-based trajectory optimization improved plan quality as defined by the composite metric. While dose differences were case specific, fluence-based optimization improved both PTV and OAR dosimetry in 80% of cases.

60 APPLIED LIFE SCIENCES↗

Three-Dimensional High-Order Spectral Volume Method for Solving Maxwell's Equations on Unstructured Grids

A three-dimensional, high-order, conservative, and efficient discontinuous spectral volume (SV) method for the solutions of Maxwell's equations on unstructured grids is presented. The concept of discontinuous 2nd high-order loca1 representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) method, but instead of using a Galerkin finite-element formulation, the SV method is based on a finite-volume approach to attain a simpler formulation. Conventional unstructured finite-volume methods require data reconstruction based on the least-squares formulation using neighboring cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In the SV method, one starts with a relatively coarse grid of triangles or tetrahedra, called spectral volumes (SVs), and partition each SV into a number of structured subcells, called control volumes (CVs), that support a polynomial expansion of a desired degree of precision. The unknowns are cell averages over CVs. If all the SVs are partitioned in a geometrically similar manner, the reconstruction becomes universal as a weighted sum of unknowns, and only a few universal coefficients need to be stored for the surface integrals over CV faces. Since the solution is discontinuous across the SV boundaries, a Riemann solver is thus necessary to maintain conservation. In the paper, multi-parameter and symmetric SV partitions, up to quartic for triangle and cubic for tetrahedron, are first presented. The corresponding weight coefficients for CV face integrals in terms of CV cell averages for each partition are analytically determined. These discretization formulas are then applied to the integral form of the Maxwell equations. All numerical procedures for outer boundary, material interface, zonal interface, and interior SV face are unified with a single characteristic formulation. The load balancing in a massive parallel computing environment is therefore easier to achieve. A parameter is introduced in the Riemann solver to control the strength of the smoothing term. Important aspects of the data structure and its effects to communication and the optimum use of cache memory are discussed. Results will be presented for plane TE and TM waves incident on a perfectly conducting cylinder for up to fifth order of accuracy, and a plane wave incident on a perfectly conducting sphere for up to fourth order of accuracy. Comparisons are made with exact solutions for these cases.

Liu, Yen↗

DEVELOPMENT OF A SCALABLE THERMOMECHANICAL PROCESS (TMP) FOR THICK BAINITIC ARMOR STEEL

A bainite phase-based alloy and associated thermomechanical process were developed to produce (2.5 to over 5 cm) thick armor-grade steel with uniform through-thickness high hardness and strength. The alloy composition and the final-critical (austenite to bainite) isothermal transformation step were specifically designed to utilize a simple and versatile air-cool/quench method to keep a low upfront capital cost and to provide the ability to continuously control the cooling rate in real time, in order to produce maximum volume fraction of bainite phase, and promote uniformly distributed strength and hardness. Final thickness of 2.5 cm and 5.7 cm steel plates were fabricated for characterization, testing and evaluation and found to possess uniform through-thickness hardness between 53 to 55 HRC and dynamic compressive strength of up to 2 GPa.

36 MATERIALS SCIENCE↗

QProR: An Efficient Framework for Quantity-of-Interest Based Progressive Retrieval with Guaranteed Error Control

Scientific applications generate an unprecedented volume of data, overwhelming the network and file systems’ bandwidth and posing challenges for efficient and scalable data retrieval and analysis. Progressive data compression offers a promising solution by enabling on-demand retrieval at reduced size. However, existing progressive methods either fail to bound the errors in essential quantities of interest (QoIs) derived from raw data or suffer from suboptimal retrieval efficiency. In this work, we propose QProR, an efficient QoI-based progressive framework that optimizes progressive retrieval for target QoIs. Our key contributions include: (1) a systematic framework that integrates error-controlled lossy compressors with bitplane encoding while decoupling the two processes for high flexibility and adaptability; (2) a novel weighted bitplane encoding method which incorperates QoI knowledge into data refactoring to enhance retrieval efficiency; (3) an optimized retrieval strategy that accounts for the varying impacts of different variables on multivariate QoIs; (4) comprehensive evaluations using six real-world datasets from multiple scientific applications and thorough comparisons against state of the arts. Experimental results demonstrate that QProR achieves up to 80.38% reduction in the retrieval size under the same requested QoI error tolerance, when compared with the best-performing existing methods. When transferring 384 GB of scientific data to remote sites, QProR delivers up to 1.68 × speedup in the end-to-end data transfer performance.

Li, Wenbo [University of Kentucky]↗

Low level CO2 effects on pulmonary function in humans

The purpose of the study was to determine whether chamber exposure to low levels of CO2 results in functional alterations in gas mixing and closing volume in humans. Four healthy volunteer subjects were exposed to 0.7% CO2 and to 1.2% CO2. Spirometry, lung volumes, single breath nitrogen washout, diffusing capacity for carbon monoxide (DLCO) by two methods, and cardiac output were measured in triplicate. Values were obtained over two non-consecutive days during the training period (control) and on days 2 or 3, 4, 6, 10, 13, and 23 of exposure to each CO2 level. Measurements were made during the same time of day. There was one day of testing after exposure, while still in the chamber but off carbon dioxide. The order of testing, up until measurements of DLCO and cardiac output, were randomized to avoid presentation effects. The consistent findings were a reduction in diffusing capacity for carbon monoxide and a fall in cardiac output, occurring to a similar degree with both exposures. For the group as a whole, there was no indication of major effects on spirometry, lung volumes, gas mixing or dead space. We conclude that small changes may occur in the function of distal gas exchanging units; however, these effects were not associated with any adverse health effects. The likelihood of pathophysiologic changes in lung function or structure with 0.7 or 1.2% CO2 exposure for this period of time, is therefore, low.

NASA Discipline Cardiopulmonary↗

WE-G-BRE-09: Targeted Radiotherapy Enhancement During Accelerated Partial Breast Irradiation (ABPI) Using Controlled Release of Gold Nanoparticles (GNPs)

Purpose: Several studies have demonstrated low rates of local recurrence with brachytherapy-based accelerated partial breast irradiation (APBI). However, long-term outcomes on toxicity (e.g. telangiectasia), and cosmesis remain a major concern. The purpose of this study is to investigate the dosimetric feasibility of using targeted non-toxic radiosensitizing gold nanoparticles (GNPs) for localized dose enhancement to the planning target volume (PTV) during APBI while reducing dose to normal tissue. Methods: Two approaches for administering the GNPs were considered. In one approach, GNPs are assumed to be incorporated in a micrometer-thick polymer film on the surface of routinely used mammosite balloon applicators, for sustained controlled in-situ release, and subsequent treatment using 50-kVp Xoft devices. In case two, GNPs are administered directly into the lumpectomy cavity e.g. via injection or using fiducials coated with the GNP-loaded polymer film. Recent studies have validated the use of fiducials for reducing the PTV margin during APBI with 6 MV beams. An experimentally determined diffusion coefficient was used to determine space-time customizable distribution of GNPs for feasible in-vivo concentrations of 43 mg/g. An analytic calculational approach from previously published work was employed to estimate the dose enhancement due to GNPs (2 and 10 nm) as a function of distance up to 1 cm from lumpectomy cavity. Results: Dose enhancement due to GNP was found to be about 130% for 50-kVp x-rays, and 110% for 6-MV external beam radiotherapy, 1 cm away from the lumpectomy cavity wall. Higher customizable dose enhancement could be achieved at other distances as a function of nanoparticle size. Conclusion: Our preliminary results suggest that significant dose enhancement can be achieved to residual tumor cells targeted with GNPs during APBI with electronic brachytherapy or external beam therapy. The findings provide a useful basis for developing nanoparticle-aided APBI, with potential to significantly reduce normal tissue/skin toxicity.

60 APPLIED LIFE SCIENCES↗

A Modified LC/MS/MS Method with Enhanced Sensitivity for the Determination of Scopolamine in Human Plasma

Intranasal scopolamine is a choice drug for the treatment of motion sickness during space flight because of its quick onset of action, short half-life and favorable sideeffects profile. The dose administered usually ranges between 0.1 and 0.4 mg. Such small doses make it difficult to detect concentrations of scopolamine in biological fluids using existing sensitive LC/MS/MS method, especially when the biological sample volumes are limited. To measure scopolamine in human plasma to facilitate pharmacokinetic evaluation of the drug, we developed a sensitive LC/MS/MS method using 96 well micro elution plates for solid phase extraction (SPE) of scopolamine in human plasma. Human plasma (100-250 micro L) were loaded onto Waters Oasis HLB 96 well micro elution plate and eluted with 50 L of organic solvent without evaporation and reconstitution. HPLC separation of the eluted sample was performed using an Agilent Zorbax SB-CN column (50 x 2.1 mm) at a flow rate of 0.2 mL/min for 3 minutes. The mobile phase for separation was 80:20 (v/v) methanol: ammonium acetate (30 mM) in water. Concentrations of scopolamine were determined using a Micromass Quattro Micro(TM) mass spectrometer with electrospray ionization (ESI). ESI mass spectra were acquired in positive ion mode with multiple reaction monitoring for the determination of scopolamine m/z = 304.2 right arrow 138.1 and internal standard hyoscyamine m/z = 290.2 right arrow 124.1. The method is rapid, reproducible, specific and has the following parameters: scopolamine and the IS are eluted at about 1.1 and 1.7 min respectively. The linear range is 25-10000 pg/mL for scopolamine in human plasma with correlation coefficients greater than 0.99 and CV less than 0.5%. The intra-day and inter-day CVs are less than 15% for quality control samples with concentrations of 75,300, and 750 pg/mL of scopolamine in human plasma. SPE using 96 well micro elution plates allows rapid sample preparation and enhanced sensitivity for the LC/MS/MS determination of scopolamine in a small volume of biological samples. The new method is also cost effective since it uses a small volume of organic solvents compared to the methods using SPE cartridges or regular 96 well SPE plates. This method can be successfully used for bioavailability and pharmacokinetic evaluations of scopolamine, especially when volumes of biological samples are limited. Further investigation to use automated SPE system with 96 well micro elution plates is planned.

Wang, Zuwei↗

Liquid and gaseous oxygen safety review, volume 1

Materials used or contained in liquid and gaseous oxygen systems are analyzed for their compatibility; and areas of possible concern in oxygen systems are outlined. Design criteria, cleaning procedures, and quality control methods are covered in detail.

Lapin, A.↗

Subplane-based Control Rod Decusping Techniques for the 2D/1D Method in MPACT

The MPACT transport code is being jointly developed by Oak Ridge National Laboratory and the University of Michigan to serve as the primary neutron transport code for the Virtual Environment for Reactor Applications Core Simulator. MPACT uses the 2D/1D method to solve the transport equation by decomposing the reactor model into a stack of 2D planes. A fine mesh flux distribution is calculated in each 2D plane using the Method of Characteristics (MOC), then the planes are coupled axially through a 1D NEM-P$_3$ calculation. This iterative calculation is then accelerated using the Coarse Mesh Finite Difference method. One problem that arises frequently when using the 2D/1D method is that of control rod cusping. This occurs when the tip of a control rod falls between the boundaries of an MOC plane, requiring that the rodded and unrodded regions be axially homogenized for the 2D MOC calculations. Performing a volume homogenization does not properly preserve the reaction rates, causing an error known as cusping. The most straightforward way of resolving this problem is by refining the axial mesh, but this can significantly increase the computational expense of the calculation. The other way of resolving the partially inserted rod is through the use of a decusping method. This paper presents new decusping methods implemented in MPACT that can dynamically correct the rod cusping behavior for a variety of problems.

Graham, Aaron M↗

Microstructure and Properties of Thermally Sprayed Functionally Graded Coatings for Polymeric Substrates

The use of polymer matrix composites (PMCs) in the gas flow path of advanced turbine engines offers significant benefits for aircraft engine performance but their useful lifetime is limited by their poor erosion resistance. High velocity oxy-fuel (HVOF) sprayed polymer/cermet functionally graded (FGM) coatings are being investigated as a method to address this technology gap by providing erosion and oxidation protection to polymer matrix composites. The FGM coating structures are based on a polyimide matrix filled with varying volume fractions of WC-Co. The graded coating architecture was produced using a combination of internal and external feedstock injection, via two computer-controlled powder feeders and controlled substrate preheating. Porosity, coating thickness and volume fraction of the WC-Co filler retained in the coatings were determined using standard metallographic techniques and computer image analysis. The pull-off strength (often referred to as the adhesive strength) of the coatings was evaluated according to the ASTM D 4541 standard test method, which measured the greatest normal tensile force that the coating could withstand. Adhesive/cohesive strengths were determined for three different types of coating structures and compared based on the maximum indicated load and the surface area loaded. The nature and locus of the fractures were characterized according to the percent of adhesive and/or cohesive failure, and the tested interfaces and layers involved were analyzed by Scanning Electron Microscopy.

Ivosevic, M.↗

Spaceflight Autonomous Multigenerational Microbial Sequencer (SAMMS) in Support of Plant-Growth Systems

As the National Aeronautics and Space Association (NASA) begins to pursue long-duration space flights, they will need to be able to provide astronauts with a nutritious and reliable food source. To meet the administration’s goal of traveling to the Moon and Mars, astronauts will need to begin to grow their own food in space. To protect their food source, extensive monitoring will occur to test for the effects of a space flight environment (e.g., radiation) as well as for early pathogen and disease detection. Genomic sequencing allows for both concerns to be tested on a regular basis. However, NASA’s current sequencer is unable to process plant tissues. Therefore, a novel method for plant DNA extraction using microneedle (MN) patches that will be able to feed into NASA’s existing system, but also require minimal human input is proposed. To support this, the design was broken down into four components (1) MN patch fabrication (2) MN patch extraction, (3) automated sampling motion control, and (4) a processing module. The MN patch is fabricated using a custom mold with conically shaped needles. The mold is filled with Polyvinyl alcohol (PVA) solution and placed in a vacuum desiccator. The mold is left in the vacuum overnight until the patch is dry and ready for use. The protocol was tested with varying pressures, drying times, volume amounts, and preparation methods to determine if highquality needles can be produced. A MN is a method of DNA extraction where the patch is applied to a leaf, the needles penetrate the leaf, breaking the rigid plant cell wall to isolate the DNA. A protocol for this method of extraction was tested to ensure the patch could produce the needed yield and purity. The tests varied by the number of patches, number of applications, and plant type. To automate the MN extraction method, motion control will utilize two separate axis tables which move in the x and y directions. The y-axis table will have an end effector that fits a MN patch and will have the ability to apply the patch to the leaf sample. This end effector will also act as a lid for a downstream processing module. The other axis will position the leaf sample and processing container so that the patch can be applied accurately. The Joint Comprehensive Sequencing System (JCSS) module integrates all the components together. The output of this module feeds into the NASA Charged Information-Storage Polymer Preparation System (CHIPPS) for genomic sequencing. The extraction module operates using a series of syringes and tubing to pump the varying reagents needed for the extraction protocol. The results of the study proved that MN patches are a viable method of DNA extraction. While fabrication of high-quality needles was unsuccessful, the protocol was able to be further developed using centrifugation. The integrated design between the motion control and the JCSS enabled the potential for automation with a complete conceptual design and prototype. Future research and development for this study would include (1) further testing for fabrication (2) expanding the range of plant species compatible with the MN patch, and (3) building a working prototype for the integrated system.

Peter Ling↗

A method to account for variation of average compressor inlet pressure during instantaneous distortion analyses

Instantaneous distortion analyses compare a time-varying value of an index (or 'surge margin used up') with a critical level (or 'available surge margin' of the compressor) to determine inlet-engine compatibility. Unless freestream conditions or propulsion system controls are changing, it is generally assumed that the available surge margin of the compressor is accurately determined from the steady-state operating point. Results are presented which show that variations of average compressor inlet pressure may occur without changes in freestream conditions or propulsion system controls. The volume dynamics of the compressor will cause these pressure variations to be attenuated and delayed by the time they reach the exit. This will cause the compressor pressure ratio (and available surge margin) to vary with time. A method is presented to calculate the available surge margin as a function of time and incorporate it into an instantaneous distortion analysis. Results show that inlet pressure variations which cause only a small change at the compressor exit can cause a significant variation in the available surge margin.

Burstadt, P. L.↗

Tumor Shrinkage Assessed by Volumetric MRI in Long-Term Follow-Up After Fractionated Stereotactic Radiotherapy of Nonfunctioning Pituitary Adenoma

Purpose: To evaluate tumor control and side effects associated with fractionated stereotactic radiotherapy (FSRT) in the management of residual or recurrent nonfunctioning pituitary adenomas (NFPAs). Methods and Materials: We assessed exact tumor volume shrinkage in 16 patients with NFPA after FSRT. All patients had previously undergone surgery. Gross tumor volume (GTV) was outlined on contrast-enhanced magnetic resonance imaging (MRI) before and median 63 months (range, 28-100 months) after FSRT. MRI was performed as an axial three-dimensional gradient echo T1-weighted sequence at 1.6-mm slice thickness without gap (3D MRI). Results: Mean tumor size of all 16 pituitary adenomas before treatment was 7.4 mL (3.3-18.9 mL). We found shrinkage of the treated pituitary adenoma in all patients. Within a median follow-up of 63 months (28-100 months) an absolute mean volume reduction of 3.8 mL (0.9-12.4 mL) was seen. The mean relative size reduction compared with the volume before radiotherapy was 51% (22%-95%). Shrinkage measured by 3D MRI was greater at longer time intervals after radiotherapy. A strong negative correlation between the initial tumor volume and the absolute volume reduction after FSRT was found. There was no correlation between tumor size reduction and patient age, sex, or number of previous surgeries. Conclusions: By using 3D MRI in all patients undergoing FSRT of an NFPA, tumor shrinkage is detected. Our data demonstrate that volumetric assessment based on 3D MRI adds additional information to routinely used radiological response measurements. After FSRT a mean relative size reduction of 51% can be expected within 5 years.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Conservative fourth-order finite-volume Vlasov–Poisson solver for axisymmetric plasmas in cylindrical (r,v r ,v θ ) phase space coordinates

In this paper, a fourth-order finite-volume Vlasov–Poisson algorithm is developed for simulating axisymmetric plasma configurations in ( r , v r , v θ ) phase space coordinates. The Vlasov equation for cylindrical phase space coordinates is cast into conservation-law form and is discretized on a structured grid. The conservative finite-volume discretization is based on fifth-order upwind reconstructions of the distribution function and a fourth-order quadrature rule that accounts for transverse variations of fluxes along control volume surfaces. High-order specular reflection boundary conditions enable high-fidelity treatment of plasma distribution functions at the axis and at wall boundaries. The numerical method is applied to simulate a confined uniform neutral gas to assess convergence properties for an equilibrium system in which effects of finite-temperature and acceleration due to centrifugal and Coriolis forces are present. The discretization is also applied to a Z-pinch configuration to study electrostatic ion confinement and the dynamics of the associated breathing mode. Simulations show that the ion distribution function exhibits non-Maxwellian features and that a temperature anisotropy develops and is sustained. Finally, the finite-volume implementation is demonstrated to converge at fourth-order for both applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FVCOM-plume – A Three-dimensional Lagrangian Outfall Plume Dilution and Transport Model for Dynamic Tidal Environments: Model Development

Although numerous steady-state effluent plume dilution models are in use, their application in tidal environments remains a challenge. Three-dimensional dynamic circulation models are also inadequate, often due to the lack of required resolution and simplifying assumptions. To overcome these limitations, we present FVCOM-plume—an outfall plume dilution and transport model that operates within the Finite Volume Coastal Ocean Model (FVCOM) framework. It provides simultaneous inclusion of near-field dilution and far-field plume transport processes. The near-field is based on UM3 model using Lagrangian Control Volume approach to compute buoyant plume trajectory and dilution from multiport diffusers. The far-field uses neutrally buoyant particles with point masses and the random walk method to solve unsteady advection-diffusion processes. A density kernel approach is used to compute concentrations at point locations and analyze far-field plume characteristics. The results demonstrate the ability of FVCOM-plume to simultaneously capture near-field and far-field effluent plume dynamics in tidal environments.

Plume Dilution Models↗