Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D radiative transfer”

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 109 records · Page 6

RELAP5-3D Capabilities for High-Temperature Gas-Cooled Reactor Analysis

RELAP5-3D is the premier systems thermal-hydraulic code in the world. The code has been used significantly for high-temperature gas-cooled reactor (HTGR) analyses over the last 10 years, but there has not been a comprehensive summary of RELAP5-3D capabilities and methods for HTGR modeling. In this paper, we provide a guide for how novice users can build models for prismatic HTGRs. This guide includes best practices for developing heat structures and setting up radiation and conduction heat transfer. We also provide a summary of relevant validation studies. These studies are divided into three categories: separate effects tests, integral effects tests, and nuclear reactor studies. We present a comprehensive review of RELAP5-3D HTGR applications for reactor modeling over the last 10 years. These applications include work for the Modular High-Temperature Gas-Cooled Reactor benchmark, design activities for a prismatic HTGR as part of the U.S. Department of Energy’s Advanced Demonstration and Test Reactor Options Study, and studies performed to design gas-cooled microreactors. Following the comprehensive literature review, we present a gap analysis highlighting capabilities of RELAP5-3D that are available and missing for a complete analysis of HTGR thermal hydraulics. We also discuss the limitations of existing code capabilities and potential improvements that could improve analyses. Current research efforts have yet to fully validate RELAP5-3D for a prismatic HTGR analysis; however, some of this is attributable to confounding factors in the available experimental data as opposed to code capabilities. In conclusion, this paper presents the most comprehensive review to date on RELAP5-3D capabilities, limitations, and applications for HTGR analysis.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis of Postulated Accident Scenarios for the Transformational Challenge Reactor

The Transformational Challenge Reactor (TCR) is a helium-cooled, yttrium-hydride-moderated reactor that was designed for the U.S. Department of Energy Office of Nuclear Energy. A key objective of the TCR was to employ advanced manufacturing techniques in a nuclear system and demonstrate their potential for revolutionizing the nuclear reactor design process. One purpose of the present work is to demonstrate the safety of the TCR under postulated accidents. Based on RELAP5-3D and COMSOL analyses, the TCR remained below all current safety limits and far below the expected failure limits for the core materials. Another purpose of this work is to provide useful insights and recommendations regarding the application of RELAP5-3D to gas-cooled or other advanced reactors. A novel approach was implemented for simultaneously modeling conduction and radiation in RELAP5-3D, which was found to provide reasonable predictions of radial core, vessel, and ex-vessel heat transfer during postulated events. A multicode approach was also applied, in which high-fidelity COMSOL calculations were used to tune the radial heat transfer parameters in RELAP5-3D. The tuned RELAP5-3D model demonstrated comparable peak temperature predictions as COMSOL, despite a coarse treatment of the core in RELAP5-3D consisting of only two lumped heat structures. This high-fidelity tuning approach enabled enhanced accuracy as well as minimal complexity within the RELAP5-3D model, even for complex fuel geometric designs as in the TCR. Finally, investigations were made into the potential for flow reversal during a pressurized loss-of-forced-flow event in the TCR. The TCR is designed with downward helium flow through the core during normal operation. The RELAP5-3D model predicted that this downward flow would persist, without flow reversal, up to several days after the circulator trip. This was attributed to natural circulation hysteresis effects as have been noted in similar thermofluidic systems. Although flow stagnation and eventual reversal did not lead to unsafe TCR conditions, interesting spatial effects were observed which may have safety relevance for other reactor system designs and coolant types that are designed for downward core flow during normal operation, warranting closer investigation of the flow reversal phenomenon.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CubeSat Active Thermal Control in Support of Advanced Payloads: The Active Thermal Architecture Project

