Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Material balance time”

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 271 records · Page 15

Comparative life cycle assessment of a modular cross-laminated timber residential building designed for disassembly and reuse versus traditional wood frame construction

There is a need for affordable housing across the U.S., with high-performance modular and prefabricated buildings providing a logical avenue for meeting some of this demand. However, there is a need to balance high performance construction – including low emissions – with affordability. To provide a proof-of-concept in meeting these goals, the Circular Home is a cross-laminated timber (CLT)-based deconstructible and reconfigurable single-family residence that meets high performance targets in moisture, energy, design, economics, and life cycle assessment (LCA). This study focuses on the LCA, presenting a cradle-to-cradle whole-building life cycle assessment (WBLCA) for the Circular Home and a functionally equivalent Baseline Home constructed with traditional materials and methods. The functional unit is 1 m 2 of gross floor area across 60 years. Revit building information models (BIM) provided material quantities and Tally LCA was utilized for impact data (inclusive of biogenic carbon sequestration), supplemented with manufacturer environmental product declarations (EPDs). The Circular Home outperforms the baseline residence in most measured impact categories, including global warming potential (GWP), producing −2.73 kgCO 2 eq/m 2 in embodied emissions, whereas the modeled baseline has an embodied GWP of 428 kgCO 2 eq/m 2 . The careful material selection and advanced building design optimizes performance, with the Circular Home containing only −0.006 times the embodied emissions and −0.02 times the operational emissions of its traditional counterpart. Finally, the unique contribution of this work is in the environmental impact comparison of a high-performance modular CLT structure that can be affordably scaled and mass produced in a U.S. market, compared to typical single family home construction.

Circularity↗

Laser heating and evaporation of a single droplet

The laser technology is being abundantly studied for controlled energy deposition for a range of applications in aerodynamic flow control, material processing, ignition, and combustion. The absorption of laser radiation by liquid droplets affects further propagation of laser in the atmosphere and causes bleaching of suspended droplets while the ignition and combustion characteristics in combustors are influenced by the evaporation rate of the sprayed fuel. In this work, we present a multi-dimensional mathematical model built on OpenFOAM for laser heating and evaporation of a single droplet in the diffusion dominated regime taking into account absorption of the laser radiation, evaporation process and vapor flow dynamics. The developed solver is validated against available experimental and numerical data for the ethanol and water droplet heating and evaporation. For continuous heating the peak temperature is established by the balance of cooling, evaporation and heating and results in high temperature for larger droplets. It has been shown that for heating by a single laser pulse the maximum temperature of droplets depends only on the peak intensity of the laser radiation. Furthermore, for the peak irradiance close to the transition to the boiling regime, temporal dynamics of the droplet temperature is independent of the droplet size. With proper normalization of time, the dynamics of the droplet shrinkage and cooling is shown to be independent of droplet sizes and peak laser intensities. The influence of cooling and evaporation processes on droplet heating was found to be controlled by the pulse repetition rate for repeated pulse operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Automated eddy current analysis of materials

This research effort focused on the use of eddy current techniques for characterizing flaws in graphite-based filament-wound cylindrical structures. A major emphasis was on incorporating artificial intelligence techniques into the signal analysis portion of the inspection process. Developing an eddy current scanning system using a commercial robot for inspecting graphite structures (and others) has been a goal in the overall concept and is essential for the final implementation for expert system interpretation. Manual scans, as performed in the preliminary work here, do not provide sufficiently reproducible eddy current signatures to be easily built into a real time expert system. The expert systems approach to eddy current signal analysis requires that a suitable knowledge base exist in which correct decisions as to the nature of the flaw can be performed. In eddy current or any other expert systems used to analyze signals in real time in a production environment, it is important to simplify computational procedures as much as possible. For that reason, we have chosen to use the measured resistance and reactance values for the preliminary aspects of this work. A simple computation, such as phase angle of the signal, is certainly within the real time processing capability of the computer system. In the work described here, there is a balance between physical measurements and finite element calculations of those measurements. The goal is to evolve into the most cost effective procedures for maintaining the correctness of the knowledge base.

