Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state process”

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,189 records · Page 66

Development of Next Generation Hierarchical Hybrid Cu-Si anode Batteries via Direct-Ink Writing Application: End of (6th) Month Report - November 2025

Sustainable renewable energy continues to be in dire need to effectively combat global warming. Emerging technology for electric vehicles/ devices remains in high demand that is not only lower in cost, more efficient, but safer in comparison to commercial materials on the market. Although first-generation lithium-ion batteries have exhibited extensive commercial application, conventional graphite no longer meets this increasing demand as an efficient anode material. Due to the fact that graphite has a subpar theoretical specific capacity (372 mAh g -1 ), thus significant limitations in rate capability (for potential faster charging at higher C-rates currently commercially available.). Alternatively, silicon has gained significant attention as a superior candidate to potentially surpass graphite. Due to silicon’s exceedingly high theoretical capacity (4,200 mAh g -1 ) in comparison to standard graphite, its abundance thus in turn it’s low-cost, in addition to exhibiting a significantly low working potential (< 0.4 V vs Li/Li + ). However, one of the main (and most detrimental) challenges is silicon’s tendency to expand in volume (> 300%) upon discharge as it begins the lithiation process. As a direct result, it causes not only for the particles to both crack and pulverize under mechanical stress as the volume continues to expand and contract during cycling. Upon assembling the cell, it needs to undergo ‘charging’ for initially discharging/ ‘activating’ the cell, otherwise commonly known as the ‘formation’ step. As a result a solid electrolyte interface (SEI) layer begins to form at the anode surface because some of the electrolyte begins to react during the formation process. However, this (SEI) layer is deemed as a ‘protective’ interlayer because in theory it prevents further reaction as the cell continues to cycle. However, due to the volume change it causes significant degradation at the interface. As cracking starts to occur at the anode surface, it results in a ‘new’ altered surface with each cycle that causes further reaction with the electrolyte as a byproduct quickly consuming active Li/ and more electrolyte. Thus, fracturing this ‘protective layer,’ causing significantly higher impedance as a result, and in turn a decline in capacity due to active Li-loss. In addition to the active material exfoliating off from the current collector, further contributing to the steady decline in capacity and overall performance. Current state of the art Si-anode batteries on the market range between a maximum content of 5-10 Si wt%. It has been previously reported Tesla has utilized SiO x -C anodes containing 5 wt% Si within their ‘Model 3/ Model X’ electric vehicles. However, more recent ‘Model 3’ vehicles have started to incorporate 10 Si wt%, in which they were able to increase their energy density upwards by approximately 30%. Recent effort has been focused on continuing to increase the wt% of Si being utilized, eventually to 100 wt% of Si, to maximize the energy density even further.

36 MATERIALS SCIENCE↗

Generalized Fluid System Simulation Program (GFSSP) - Version 6

The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, flow control valves and external body forces such as gravity and centrifugal. The thermo-fluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids, and 24 different resistance/source options are provided for modeling momentum sources or sinks in the branches. Users can introduce new physics, non-linear and time-dependent boundary conditions through user-subroutine.

Majumdar, Alok↗

Case Study: Calorimetry

This chapter presents the use of adsorption microcalorimetry to study acid-base catalysis over mixed metal oxides. Understanding the acid-base properties of catalysts is of paramount importance for steering product selectivity in a wide variety of catalytic industrial processes. In particular, catalytic conversion of biomass into biofuels and chemicals demands comprehensive understanding of the acid-base properties of the catalyst to selectively convert functional groups. Herein, the use of probe molecules to measure the strength and density of acid and basic sites on mixed metal oxide surfaces is introduced. Further, the challenges when characterizing mixed metal oxides, presenting surface reconstruction, are explained. An example displays how combining adsorption calorimetry with ab initio calculations explains the activation of a reactant on the surface of a mixed metal oxide. Acid-base catalysis over nanoshapes is studied to unveil that different nanoshapes lead to different extents of surface reconstruction, which ultimately affect acid-base catalysis. The use of steady-state isotopic transient kinetic analysis (SSITKA), in combination with adsorption microcalorimetry, unveils that surface reconstruction of a set of perovskites impacts the density of active sites to a greater extent than the intrinsic reactivity of the active site. At the end of this chapter, fundamental insights into catalyst transient behavior during isomerization of n-butane, and the reaction mechanism for conversion of 4-methylpentan-2-ol, are presented.