The Active Thermal Architecture (ATA) is an advanced sub-1U Active Thermal Control technology (ATC) for high power payload support in 6U CubeSat form factors and above. The design utilizes a two-stage, single-phase mechanically pumped fluid loop coupled through a two-axis flexible rotary fluid hinge, to reject thermal power to a deployable tracking radiator. A COTS Ricor K508N cryocooler forms the second stage and provides cryogenic cooling to a custom Kevlar detector mount through a TMT pyrolytic graphene thermal strap. Passive vibration isolation and damping technologies prevent the transfer of jitter to the satellite systems. The ATA design utilizes state-of-the-art 3D fabrication techniques such as Ultrasonic Additive Manufacturing (UAM) to directly embed the working fluid channels into the HX, radiator, and CubeSat chassis allowing for the miniaturization and simplification of the ATA system into an integrated thermal control solution. This paper will focus on the design and ground-based characterization and qualification of the ATA system and provide performance metrics for its use as a thermal support subsystem for advanced infrared electro-optical CubeSat payloads. The ATA project is funded through a NASA Small Satellite Technology Program (SSTP) and is a partnership between the Center for Space Engineering at Utah State University and the Jet Propulsions Laboratory. The ATA active thermal control system has been raised to a TRL of 6 and hopes to provide payload support to advanced missions such as the SABER-Lite and JPL CIRAS projects.

Mok, Mason↗

Characterization of radiation damage in 3D printed SiC

The SiC fuel matrix for advanced gas-cooled high temperature reactors as part of the Transformational Challenge Reactor program serves as the fuel particle container structure, a barrier to fission gas release, and a heat transfer medium. Its performance is particularly important because the fuel matrix must demonstrate good structural stability and thermal behaviors. An additive manufacturing methodology combining a binder jet 3D printing process with chemical vapor infiltration (CVI) for the production of SiC was recently developed. In this study, post irradiation examination by transmission electron microscopy shows that defect accumulation within the printed particles is very similar to other forms of high-purity SiC. However, damage accumulation was not directly observed in the CVI matrix because black spot damage and dislocation loops are difficult to image within the nanoscale highly faulted CVI matrix and because interstitial defects may rapidly annihilate at the stacking faults. Therefore, electron energy loss spectroscopy (EELS) analysis was used to analyze defect swelling in both the printed particles and the CVI matrix. The EELS analysis helped reveal that the radiation-induced swelling in the CVI matrix is similar to that of the printed SiC particles. This work shows that 3D printed SiC has behavior that is comparable to SiC processed by other means and that 3D printing could serve as a suitable processing technique for high-purity SiC for nuclear applications.

36 MATERIALS SCIENCE↗

Modeling cross-beam energy transfer with sector ray tracing

Ray-based cross-beam energy transfer (CBET) models are an essential feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion implosions, but full 3D ray-based CBET calculations can have a prohibitively high computational cost. Sector ray tracing can be used to reduce the cost by orders of magnitude in cases where the coronal plasma and laser drive can be approximated as spherically symmetric. An extension of sector ray tracing (section ray tracing) can be used to relax the assumption of a spherically symmetric drive while still retaining most of the computational savings of a sector ray trace. We discuss the foundations of sector and section ray tracing and compare them to full ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

3D simulations of inertial confinement fusion implosions part 2: systematic flow anomalies and impact of low modes on performances in OMEGA experiments

We present the first 3D radiation-hydrodynamic simulations of directly driven inertial confinement fusion implosions with an inline package for polarized crossed-beam energy transfer, which were used to assess the impact of the current distributed polarization rotators (DPRs) on OMEGA as well as other known sources of asymmetry. Applied to OMEGA implosions, the simulations predict bang times with no need for ad hoc multipliers, as well as yields—if you separately account for the impacts of imprint and fuel age. The magnitude of the flow is well reproduced when the low mode sources are large, whereas the modeling of the stalk is thought to be required to match the flow magnitude in the remaining cases. For the cases explored in more detail, polarized cross-beam energy transfer (CBET)—the only known systematic drive asymmetry, brought the results closest to the measured flow vectors. The remaining discrepancies are shown to be stemming from the limited knowledge of the laser pointing modes. For typical current levels of beam mispointing, power imbalance, target offset, and asymmetry caused by polarized CBET, low modes degrade the yield by more than 40%. The current strategy of attempting to compensate the mode-1 asymmetry with a preimposed target offset recovers only about one-third of the losses caused by the low modes due to the dynamic nature of the multiple asymmetries and the presence of low modes other than l = 1. Furthermore, addressing the root causes of the drive asymmetries is apt to be more beneficial. To that end, one possible solution to the specific issue of polarized CBET (10 µm DPRs) is shown to work well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Space-Based Lidar Observations of the 3D Structure of the Earth System

