Search NASA⌕ Search

SEARCH · Search NASA

Results for “Compositional simulation”

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,549 records · Page 86

Structure–property relations of sodium iron phosphate nuclear waste glasses: Effects of iron redox ratio and glass composition

Iron phosphate glasses, known for their exceptional chemical durability and potential applicability in nuclear waste management, have gained significant attention over the years. The structures of these glasses are complicated by the coexistence of Fe 3+ and Fe 2+ , which plays a crucial role in determining their structures and properties. Here, this work uses molecular dynamics simulations to study the structural changes in Na 2 O–Fe 2 O 3 –P 2 O 5 glasses with varying glass composition and Fe 2+ /Fe 3+ redox ratio. It was found that the redox ratio and modifier contents significantly affected the short-range and medium-range orders in the glasses. Significant changes in the local environments around P 5+ and Fe 3+ were observed, as reflected by the bond distances and coordination numbers. Na + cations are found to preferentially associate with Fe 3+ (rather than Fe 2+ ), whereas Fe 2+ has stronger association with P 5+ than Na + , confirming the structural role of Fe 2+ as a glass modifier. The disruptions in P–O–P linkages upon increasing FeO suggest that FeO causes glass depolymerization. These glasses achieved higher connectivity with increasing Fe 3+ / (Fe 3+ + Fe 2+ ) ratios, conerting phosphorous Q 2 to Q 3 units and iron Q 5 units to Q 4 units. The decrease of nonbridging oxygen fractions with increasing Fe 3+ / (Fe 3+ + Fe 2+ ) ratios, through creating P–O–Fe linkages, is the main reason of enhanced network connectivity. Quantitative structure–property relationship analyses with different structural descriptors were used to correlate with measured properties. The analyses provided valuable insights into structure–property relationships, emphasizing the importance of choosing relevant energy parameters and defining glass network connectivity, particularly in F net descriptors. It was found the Fe–O–P linkage density exhibits strong correlations to measured dissolution rates, supporting the importance of these linkages in improving the chemical durability in iron phosphate glasses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A comparative study of the effectiveness of empirical potentials for molecular dynamics simulations of borosilicate glasses

Despite their practical importance, atomistic modeling of B 2 O 3 containing glasses have been challenging due to the lack of reliable empirical potentials. Fortunately, a few recent developments have shown promises to simulate these glasses where the boron coordination has complex and non-linear dependence on glass composition. This work aims to provide an evaluation of the effectiveness of three recently developed potentials by a systematic study of a series (~20) of sodium borosilicate glasses with constant K (ratio of [SiO 2 ]/[B 2 O 3 ]) =2 and varying R (ratio of [Na 2 O]/[B 2 O 3 ] ranging from 0.1 to 4) values and several sodium boroaluminosilicate glasses. A comparative assessment was established on the basis of the short- and medium-range structure features, such as boron N 4 values, total correlation functions, bond angle distribution, oxygen speciation, and mechanical properties using experimental or well-established models as criteria. Furthermore, this work provides insights on the choices of empirical potentials for MD simulations of borosilicate glasses and gives directions of future potential development and refining.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Characterization of the Sulfur-Saturated Melt Versions of the HFG1 Study Glasses

The U.S. Department of Energy (DOE) is responsible for building the Hanford Tank Waste Treatment and Immobilization Plant (WTP) at the Hanford Site in Washington to remediate 55 million gallons of radioactive waste that is temporarily stored in 177 underground tanks. The Office of River Protection (ORP) has requested that the Savannah River National Laboratory (SRNL) contribute to current glass studies to support successful startup of the WTP, due to SRNL’s recognized capabilities and expertise for glass waste form development. As part of this effort, the Pacific Northwest National Laboratory (PNNL) is developing, batching, and fabricating simulated high-level waste (HLW) glasses to generate the associated property data needed to validate or identify areas of potential characterization improvements in the current glass property/ composition models. Currently, the models cover only a small fraction of the waste compositions projected in the Hanford tank farm; therefore, the models must be expanded to include high fluoride compositions to successfully complete the WTP mission. SRNL support of this work is defined in the Task Technical and Quality Assurance Plan (TTQAP). This report provides results from the chemical analyses of a series of sulfur-saturated melt (SSM) versions of simulated nuclear waste glasses fabricated at PNNL. The glasses were selected as part of a broader study of the influence of glass composition on chemical durability, sulfur retention, and other properties. The glasses were designated the High Fluoride Glasses-1 (HLG1).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Polymer-cement composites with adhesion and re-adhesion (healing) to casing capability for geothermal wellbore applications

