SEARCH · Search NASA

Results for “generalized spherical functions”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2021

A family of ionic supersalts with covalent-like directionality and unconventional multiferroicity

Ionic crystals composed of elemental ions such as NaCl are non-polar due to directionless ionic bonding interactions. Here, we show that these can develop polarity by changing their building blocks from elemental ions to superalkali and superhalogen cluster-ions, which mimic the chemistry of alkali and halogen atoms, respectively. Due to the non-spherical geometries of these cluster ions, corresponding supersalts form anisotropic polar structures with ionic bonding, yet covalent-like directionality, akin to sp 3 hybridized systems. Using density functional theory and extensive structure searches, we predict a series of stable ferroelectric/ferroelastic supersalts, PnH 4 MX 4 (Pn = N, P; M = B, Al, Fe; X = Cl, Br) composed of superalkali PnH 4 and superhalogen MX 4 ions. Unlike traditional ferroelectric/ferroelastic materials, the cluster-ion based supersalts possess ultra-low switching barrier and can endure large ion displacements and reversible strain. In particular, PH 4 FeBr 4 exhibits triferroic coupling of ferroelectricity, ferroelasticity, and antiferromagnetism with controllable spin directions via either ferroelastic or 90-degree ferroelectric switching.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Linear Velocity Problems Part 2: Thermodynamic Considerations

A compressible polytropic gas, in one dimensional planar, cylindrical, or spherical coordinate systems, is examined through four thermodynamic properties: mass density, pressure, specific internal energy (SIE), and entropy, under the assumptions that the fluid velocity is linearly proportional to the radial coordinate and that the Euler gas dynamics equations assert a self similar solution class. The behavior of the thermodynamics is almost entirely determined by an arbitrary function that results from the derivation of the density distribution. Through physical and other arguments, the arbitrary function can be plotted spatially with time profiles, providing a deeper understanding of the system and how the arbitrary function affects the behavior of the compressible gas.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Group structure selection with random forests

Choosing an appropriate group structure for multigroup transport is far from an exact science. For some applications, one blindly uses a group structure developed years ago by forgotten methods. Furthermore, one sometimes uses the same group structure for a variety of problems, even if the group structure was originally developed with a certain application in mind. In this work, we create optimized group structures with simulated annealing for critical assembly test problems and apply a random forest regressor with bagging to choose the best group structure based on parameters of the different test problems. The optimized group structures were generated using a simulated annealing optimizer for several simple, spherical, and unreflected problems. The optimization was performed to minimize a cost function that included fission rate, absorption rate, leakage, and k{sub eff}. A random forest regressor was then trained on a set of International Criticality Safety Benchmark Evaluation Project inputs and used to select one of these six group structures. The trained machine learning model chose the best group structure 65% of the time, and one of the three best 89% of the time. Furthermore, it decreased the L2 error over all test problems by a factor of 25 when compared to the standard Los Alamos 70-group structure. In other words, the model chose group structure that were far more appropriate for the test problems than the traditional LANL group structure. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Parametric raytracing modeling for NSTX-U scenario development with high harmonic fast waves and neutral beam injection

High harmonic fast waves (HHFW) are a versatile heating and current drive tool for scenario development. Extensive modeling scans were performed to find optimal parameters for different uses of HHFW in National Spherical Tokamak Experiment (NSTX-U). Scans of plasma density, temperature, magnetic field, and antenna phasing were performed both with and without neutral beam injection. For speed of calculation, the ray-tracing code GENRAY coupled to the quasilinear Fokker–Planck code CQL3D was used. CQL3D allows for a more accurate description of the fast ion population, as well as for quasilinear effects such as HHFW-induced modifications of the distribution function. Best current drive results are obtained at elevated electron temperatures and with the lowest k φ phasing. Adding neutral beams however typically strongly reduces the HHFW current drive efficiency at the low density cases due to HHFW absorption on beam ions. Results of this parametric study will feed into scenario development and predict-first whole-shot modeling of NSTX-U discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Semi-Analytic Functions to Calculate the Deposition Coefficients for Ice Crystal Vapor Growth in Bin and Bulk Microphysical Models