Polo Garzon, Felipe↗

Non-Hermitian approach for quantum plasmonics

Herein, we examine the limits of applicability of a simple non-Hermitian model for exciton/plasmon interactions in the presence of dissipation and dephasing. The model can be used as an alternative to the more complete Lindblad density matrix approach and is computationally and conceptually simpler. We find that optical spectra in the linear regime can be adequately described by this approach. The model can fail, however, under continuous optical driving in some circumstances. In the case of two quantum dots or excitons interacting with a plasmon, the model can also describe coherences and entanglement qualitatively when both dissipation and dephasing are present and quantitatively in the limit with no dephasing. The approach, within a single excitation manifold, is also applied to assess the role of disorder for 50 quantum dots interacting with a plasmon, where we find that, on average, large enough disorder can help stabilize the ensemble average of the open quantum system toward a dark quasi-steady-state much faster than without disorder. Furthermore, while such single excitation manifold calculations in this size limit can readily be done with either the non-Hermitian or Lindblad forms, as one goes to larger Hilbert space sizes, the computational and storage advantages of the non-Hermitian approach can become more useful. Published under license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electrical monitoring of saline tracers to reveal subsurface flow pathways in a flat ditch-drained field

Time-lapse electrical resistivity imaging (ERI) combined with salt tracers can enhance our understanding of subsurface flow and transport in ditch-drained agroecosystems. In a nearly level field adjacent to a drainage ditch, we installed 192 electrodes in a 72-m2 plot along with five hydrometric stations to enable time-lapse ERI monitoring during storms. Immediately prior to a wintertime rainstorm, we applied 5 kg of KBr to a shallow trench located roughly 11 m upgradient from the ditch. Antecedent conditions prior to the storm were wet, assuring a hydrologically connected state that favored tracer dissolution and movement in the subsurface. We monitored the evolution of the tracer plume with ERI over a seven-day period. Four-dimensional (4D) delineation of precipitation and tracer transport using time-lapse ERI showed that the tracer plume moved preferentially downward through a discontinuity in the argillic horizon within the first 57 hrs. Over the following 79 hrs, the tracer plume exhibited steady, lateral movement towards the drainage ditch within a more permeable sandy horizon situated just below the argillic. A spatial moment analysis of logarithmic changes in conductivity in the sand layer revealed seepage velocities ranging from 0.0027 to 0.0482 m d-1; these values were positively correlated with hydraulic gradients, suggesting laminar flow processes typical of sandy aquifers. Increasing groundwater velocities with hydraulic gradients suggest that solute transfers from groundwater to ditch are likely to be faster in storms, especially during wet periods. Even so, predominantly slow lateral transport within the sandy horizon, even during storms, indicates that critical source areas of pollutant loss are likely within a few meters of the ditch.

Robinson, Judith L.↗

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗

Ducted Assembly Steady-State Heat Transfer Software (DASSH): Theory Manual

The Ducted Assembly Steady-State Heat transfer software (DASSH) is being developed at Argonne National Laboratory to obtain a steady-state thermal fluids result for the coolant flow and temperature distribution in a hexagonally gridded reactor core with ducted assemblies. In DASSH, each assembly is assumed to contain a hexagonal bundle of wire-wrapped pins. DASSH divides the coolant in the pin bundle into subchannels on which it applies energy balances to determine coolant and duct temperatures. DASSH offers multiple models for inter assembly heat transfer and features the ability to calculate pin temperatures. This document describes the models and equations implemented in DASSH. It includes a description of the geometric system to which the models are applied and a derivation of the subchannel energy balance formulation is presented. It also reviews the available correlations and provides an overview of the process by which power distributions are generated. Finally, it covers the code execution procedure, discussing the main code objects and their interactions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Supraglacial River Forcing of Subglacial Water Storage and Diurnal Ice Sheet Motion