Deterioration of cement/casing adhesion in wellbore scenarios can result in unwanted and potentially harmful leakage with the potential of serious repair costs. In this work, we explore the use of self-healing polymers added to conventional wellbore cements as a way to bring about self-healing and readhering (to steel casing) properties to the composite material. The polymers are pH resistant and seem to improve the cement integrity after exposure to typical chemical and thermal stresses encountered under geothermal wellbore conditions. We find that addition of about 10-15 wt% of polymer to the cement visually increases its resistance to fracturing from exposure to geothermal conditions, while the adhesive strength of cement/stainless steel increases with curing time for a period of about 10 days. Self-healing capability was demonstrated by permeability analysis showing that polymer-cement composites reduce flow by 50-70% at cement bulk and at the cement/steel interface. Use of atomistic simulations imply that these polymers have good wetting properties on the steel surfaces. Analysis of the interactions between steel/polymer and cement/polymer show that they are complementary, resulting in a wider range of bonding patterns. Cracks are likely to expose under-coordinated sites that result in more bonding interactions, which agrees well with the permeability measurements showing high degree of healed cracks and healed (cement-steel) interfacial gaps together with an overall increased in structural integrity of these advanced polymer-cement composite materials.

Rod, Kenton A.↗

Evaluating Heat Extraction Performance of Closed-Loop Geothermal Systems with Thermally Conductive Enhancements in Conduction-Only Reservoirs: Preprint

We investigated the impact of disc-shaped and linear thermal enhancements on performance of closed-loop geothermal systems using steady-state and transient COMSOL numerical simulations. Thermal enhancements refer to thermally conductive material introduced in the rocks surrounding a wellbore to compensate for the relatively low rock thermal conductivity and increase heat extraction. Materials proposed include composite cements with thermally conductive additives such as metals or graphite to obtain overall cement thermal conductivities of 10 to 100x the rock thermal conductivity. Our simulation approach was modeling a subsegment (e.g., 100 m long) of a wellbore with idealized (e.g., perfect disc-shaped) thermal enhancements, and we report results as relative increase in heat extraction with respect to systems without thermal enhancements. We did not evaluate technical feasibility or cost of installing thermal enhancements. Simulations indicate that for 5-mm thick disc-shaped thermal enhancements with 5-m radius, repeated every 1 m along the wellbore and with thermal conductivity of 100x the rock thermal conductivity, the thermal output increases roughly 20% with respect to a closed-loop system without thermal enhancements. For fishbone structure thermal enhancements, we estimate a roughly 7% increase in thermal output for 5-m long radially outward pipes with half the radius of the main wellbore, repeated every 2.5 m, with thermal conductivity of 100x the rock thermal conductivity.

advanced geothermal systems↗

Evolution of A-Type Macrosegregation in Large Size Steel Ingot After Multistep Forging and Heat Treatment

A-type macrosegregation refers to the channel chemical heterogeneities that can be formed during solidification in large size steel ingots. In this research, a combination of experiment and simulation was used to study the influence of open die forging parameters on the evolution of A-type macrosegregation patterns during a multistep forging of a 40 metric ton (MT) cast, high-strength steel ingot. Macrosegregation patterns were determined experimentally by macroetch along the longitudinal axis of the forged and heat-treated ingot. Mass spectroscopy, on more than 900 samples, was used to determine the chemical composition map of the entire longitudinal sectioned surface. FORGE NxT 1.1 finite element modeling code was used to predict the effect of forging sequences on the morphology evolution of A-type macrosegregation patterns. For this purpose, grain flow variables were defined and implemented in a large scale finite element modeling code to describe oriented grains and A-type segregation patterns. Examination of the A-type macrosegregation showed four to five parallel continuous channels located nearly symmetrical to the axis of the forged ingot. In some regions, the A-type patterns became curved or obtained a wavy form in contrast to their straight shape in the as-cast state. Mass spectrometry analysis of the main alloying elements (C, Mn, Ni, Cr, Mo, Cu, P, and S) revealed that carbon, manganese, and chromium were the most segregated alloying elements in A-type macrosegregation patterns. The observed differences were analyzed using thermodynamic calculations, which indicated that changes in the chemical composition of the liquid metal can affect the primary solidification mode and the segregation intensity of the alloying elements. Finite element modeling simulation results showed very good agreement with the experimental observations, thereby allowing for the quantification of the influence of temperature and deformation on the evolution of the shape of the macrosegregation channels during the open die forging process.