Numerical cloud models require estimates of the vapor growth rate for ice crystals. Current bulk and bin microphysical parameterizations generally assume that vapor growth is diffusion limited, though some parameterizations include the influence of surface attachment kinetics through a constant deposition coefficient. A parameterization for variable deposition coefficients is provided herein. The parameterization is an explicit function of the ambient ice supersaturation and temperature, and an implicit function of crystal dimensions and pressure. The parameterization is valid for variable surface types including growth by dislocations and growth by step nucleation. Deposition coefficients are predicted for the two primary growth directions of crystals, allowing for the evolution of the primary habits. Comparisons with benchmark calculations of instantaneous mass growth indicate that the parameterization is accurate to within a relative error of 1%. Parcel model simulations using Lagrangian microphysics as a benchmark indicate that the bulk parameterization captures the evolution of mass mixing ratio and fall speed with typical relative errors of less than 10%, whereas the average axis lengths can have errors of up to 20%. The bin model produces greater accuracy with relative errors often less that 10%. The deposition coefficient parameterization can be used in any bulk and bin scheme, with low error, if an equivalent volume spherical radius is provided.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A coarse-grain force field based on quantum mechanics (CGq FF) for molecular dynamics simulation of poly(ethylene glycol)- block -poly(ε-caprolactone) (PEG- b -PCL) micelles

In order to provide the means to predict from molecular dynamics (MD) simulations the structures of copolymer-based micelles in solution, we developed coarse grain force field (CGq FF) parameters for poly(ethylene glycol) (PEG) and for poly(ε-caprolactone) (PCL). A key advance here is the use of quantum mechanics to train the parameters describing the non-bonded (NB) interactions between the CG beads. The functional forms are the same as the MARTINI CG FF so standard MD codes can be used. Our CGq FF describes well the experimentally observed properties for the polymer–air and polymer–water interfaces, indicating the accuracy of the NB interactions. Additionally, the structural properties (density, radius of gyration ( R g ), and end-to-end distance ( h )) match both experiment and all atom (AA) simulations. We illustrate the application of this CGq FF by following the formation of a spherical micelle from 250 chains of PEG 23 - b -PCL 9 diblock copolymer, each block with molecular weight of 1000 Daltons (10 500 beads, corresponding to 123 250 atoms), in a water box with 119 139 water beads (426 553 water molecules).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Near-field radiative heat transfer between irregularly shaped dielectric particles modeled with the discrete system Green's function method

Near-field radiative heat transfer (NFRHT) between irregularly shaped dielectric particles made of SiO 2 and morphology characterized by Gaussian random spheres is studied. Particles are modeled using the discrete system Green's function (DSGF) approach, which is a volume integral numerical method based on fluctuational electrodynamics. This method is applicable to finite, three-dimensional objects, and all system interactions are defined independent of thermal excitation by a generalized system Green's function. The DSGF method is deemed suitable to model NFRHT between irregularly shaped particles after verification against the analytical solution for chains of two and three SiO 2 spheres. The NFRHT results reveal that geometric irregularity in particles leads to a reduction of the total conductance from that of comparable perfect spheres at vacuum separation distances smaller than the particle size, a regime in which NFRHT is a surface phenomenon. At vacuum separation distances larger than the particle size, NFRHT becomes a volumetric process, and the total conductance between irregularly shaped particles converges to that of comparable perfect spheres. Spectral analysis reveals, however, that particle irregularity leads to damping and broadening of resonances at all separation distances, thereby highlighting the importance of the DSGF method for spectral engineering in the near field. The reduced spectral coherence when particle size is larger than the vacuum separation distance is attributed to coupling of surface phonon-polaritons within the randomly generated, distorted particle features. For particle size smaller than the vacuum separation distance, resonance broadening and damping are linked with the multiple localized surface phonon modes supported by the composite spherical harmonic morphologies of the Gaussian random spheres. In conclusion, this paper has direct implications for thermal management of packed particle systems, with applications in radiative property control, electronics, energy conversion, and nanomanufacturing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1996