Workman, Gary L.↗

Low-phase-noise surface-acoustic-wave oscillator using an edge mode of a phononic band gap

Low-phase-noise microwave-frequency integrated oscillators provide compact solutions for various applications in signal processing, communications, and sensing. Surface acoustic waves (SAWs), featuring orders-of-magnitude shorter wavelength than electromagnetic waves at the same frequency, enable integrated microwave-frequency systems with much smaller footprint on chip. SAW devices also allow higher-quality (Q) factors than electronic components at room temperature. Here, we demonstrate a low-phase-noise gigahertz-frequency SAW oscillator on 128°Y-cut lithium niobate, where the SAW resonator occupies a footprint of 0.05 mm 2 . Leveraging phononic crystal band-gap-edge modes to balance between Q factors and insertion losses, our 1-GHz SAW oscillator features a low phase noise of -132.5 dBc/Hz at a 10-kHz offset frequency and an overlapping Hadamard deviation of 6.5 × 10 -10 at an analysis time of 64 ms. The SAW resonator-based oscillator holds high potential in developing low-noise sensors and acousto-optic integrated circuits.

42 ENGINEERING↗

Capability in Theory, Modeling, and Validation for a Range of Innovative Fusion Concepts using High-Fidelity Moment-Kinetic Models

A computational modeling capability is created and available to the fusion community to understand and design lower-cost and innovative fusion concepts. The approach uses high- fidelity kinetic, moment-kinetic, and moment models and includes sophisticated plasma- boundary interactions. A majority of fusion-relevant simulations are performed with magnetohydrodynamic models and hybrid particle-in-cell codes, with limited-fidelity electron and kinetic physics. However, in fusion configurations like Z-pinches, field-reversed- configurations, plasma jet magneto-inertial fusion, spinning mirrors, and others, kinetic effects (both electron and ions) are critical to understand the physics and design scaling into the highly kinetic regime of a burning fusion plasma. Furthermore, as present fusion machines move towards a burning plasma regime, liquid-metal blankets are needed to handle first-wall heat- flux, reduce erosion, and eventually for energy conversion and fuel breeding. The work performed under this ARPA-E BETHE Capability Team advances the state-of-the-art in modeling and understanding plasma dynamics in fusion devices and its coupling with liquid-metal dynamics. These are critical areas of research for fusion energy to become realizable. To address these complex problems, we have leveraged and extended computational capabilities through the code, Gkeyll (developed jointly with Princeton Plasma Physics Laboratory and academic partners), for kinetic and moment modeling of fusion plasmas. The Concept Teams supported by this Capability Team include the Wisconsin High-field Axisymmetric Mirror (WHAM), Centrifugal Mirror Experiment (CFME), Plasma-Jet Magneto- Inertial Fusion (PJMIF), and solid and liquid wall plasma-material interaction studies relevant to a number of fusion concepts including Zap Energy’s Z-pinch. This software is open-source and available to the fusion community as a high-fidelity tool for the design of lower-cost fusion experiments. 3D gyrokinetic simulations of WHAM are now possible for long enough time scales to understand the evolution of interchange instabilities. 3D multi-fluid simulations of CMFE at higher Mach numbers are now possible for detailed design iterations with the goal of stability. The state-of-the-art in understanding shock formation and shock mitigation regimes in merging liners for PJMIF have been furthered by our kinetic simulations. Our novel models and frameworks studying plasma-material interaction by incorporating wall emission for various solid wall materials of relevance to pulsed and steady fusion concepts have advanced the state-of-the-art in our understanding of particle fluxes, heat fluxes, and other quantities at cathodes and anodes. The results from this work may explain discrepancies between experimental and theoretical predictions of achieved current densities in pulsed concepts such as Z-pinches. Another significant contribution of this Capability Team is the development and deployment of a novel experimental platform, LEX (Liquid Electrode eXperiment), at Virginia Tech to understand liquid metal free-surface response to electromagnetic pulses. The novel experiments along with model validation quantified the effect of different materials and sizes of liquid metal droplets on the radiative power balance of fusion plasmas for pulsed concepts. Furthermore, these experiments provided mitigation strategies for violent liquid metal response for high current pulses as would be expected in fusion regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resolving Salt-Induced Agglomeration of Laponite Suspensions Using X-ray Photon Correlation Spectroscopy and Molecular Dynamics Simulations