Surface melting impacts ice sheet sliding by supplying water to the bed, but subglacial processes driving ice accelerations are complex. We examine linkages between surface runoff, transient subglacial water storage, and short-term ice motion from 168 consecutive hourly measurements of meltwater discharge (moulin input) and GPS-derived ice surface motion for Rio Behar, a ∼60 sq. km moulin-terminating supraglacial river catchment on the southwest Greenland Ice Sheet. Short-term accelerations in ice speed correlate strongly with lag-corrected measures of supraglacial river discharge (r = 0.9, τ = 0.7, p < 0.01). Though our 7 days record cannot address seasonal-scale forcing, diurnal ice accelerations align with normalized differenced supraglacial and proglacial discharge, a proxy for subglacial storage change, better than GPS-derived ice surface uplift. These observations counter theoretical steady state basal sliding laws and suggest that moulin and proglacially induced fluctuations in subglacial water storage, rather than absolute subglacial water storage, drive short-term ice accelerations.

Glacial↗

Coupling of Pronghorn and RELAP-7 for a Pebble Bed Reactor

High temperature gas cooled reactors (HTGR) are a candidate for timely Gen-IV reactor technology deployment because of high technology readiness and walk-away safety. Among HTGRs, pebble bed reactors (PBRs) have attractive features such as low excess reactivity and online refueling. Pebble bed reactors pose unique challenges to analysts and reactor designers such as continuous burnup distribution depending on pebble motion and recirculation, radiative heat transfer across a variety of gas-filled gaps, and long design basis transients such as pressurized and depressurized loss of forced circulation. Modeling and simulation is essential for both the PBR’s safety case and design process. In order to verify and validate the new generation codes the Nuclear Energy Agency (NEA) Data bank provide a set of benchmarks data together with solutions calculated by the participants using the state of the art codes of that time. An important milestone to test the new PBR simulation codes is the OECD NEA PBMR-400 benchmark which includes thermal hydraulic and neutron kinetic standalone exercises as well as coupled exercises and transients scenarios. In this work, the reactor multiphysics code MAMMOTH and the thermal hydraulics code Pronghorn, both developed by the Idaho National Laboratory (INL) within the multiphysics object-oriented simulation environment (MOOSE), have been used to solve Phase 1 exercises 1 and 2 of the PBMR-400 benchmark. The steady state results are in agreement with the other participants’ solutions demonstrating the adequacy of MAMMOTH and Pronghorn for simulating PBRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Validation of a detonation product equation of state for an insensitive high explosive via slab geometry expansion tests

Slab expansion (SLABEX) tests are conducted to validate a process for calibrating the detonation products equation of state (EOS) of a high explosive (HE). The SLABEX tests use rectangular slabs of PBX 9502, a polymer-bonded HE formulation consisting of 95 wt.% 1,3,5-triamino-2,4,6-trinitrobenzene bound with Kel F-800, a co-polymer of chlorotrifluoroethylene and vinylidene-fluoride. Three PBX 9502 slab thicknesses are examined, each confined symmetrically by two rectangular copper (Cu) plates approximately one-tenth the thickness of the HE slab. For the duration of each experiment, the detonation flow along the central axis of the PBX 9502 slab remains two-dimensional. The lateral flow velocity component of the outer surfaces of the expanding Cu plates is measured, along with the steady axial detonation speed along the central axis of the SLABEX. Hydrodynamic simulations of the Cu plate expansion in the SLABEX geometry, driven by the energy stored in the detonation products by the detonation combustion event, are conducted using a Jones–Wilkins–Lee EOS for the detonation products. This EOS form was recently parameterized for PBX 9502 in the cylinder expansion test geometry using a newly developed calibration technique [Voelkel et al., Combust. Flame 246, 112373 (2022)]. Good agreement between the experiment and prediction is found in each SLABEX test, demonstrating that the detonation product EOS calibration technique produces EOSs that are predictive when applied to other geometries.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