Waveform Perturbations os Spherical Transiet Waves Propagating in a Random Medium

Of those aspects of sonic-boom propagation that are not yet fully understood, one of the more important is that which relates to the perturbations of the waveform. This phenomenon, which arises also in connection with the propagation of other types of transient waves in real media, generally take the form, in the case of sonic-boom N-waves, of a random high-frequency structure (sometimes called fine structure) that is most prominent in the regions immediately behind each of the shocks. The perturbations in those regions can be large, occasionally attaining magnitudes comparable to that of th incident wave itself. Such magnitudes, in combination with the high-frequency character of the perturbations, can lead to a considerable increase in the perceived noisiness of the sonic boom. Waveform perturbations are consequently an important factor as regards the question of sonic-boom acceptability. On the basis of observations, and some early theoretical studies, it is now generally accepted that perturbations of sonic-boom waveforms are a manifestation of the effect on the propagating wave of relatively small-scale variations in the acoustic properties of the atmosphere- variations that are usually associated with turbulence. Although the mechanism underlying perturbations of sonic-boom waveforms seems thus to be well understood, no fully-satisfactory theory of such perturbations has emerged. Indeed, even for the relatively simple case of an incident step-function pulse, no theory of waveform perturbations, formulated in a realistic three-dimensional context, has been advanced that is valid in the region of strongest perturbations; viz., the region immediately behind, and including, the wave front. The work reported herein represents an attempt to develop such a theory.

Wenzel, Alan R.↗

NASA NTRS2001

Satellite Remote Sensing of Tropical Precipitation and Ice Clouds for GCM Verification