Lidar provides precise measurements of the three-dimensional structure of the clouds, aerosols, ocean/land/snow/ice surfaces, as well as ocean subsurface. Lidar also provides unique information about physical propertiesof particulates in the atmosphere for both radiative transfer and air quality applications.In this talk, I will present an overview of our recent studies of aerosols, clouds, ocean and snow using space-based lidar measurements (e.g., LITE, CALIPSO and ICESat-2), such as classifications of aerosols and thermodynamics phase of clouds, cloud microphysical properties, snow depths and phytoplankton biomass. I will also introduce a new concept of 3D Earth system observations with data fusion though combined active/passive remote sensing and machine learning. The new concept aims to reveal vertical structure of the aerosols/clouds/surfaces/subsurface from passive sensors by taking advantage of lidar measurements to effectively resolve the vertical structureby unscrambling the highly convoluted multi-angle, spectral and polarization information from passive sensors and apply the knowledge to a large swath where lidar measurements are not available.

Ali Omar↗

Assessing impacts of canopy 3D structure on chlorophyll fluorescence radiance and radiative budget of deciduous forest stands using DART

Although remote sensing (RS) of solar-induced chlorophyll fluorescence (SIF) is increasingly used as a valuable source of information about vegetation photosynthetic activity, the RS SIF observations are significantly influenced by canopy-specific structural features (i.e., canopy architecture including leaf area index and presence of woody components), atmospheric conditions during their acquisition (e.g., proportion of direct and diffuse irradiance) and observational geometric configurations (e.g., sun and viewing directions). Radiative transfer (RT) models have the potential to provide a better understanding of the canopy structural effects on the SIF emission and RS signals. Here, we used the DART model to assess the daily influence, from morning to evening, of forest 3D architecture on SIF nadir radiance, emission, escape factor and nadir yield of eight 100 m × 100 m forest study plots established in a temperate deciduous forest of the Smithsonian Environmental Research Center (Edgewater, MD, USA). The 3D architecture of each plot was derived from airborne LiDAR. DART simulations of these 3D forest plots and their 1D (i.e., vertical profile of sun-adapted and shade-adapted leaves) and 0D (i.e., homogeneous layer of sun-adapted leaves above an homogeneous layer of shade-adapted leaves) abstractions were compared to assess the relative errors (ε1D−3D and ε0D−3D) associated with horizontal and vertical structural heterogeneity, respectively. Forest 3D structure, especially horizontal heterogeneity, had a great influence on forest nadir SIF radiance, resulting in ε1D−3D up to 55% at 8:00 and 18:00 (i.e., for oblique sun directions). The key indicators of this impact, in the descending order of importance, were the SIF escape factor (ε1D−3D up to 40%), the attenuation of incident photosynthetically active radiation (ε1D−3D less than 5%), and the SIF emission yield (ε1D−3D less than 2%). The influence of forest architecture on the nadir SIF escape factor and SIF yield (ε1D−3D up to 40%) varied over time, with differences in forest stand structure, and per spectral domain, being always larger between 640 and 700 nm than between 700 and 850 nm. In addition, woody elements demonstrated a large influence on forest SIF radiance due to their “shading” effect (ε up to 17%) and their “blocking” effect (ε ≈ 10%), both of them higher for far-red than for red SIF. These results underline the importance of 3D forest canopy architecture, especially 2D heterogeneity, and inclusion of woody elements in RT modeling used for interpretation of the RS SIF signal, and subsequently for the estimation of gross primary production and detection of vegetation stress.

Omar Regaieg↗

Assessing the Impact of Effective Thermal Conductivity on Gas-cooled Reactor Transients in RELAP5-3D