NASTRAN solutions of problems described by simultaneous parabolic differential equations

NASTRAN solution techniques are shown for a numerical analysis of a class of coupled vector flow processes described by simultaneous parabolic differential equations. To define one physical problem type where equations of this form arise, the differential equations describing the coupled transfers of heat and mass in mechanical equilibrium with negligible mass average velocity are presented and discussed. Also shown are the equations describing seepage when both electrokinetic and hydrodynamic forces occur. Based on a variational statement of the general problem type, the concepts of scalar transfer elements and parallel element systems are introduced. It is shown that adoptation of these concepts allows the direct use of NASTRAN's existing Laplace type elements for uncoupled flow (the heat transfer elements) for treating multicomponent coupled transfer. Sample problems are included which demonstrate the application of these techniques for both steady-state and transient problems.

Mason, J. B.↗

Validation of the SINDA/FLUINT code using several analytical solutions

The Systems Improved Numerical Differencing Analyzer and Fluid Integrator (SINDA/FLUINT) code has often been used to determine the transient and steady-state response of various thermal and fluid flow networks. While this code is an often used design and analysis tool, the validation of this program has been limited to a few simple studies. For the current study, the SINDA/FLUINT code was compared to four different analytical solutions. The thermal analyzer portion of the code (conduction and radiative heat transfer, SINDA portion) was first compared to two separate solutions. The first comparison examined a semi-infinite slab with a periodic surface temperature boundary condition. Next, a small, uniform temperature object (lumped capacitance) was allowed to radiate to a fixed temperature sink. The fluid portion of the code (FLUINT) was also compared to two different analytical solutions. The first study examined a tank filling process by an ideal gas in which there is both control volume work and heat transfer. The final comparison considered the flow in a pipe joining two infinite reservoirs of pressure. The results of all these studies showed that for the situations examined here, the SINDA/FLUINT code was able to match the results of the analytical solutions.

Keller, John R.↗

Automated Propulsion Data Screening demonstration system

A fully-instrumented firing of a propulsion system typically generates a very large quantity of data. In the case of the Space Shuttle Main Engine (SSME), data analysis from ground tests and flights is currently a labor-intensive process. Human experts spend a great deal of time examining the large volume of sensor data generated by each engine firing. These experts look for any anomalies in the data which might indicate engine conditions warranting further investigation. The contract effort was to develop a 'first-cut' screening system for application to SSME engine firings that would identify the relatively small volume of data which is unusual or anomalous in some way. With such a system, limited and expensive human resources could focus on this small volume of unusual data for thorough analysis. The overall project objective was to develop a fully operational Automated Propulsion Data Screening (APDS) system with the capability of detecting significant trends and anomalies in transient and steady-state data. However, the effort limited screening of transient data to ground test data for throttle-down cases typical of the 3-g acceleration, and for engine throttling required to reach the maximum dynamic pressure limits imposed on the Space Shuttle. This APDS is based on neural networks designed to detect anomalies in propulsion system data that are not part of the data used for neural network training. The delivered system allows engineers to build their own screening sets for application to completed or planned firings of the SSME. ERC developers also built some generic screening sets that NASA engineers could apply immediately to their data analysis efforts.

Hoyt, W. Andes↗

Modified Penning-Malmberg Trap for Storing Antiprotons

