Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Voids”

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 415 records · Page 23

Photonic Crystal Surface Emitting GaSb-Based Type-I Quantum Well Diode Lasers

The GaSb-based epitaxially regrown monolithic diode PCSELs operating near 2 μm at room temperature in continuous wave regime and generating 30 mW of output power from 200 μm diameter aperture have been designed and fabricated. Here, the devices demonstrated CW threshold current density of about 500 A/cm 2 . The laser output power was enhanced thanks to increased buried void area fill-factor in the photonic crystal layer with multiple voids per unit cell. The PCSEL generated ultra-low divergence donut shape beams at the currents near threshold. At higher injection currents, the device brightness was limited by excitation of the higher order lateral modes. Generation of the vector-vortex beams of different types by different band edge states of the buried photonic crystal was observed.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Oxidation and microstructure of SiC f /SiC composites in moist air up to 1600°C by X‐ray tomographic characterization

Abstract SiC f /SiC composites that possess PyC or BN interface layers were fabricated and then oxidized in moist air at 1000, 1200, 1400, and 1600°C. High‐resolution CT was used for capturing 3D images and quantifying the SiC phase, mesophase, and voids. The oxidation behavior and microstructural evolution of SiC f /SiC with PyC or BN interface are discussed in this study. The microstructure of the SiC f /SiC with a PyC layer was seriously damaged in moist air at high temperature, whereas the BN interface layer enhanced the oxidation resistance of the SiC f /SiC. These results are also confirmed by using XRD, oxidation mass gain, tensile testing, and SEM measurements. The results of the oxidation behavior and microstructural evolution for SiC f /SiC oxidized in dry air are also compared with the results of this study. Comparing the SiC f /SiC with a PyC interface layer, the composite with a BN interface layer oxidized in moist air exhibits a high void growth rate and a low SiO 2 grain growth rate from 1000 to 1600°C. This work will provide guidance for predicting the service life of SiC f /SiC for multiscale damage rate models of materials at a local scale and will also provide guidance on the life service design of SiC f /SiC materials.

Lu, Zilong↗

Data-Efficient Dimensionality Reduction and Surrogate Modeling of High-Dimensional Stress Fields

Tensor datatypes representing field variables like stress, displacement, velocity, etc., have increasingly become a common occurrence in data-driven modeling and analysis of simulations. Numerous methods [such as convolutional neural networks (CNNs)] exist to address the meta-modeling of field data from simulations. As the complexity of the simulation increases, so does the cost of acquisition, leading to limited data scenarios. Modeling of tensor datatypes under limited data scenarios remains a hindrance for engineering applications. Here, in this article, we introduce a direct image-to-image modeling framework of convolutional autoencoders enhanced by information bottleneck loss function to tackle the tensor data types with limited data. The information bottleneck method penalizes the nuisance information in the latent space while maximizing relevant information making it robust for limited data scenarios. The entire neural network framework is further combined with robust hyperparameter optimization. We perform numerical studies to compare the predictive performance of the proposed method with a dimensionality reduction-based surrogate modeling framework on a representative linear elastic ellipsoidal void problem with uniaxial loading. The data structure focuses on the low-data regime (fewer than 100 data points) and includes the parameterized geometry of the ellipsoidal void as the input and the predicted stress field as the output. The results of the numerical studies show that the information bottleneck approach yields improved overall accuracy and more precise prediction of the extremes of the stress field. Additionally, an in-depth analysis is carried out to elucidate the information compression behavior of the proposed framework.

artificial intelligence↗

Convective Heat Transfer Potential of Particles/Airflow Through Single Cell Thick Additively Manufactured Octet-Shaped Lattice Frame Material