36 MATERIALS SCIENCE↗

Composition Driven, Statistical Review of Cast Stone Screening Matrix Results

Recent studies in support of Cast Stone development for the Hanford Site were designed to screen the effects of mix components on properties and performance of the waste form. Dry blend components, including portland cement, blast furnace slag, and Class F fly ash, sourced from both southeastern and northwestern suppliers, were included in the screening matrix. Multiple decontaminated salt solutions (simulated) were combined with the dry blend components at various liquid to solid ratios following a statistically designed matrix of compositions. Properties of the resulting mixes and performance of the cured materials were then measured and related to the dry blend components and salt solutions that were used for each mix. At the conclusion of the screening matrix work, SRNL felt that additional insight into the study might be obtained by relating the properties and performance of the materials to the composition of each Cast Stone mix, as opposed to just the dry blend – salt solution combination. This memorandum provides the results of a statistically driven review of the screening matrix data as a function of the composition of each mix. It is meant to augment the results obtained by reviewing the screening matrix data as a function of the dry blend and salt solution combination.

Edwards, T. B.↗

High fidelity polycrystalline CdTe/CdS heterostructures via molecular dynamics

Molecular dynamics simulations of polycrystalline growth of CdTe/CdS heterostructures have been performed. First, CdS was deposited on an amorphous CdS substrate, forming a polycrystalline film. Subsequently, CdTe was deposited on top of the polycrystalline CdS film. Cross-sectional images show grain formation at early stages of the CdS growth. During CdTe deposition, the CdS structure remains almost unchanged. Concurrently, CdTe grain boundary motion was detected after the first 24.4 nanoseconds of CdTe deposition. With the elapse of time, this grain boundary pins along the CdS/CdTe interface, leaving only a small region of epitaxial growth. CdTe grains are larger than CdS grains in agreement with experimental observations in the literature. Crystal phase analysis shows that zinc blende structure dominates over the wurtzite structure inside both CdS and CdTe grains. Composition analysis shows Te and S diffusion to the CdS and CdTe films, respectively. Lastly, these simulated results may stimulate new ideas for studying and improving CdTe solar cell efficiency.

36 MATERIALS SCIENCE↗

Thermodynamic modeling of the Ti-Nb-Sn system and its preliminary applications in aging precipitation and solidification simulations

Abstract Gasping the thermodynamic and phase diagram of the materials being studied is essential for further guiding the development process, like composition design, heat treatment optimization and so on. Nb and Sn are common alloying elements in Ti-based biomedical alloys. In this work, the Ti-Nb-Sn system has been thermodynamically assessed based on experimental phase equilibria. A self-consistent thermodynamic description for the Ti-Nb-Sn system including two ternary compounds Ti 3 Nb 3 Sn 2 and Ti 3 NbSn has been obtained. Based on the comparison between the calculated and experimental phase equilibria, the reliability of the presently obtained thermodynamic description has been verified. The influence of Nb and Sn contents on the precipitation temperature and amount of hcp α -Ti phase in bcc β -(Ti, Nb) has been investigated based on thermodynamic calculations. Results show that the increase of Nb or Sn content in Ti-based alloys would significantly reduce the precipitation amount and transformation temperature of hcp α -Ti phase in bcc β -(Ti, Nb) and the influence of Nb content is more effective than that of Sn. The solidified phases in several as-cast alloys have been thermodynamically predicted based on Scheil solidification simulations. The presently developed thermodynamic description of the Ti-Nb-Sn system can serve as a sub-ternary one in the thermodynamic database for multicomponent Ti-based biomedical alloys.