This project, supported by the NASA New Investigator Program, has primarily been funding a graduate student, Darren McKague. Since August 1999 Darren has been working part time at Raytheon, while continuing his PhD research. Darren is planning to finish his thesis work in May 2001, thus some of the work described here is ongoing. The proposed research was to use GOES visible and infrared imager data and SSM/I microwave data to obtain joint distributions of cirrus cloud ice mass and precipitation for a study region in the Eastern Tropical Pacific. These joint distributions of cirrus cloud and rainfall were to be compared to those from the CSU general circulation model to evaluate the cloud microphysical amd cumulus parameterizations in the GCM. Existing algorithms were to be used for the retrieval of cloud ice water path from GOES (Minnis) and rainfall from SSM/I (Wilheit). A theoretical study using radiative transfer models and realistic variations in cloud and precipitation profiles was to be used to estimate the retrieval errors. Due to the unavailability of the GOES satellite cloud retrieval algorithm from Dr. Minnis (a co-PI), there was a change in the approach and emphasis of the project. The new approach was to develop a completely new type of remote sensing algorithm - one to directly retrieve joint probability density functions (pdf's) of cloud properties from multi-dimensional histograms of satellite radiances. The usual approach is to retrieve individual pixels of variables (i.e. cloud optical depth), and then aggregate the information. Only statistical information is actually needed, however, and so a more direct method is desirable. We developed forward radiative transfer models for the SSM/I and GOES channels, originally for testing the retrieval algorithms. The visible and near infrared ice scattering information is obtained from geometric ray tracing of fractal ice crystals (Andreas Macke), while the mid-infrared and microwave scattering is computed with Mie scattering. The radiative transfer is performed with the Spherical Harmonic Discrete Ordinate Method (developed by the PI), and infrared molecular absorption is included with the correlated k-distribution method. The SHDOM radiances have been validated by comparison to version 2 of DISORT (the community "standard" discrete-ordinates radiative transfer model), however we use SHDOM since it is computationally more efficient.

Evans, K. Franklin↗

NASA NTRS1993

A Regional Analysis of Cloudy Mean Spherical Albedo over the Marine Stratocumulus Region and the Tropical Atlantic Ocean

Since clouds are the largest variable in Earth's radiation budget, it is critical to determine both the spatial and temporal characteristics of their radiative properties. This study examines the relationships between cloud properties and cloud fraction in order to supplement grid scale parameterizations. The satellite data used in this study is from three hourly ISCCP (International Satellite Cloud Climatology Project) and monthly ERBE (Earth Radiation Budget Experiment) data on a 2.50 x 2.50 latitude-longitude grid. Mean cloud spherical albedo, the mean optical depth distribution and cloud fraction are examined and compared off the coast of California and the mid-tropical Atlantic for July 1987 and 1988. Individual grid boxes and spatial averages over several grid boxes are correlated to Coakleys (1991) theory of reflection for uniform and broken layered cloud and to Kedem, et al.(1990) findings that rainfall volume and fractional area of rain in convective systems is linear. Kedem's hypothesis can be expressed in terms of cloud properties. That is, the total volume of liquid in a box is a linear function of cloud fraction. Results for the marine stratocumulus regime indicate that albedo is often invariant for cloud fractions of 20% to 80%. Coakley's satellite model of small and large clouds with cores (1 km) and edges (100 in) is consistent with this observation. The cores maintain high liquid water concentrations and large droplets while the edges contain low liquid water concentrations and small droplets. Large clouds are just a collection of cores. The mean optical depth (TAU) distributions support the above observation with TAU values of 3.55 to 9.38 favored across all cloud fractions. From these results, a method based upon Kedem, et al. theory is proposed to separate the cloud fraction and liquid water path (LWP) calculations in a general circulation model (GCM). In terms of spatial averaging, a linear relationship between albedo and cloud fraction is observed. For tropical locations outside the Intertropical Convergence Zone (ITCZ), results of cloud fraction and albedo spatial averaging followed that of the stratus boxes containing few overcast scenes. Both the ideas of Coakley and Kedem, et al. apply. Within the ITCZ, the grid boxes tended to have the same statistical properties as stratus boxes containing many overcast scenes. Because different dynamical forcing mechanisms are present, it is difficult to devise a method for determining subgrid scale variations. Neither of the theories proposed by Kedem, et al. or Coakley works well for the boxes with numerous overcast scenes.

Ginger, Kathryn M.↗

NASA NTRS1993

A regional analysis of cloudy mean spherical albedo over the marine stratocumulus region and the tropical Atlantic Ocean

Since clouds are the largest variable in Earth's radiation budget, it is critical to determine both the spatial and temporal characteristics of their radiative properties. The relationships between cloud properties and cloud fraction are studied in order to supplement grid scale parameterizations. The satellite data used is from three hourly ISCCP (International Satellite Cloud Climatology Project) and monthly ERBE (Earth Radiation Budget Experiment) data on a 2.5 deg x 2.5 deg latitude-longitude grid. Mean cloud spherical albedo, the mean optical depth distribution, and cloud fraction are examined and compared off the coast of California and the mid-tropical Atlantic for July 1987 and 1988. Individual grid boxes and spatial averages over several grid boxes are correlated to Coakley's theory of reflection for uniform and broken layered cloud and to Kedem, et al.'s findings that rainfall volume and fractional area of rain in convective systems is linear. Kedem's hypothesis can be expressed in terms of cloud properties. That is, the total volume of liquid in a box is a linear function of cloud fraction. Results for the marine stratocumulus regime indicate that albedo is often invariant for cloud fractions of 20% to 80%. Coakley's satellite model of small and large clouds with cores (1 km) and edges (100 m) is consistent with this observation. The cores maintain high liquid water concentrations and large droplets while the edges contain low liquid water concentrations and small droplets. Large clouds are just a collection of cores. The mean optical depth (TAU) distributions support the above observation with TAU values of 3.55 to 9.38 favored across all cloud fractions. From these results, a method based upon Kedem, et al's theory is proposed to separate the cloud fraction and liquid water path (LWP) calculations in a general circulation model (GCM). In terms of spatial averaging, a linear relationship between albedo and cloud fraction is observed. For tropical locations outside the Intertropical Convergence Zone (ITCZ), results of cloud fraction and albedo spatial averaging followed that of the stratus boxes containing few overcast scenes. Both the ideas of Coakley and Kedem, et al. apply. Within the ITCZ, the grid boxes tended to have the same statistical properties as stratus boxes containing many overcast scenes. Because different dynamical forcing mechanisms are present, it is difficult to devise a method for determining subgrid scale variations. Neither of the theories proposed by Kedem, et al. or Coakley works well for the boxes with numerous overcast scenes.

Ginger, Kathryn M.↗

NASA NTRS2020

Analysis of Sample Acquisition Dynamics Using Discrete Element Method

The analysis presented in this paper is conducted in the framework of the Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) project, currently under development at NASA Ames Research Center. OceanWATERS aims at designing a simulation environment which allows for testing autonomy of scientific lander missions to the icy moons of our solar system. Mainly focused on reproducing the end effector interaction with the inherent terrain, this paper introduces a novel discrete element method (DEM)-based approach to determine forces and torques acting on the lander’s scoop during the sample acquisition process. An accurate force feedback from the terrain on the scoop is required by fault-detection and autonomous decision-making algorithms to identify when the requested torque on the robotic arm’s joints exceeds the maximum available torque. Knowledge of the terrain force feedback significantly helps evaluating the arm’s links structural properties and properly selecting actuators for the joints. Models available in literature constitute a partial representation of the dynamics of the interaction. As an example, Balovnev derived an analytical expression of the vertical and horizontal force acting on a bucket while collecting a sample as a function of its geometry and velocity, soil parameters and reached depth. Although the model represents an adequate approximation of the two force components, it ignores the direction orthogonal to the scoop motion and neglects the torque. This work relies on DEM analysis to compensate for analytical models’ deficiencies and inaccuracies, i. e. provide force and torque 3D vectors, defined in the moving reference (body) frame attached to the scoop, at each instant of the sample collection process. Results from the first presented analysis relate to the specific OceanWATERS sampling strategy, which consists of collecting the sample through five consecutive passes with increasing depth, each pass following the same circularlinear- circular trajectory. Data is collected given a specific scoop design interacting with two types of bulk materials, which may characterize the surface of icy planetary bodies: snow and ice. Although specifically concerned with the OceanWATERS design, this first analysis provides the expected force trends for similar sampling strategies and allows to deduce phenomenological information about the general scooping process. In order to further instruct the community on the use of DEM tools as a solution to the sampling collection problem, two more analyses have been carried out, mainly focused on reducing the DEM computation time, which increases with a decrease in particle size. After running a set of identical simulations, where the only changing parameter is the size of the spherical particle, it is observed that the resulting force trajectories, starting from a given particle size, converge to the true trend. It is deducible that a further decrease in size yields negligible improvements in the accuracy, while it sensibly increases computation time. A final analysis aims at discussing limitations of approximating bulk material particles having a complex shape, e. g. ice fragments, with spheres, by comparing force trends resulting in the two cases for the same simulation scenario.