Abstract Particle-to-Supercritical Carbon Dioxide (sCO2) heat exchangers are one of the most critical components of the next-generation Concentrating Solar Power (CSP) plants. There have been several efforts to enhance the overall heat exchanger performance which essentially comprises of thermal resistances offered by sCO2 channel, wall (separating sCO2 with particles) thickness, particle-wall contact resistance and particle-side effective heat transfer coefficient. This study is focused towards reducing the particle side thermal resistance by incorporating single unit cell thick reticulated Octet lattice frame structures on the falling particle side to enhance the effective thermal conductivity of the particle channel and to enhance convective heat transfer between the falling particles and the solid phase of the falling particle channel (endwalls and fibers). Steady-state experiments were conducted to measure the effective thermal conductivity of lattice frame material for two cases: a) when void space was occupied by air, b) when void space was occupied with particles. Further, convective heat transfer experiments have been conducted with both air (steady-state) and particles (quasi steady-state) as “working fluid” for panels sandwiching the Octet array. Three different lattice porosities ranging from 0.75 to 0.9 have been tested for a wide range of air flow rates and a fixed particle flow rate (highest potential).

Aider, Youssef↗

In situ measurement of damage evolution in shocked magnesium as a function of microstructure

Accurate modeling and prediction of damage induced by dynamic loading in materials have long proved to be a difficult task. Examination of postmortem recovered samples cannot capture the time-dependent evolution of void nucleation and growth, and attempts at analytical models are hindered by the necessity to make simplifying assumptions, because of the lack of high-resolution, in situ, time-resolved experimental data. We use absorption contrast imaging to directly image the time evolution of spall damage in metals at ~1.6-μm spatial resolution. We observe a dependence of void distribution and size on time and microstructure. The insights gained from these data can be used to validate and improve dynamic damage prediction models, which have the potential to lead to the design of superior damage-resistant materials.

36 MATERIALS SCIENCE↗

Stabilizing perovskite-substrate interfaces for high-performance perovskite modules

The interfaces of perovskite solar cells (PSCs) are important in determining their efficiency and stability, but the morphology and stability of imbedded perovskite-substrate interfaces have received less attention than have top interfaces. We found that dimethyl sulfoxide (DMSO), which is a liquid additive broadly applied to enhance perovskite film morphology, was trapped during film formation and led to voids at perovskite-substrate interfaces that accelerated the film degradation under illumination. Partial replacement of DMSO with solid-state carbohydrazide reduces interfacial voids. A maximum stabilized power conversion efficiency (PCE) of 23.6% was realized for blade-coated p-type/intrinsic/n-type (p-i-n) structure PSCs with no efficiency loss after 550-hour operational stability tests at 60°C. The perovskite mini-modules showed certified PCEs of 19.3 and 19.2%, with aperture areas of 18.1 and 50.0 square centimeters, respectively.

Chen, Shangshang↗

Surface Film Formation on Al-V Alloys with Far-From-Equilibrium Microstructure

Nanocrystalline supersaturated Al-V alloys produced by high-energy ball milling have been reported to exhibit enhanced corrosion resistance and mechanical properties compared to commercial Al alloys. Corrosion of passive alloys such as Al-V alloy relies on the characteristics of the surface film, which is studied using scanning/transmission electron microscopy and time-of-flight secondary ion mass spectrometry. The effect of microstructure and composition on the surface film has been investigated after different immersion periods (30 min, 2 h, and 1 day) in 0.1 M NaCl. The surface film was complex and composed of oxidized Al and V. The heterogeneous surface film was observed due to the presence of secondary phases and initiation of localized corrosion. The void formation was observed beneath the surface film that would potentially cause pitting corrosion. The generation of nano-sized voids was dependent on grain orientation. Compared to pure Al, the chloride penetration is suppressed in Al-V alloys. The effect of composition and microstructure on surface film formation and attendant corrosion behavior is discussed herein.

36 MATERIALS SCIENCE↗

Tomo2Mesh: Fast Reconstruction and Visualization of Tomography Data in Mesh Format

Tomo2Mesh is an open-source project targeted towards real-time reconstruction, segmentation, and visualization of computed tomography (CT) data in mesh format. The CT reconstruction scheme is based on filtered back-projection of voxel subsets. The segmentation scheme uses a 3D convolutional neural network. To allow for fast, real-time reconstruction, voxel subsets are first identified by coarse reconstruction. The detail in specific regions of interest is improved through full reconstruction of voxel subsets in that region. Data structures are implemented to store and process voxel subsets. Finally, voxel data is labeled using connected components to detect disconnected regions such as voids whose morphological attributes (e.g., Feret diameter, principal axis orientation, local number density, size, etc.) can be measured also in real-time. Finally, a fast marching cubes implementation processes labeled voxel data into a triangular face mesh (vertices and faces) in .ply format for visualization in Paraview or other mesh visualization tools. The code provides a simple programming interface for detecting, classifying, and visualizing regions of interest based on morphology. For example, detected voids can be classified as round pores or extended cracks. Highly porous neighborhoods can be identified based on local number density. The mesh texture (or color) is assigned based these morphological attributes to allow smart visualization scenarios in real-time (e.g., show only long cracks). At the time of first release (July 2022), extraction of face mesh for visualization for raw CT data from a 2 megapixel camera would take between 1-5 minutes for most scenarios.