A modified Penning-Malmberg trap that could store a small cloud of antiprotons for a relatively long time (weeks) has been developed. This trap is intended for use in research on the feasibility of contemplated future matter/antimatter-annihilation systems as propulsion sources for spacecraft on long missions. This trap is also of interest in its own right as a means of storing and manipulating antiprotons for terrestrial scientific experimentation. The use of Penning-Malmberg traps to store antiprotons is not new. What is new here is the modified trap design, which utilizes state-of-the-art radiofrequency (RF) techniques, including ones that, heretofore, have been used in radio-communication applications but not in iontrap applications. A basic Penning-Malmberg trap includes an evacuated round tube that contains or is surrounded by three or more collinear tube electrodes. A steady axial magnetic field that reaches a maximum at the geometric center of the tube is applied by an external source, and DC bias voltages that give rise to an electrostatic potential that reaches a minimum at the center are applied to the electrodes. The combination of electric and magnetic fields confines the charged particles (ions or electrons) for which it was designed to a prolate spheroidal central region. However, geometric misalignments and the diffusive cooling process prevent the steady fields of a basic Penning- Malmberg trap from confining the particles indefinitely. In the modified Penning-Malmberg trap, the loss of antiprotons is reduced or eliminated by use of a "rotating-wall" RF stabilization scheme that also heats the antiproton cloud to minimize loss by matter/antimatter annihilation. The scheme involves the superposition of a quadrupole electric field that rotates about the cylindrical axis at a suitably chosen radio frequency. The modified Penning-Malmberg trap (see Figure 1) includes several collinear sets of electrodes inside a tubular vacuum chamber. Each set comprises either a single metal tube or else a tube that is segmented into four electrodes that subtend equal angles about the cylindrical axis. The output of an RF signal generator is fed through a 90 hybrid coupler and then through two baluns to generate four replicas of the signal at relative phase shifts of 0 , 90 , 180 , and 270 (see Figure 2). These signal replicas are fed through 6-dB directional couplers, then via coaxial cables to the vacuum chamber. The signal is then routed to a phase cancellation network, which filters out the drive signal with the difference representing the plasma interaction. Inside the vacuum chamber, twisted-pair wires feed the signals from the coaxial cables to the four electrodes of each segmented electrode tube.

Sims, William H.↗

Modeling of Laser Vaporization and Plume Chemistry in a Boron Nitride Nanotube Production Rig

Flow in a pressurized, vapor condensation (PVC) boron nitride nanotube (BNNT) production rig is modeled. A laser provides a thermal energy source to the tip of a boron ber bundle in a high pressure nitrogen chamber causing a plume of boron-rich gas to rise. The buoyancy driven flow is modeled as a mixture of thermally perfect gases (B, B2, N, N2, BN) in either thermochemical equilibrium or chemical nonequilibrium assuming steady-state melt and vaporization from a 1 mm radius spot at the axis of an axisymmetric chamber. The simulation is intended to define the macroscopic thermochemical environment from which boron-rich species, including nanotubes, condense out of the plume. Simulations indicate a high temperature environment (T > 4400K) for elevated pressures within 1 mm of the surface sufficient to dissociate molecular nitrogen and form BN at the base of the plume. Modifications to Program LAURA, a finite-volume based solver for hypersonic flows including coupled radiation and ablation, are described to enable this simulation. Simulations indicate that high pressure synthesis conditions enable formation of BN vapor in the plume that may serve to enhance formation of exceptionally long nanotubes in the PVC process.

Gnoffo, Peter A.↗

Review of European microgravity measurements