Linking the physics of the relaxation behavior of viscoelastic fluids as they form arrested gel states to the underlying chemical changes is essential for developing predictive controls on the properties of the suspensions. In this study, 3 wt.% laponite suspensions are studied as model systems to probe the influence of salt-induced relaxation behavior arising from the assembly of laponite nanodisks. X-ray Photon Correlation Spectroscopy (XPCS) measurements show that laponite suspensions prepared in the presence of 5 mM concentrations of CaCl2, MgCl2 and CsCl salts accelerate the formation of arrested gel states, with CaCl2 having a significant impact followed by CsCl and MgCl2 salts. The competing effects of ion size and charge on relaxation behavior are noted. For example, the relaxation times of laponite suspensions in the presence of Mg2+ ions are slower compared to Cs+ ions despite the higher charge, suggesting that cation size dominates in this scenario. The faster relaxation behavior of laponite suspensions in the presence of Ca2+ ions compared to Cs+ ions shows that a higher charge dominates the size of the ion. The trends in relaxation behavior are consistent with the cluster formation behavior of laponite suspensions and the electrostatic interactions predicted from MD simulations. Charge balance is achieved by the intercalation of the cations at the negatively charged surfaces of laponite suspensions. These studies show that the arrested gel state of laponite suspensions is accelerated in the presence of salts, with ion sizes and charges having a competing effect on relaxation behavior.

36 MATERIALS SCIENCE↗

Optimization of the Second Target Station cold source moderators using an automated workflow

The Second Target Station (STS) at the US Department of Energy’s Oak Ridge National Laboratory is designed to become the world’s highest peak-brightness spallation source of cold neutrons. Successful completion of the STS, which is currently in the preliminary design phase, will provide transformative new capabilities to examine novel materials for future technologies. At STS, neutrons will be generated by spallation reactions in a solid tungsten target. They will be moderated and thermalized in two cold (20 K) para-hydrogen moderators. Careful optimization of these moderators is essential to the project’s success. To find optimal moderator designs, an advanced optimization workflow integrates high-fidelity neutronics calculations using the Monte Carlo N-Particle (MCNP) transport code MCNP6.2 with state-of-the-art optimization algorithms in the Dakota optimization toolkit. For each design iteration, a parametrized solid CAD geometry is generated in Creo and automatically converted into an unstructured mesh geometry by Attila 4MC for the neutronics calculation with MCNP. Iterations repeat until optimal designs are found. Herein this paper presents the results of a sensitivity and optimization study for the cylindrical and tube moderators. Both moderators can be optimized for maximum peak brightness, maximum time-integrated brightness, or any combination between these extremes. Maximum peak brightness is achieved by using smaller optimal dimensions of the moderators, whereas maximum time-integrated brightness is achieved by using larger dimensions. A Pareto front details the designs that optimally balance both brightness metrics. The Pareto front can be found in only 40–110 iterations with 4–5 design parameters when using the efficient global and Pareto-set optimization algorithms in Dakota. Additionally, important engineering constraints can be taken into account, such as the coupling between the cylindrical moderator radius and aluminum vessel wall thicknesses required to ensure structural integrity of the vessels. This interaction has a significant impact on the resulting optimal designs. Our new, highly efficient, fully automated optimization workflow will be used to optimize additional STS components in the future and can be adopted for design and optimization studies at other experimental neutron and accelerator facilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Laboratory Spectroscopy of Large Carbon Molecules and Ions in Support of Space Missions. A New Generation of Laboratory & Space Studies