TEKAWADE, ANIKET↗

Covariate Dependent Sparse Functional Data Analysis

This study proposes a method to incorporate covariate information into sparse functional data analysis. The method aims at cases where each subject has a limited number of longitudinal measurements and is associated with static covariates. This research is motivated by several use cases in practice. One representative example is void swelling, a nuclear-specific material degradation mechanism. Void swelling is affected by many covariates, including alloy composition and irradiation type. How to accurately model the complicated joint effects of such covariates on the swelling process is the key to mitigating the effect of swelling and ensuring safe operation. Unlike most of the existing methods, the proposed method can handle high-dimensional covariates with the informative covariate identification procedure and sparse and irregularly spaced measurements, that is, does not require complete or dense observations. The main innovation of the proposed method is that we model the variation coming from covariates and the variation left conditioned on covariates, such that the functional principal component analysis and Gaussian process can be conducted in a unified manner. Further, we also propose a systematic approach to identify important covariates in the hypothesis testing context. The methodology is demonstrated on applications in nuclear engineering and healthcare and simulation studies.

42 ENGINEERING↗

Improved fabrication of scalable hybrid silicon nitride–thin-film lithium niobate electro-optic modulators

Strategies for improved fabrication of integrated thin-film lithium niobate electro-optic (EO) Mach–Zehnder modulators (MZMs) are reported using scalable processes and designs. The MZM devices utilize direct bonding of unetched and unpatterned thin-film lithium niobate to patterned and planarized silicon photonic microchips. The latter contains silicon nitride waveguide structures of various widths that are used to form hybrid modes that are suitable for high-bandwidth low-voltage EO modulators based on Pockels effect. We report that the incorporation of appropriately designed outgassing channels and certain modifications to key processing steps helped achieve a greater than 99% reduction in void density during bonding. Void reduction is critically important for these traveling-wave hybrid MZM devices in which the optical mode is controllably distributed between multiple thin layers and propagates over millimeter-scale lengths.

Mere, Viphretuo↗

Bragg edge imaging (BEI) of B12W and M8N socket sections of the Arecibo telescope

From neutron user principal investigator: We kindly request the public release of three neutron imaging datasets through ONCat. All datasets were collected from two forensic specimens, B12W and M8N, sectioned from zinc-filled steel-wire sockets recovered from the collapsed Arecibo Telescope. The dataset titled “Neutron radiographs of B12W and M8N socket sections of the Arecibo telescope” contains normalized two-dimensional (2D) neutron radiographs of the specimens, showing the geometry and spatial distribution of the steel wires embedded within the zinc matrix, as well as internal features such as voids and cracks. The dataset titled “Neutron computed tomography of B12W and M8N socket sections of the Arecibo telescope” contains normalized 2D neutron projection images acquired over a range of specimen rotation angles for one selected region of each specimen. These projection images were used to reconstruct three-dimensional (3D) tomographic volumes that reveal the embedded-wire geometry and internal defects. The dataset titled “Bragg edge imaging (BEI) of B12W and M8N socket sections of the Arecibo telescope” contains six time-of-flight (TOF) neutron imaging datasets, three from each specimen, acquired at regions of interest selected based on the radiographs. The spatially resolved 2D TOF images show the zinc matrix and embedded steel wires, and the wavelength-dependent neutron transmission data were used to characterize crystallographic texture within the zinc. All components and their condition are in the public domain as they are the property of the National Science Foundation (NSF). The neutron imaging data, part geometries, and detailed forensic information have been widely published in the Arecibo Telescope Collapse Forensic Report by Thornton Tomasetti Engineers and others (NASA report and NASEM report).

Bilheux, Hassina↗