AA In a French/Russion cooperation, CNES developed a microgravity detection system for analyzing the Mir space station micro-g-environment for the first time. European efforts to characterize the microgravity (1/9) environment within a space laboratory began in the late seventies with the design of the First Spacelab Mission SL-1. Its Material Science Double Rack was the first payload element to carry its own tri-axial acceleration package. Even though incapable for any frequency analysis, the data provided a wealth of novel information for optimal experiment and hardware design and operations for missions to come. Theoretical investigations under ESA contract demonstrated the significance of the detailed knowledge of micro-g data for a thorough experiment analysis. They especially revealed the high sensitivity of numerous phenomena to low frequency acceleration. Accordingly, the payloads of the Spacelab missions D-1 and D-2 were furnished with state-of-the-art detection systems to ensure frequency analysis between 0.1 and 100 Hz. The Microgravity Measurement Assembly (MMA) of D-2 was a centralized system comprising fixed installed as well as mobile tri-axial packages showing real-time data processing and transmission to ground. ESA's free flyer EURECA carried a system for continuous measurement over the entire mission. All EURECA subsystems and experimental facilities had to meet tough requirements defining the upper acceleration limits. In a French/Russion cooperation, CNES developed a mi crogravity detection system for analyzing the Mir space station micro-g-environment for the first time. An approach to get access to low frequency acceleration between 0 and 0.02 Hz will be realized by QSAM (Quasi-steady Acceleration Measurement) on IML-2, complementary to the NASA system Spacelab Acceleration Measurement System SAMS. A second flight of QSAM is planned for the Russian free flyer FOTON.

Hamacher, Hans↗

Modeling the effect of water on mantle rheology

To study the thermal history of the Earth we use a parameterized model of mantle convection. This model includes a mathematical description of de- and regassing processes of water from the Earth's mantle. The rates of this processes are considered to be directly proportional to the seafloor spreading rate. The kinematic viscosity of the mantle depends on the temperature/pressure as well as on the volatile content. Dissolved volatiles such as water weaken the minerals by reducing their activation energy for solid state creep. Karato and Toriumi showed a power law dependence between creep rate and water fugacity derived from experimental results. Therefore, we use such flow parameters of diffusion creep in olivine under wet and dry conditions to calculate the mantle viscosity as a function of the water content. Because the creep rate is proportional to the concentration of water-related point deflects we assume that the water fugacity is proportional to the water weight fraction. An equation for the steady-state strain rate under wet conditions is established. To assess the unknown constant K in this equation, we use flow law parameters given by Karato and Wu as well as the results of McGovern and Schubert.

Bounama, CH.↗

Modeling of Melt Growth During Carbothermal Processing of Lunar Regolith

The carbothermal processing of lunar regolith has been proposed as a means to produce carbon monoxide and ultimately oxygen to support human exploration of the moon. In this process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone of lunar regolith and is converted to carbon and hydrogen. Carbon gets deposited on the surface of the melt, and mixes and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. Carbon monoxide is further processed in other reactors downstream to ultimately produce oxygen. The amount of oxygen produced crucially depends on the amount of regolith that is molten. In this paper we develop a model of the heat transfer in carbothermal processing. Regolith in a suitable container is heated by a heat flux at its surface such as by continuously shining a beam of solar energy or a laser on it. The regolith on the surface absorbs the energy and its temperature rises until it attains the melting point. The energy from the heat flux is then used for the latent heat necessary to change phase from solid to liquid, after which the temperature continues to rise. Thus a small melt pool appears under the heated zone shortly after the heat flux is turned on. As time progresses, the pool absorbs more heat and supplies the energy required to melt more of the regolith, and the size of the molten zone increases. Ultimately, a steady-state is achieved when the heat flux absorbed by the melt is balanced by radiative losses from the surface. In this paper, we model the melting and the growth of the melt zone with time in a bed of regolith when a portion of its surface is subjected to a constant heat flux. The heat flux is assumed to impinge on a circular area. Our model is based on an axisymmetric three-dimensional variation of the temperature field in the domain. Heat transfer occurs only by conduction, and effects of convective heat transport are assumed negligible. Radiative heat loss from the surface of the melt and the regolith to the surroundings is permitted. We perform numerical computations to determine the shape and the mass of the melt at steady state and its time evolution. We first neglect the volume change upon melting, and subsequently perform calculations including it. Predictions from our model are compared to test data to determine the effective thermal conductivities of the regolith and the melt that are compatible with the data

Balasubramaniam, R.↗