Polycyclic Aromatic Hydrocarbons (PAHs) are an important and ubiquitous component of carbon-bearing materials in space. A long-standing and major challenge for laboratory astrophysics has been to measure the spectra of large carbon molecules in laboratory environments that mimic (in a realistic way) the physical conditions that are associated with the interstellar emission and absorption regions [1]. This objective has been identified as one of the critical Laboratory Astrophysics objectives to optimize the data return from space missions [2]. An extensive laboratory program has been developed to assess the properties of PAHs in such environments and to describe how they influence the radiation and energy balance in space. We present and discuss the gas-phase electronic absorption spectra of neutral and ionized PAHs measured in the UV-Visible-NIR range in astrophysically relevant environments and discuss the implications for astrophysics [1]. The harsh physical conditions of the interstellar medium characterized by a low temperature, an absence of collisions and strong VUV radiation fields - have been simulated in the laboratory by associating a pulsed cavity ringdown spectrometer (CRDS) with a supersonic slit jet seeded with PAHs and an ionizing, penning-type, electronic discharge. We have measured for the {\it first time} the spectra of a series of neutral [3,4] and ionized [5,6] interstellar PAHs analogs in the laboratory. An effort has also been attempted to quantify the mechanisms of ion and carbon nanoparticles production in the free jet expansion and to model our simulation of the diffuse interstellar medium in the laboratory [7]. These experiments provide {\it unique} information on the spectra of free, large carbon-containing molecules and ions in the gas phase. We are now, for the first time, in the position to directly compare laboratory spectral data on free, cold, PAH ions and carbon nano-sized carbon particles with astronomical observations in the UV-NIR range (interstellar UV extinction, DIBs in the NUV-NIR range). This new phase offers tremendous opportunities for the data analysis of current and upcoming space missions geared toward the detection of large aromatic systems Le., the "new frontier space missions" (Spitzer, HST, COS, JWST, SOFIA,...).

Salama, Farid↗

Thermal Testing and Model Correlation of the Magnetospheric Multiscale (MMS) Observatories