Carbide Particle Redistribution in U-10Mo Alloy during Hot and Cold Rolling Processes: Synthetic Microstructure-Based Finite-Element Approach

The influence of initial carbide particle distribution features on particle fracture and the evolution of stringer distribution within an alloy of uranium with 10 weight percent molybdenum during cold rolling was investigated in this study. This work will provide guidance to casting optimization work and identify actionable process changes to eliminate the surface features directly attributed to carbides. In a previous study, various synthetic microstructure-based finite-element (FE) models were generated with consideration of different parameters (i.e., particle shape, particle volume fraction, and particle distribution), and used to examine the effects of hot rolling on the microstructure evolution. In this study, microstructures after ~80% hot rolling were used as the initial microstructures for the subsequent cold rolling simulation. Plane-strain compression loading was employed in the FE models to simulate cold rolling with reductions from ~80% to ~94%. Stringer analyses with the consideration of particle fracture were then performed with the simulated cold rolling results. The results of void evolution analyses show that microstructures with elliptical particles incorporate a much larger (almost double) void volume fraction than those with circular particles, indicating that the elliptical particles are fractured more easily than the circular ones. The results from the fracture mode analyses show that, under the adopted simple compression loading assumption, almost all the particles (i.e., >99%) were fractured by tensile separation. The results of stringer analyses show that, during cold rolling, the stringer volume fractions increase substantially due to the particle fracture, especially for the microstructures with the elliptical particles, and that the stringer volume fractions from stringers comprising many particles tend to increase faster. The maximum stringer lengths for the microstructures with circular particles tend to remain constant or decrease during cold rolling, whereas those for elliptical particles tend to slightly increase after a slight decrease. The number of particles at the surface increases very rapidly due to particle fracture during cold rolling.

36 MATERIALS SCIENCE↗

Preparation of Mesophase Pitch Feedstock for Carbon Fiber

A filtration process of the low cost mesophase pitch developed by Advanced Carbon Products, LLC (ACP) was developed. The mesophase pitch was melt spun into fiber at the University of Kentucky Center for Applied Energy Research (UKY), and converted to carbon fiber. The melt spinning process was moderately stable, but resulted in some void structure in the green fibers. The tensile properties of the resulting carbon fibers, after graphitization heat treatment, were tested yielding a modulus value of approximately 560 GPa or 81 MSI. This demonstrated that a high modulus graphitic carbon fiber is certainly possible from the ACP 20 mesophase pitch. For future work, a better balance between the softening point temperature, the spinning temperature, and filtering must be met to reduce the volatiles generated in the pitch during spinning. This would reduce voids in the fiber, facilitate spinning, and generate finer, higher performance graphitic fiber.

36 MATERIALS SCIENCE↗

Development and Validation of a Two-Phase Thermal-Hydraulic CFD Code NEK-2P

A project is underway to develop, verify and validate an advanced two-phase flow modeling capability for the highly-scalable, high-performance Computational Fluid Dynamics (CFD) code NEK5000. The goal of this work is to verify and validate the two-phase version of the NEK5000 code, named NEK-2P, to simulate the two-phase flow and heat transfer phenomena that occur in a Boiling Water Reactor (BWR) fuel bundle under various operating conditions. The NEK-2P two-phase flow models follow the approach used for the Extended Boiling Framework (EBF) previously developed at Argonne but include more fundamental physical models of boiling phenomena and advanced numerical algorithms for improved computational accuracy, robustness, and computational speed. The development of the NEK-2P two-phase solver and the implementation of the Extended Boiling Framework two-phase models were initially supported by Argonne National Laboratory (Argonne) through a Laboratory Directed Research and Development (LDRD) project during FY2014-2016. The development and validation of the two-phase models through analyses of selected two-phase boiling flow experiments was supported by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY2017-2020. This report focuses on verification and validation of the water-steam boiling model NEK-2P Two-Phase, CFD code. The NEK-2P was validated with Nuclear Power Engineering Corporation (NUPEC) Pressurized Water Reactor (PWR) Sub-channel and Bundle Test (PSBT) void distribution benchmark. Three different simulations were performed and analyzed for various operating conditions such as wall-heat flux and sub-cooled inlet temperatures. Reasonably good agreement with measured data was obtained in predicting the measured void distributions. Simulations were performed for Virginia Tech. (VT) 3x3 rod bundle geometry with and without spacers. The preliminary results were presented for Simplified Spacer Grid (SSG). In addition, the implementation of interface reconstruction model was tested with one of the Becker benchmark Critical Heat Flux (CHF) experiments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development of CTF modeling of interfacial drag, wall shear, and interfacial heat transfer for bubbly and annular-mist flow regimes