Peng, Yingbiao↗

Persistence of systematic errors in the Asian-Australian monsoon precipitation basic states in climate models: a way forward

The annual cycle of the Asian-Australian monsoon (AAM) system can be regarded as the seasonal displacement of the large-scale Intertropical Convergence Zone (ITCZ), which is anchored by the north-south migration of the Indo-Pacific warm pool. In the respective hemispheres, intense solar heating over land during spring and early summer provide the necessary thermodynamic conditions for the occurrence of deep convection off the equator. The rainfall and diabatic heating associated with the AAM is perhaps the most vigorous of all the regional monsoon components in the globe. Yet, skill in monsoon prediction (days to seasons) by dynamical models remains low, partly due to our lack of understanding of the entirety of the monsoon system and our inability to model the interactive processes that govern it. Simulating the monsoon precipitation climatology remains a grand challenge. The multi-model mean (MMM) errors for summertime precipitation relative to GPCP observations have shown little improvement in CMIP5 as compared to CMIP3. For the Asian Summer Monsoon, the MMM monsoon rainfall is weaker over South Asia, the western equatorial Indian Ocean and tropical West Pacific. For the Australian summer monsoon, it is excessive over most of the Maritime Continent, and deficient over northern Australia. Throughout the year, excessive rainfall is simulated over South China Sea – Maritime Continent regions, and CMIP5 models do not capture the annual cycle of the AAM. One implication is that uncertainties in future projections (e.g., IPCC 2013) of AAM mean rainfall may not have reduced from CMIP3 to CMIP5. Solutions from an intermediate model show that diabatic heating (Q) associated with the AAM influences the global circulation. Compared to the reanalysis, many models tend to have maxima at the mid-troposphere but their simulated amplitude is overestimated in the lower troposphere (900-700 hPa) and underestimated from 600-300 hPa, a feature readily apparent in the MMM composite (dashed black line). The lower troposphere peak may be attributed to misrepresentations in shallow convection. Some outliers, such as CSIROMk3- 6-0 and ACCESS1-3, do not show any appreciable vertical structure, and the simulated monsoon over South Asia is virtually absent in these models. TRMM observations indicate that over the monsoon region the contribution from stratiform rainfall is about 40% to the Q intensity. In contrast, most of the CMIP3 models produce too much convective rainfall (95% of the total) and too little stratiform precipitation. Given the persistence of errors, we speculate that errors in the partitioning of total rainfall into convective-stratiform still persist in CMIP5, and may be one of the reasons for underestimation (overestimation) of Q in the layer 600-300 (900-700) hPa. In summary, the rather slow progress in modeling in the last decade or so led us to wonder: Has the scientific community reached a “plateau” in modeling mean monsoon precipitation?

54 ENVIRONMENTAL SCIENCES↗

Improving Stability of an Optical Fiber pH Sensor with a Calcined Polyethylenimine-Coating at High Pressures and Temperatures

With an increased interest in subsurface gas storage technology for various energy applications, monitoring wellbore structural stability and subsurface geochemistry has become more pressing, and pH is a key parameter to measure. As high pressure and elevated temperatures in subsurface conditions are comparatively harsh relative to that expected for most standard pH sensor designs, any pH monitoring hardware must be designed for extended exposure to high pressures and temperatures. We previously reported that an optical fiber pH sensor functionalized with a calcined polyethylenimine coating had shown some promise as a high temperature and pressure pH sensor but with some drifting when operating for longer than 8 hours. In this paper, we investigated the coating composition and potential cause of the drifting and improved the stability of the prepared coating to minimize sensor drift under simulated wellbore conditions. Scanning electron microscopy (SEM) had previously shown moderate cracking at high pressures over short tests. By applying X-ray photoelectron spectroscopy (XPS) to characterize the sensor coating before and after one week of testing in an H2/CH4 gas blend at 80 °C and 900 psi, a compositional change in the coating was observable, which may indicate susceptibility to alteration by subsurface gas storage conditions. Non-reducing (CH4, N2) environments were also tested, and confirmed that both temperature and pressure were also contributing to the drift.

energy infrastructure↗

Improving Stability of an Optical Fiber pH Sensor with a Calcined Polyethylenimine-Coating at High Pressures and Temperatures