High-temperature gas-cooled reactors (HTGRs) are a relatively mature advanced reactor concept that are of interest to industry and government for near-term deployment. These systems feature passive safety through low power density, coated particle fuels, and large graphite volumes that lead to slow heatup. In loss of flow transients, cooldown is achieved through conduction and radiation heat transfer. In prismatic block-type reactors, such as the mHTGR-350, the presence of coolant holes and fuel compacts alters the flow of heat through the blocks. To capture this effect in systems codes like RELAP5-3D, relationships for effective thermal conductivity (ETC) must be used to appropriately degrade the thermal conductivity. This work presents an assessment of the impact of ETC on performance of block-type gas-cooled reactors in a pressurized conduction cooldown (PCC) and depressurized conduction cooldown (DCC) transients. Models for PCC and DCC with bulk material thermal conductivity and ETC were run for both the mHTGR-350 and the High-Temperature Test Facility (HTTF) to determine the impact of ETC on transient performance.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

General relativistic MHD simulations of non-thermal flaring in Sagittarius A*

Sgr A* exhibits regular variability in its multiwavelength emission, including daily X-ray flares and roughly continuous near-infrared (NIR) flickering. The origin of this variability is still ambiguous since both inverse Compton and synchrotron emission are possible radiative mechanisms. The underlying particle distributions are also not well constrained, particularly the non-thermal contribution. In this work, we employ the GPU-accelerated general relativistic magnetohydrodynamics code H-AMR to perform a study of flare flux distributions, including the effect of particle acceleration for the first time in high-resolution 3D simulations of Sgr A*. For the particle acceleration, we use the general relativistic ray-tracing code bhoss to perform the radiative transfer, assuming a hybrid thermal+non-thermal electron energy distribution. We extract ~60 h light curves in the sub-millimetre, NIR and X-ray wavebands, and compare the power spectra and the cumulative flux distributions of the light curves to statistical descriptions for Sgr A* flares. Our results indicate that non-thermal populations of electrons arising from turbulence-driven reconnection in weakly magnetized accretion flows lead to moderate NIR and X-ray flares and reasonably describe the X-ray flux distribution while fulfilling multiwavelength flux constraints. These models exhibit high rms per cent amplitudes, $\gtrsim 150{{\ \rm per\ cent}}$ both in the NIR and the X-rays, with changes in the accretion rate driving the 230 GHz flux variability, in agreement with Sgr A* observations.

79 ASTRONOMY AND ASTROPHYSICS↗

The Use of Vegetation For Estimating Broken-Cloud Optical Properties From Surface Measurements

It is known that because of complex three-dimensional (3D) radiative effects of broken clouds, the retrieval of cloud optical properties from upward measurements based on a one-dimensional (1D) inversion technique almost surely fails. To remove radiative effects of 3D cloud structure, we have developed a new technique that retrieves cloud optical thickness for broken clouds above green vegetation from simultaneous surface measurements in the VIS and Near Infrared (NIR) spectral regions. The theoretical basis of the method is the very different spectral behavior of cloud liquid water drops and green vegetation. For example, cloud optical properties, and hence cloud reflectivities, change little between 650 and 860 nm, while the vegetated surface albedo changes from 0.05 to 0.5 between the same two wavelengths. This spectral contrast in surface albedo suggests using ground measurements at both wavelengths not independently, but as an algebraic combination (a spectral index). For a spectral band in the NIR region, the green vegetation acts as a powerful reflector that "illuminates" horizontally inhomogeneous clouds from below. This provides the extra information needed to largely remove the 3D radiative effects, especially in the case of broken clouds; this in turn allows the retrieval of cloud optical depth using traditional 1D radiative transfer theory. This approach is similar to the so-called Green's function problem for radiative transfer where a laser beam illuminates clouds and the resulting "spot-size" of the reflected light around the beam characterizes cloud properties. We generalize Green's function theory to surf ace-cloud interaction and develop new spectral indices from which broken-cloud optical depth can be retrieved.