This milestone report discusses improvements to CTF closure models for modeling of bubbly and annular mist two-phase flow conditions, which are important for accurate modeling and simulation for boiling water reactor (BWR) conditions. Models that were improved include the two-phase wall shear, wall boiling, and phase mass and heat transfer. A new flow regime map was implemented to improve on the prediction of transition to annular flow, which was generally largely over-predicted by the legacy CTF flow regime map. The annular mist interfacial drag models were also reviewed and new models were proposed for implementation into the code. In addition to implementing new models, an activity was performed to calibrate these new models to experimental data using a Markov Chain Monte Carlo technique, thus further improving their accuracy. Finally, new experiments were implemented into the CTF validation matrix which add more two-phase pressure drop and void measurements for the assessment of CTF for modeling of BWRs. The results show that the Chisholm model, which was implemented and then calibrated, results in significantly improved comparison against the Riso two-phase pressure drop data. Furthermore, improvements in the implicit heat transfer and wall boiling models resulted in good agreement with the void measurement data.

42 ENGINEERING↗

IACMI Project 4.2: Thermoplastic Composite Development for Wind Turbine Blades

(Section 5.1) Composites made from Arkema’s Elium® thermoplastic resin and Johns Manville fiberglass were researched during this project for applications in wind blade manufacturing. A techno-economic model was developed to model this wind blade manufacturing process using these materials in place of traditional composites made with thermoset resin. This model was based on manufacturing a 61.5-meter wind blade, which showed a 4.7% reduction in wind blade cost as compared traditional thermoset materials. These cost savings were not from the thermoplastic material costing less than traditional thermoset materials, but rather from decreased capital costs, faster cycle times and reduced energy requirements and labor costs. (Section 5.2) An infusion and curing model was developed for thermoplastic composite wind blades using PAM-RTM. The primary goal was to demonstrate the infusion simulation for the Elium® resin system on a 13-meter wind blade. Additionally, the exotherm temperature was predicted and compared to measurements, which showed model results within 10% of actual measurements. (Section 5.3) Composite laminate panels and composite sandwich panels with a balsa core were produced; specimens were cut and characterized. Similar composite specimens were made with Elium® thermoplastic resin and Hexion thermoset epoxy (RIMR135/RIMH1366) to enable comparisons between these resin systems. The static test methods included: tensile, compression, in-plane shear, interlaminar shear, flexural, sandwich core shear flexure, and single cantilever beam tests for sandwich beams. Fatigue testing at room temperature was completed to composite laminate panels at a stress ratio of R=0.1 and R=10. In addition, fatigue testing to laminate panels was completed at -30°C, and at room temperature after conditioning specimens at 70°C and 90% relative humidity. Overall, mechanical test results from Elium® composites are similar to epoxy composites. (Section 5.4) Elium composite panels were produced with intentional defects such as voids and nonwetting of fibers to begin to understand performance sensitivity to defects. A thermal digital image correlation (TDIC) method provides high spatial resolution strain field at elevated temperatures and can be used to identify defective regions within composite panels. Flexural modulus differences of 21% were seen between defect and non-defect panels. Other Elium® composite panels were forced to be defective by boiling the resin after infusion, which created voids throughout the composite laminate. X-ray computed tomography scanning was used to view the internal structure of the defect panels. Defect panels had a significant reduction in fatigue life as compared to baseline panels produced without intentional defects. (Section 5.5) Lap shear specimens were fabricated to compare the lap shear strength of an off-the-shelf adhesive (Plexus MA590) and two new adhesives developed by Arkema (Bostik SAF30 90 and Bostik SAF30 120). ISO standard 4587:2003 was used to standardize the testing method and sample fabrication. Lap shear specimens were made at 1mm, 3mm, and 10mm thicknesses. The Bostik adhesive lap shear test results were similar to Plexus for all thicknesses. (Section 5.6) Fiber-reinforced polymer (FRP) composites are typically used in high-performance applications (e.g., aerospace), and their expansion into high-volume industries (e.g. consumer automotive and wind turbine blade manufacturer or similar) is hindered by their cost and a lack of efficient manufacturing techniques. Monitoring the curing process of these composites during manufacturing can improve the efficiency of the process, and therefore reduce the manufacturing cost. Cure monitoring techniques were developed that use probabilistic estimation methods and surface temperature measurements made using infrared cameras. These techniques enable real-time monitoring of the infusion process to locate manufacturing flaws, and they can, potentially, estimate residual stresses in the part. Their commercialization will help facilitate expansion of FRP composites in high-volume industries. (Section 5.7) A 13-meter composite wind blade was produced with Elium® resin and Johns Manville fiberglass; this blade was made with VARTM processing similar to how megawatt-scale wind blades are currently manufactured, but no post-mold heating was used for this thermoplastic composite blade. The wind blade underwent full-scale validation for static loading (4-different load orientations) and flapwise fatigue loading to simulate 20-years of operational loads. The thermoplastic composite wind blade withstood the loading without any noted issues and performed similar to results from a previous full-scale validation to an equivalent epoxy composite wind blade produced with the same blade molds. (Section 5.8) A study was conducted to determine the feasibility of recycling composite wind turbine blade components fabricated with glass fiber reinforced Elium® thermoplastic resin. Dissolution, which is a process unique to thermoplastic matrices, allows recovery of both the polymer matrix and full-length glass fibers, while maintaining their stiffness and strength throughout the recovery process. The economics of recycling is favorable if 50% of the glass fiber is recovered and resold for a process of $\$$ 0.28/kg, and 90% of the resin is recovered and resold at a price of $\$$ 2.50/kg.(Section 10) Recommendations are outlined for commercializing thermoplastic resin for composite wind blade production, in addition to recommended areas for future research.