Catanoso, Damiana↗

DOE OSTI2023

Kinetic model describing self-limiting CO 2 diffusion in supported amine adsorbents

A reaction–diffusion shrinking core model describing the decay in diffusivity of supported amine sorbents upon CO 2 sorption under both simulated direct air capture and point source capture conditions is described. The decay in CO 2 diffusivity is associated with crosslinking in the aminopolymer samples and general pore blockage in the amino-silane derived samples, which occurs as CO 2 is adsorbed. The model is used to extract four kinetic parameters that govern the CO 2 uptake kinetics and working capacity: an apparent reaction rate constant, an initial effective diffusivity, and two dimensionless decay parameters. Ideally, an initially reaction limited system would allow for direct determination of the intrinsic reaction rate constant; however, sorption experiments suggest mass transfer resistances related to gas mixing, external boundary layers and intraparticle diffusion are present. Reaction rate constants are determined and agree well with theoretical values predicted with the Eyring equation parameterized using density functional theory energies from literature sources. The kinetic performance is expressed as the average effective diffusivity as a function of average conversion, which can be correlated to the dispersion of sorption sites on the support and the morphology of the active sorbent phase. Four supports are impregnated or grafted with amines, SBA-15, single-walled zeolite nanotubes (ZNT), Syloid SiO 2 , and γ-Al 2 O 3 . Due to its pore structure, γ-Al 2 O 3 supported amines can be modeled at the μm scale or at the nm scale, where the shell balance is on the μm-sized macroporous particle aggregate or on the nm-sized amine film on the surface of the Al 2 O 3 nanoparticles, which comprise the spherical particle aggregates. Faster diffusion rates are maintained under 400 ppm rather than 10% CO 2 due to a slower reaction rate giving a slower decay in diffusivity. In conclusion, this work provides a first principles kinetic analysis of CO 2 sorption where previous models are semi-empirical and use arbitrary kinetic parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Isoscalar and isovector giant resonances in 44 Ca, 54 Fe, 64,68 Zn and 56,58,60,68 Ni