Marshak, A.↗

TLSLEAF: Automatic Leaf Angle Estimates From Single-Scan Terrestrial Laser Scanning

Leaf angle distribution (LAD) in forest canopies affects estimates of leaf area, light interception, and global-scale photosynthesis, but is often simplified to a single theoretical value. Here, we present TLSLeAF (Terrestrial Laser Scanning Leaf Angle Function), an automated open-source method of deriving LADs from terrestrial laser scanning. TLSLeAF produces canopy-scale leaf angle and LADs by relying on gridded laser scanning data. The approach increases processing speed, improves angle estimates, and requires minimal user input. Key features are automation, leaf–wood classification, beta parameter output, and implementation in R to increase accessibility for the ecology community. TLSLeAF precisely estimates leaf angle with minimal distance effects on angular estimates while rapidly producing LADs on a consumer-grade machine. We challenge the popular spherical LAD assumption, showing sensitivity to ecosystem type in plant area index and foliage profile estimates that translate to c. 25% and c. 11% increases in canopy net photosynthesis (c. 25%) and solar-induced chlorophyll fluorescence (c. 11%). TLSLeAF can now be applied to the vast catalog of laser scanning data already available from ecosystems around the globe. The ease of use will enable widespread adoption of the method outside of remote-sensing experts, allowing greater accessibility for addressing ecological hypotheses and large-scale ecosystem modeling efforts.

3D↗

NRC Multiphysics Analysis Capability Deployment FY21: Part 3

This report details the progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” The deliverables completed for this report are: Deliverable 1c: the capability to model gas mixtures was added to Pronghorn. A test problem mimicking the conditions achieved in a depressurized loss of forced cooling (DLOFC) event was solved with both RELAP-5 and Pronghorn. Pronghorn employed a finite vol ume method with the Kurganov-Tadmor discretization. The comparison between the mass fraction spatial profiles computed with RELAP-5 and Pronghorn clearly shows the presence of numerical artifacts (i.e., overly diffusive behavior at low Mach numbers). We confirmed that the problem disappears at higher Mach numbers. We recommend future work on the implementation of a low Mach finite volume formulation to better treat low Mach number problems. Deliverable 2a: we demonstrated two approaches to model the radiation/conduction/natural convection heat transfer across a stagnant gas for the PBMR-400 design using Pronghorn. The first approach is based on the net radiation method, which relies on the computation of view factors with the Multiphysics Object-Oriented Simulation Environment (MOOSE) ray tracing capability. The second method is a traditional thermal resistance approach. The test problems include both 2D and 3D geometries. In all cases, the results show very good agreement during a DLOFC transient. This confirms that the faster thermal resistance method produces solutions that are equivalent to the net radiation method for this geometry. Deliverable 3d: we demonstrated the use of the advection kernel for the delayed neutron precursor equation in Griffin with a 2D MSFR model. The results appear physical but further verification is recommended. We also recommend the addition of conjugate heat transfer to compute the temperatures and model the thermomechanic behavior of the reflectors and other structures. Significant memory and performance issues were encountered in the 3D axisymmetric model. Future work is recommended in this area. Task 8g: this task allows multidimensional MOOSE applications to be coupled to system codes (RELAP-7 and SAM). We implemented a faster multiphysics iteration coupling algorithm, which provides an overall 6× acceleration of the 3D-1D coupling of the core multidi- mensional fluid flow solver and the 1D primary and secondary loop model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Modeling Radiation-Induced Epithelial Cell Injury in Murine Three-Dimensional Esophageal Organoids