With an increased interest in subsurface gas storage technology for various energy applications, monitoring wellbore structural stability and subsurface geochemistry has become more pressing, and pH is a key parameter to measure. As high pressure and elevated temperatures in subsurface conditions are comparatively harsh relative to that expected for most standard pH sensor designs, any pH monitoring hardware must be designed for extended exposure to high pressures and temperatures. We previously reported that an optical fiber pH sensor functionalized with a calcined polyethylenimine coating had shown some promise as a high temperature and pressure pH sensor but with some drifting when operating for longer than 8 hours. In this paper, we investigated the coating composition and potential cause of the drifting and improved the stability of the prepared coating to minimize sensor drift under simulated wellbore conditions. Scanning electron microscopy (SEM) had previously shown moderate cracking at high pressures over short tests. By applying X-ray photoelectron spectroscopy (XPS) to characterize the sensor coating before and after one week of testing in an H2/CH4 gas blend at 80°C and 900 psi, a compositional change in the coating was observable, which may indicate susceptibility to alteration by subsurface gas storage conditions. Non-reducing (CH4, N2) environments were also tested, and confirmed that both temperature and pressure were also contributing to the drift.

energy infrastructure↗

Black Hole Supernovae, Their Equation of State Dependence, and Ejecta Composition

Abstract Recent literature on core-collapse supernovae suggests that a black hole (BH) can form within ∼1 s of shock revival, while still culminating in a successful supernova. We refer to these as BH supernovae, as they are distinct from other BH formation channels in both timescale and impact on the explosion. We simulate these events self-consistently from core collapse until 20–50 days after collapse using three axisymmetric models of a 60 M ⊙ zero-age main-sequence progenitor star and investigate how the composition of the ejecta is impacted by the BH formation. We employ Skyrme-type equations of state (EOSs) and vary the uncertain nucleonic effective mass, which affects the pressure inside the proto–neutron star through the thermal part of the EOS. This results in different BH formation times and explosion energies at BH formation, yielding final explosion energies between 0.06 and 0.72 × 10 51 erg with 21.8–23.3 M ⊙ of ejecta, of which 0–0.018 M ⊙ is 56 Ni. Compared to expectations from 1D simulations, we find more nuanced EOS dependences of the explosion dynamics, the mass of the BH remnant, and the elemental composition of the ejecta. We investigate why the explosions survive despite the massive overburden and link the shape of the diagnostic energy curve and character of the ejecta evolution to the progenitor structure.