We have studied the uncharacteristic behavior of the measured values of the isoscalar and isovector centroid energies, E CEN , of nuclear giant resonances of multipolarity from L = 0 to L = 3 in 44 Ca, 54 Fe, 64,68 Zn and 56,58,60,68 Ni. For this purpose, we carried out calculations of E CEN within the spherical Hartree-Fock (HF)-based random phase approximation (RPA) theory with 33 distinct Skyrme-like effective nucleon-nucleon interactions. We have also determined the Pearson linear correlation coefficients between centroid energies, obtained from the HF-RPA, and the various properties of nuclear matter (NM) of each interaction and determined the sensitivity of E CEN to NM properties. We compared the theoretical values of E CEN obtained from the HF-RPA calculations with experimental data and discuss the results, pointing out significant disagreements between theoretical and experimental values. We note in particular, that we obtain good agreement for the theoretical E CEN of the isovector giant dipole resonance and the available experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗

NASA NTRS2010

Multiphysics Modeling for Dimensional Analysis of a Self-Heated Molten Regolith Electrolysis Reactor for Oxygen and Metals Production on the Moon and Mars

The technology of direct electrolysis of molten lunar regolith to produce oxygen and molten metal alloys has progressed greatly in the last few years. The development of long-lasting inert anodes and cathode designs as well as techniques for the removal of molten products from the reactor has been demonstrated. The containment of chemically aggressive oxide and metal melts is very difficult at the operating temperatures ca 1600 C. Containing the molten oxides in a regolith shell can solve this technical issue and can be achieved by designing a self-heating reactor in which the electrolytic currents generate enough Joule heat to create a molten bath. In a first phase, a thermal analysis model was built to study the formation of a melt of lunar basaltic regolith irradiated by a focused solar beam This mode of heating was selected because it relies on radiative heat transfer, which is the dominant mode of transfer of energy in melts at 1600 C. Knowing and setting the Gaussian-type heat flux from the concentrated solar beam and the phase and temperature dependent thermal properties, the model predicts the dimensions and temperature profile of the melt. A validation of the model is presented in this paper through the experimental formation of a spherical cap melt realized by others. The Orbitec/PSI experimental setup uses an 3.6-cm diameter concentrated solar beam to create a hemispheric melt in a bed of lunar regolith simulant contained in a large pot. Upon cooling, the dimensions of the vitrified melt are measured to validate the thermal model. In a second phase, the model is augmented by multiphysics components to compute the passage of electrical currents between electrodes inserted in the molten regolith. The current through the melt generates Joule heating due to the high resistivity of the medium and this energy is transferred into the melt by conduction, convection and primarily by radiation. The model faces challenges in two major areas, the change of phase as temperature increases, and the dominance of radiative heat flux as heat transfer mechanism within the melt the change of phase concerns the regolith itself which is present in states ranging from a fine grain regolith with low thermal conductivity and low density to a vitrified melt with much higher thermal conductivity, and higher density. As the regolith is heated, it starts to soften around 1300 C the melt iS very viscous and evolving gas bubbles out in thick, lava-like fashion. By 1600 C the regolith is completely melted and the viscosity is low The second challenge resides in the proper modeling of the radiative heat flux requiring the addition of the computing-demanding radiative-heat-transfer function to the general heat transfer equation. The model Includes temperature-dependent properties (density, thermal conductivity, heat capacity, and viscosity, and absorption coefficients) and solves the radiative heat flux equation assuming gray (fine grains) and semi-transparent (melt) media and using an absorption coefficient spectral found in the literature for terrestrial minerals similar in composition to those of lunar regolith simulant