17 WIND ENERGY↗

HyMARC Seedling: Optimized Hydrogen Adsorbents via Machine Learning and Crystal Engineering

This final report is presented in two parts: In Part 1, the H 2 capacities of a diverse set of 918,734 metal-organic frameworks (MOFs) sourced from 19 databases is predicted via machine learning (ML). Using only 7 structural features as input, ML identifies 8,282 MOFs with the potential to exceed the capacities of state-of-the-art materials. The identified MOFs are predominantly hypothetical compounds having low densities (<0.31 g/cm3) in combination with high surface areas (>5,300 m2/g), void fractions ($0.90), and pore volumes (>3.3 cm3/g). The relative importance of the input features are characterized, and dependencies on the ML algorithm and training set size are quantified. The most important features for predicting H 2 uptake are pore volume (for gravimetric capacity) and void fraction (for volumetric capacity). The ML models are available on the web, allowing for rapid and accurate predictions of the hydrogen capacities of MOFs from limited structural data; the simplest models require only a single crystallographic feature. In part 2, ways to improve the poor powder packing density of MOFs is discussed. More specifically, a strategy that improves packing efficiency and volumetric hydrogen gas storage density dramatically through engineered morphologies and controlled-crystal size distributions is presented that holds promise for maximizing storage capacity for a given MOF. The packing density improvement, demonstrated for the benchmark sorbent MOF-5, leads to a significant enhancement of volumetric hydrogen storage performance relative to commercial MOF-5. System model projections demonstrate that engineering of crystal morphology/size or use of a bimodal distribution of cubic crystal sizes in tandem with system optimization can surpass the 25 g/L volumetric capacity of a typical 700 bar compressed storage system and exceed the DOE targets 2020 volumetric capacity (30 g/L). Finally, a critical link between improved powder packing density and reduced damage upon compaction is revealed leading to sorbents with both high surface area and high density.

08 HYDROGEN↗

Truchas-Flow

The truchas flow algorithm is designed for simulating the behavior of incompressible (but not necessarily constant density) fluids on unstructured meshes while allowing for solid and void phases to serve as moving boundaries. The interface between fluids, solids and void is tracked via a geometric volume of fluid approach. This allows for a single fluid formulation of the governing equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