International Conference on Envronmental Systems (ICES), Seattle WA NCTS 20964-15. The Magnetospheric Multiscale (MMS) mission is a Solar Terrestrial Probes mission comprising four identically instrumented spacecraft that will use Earths magnetosphere as a laboratory tostudy the microphysics of three fundamental plasma processes: magnetic reconnection, energetic particle acceleration, and turbulence. This paper presents the complete thermal balance (TB) test performed on the first of four observatories to go through thermal vacuum (TV) and the minibalance testing that was performed on the subsequent observatories to provide a comparison of all four. The TV and TB tests were conducted in a thermal vacuum chamber at the Naval Research Laboratory (NRL) in Washington, D.C. with the vacuum level higher than 1.3 x 10-4 Pa (10-6 torr)and the surrounding temperature achieving -180 C. Three TB test cases were performed that included hot operational science, cold operational science and a cold survival case. In addition to the three balance cases a two hour eclipse and a four hour eclipse simulation was performed during the TV test to provide additional transient data points that represent the orbit in eclipse (or Earth's shadow) The goal was to perform testing such that the flight orbital environments could be simulated as closely as possible. A thermal model correlation between the thermal analysis and the test results was completed. Over 400 1-Wire temperature sensors, 200 thermocouples and 125 flight thermistor temperature sensors recorded data during TV and TB testing. These temperatureversus time profiles and their agreements with the analytical results obtained using Thermal Desktop and SINDAFLUINT are discussed. The model correlation for the thermal mathematical model (TMM) is conducted based on the numerical analysis results and the test data. The philosophy of model correlation was to correlate the model to within 3 C of the test data using the standard deviation and mean deviation error calculation. Individual temperature error goal is to be within 5 C and the heater power goal is to be within 5 of test data. The results of the model correlation are discussed and the effect of some material and interface parameters on the temperature profiles are presented.

Thermal↗

Sensitivity of magnetohydrodynamic simulations of Joule-heated conductors to the vaporization curve in equations of state

Magnetohydrodynamic (MHD) simulations of electrically exploded aluminum and copper rods demonstrate a technique to validate equations of state (EOS) for rapidly Joule-heated conductors. The balance of internal and magnetic forces at the conductor-insulator interface drives the metal there along the vaporization phase boundary. Variations between critical points and vaporization curves in existing models predict differing densities and temperatures in MHD simulations for these models. Here, the inclusion of Maxwell constructs in the liquid-vapor biphase region of the EOS caused the rod surface to vaporize earlier in time than unmodified tables with van der Waals loops. Velocimetry of recent experiments is used to validate the location of the vaporization curve in existing EOS models and differentiate between the vapor dome treatments. Dielectric coatings applied to the metal surface restricted the conductor’s expansion and diverted the metal into the warm dense matter regime.

36 MATERIALS SCIENCE↗

Investigating the Electromechanical Sensitivity of Carbon-Nanotube-Coated Microfibers

The piezoresistance of carbon nanotube (CNT)-coated microfibers is examined using diametric compression. Diverse CNT forest morphologies were studied by changing the CNT length, diameter, and areal density via synthesis time and fiber surface treatment prior to CNT synthesis. Large-diameter (30–60 nm) and relatively low-density CNTs were synthesized on as-received glass fibers. Small-diameter (5–30 nm) and-high density CNTs were synthesized on glass fibers coated with 10 nm of alumina. The CNT length was controlled by adjusting synthesis time. Electromechanical compression was performed by measuring the electrical resistance in the axial direction during diametric compression. Gauge factors exceeding three were measured for small-diameter (<25 μm) coated fibers, corresponding to as much as 35% resistance change per micrometer of compression. The gauge factor for high-density, small-diameter CNT forests was generally greater than those for low-density, large-diameter forests. A finite element simulation shows that the piezoresistive response originates from both the contact resistance and intrinsic resistance of the forest itself. The change in contact and intrinsic resistance are balanced for relatively short CNT forests, while the response is dominated by CNT electrode contact resistance for taller CNT forests. These results are expected to guide the design of piezoresistive flow and tactile sensors.

36 MATERIALS SCIENCE↗

Update to NASA's Clinical Practice Guideline for Renal Stone Surveillance in U.S. Astronauts

BACKGROUND: For the past 9 years all NASA astronauts received specialized renal ultrasound screening for precursors to renal stones, our so called “mineralized renal material” (MRM), to reduce the possibility of sending a person to space with a renal stone that could cause morbidity and threaten mission objectives. A recent analysis of this data shows no increased incidence of clinical or sub-clinical renal stones associated with spaceflight; thus, our clinical practice guideline has been updated. OVERVIEW: Astronauts previously received a renal ultrasound screening at selection and at every annual recertification exam. In total, 73 individuals have been scanned and 282 total studies have been performed as of April 2023.This data shows no observed statistical change seen in MRM burden pre-to postflight, therefore a reduction in the frequency of the screening exams is appropriate. DISCUSSION: Data presented in this panel have shown no direct correlation between spaceflight and development of MRM or renal stones in astronaut kidneys out to at least 6 months of spaceflight. The MRM burden does not seem to change appreciably over time and so far, has no demonstrated correlation to external factors. Therefore, we find it appropriate, in lieu of yearly exams, to reduce the frequency of screening to: initial selection, first annual exam after selection, and pre/postflight. This balances the screening frequency of our astronauts, against the potential for unnecessary medical procedures, while still serving the purpose of risk reduction.

D Reyes↗

Covalent Functionalization of Silicon with Plasma-Grown “Fuzzy” Graphene: Robust Aqueous Photoelectrodes for CO 2 Reduction by Molecular Catalysts

Carbon electrodes are ideal for electrochemistry with molecular catalysts, exhibiting facile charge transfer and good stability. Yet for solar-driven catalysis with semiconductor light absorbers, stable semiconductor/carbon interfaces can be difficult to achieve, and carbon’s high optical extinction means it can only be used in ultrathin layers. Here, we demonstrate a plasma-enhanced chemical vapor deposition process that achieves well-controlled deposition of out-of-plane “fuzzy” graphene (FG) on thermally oxidized Si substrates. The resulting Si|FG interfaces possess a silicon oxycarbide (SiOC) interfacial layer, implying covalent bonding between Si and the FG film that is consistent with the mechanical robustness observed from the films. The FG layer is uniform and tunable in thickness and optical transparency by deposition time. Using p-type Si|FG substrates, noncovalent immobilization of cobalt phthalocyanine (CoPc) molecular catalysts was employed for the photoelectrochemical reduction of CO 2 in aqueous solution. The Si|FG|CoPc photocathodes exhibited good catalytic activity, yielding a current density of ∼1 mA/cm 2 , Faradaic efficiency for CO of ∼70% (balance H 2 ), and stable photocurrent for at least 30 h at −1.5 V vs Ag/AgCl under 1-sun illumination. Furthermore, the results suggest that plasma-deposited FG is a robust carbon electrode for molecular catalysts and suitable for further development of aqueous-stable Si photocathodes for CO 2 reduction.

CO2 reduction↗

Electron dynamics in extended systems within real-time time-dependent density-functional theory

Due to a beneficial balance of computational cost and accuracy, real-time time-dependent density-functional theory has emerged as a promising first-principles framework to describe electron real-time dynamics. Here we discuss recent implementations around this approach, in particular in the context of complex, extended systems. Results include an analysis of the computational cost associated with numerical propagation and when using absorbing boundary conditions. Here we extensively explore the shortcomings for describing electron–electron scattering in real time and compare to many-body perturbation theory. Modern improvements of the description of exchange and correlation are reviewed. In this work, we specifically focus on the Qb@ll code, which we have mainly used for these types of simulations over the last years, and we conclude by pointing to further progress needed going forward.

2D materials↗

Thermal and Mechanical Performance of a Carbon/Carbon Composite Spacecraft Radiator

Carbon-carbon composite materials offer greater thermal efficiency, stiffness to weight ratio, tailorability, and dimensional stability than aluminum. These lightweight thermal materials could significantly reduce the overall costs associated with satellite thermal control and weight. However, the high cost and long lead-time for carbon-carbon manufacture have limited their widespread usage. Consequently, an informal partnership between government and industrial personnel called the Carbon-Carbon Spacecraft Radiator Partnership (CSRP) was created to foster carbon-carbon composite use for thermally and structurally demanding space radiator applications. The first CSRP flight opportunity is on the New Millennium Program (NMP) Earth Orbiter-1 (EO-1) spacecraft, scheduled for launch in late 1999. For EO-1, the CSRP designed and fabricated a Carbon-Carbon Radiator (CCR) with carbon-carbon facesheets and aluminum honeycomb core, which will also serve as a structural shear panel. While carbon-carbon is an ideal thermal candidate for spacecraft radiators, in practice there are technical challenges that may compromise performance. In this work, the thermal and mechanical performance of the EO-1 CCR is assessed by analysis and testing. Both then-nal and mechanical analyses were conducted to predict the radiator response to anticipated launch and on-orbit loads. The thermal model developed was based on thermal balance test conditions. The thermal analysis was performed using SINDA version 4.0. Structural finite element modeling and analysis were performed using SDRC/1-DEAS and UAI/NASTRAN, respectively. In addition, the CCR was subjected to flight qualification thermal/vacuum and vibration tests. The panel meets or exceeds the requirements for space flight and demonstrates promise for future satellite missions.

Kuhn, Jonathan↗

Thermal Testing and Model Correlation of the Magnetospheric Multiscale (MMS) Observatories

The Magnetospheric Multiscale (MMS) mission is a Solar Terrestrial Probes mission comprising four identically instrumented spacecraft that will use Earth's magnetosphere as a laboratory to study the microphysics of three fundamental plasma processes: magnetic reconnection, energetic particle acceleration, and turbulence. This paper presents the complete thermal balance (TB) test performed on the first of four observatories to go through thermal vacuum (TV) and the minibalance testing that was performed on the subsequent observatories to provide a comparison of all four. The TV and TB tests were conducted in a thermal vacuum chamber at the Naval Research Laboratory (NRL) in Washington, D.C. with the vacuum level higher than 1.3 x 10 (sup -4) pascals (10 (sup -6) torr) and the surrounding temperature achieving -180 degrees Centigrade. Three TB test cases were performed that included hot operational science, cold operational science and a cold survival case. In addition to the three balance cases a two hour eclipse and a four hour eclipse simulation was performed during the TV test to provide additional transient data points that represent the orbit in eclipse (or Earth's shadow) The goal was to perform testing such that the flight orbital environments could be simulated as closely as possible. A thermal model correlation between the thermal analysis and the test results was completed. Over 400 1-Wire temperature sensors, 200 thermocouples and 125 flight thermistor temperature sensors recorded data during TV and TB testing. These temperature versus time profiles and their agreements with the analytical results obtained using Thermal Desktop and SINDA/FLUINT are discussed. The model correlation for the thermal mathematical model (TMM) is conducted based on the numerical analysis results and the test data. The philosophy of model correlation was to correlate the model to within 3 degrees Centigrade of the test data using the standard deviation and mean deviation error calculation. Individual temperature error goal is to be within 5 degrees Centigrade and the heater power goal is to be within 5 percent of test data. The results of the model correlation are discussed and the effect of some material and interface parameters on the temperature profiles are presented.