Dominguez, Jesus A.↗

NASA NTRS1994

A spectral formalism for computing three-dimensional deformations due to surface loads. 1: Theory

We outline a complete spectral formalism for computing high spatial resolution three-dimensional deformations arising from the surface mass loading of a spherically symmetric planet. The main advantages of the formalism are that all surface mass loads are always described using a consistent mathematical representation and that calculations of deformation fields for various spatial resolutions can be performed by simpley altering the spherical harmonic degree truncation level of the procedure. The latter may be important when incorporating improved observational constraints on a particular surface mass load, when considering potential errors in the computed field associated with mass loading having a spatial scale unresolved by the observational constraints, or when treating a number of global surface mass loads constrained with different spatial resolutions. The advantages do not extend to traditional 'Green's function' approaches which involve surface element discretizations of the global mass loads. Another advantage of the spectral formalism, over the Green's function approach, is that a posteriori analyses of the computed deformation fields are easily performed. In developing the spectral formalism, we consider specific cases where the Earth's mantle is assumed to respond as an elastic, slightly anelastic, or linear viscoelastic medium. In the case of an elastic or slightly anelastic mantle rheology the spectral response equations incorporate frequency dependent Love numbers. The formalism can therefore be used, for example, to compute the potentially resonant deformational response associated with the free core nutation and Chandler wobble eigenfunctions. For completeness, the spectral response equations include both body forces, as arise from the gravitational attraction of the Sun and the Moon, and surface mass loads. In either case, and for both elastic and anelastic mantle rheologies, we outline a pseudo-spectral technique for computing the ocean adjustment associated with the total gravitational perturbation induced by the external forcing. Three-dimensional deformations computed using the usual Love number approach are generally referenced to an origin at the center of mass of the undeformed planet. We derive a spectral technique for transforming the results to an origin located at the center of mass of the deformed planet.

Mitrovica, J. X.↗

NASA NTRS2005

On Modeling Air/Space-Borne Radar Returns in the Melting Layer

The bright band is the enhanced radar echo associated with the melting of hydrometeors in stratiform rain where the melting process usually occurs below 0 C isotherm over a distance of about 500m. To simulate this radar signature, a scattering model of melting snow is proposed in which the fractional water content is prescribed as a function of the radius of a spherical mixed- phase particle consisting of air, ice and water. The model is based on the observation that melting starts at the surface of the particle and then gradually develops towards the center. To compute the scattering parameters of a non-uniform melting particle, the particle is modeled as a sphere represented by a collection of 64(exp 3) cubic cells of identical size where the probability of water at any cell is prescribed as a function of the radius. The internal field of the particle, used for deriving the effective dielectric constant, is computed by the Conjugate Gradient and Fast Fourier Transform (CGFFT) numerical methods. To make computations of the scattering parameters more efficient, a multi-layer stratified-sphere scattering model is introduced after demonstrating that the scattering parameters of the non-uniformly melting particle can be accurately reproduced by the stratified sphere. In conjunction with a melting layer model that describes the melting fractions and fall velocities of hydrometeors as a function of the distance from the 0 C isotherm, the stratified-sphere model is used to simulate the radar bright band profiles. These simulated profiles are shown to compare well with measurements from the Precipitation Radar (PR) aboard the Tropical Rainfall Measuring Mission (TRMM) satellite and a dual-wavelength airborne radar. The results suggest that the proposed model of a melting snow particle may be useful in studying the characteristics of the bright-band in particular and mixed- phase hydrometeors in general.

Liao, Liang↗