Eggenberger Andersen, Oliver (ORCID:00000002966079↗

Modeling, analysis, and code/data validation of DIII-D tokamak divertor experiments on ELM and non-ELM plasma tungsten sputtering erosion

Abstract We analyzed recent DIII-D tokamak tungsten divertor probe experiments using advanced, coupled, sputter erosion/redeposition, plasma, and surface response code packages. Modeling is done for ELMing H-mode, and L-mode plasmas, impinging on various size tungsten deposits on Divertor Material Evaluation System (DiMES) carbon probes. The simulations compute 3D, full kinetic, sub-gyromotion, impurity sputtering and transport, including changes in tungsten surface composition and response due to mixed deuterium and carbon ions irradiation. Per our analysis, ELM (edge localized mode) plasma sputtering in DIII-D mostly involves free-streaming high energy (∼500–1000 eV) D + and C +6 ions, with high near-surface plasma density. L-Mode sputtering is due to impurity sputtering (C, W) only, with lower density. All cases show complete redeposition of tungsten on the divertor, with significant redeposition on the tungsten spots themselves, and low self-sputtering. Comparison of ELM plasma gross tungsten erosion simulation results with in-situ spectroscopic data is good, as are code/data comparisons of net erosion using post-exposure Rutherford backscattering (RBS) data for the L-mode probes. The analysis, extrapolated to a full tungsten divertor, implies low net erosion and negligible plasma contamination from sputtering. These results support the use of high-Z plasma facing surfaces in ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulated changes in biogenic VOC emissions and ozone formation from habitat expansion of Acer Rubrum (red maple)

A new vegetation trend is emerging in northeastern forests of the United States, characterized by an expansion of red maple at the expense of oak. This has changed emissions of biogenic volatile organic compounds (BVOCs), primarily isoprene and monoterpenes. Oaks strongly emit isoprene while red maple emits a negligible amount. This species shift may impact nearby urban centers because the interaction of isoprene with anthropogenic nitrogen oxides can lead to tropospheric ozone formation and monoterpenes can lead to the formation of particulate matter. Here in this study the Global Biosphere Emissions and Interactions System was used to estimate the spatial changes in BVOC emission fluxes resulting from a shift in forest composition between oak and maple. A 70% reduction in isoprene emissions occurred when oak was replaced with maple. Ozone simulations with a chemical box model at two rural and two urban sites showed modest reductions in ozone concentrations of up to 5–6 ppb resulting from a transition from oak to red maple, thus suggesting that the observed change in forest composition may benefit urban air quality. This study illustrates the importance of monitoring and representing changes in forest composition and the impacts to human health indirectly through changes in BVOCs.

54 ENVIRONMENTAL SCIENCES↗

Bridging adsorption behavior of confined CH 4 -CO 2 binary mixtures across scales

An accurate understanding of the competitive adsorption of CH 4 -CO 2 binary mixtures in nano-confined systems is critical for engineering CO 2 storage in shale gas reservoirs. Due to difficulties in making reliable experimental observations in nano-scale, atomistic simulations (ASs), such as the Grand Canonical Monte Carlo (GCMC) method, provide a viable approach to studying the adsorption behavior of confined fluids. ASs are, however, limited in the size of the compositional domain due to the high computational cost. This work proposes a framework that combines AS and the lattice Boltzmann (LB) method to bridge the physics of confined fluids across scales. The Peng–Robinson equation of state (PR-EoS) produces fugacity coefficients, which serve as input for conducting multi-component GCMC simulations. These GCMC simulations explore the competitive adsorption behavior of CH 4 -CO 2 in nano-slits at various composition, pressure, and channel-width conditions. Both components generate adsorption layers with high densities near the walls with CO 2 preferentially adsorbing compared to CH 4 on the organic walls of carbon sheets. At the mesoscale, a pseudopotential model represents the intermolecular forces in multi-component, multiple-relaxation-time LB simulations. The LB simulations are in good agreement with the GCMC results, allowing us to obtain values for tunable LB parameters. We then extend the use of LB to simulate adsorption behavior in complex networks with nano-sized channels. The phase behavior and fluid properties in the complex geometries of nano-channels differ from nano-slits and bulk systems. Furthermore, the bridging of physics from GCMC (microscale) to LB (mesoscale) via the macroscale PR-EoS connects the adsorption behavior of binary systems across scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation of macromolecule self-assembly in solution: A multiscale approach

One of the most common processes to produce polymer nanoparticles is to induce self-assembly by using the solvent-displacement method, in which the polymer is dissolved in a “good” solvent and the solution is then mixed with an “anti-solvent”. The polymer ability to self-assemble in solution is therefore determined by its structural and transport properties in solutions of the pure solvents and at the intermediate compositions. In this work, we focus on poly-ε-caprolactone (PCL) which is a biocompatible polymer that finds widespread application in the pharmaceutical and biomedical fields, performing simulation at three different scales using three different computational tools: full atomistic molecular dynamics (MD), population balance modeling (PBM) and computational fluid dynamics (CFD). Simulations consider PCL chains of different molecular weight in solution of pure acetone (good solvent), of pure water (anti-solvent) and their mixtures, and mixing at different rates and initial concentrations in a confined impinging jets mixer (CIJM). Our MD simulations reveal that the nano-structuring of one of the solvents in the mixture leads to an unexpected identical polymer structure irrespectively of the concentration of the two solvents. In particular, although in pure solvents the behavior of the polymer is, as expected, very different, at intermediate compositions, the PCL chain shows properties very similar to those found in pure acetone as a result of the clustering of the acetone molecules in the vicinity of the polymer chain. We derive an analytical expression to predict the polymer structural properties in solution at different solvent compositions and use it to formulate an aggregation kernel to describe the self-assembly in the CIJM via PBM and CFD. Simulations are eventually validated against experiments.

36 MATERIALS SCIENCE↗