Correlation↗

Three-Dimensional High-Order Spectral Volume Method for Solving Maxwell's Equations on Unstructured Grids

A three-dimensional, high-order, conservative, and efficient discontinuous spectral volume (SV) method for the solutions of Maxwell's equations on unstructured grids is presented. The concept of discontinuous 2nd high-order loca1 representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) method, but instead of using a Galerkin finite-element formulation, the SV method is based on a finite-volume approach to attain a simpler formulation. Conventional unstructured finite-volume methods require data reconstruction based on the least-squares formulation using neighboring cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In the SV method, one starts with a relatively coarse grid of triangles or tetrahedra, called spectral volumes (SVs), and partition each SV into a number of structured subcells, called control volumes (CVs), that support a polynomial expansion of a desired degree of precision. The unknowns are cell averages over CVs. If all the SVs are partitioned in a geometrically similar manner, the reconstruction becomes universal as a weighted sum of unknowns, and only a few universal coefficients need to be stored for the surface integrals over CV faces. Since the solution is discontinuous across the SV boundaries, a Riemann solver is thus necessary to maintain conservation. In the paper, multi-parameter and symmetric SV partitions, up to quartic for triangle and cubic for tetrahedron, are first presented. The corresponding weight coefficients for CV face integrals in terms of CV cell averages for each partition are analytically determined. These discretization formulas are then applied to the integral form of the Maxwell equations. All numerical procedures for outer boundary, material interface, zonal interface, and interior SV face are unified with a single characteristic formulation. The load balancing in a massive parallel computing environment is therefore easier to achieve. A parameter is introduced in the Riemann solver to control the strength of the smoothing term. Important aspects of the data structure and its effects to communication and the optimum use of cache memory are discussed. Results will be presented for plane TE and TM waves incident on a perfectly conducting cylinder for up to fifth order of accuracy, and a plane wave incident on a perfectly conducting sphere for up to fourth order of accuracy. Comparisons are made with exact solutions for these cases.

Liu, Yen↗

A Parallelized Oxidation-Driven Surface Recession Framework in DSMC Code, SPARTA

Spacecrafts rely on ablative thermal protection systems (TPS) made of composites consisting of a carbon-based reinforcement and a polymeric matrix. These materials are designed to withstand high-temperature oxidation and surface recession during re-entry into the Earth's atmosphere. However, ablation occurs due to a complex interplay of thermal, mechanical, and chemical factors, making it challenging to determine the individual impact of each on the TPS's overall degradation. In this study, we have developed an ablation model that can leverage a finite rate carbon oxidation model to predict material recession and surface states more accurately. Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al can perform both gas-surface and pure-surface reactions and is now adapted to tally surface data on a per grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also briefly discuss some unique challenges associated with parallelizing this dynamic surface state and geometry. Finally, we analyze the performance of this parallelized implicit chemistry model with simple 2D and 3D benchmark cases by producing surface state statistics, area changes over time, and visualization across a range of surface temperatures and processors with and without load-balancing.

DSMC↗