Esophageal squamous cell carcinoma (ESCC) is a deadly consequence of radiation exposure to the esophagus. ESCC arises from esophageal epithelial cells that undergo malignant transformation and features a perturbed squamous cell differentiation program. Understanding the dose- and radiation quality-dependence of the esophageal epithelium response to radiation may provide insights into the ability of radiation to promote ESCC. We have explored factors that may play a role in esophageal epithelial radiosensitivity and their potential relationship to ESCC risk. We have utilized a murine three-dimensional (3D) organoid model that recapitulates the morphology and functions of the stratified squamous epithelium of the esophagus to study persistent dose- and radiation quality-dependent changes. Interestingly, although high-linear energy transfer (LET) Fe ion exposure induced a more intense and persistent alteration of squamous differentiation and 53BP1 DNA damage foci levels as compared to Cs, the MAPK/SAPK stress pathway signaling showed similar altered levels for most phospho-proteins with both radiation qualities. In addition, the lower dose of high-LET exposure also revealed nearly the same degree of morphological changes, even though only ~36% of the cells were predicted to be hit at the lower 0.1 Gy dose, suggesting that a bystander effect may be induced. Although p38 and ERK/MAPK revealed the highest levels following high-LET exposure, the findings reveal that even a low dose (0.1 Gy) of both radiation qualities can elicit a persistent stress signaling response that may critically impact the differentiation gradient of the esophageal epithelium, providing novel insights into the pathogenesis of radiation-induced esophageal injury and early stage esophageal carcinogenesis.

3D organoid culture↗

The extended atmosphere and circumstellar environment of the cool evolved star VX Sagittarii as seen by MATISSE

Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3–4 μm), M (4.5–5 μm), and N (8–13 μm) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes OPTIM3D and RADMC3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N-band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of VX Sgr’s stellar surface and circumstellar environment across a very large spectral domain for the first time.

stars: atmospheres↗

Evidence for Three-Dimensional Radiative Effects in MODIS Cloud Optical Depths Retrieved at Back Scattering View Angles

This study addresses the question whether 1D radiative transfer theory describes well the angular distribution of shortwave cloud reflection. The statistical analysis of a large set of MODIS observations indicates that in oblique backward scattering directions, cloud reflection is stronger than 1D theory would predict. After considering a variety of possible causes, the paper concludes that the most likely reason for the increase lies in 3D radiative interactions. The results' main implication is that cloud optical depths retrieved at back scattering view angles larger than about 50 degrees tend to be overestimated and should be used only with great caution.

Varnai, Tamas↗

Simulation-Aided Design and Interpretation of Results of Microgravity-Based Solidification Experiments

The preparation and interpretation of microgravity experiments conducted on the International Space Station (ISS) can be challenging because of the remote location of the equipment. Although a “ground unit” is typically available on Earth, it is not located in the laboratory of the principal investigator. Furthermore, test experiments in the ground unit cannot fully duplicate heat transfer conditions of the weightless environment, and thus be a reliable reference for experiment design. Under such conditions, numerical simulations are a valuable tool both for preparing the experiments and interpreting the results obtained under terrestrial and micro-gravity conditions. In this work, we developed a detailed 3D model of the “SUBSA” furnace located at the ISS, and demonstrate its use in our present investigation which includes melt growth and vapor growth of InI crystals on the ISS. In 2002, the SUBSA furnace was designed and built for the growth of InSb crystals in the Microgravity Science Glovebox (MSG) at the ISS. The SUBSA furnace features a transparent section without heating elements, where the axial temperature gradient is high. In our present investigation, we have used the same furnace to grow four InI crystals from the melt, and two InI crystals from the vapor phase (2018-2021). InI crystals have ~20 times lower thermal conductivity than InSb crystals. As a result, during melt growth, the InI crystals were exposed to extremely high temperature gradients. To reduce the axial temperature gradient, we added external jackets made from high thermal conductivity material to the growth ampoules. The SUBSA furnace could not be modeled as asymmetric because of its rectangular shape and right-angled transparent section, which contains a convoluted geometry. The rectangular shape of the transparent section further complicates modelling of radiative heat transfer, which dominates in the zone. At present, numerical simulations are being conducted to resolve an apparent discrepancy between the readings of the 5 thermocouples located in the SUBSA furnace, and the observed position of the solid and molten material. The 3D simulation will be presented, along with the calibration and growth experiment results conducted in microgravity and in the SUBSA ground unit. The developed model should be useful for upcoming microgravity investigations planned to be conducted in the SUBSA facility.

V Riabov↗