Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point Kinetics”

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 847 records · Page 47

Mars Ground Level Enhancements in the Context of the Solar Energetic Particle Clock

In this work we discuss the growing ensemble of solar particle events registered on the Martian surface, including their temporal appearance and solar sources. Solar energetic particle events from the surface of Mars have been observed starting soon after the August 2012 landing of the Radiation Assessment Detector onboard Curiosity. The Martian atmosphere prevents protons and heavy ions of up to 180 MeV/n kinetic energy from directly reaching the Martian surface. This cut-off is high enough to limit the number of solar energetic particle events measured on the surface to only 15 in ~12 ½ years. Yet we find in this analysis that Mars ground level enhancements follow the distinct SEP clock pattern as proton events observed at lower energies as reported in Posner, Richardson and Strauss (2024). Proton acceleration occurs predominantly near the solar surface, while transport to Mars incurs a delay in onset, and, as we show here, peak, that is a function of the longitudinal magnetic connection distance. between the foot point of solar wind magnetic field lines that intersect the Mars environment with the source longitude of the solar magnetic eruption. A distinct clustering of relative solar source locations at or near the Mars foot points at the Sun’s western limb is apparent, indicating lower flux thresholds from such preferred locations. Our findings have implications for astronaut safety at Mars.

Mars Ground Level Enhancements↗

Comprehensive uncertainty quantification (UQ) for full engineering models by solving probability density function (PDF) equation

This report details a new method for propagating parameter uncertainty (forward uncertainty quantification) in partial differential equations (PDE) based computational mechanics applications. The method provides full-field quantities of interest by solving for the joint probability density function (PDF) equations which are implied by the PDEs with uncertain parameters. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known apriori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the driving (dynamics) term of the PDEs based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the dynamics term of the governing PDEs. The efficacy of the method is demonstrated by applying it to two target problems: solid mechanics explicit dynamics with uncertain material model parameters, and reacting hypersonic fluid mechanics with uncertain chemical kinetic rate parameters. A minimally invasive implementation of the method for representative codes SPARC (reacting hypersonics) and NimbleSM (finite- element solid mechanics) and associated software details are described. For solid mechanics demonstration problems the method shows order of magnitudes improvement in accuracy over traditional stochastic collocation. For the reacting hypersonics problem, the method is implemented as a streamline integration and results show very good accuracy for the approximate sample solutions of re-entry flow past the Apollo capsule geometry at Mach 30.

42 ENGINEERING↗

Optimal aeroassisted coplanar orbital transfer using an energy model

The atmospheric portion of the trajectories for the aeroassisted coplanar orbit transfer was investigated. The equations of motion for the problem are expressed using reduced order model and total vehicle energy, kinetic plus potential, as the independent variable rather than time. The order reduction is achieved analytically without an approximation of the vehicle dynamics. In this model, the problem of coplanar orbit transfer is seen as one in which a given amount of energy must be transferred from the vehicle to the atmosphere during the trajectory without overheating the vehicle. An optimal control problem is posed where a linear combination of the integrated square of the heating rate and the vehicle drag is the cost function to be minimized. The necessary conditions for optimality are obtained. These result in a 4th order two-point-boundary-value problem. A parametric study of the optimal guidance trajectory in which the proportion of the heating rate term versus the drag varies is made. Simulations of the guidance trajectories are presented.

Halyo, Nesim↗

Development of Thermostable Carbonic Anhydrases Using Structure-Guided Recombination for Use in CO2 Removal Systems on Spacecraft

Carbon capture and storage has been a research area of great interest in recent years as a method for mitigation of CO2 emissions, due to the effects of climate change. The development of technologies for the efficient capture of CO2 are also of great interest for human spaceflight applications. One of the most promising technologies in this area is CO2 scrubbing using liquid amines, unfortunately, liquid amines with low heats of desorption tend to have slow CO2 binding kinetics. One potential solution to this problem is to use the enzyme carbonic anhydrase (CA) to enhance the kinetics of CO2 binding to liquid amines, allowing the overall process to be more energy efficient. Interest in using carbonic anhydrase as a biocatalyst has led to a number of efforts to improve the thermostability and solvent tolerance of several distinct carbonic anhydrase enzymes. In the work described here, we screened through a diverse set of natural carbonic anhydrases to identify candidates for protein engineering aimed at increased stability and activity in various liquid amines. We then used SCHEMA structure-guided recombination to develop a set of chimeric carbonic anhydrases with high thermostability and activity. These chimeras were used as the starting points for further protein engineering work targeting activity in liquid amine systems. Our ultimate goal is to test the engineered enzymes in a liquid amine system for cabin air revitalization on ISS or other spacecraft.

Life Support, Carbon dioxide, Carbonic anhydrase, ↗

Acoustic misalignment mechanism for axion dark matter

A rotation in the field space of a complex scalar field corresponds to a Bose-Einstein condensation of U(1) charges. We point out that fluctuations in this rotating condensate exhibit sound-wave modes, which can be excited by cosmic perturbations and identified with axion fluctuations once the U(1) charge condensate has been sufficiently diluted by cosmic expansion. We consider the possibility that these axion fluctuations constitute dark matter and develop a formalism to compute its abundance. We carefully account for the growth of fluctuations during the epoch where the complex scalar field rotates on the body of the potential and possible nonlinear evolution when the fluctuations become non-relativistic. We find that the resultant dark matter abundance can exceed the conventional and kinetic misalignment contributions if the radial direction of the complex scalar field is sufficiently heavy. The axion dark matter may also be warm enough to leave imprints on structure formation. We discuss the implications of this novel dark matter production mechanism — acoustic misalignment mechanism — for the axion rotation cosmology, including kination domination and baryogenesis from axion rotation, as well as for axion searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Sub-Ambient Performance of Potassium Sarcosinate for Direct Air Capture Applications: CO 2 Flux and Viscosity Measurements

Absorption based direct air capture (DAC) technologies have garnered significant interest in recent years due to their scalability, competitive regeneration energy requirements, and low susceptibility to degradation. One of the key advantages of DAC lies in its flexible siting options and the potential to utilize low-value land. However, most of the research in this field has been focused on ambient climate zones (T > 20 °C), overlooking sub-ambient (–30 °C < T < 20 °C) regions, which comprise approximately 70 % of the Earth’s surface. To fully realize the potential of DAC, it is essential to understand how DAC solvents perform in these sub-ambient conditions before any large-scale deployment can be considered. Among DAC solvents of interest, potassium sarcosinate (K-SAR) has emerged as a promising candidate due to its high CO 2 capacity, fast uptake kinetics, compatibility with contactor packing materials, low volatility, good thermal and oxidative stability, and competitive regeneration energy requirements compared to current industry standards. This paper characterizes the CO 2 flux and viscosity of K-SAR at sub-ambient conditions and explores the potential of using ethylene glycol and triethylene glycol as additives to prevent solvent freezing in DAC applications. For 1 M K-SAR, the CO 2 flux ranges between 1.3 × 10 -5 and 8.0 × 10 -5 mol m –2 s –1 across a temperature range of –5 °C to 45 °C. Ethylene glycol is shown to effectively suppress the freezing point of K-SAR below –30 °C with volumetric loadings of the additive as low as 0.1. Here, a reaction model was developed to predict the CO 2 flux for 1 M K-SAR at different temperatures, demonstrating good agreement between experimental and theoretical fluxes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamics and its prediction and CALPHAD modeling: Review, state of the art, and perspectives

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable, or unstable, when interacting with its surroundings. The combined law of thermodynamics derived by Gibbs about 150 years ago laid the foundation of thermodynamics. In Gibbs combined law, the entropy production due to internal processes was not included, and the 2nd law was thus practically removed from the Gibbs combined law, so it is only applicable to systems under equilibrium, thus commonly termed as equilibrium or Gibbs thermodynamics. Gibbs further derived the classical statistical thermodynamics in terms of the probability of configurations in a system in the later 1800's and early 1900's. With the quantum mechanics (QM) developed in 1920's, the QM-based statistical thermodynamics was established and connected to classical statistical thermodynamics at the classical limit as shown by Landau in the 1940's. In 1960's the development of density functional theory (DFT) by Kohn and co-workers enabled the QM prediction of properties of the ground state of a system. On the other hand, the entropy production due to internal processes in non-equilibrium systems was studied separately by Onsager in 1930's and Prigogine and co-workers in the 1950's. In 1960's to 1970's the digitization of thermodynamics was developed by Kaufman in the framework of the CALculation of PHAse Diagrams (CALPHAD) modeling of individual phases with internal degrees of freedom. CALPHAD modeling of thermodynamics and atomic transport properties has enabled computational design of complex materials in the last 50 years. Our recently termed zentropy theory integrates DFT and statistical mechanics through the replacement of the internal energy of each individual configuration by its DFT-predicted free energy. The zentropy theory is capable of accurately predicting the free energy of individual phases, transition temperatures and properties of magnetic and ferroelectric materials with free energies of individual configurations solely from DFT-based calculations and without fitting parameters, and is being tested for other phenomena including superconductivity, quantum criticality, and black holes. Those predictions include the singularity at critical points with divergence of physical properties, negative thermal expansion, and the strongly correlated physics. Furthermore, those individual configurations may thus be considered as the genomic building blocks of individual phases in the spirit of the materials genome®. This has the potential to shift the paradigm of CALPHAD modeling from being heavily dependent on experimental inputs to becoming fully predictive with inputs solely from DFT-based calculations and machine learning models built on those calculations and existing experimental data through newly developed and future open-source tools. Furthermore, through the combined law of thermodynamics including the internal entropy production, it is shown that the kinetic coefficient matrix of independent internal processes is diagonal with respect to the conjugate potentials in the combined law, and the cross phenomena that the phenomenological Onsager flux and reciprocal relationships are due to the dependence of the conjugate potential of a molar quantity on nonconjugate molar quantities and other potentials, which can be predicted by the zentropy theory and CALPHAD modeling.

42 ENGINEERING↗

Contributions to the mechanistic understanding of the microstructural evolution in irradiated U-Mo dispersion fuel

Here, advanced microstructural characterization techniques, such as scanning electron microscopy (SEM) and scanning transmission electron microscopy - energy dispersive x-ray spectroscopy (STEM-EDS), were used to interpret the fuel microstructure evolution and fission products behavior in U-Mo dispersion fuel irradiated in the Advanced Test Reactor (ATR) as part of the European Mini-Plate Irradiation Experiment (EMPIrE) test. The larger as-fabricated fuel grain size achieved by heat-treating the U-Mo powder resulted in slower high burnup structure (HBS) development and reduced fission gas porosity. Slower HBS kinetics was observed at the fuel kernels’ periphery, which contained smaller and less fission gas bubbles at all fission densities (FDs) investigated and was attributed to a locally reduced damage density and fission products concentration, as corroborated with Monte Carlo simulations. The non-refined grains at the fuel kernel periphery hosted a perfectly ordered fission Gas Bubble Superlattice (GBS) up to 6.3 × 10 21 fissions/cm 3 . Nano-scale STEM-EDS analysis presented in this study provided useful information on the GBS characteristic morphology and evolution in U-Mo fuel. The concentration of fission gas in the GBS progressively increased with FD, pointing to an evolution of the nanobubble pressure status with irradiation. A possible connection between the GBS collapse and HBS onset is proposed for which there exists a threshold in the misorientation of the refined sub-grains above which the GBS stability during irradiation is no longer preserved, resulting in the GBS collapse.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Radial electric field and density fluctuations measured by Doppler reflectometry during the post-pellet enhanced confinement phase in W7-X

Radial profiles of density fluctuations and the radial electric field, Er, have been measured using Doppler reflectometry during the post-pellet enhanced confinement phase achieved, under different heating power levels and magnetic configurations, during the 2018 W7-X experimental campaign. A pronounced Er-well is measured with local values as high as -40 kV m -1 in the radial range ρ ~ 0.7–0.8 during the post-pellet enhanced confinement phase. The maximum Er intensity scales with both the plasma density and electron cyclotron heating power level, following a similar trend to the plasma energy content. A good agreement is found when the experimental Er profiles are compared to simulations carried out using the neoclassical codes, the drift kinetic equation solver (DKES) and kinetic orbit-averaging solver for stellarators (KNOSOS). The density fluctuation level decreases from the plasma edge toward the plasma core and the drop is more pronounced in the post-pellet enhanced confinement phase than in reference gas-fuelled plasmas. Besides, in the post-pellet phase, the density fluctuation level is lower in the high iota magnetic configuration than in the standard one. Therefore, to determine whether this difference is related to the differences in the plasma profiles or to the stability properties of the two configurations, gyrokinetic simulations have been carried out using the codes stella and EUTERPE. The simulation results point to the plasma profile evolution after the pellet injection and the stabilization effect of the radial electric field profile as the dominant players in the stabilization of the plasma turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Onset of Magnetic Reconnection in Tail-Like Equilibria

Magnetic reconnection is a fundamental mode of dynamics in the magnetotail, and is recognized as the basic mechanisms converting stored magnetic energy into kinetic energy of plasma particles. The effects of the reconnection process are well documented by spacecraft observations of plasmoids in the distant magnetotail, or bursty bulk flows, and magnetic field dipolarizations in the near Earth region. Theoretical and numerical analyses have, in recent years, shed new light on the way reconnection operates, and, in particular, which microscopic mechanism supports the dissipative electric field in the associated diffusion region. Despite this progress, however. the question of how magnetic reconnection initiates in a tail-like magnetic field with finite flux threading the current i.sheet remains unanswered. Instead, theoretical studies supported by numerical simulations support the point-of-view that such plasma and current sheets are stable with respect to collisionless tearing mode. In this paper, we will further investigate this conclusion, with emphasis on the question whether it remains valid in plasma sheets with embedded thin current sheets. For this purpose, we perform particle-in-cell simulations of the driven formation of thin current sheets, and their subsequent evolution either to equilibrium or to instability of a tearing-type mode. In the latter case we will pay particular attention to the nature of the electric field contribution which unmagnetizes the electrons.

Hesse, Michael↗

Time-dependent numerical modeling of dust halo formation at comets

The evolution of gas and dust distributions following a spatially and temporally localized comet outburst was calculated using a hybrid kinetic-hydrodynamic method. In the inner coma the time-dependent continuity, momentum, and energy equations of the dusty gas flow were solved simultaneously using 12 dust sizes. Beyond 300 km a three-dimensional kinetic model was used to calculate the trajectory of each individual dust grain. It was found that following the onset of the comet outburst a gas-dust blast wave propagates outward in the inner coma. About 15 minutes after the increased gas and dust production was initiated at the nucleus, a new equilibrium was reached in the inner coma. The most important feature of this new steady state was the significant increase of the dust terminal velocities. These higher terminal velocity values resulted in larger apex distances for dust particles emitted during the outburst. The dust particles spend a relatively long time near their apex points; therefore, the outburst generates long-lasting distinct dust envelopes in front of the regular dust coma.

Gombosi, T. I.↗

Assessment of two-temperature kinetic model for ionizing air

A two-temperature chemical-kinetic model for air is assessed by comparing theoretical results with existing experimental data obtained in shock-tubes, ballistic ranges, and flight experiments. In the model, named the TTv model, one temperature (T) is assumed to characterize the heavy-particle translational and molecular rotational energies, and another temperature (Tv) to characterize the molecular vibrational, electron translational, and electronic excitation energies. The theoretical results for nonequilibrium air flow in shock tubes are obtained using the computer code STRAP (Shock-Tube Radiation Program), and for flow along the stagnation streamline in the shock layer over spherical bodies using the newly developed code STRAP (Stagnation-Point Radiation Program). Substantial agreement is shown between the theoretical and experimental results for relaxation times and radiative heat fluxes. At very high temperatures the spectral calculations need further improvement. The present agreement provides strong evidence that the two-temperature model characterizes principal features of nonequilibrium air flow. New theoretical results using the model are presented for the radiative heat fluxes at the stagnation point of a 6-m-radius sphere, representing an aeroassisted orbital transfer vehicle, over a range of free-stream conditions. Assumptions, approximations, and limitations of the model are discussed.

Park, Chul↗

Consequences of impacts of small asteroids and comets with Earth

The fragmentation of a small asteroid in the atmosphere greatly increases its cross sections for aerodynamic braking and energy dissipation. At a typical impact velocity of 22 km/s, the atmosphere absorbs more than half the kinetic energy of stony meteoroids with diameters, D(sub m), less than 220 m and iron meteoroids with D(sub m) less than 80 m. The corresponding diameter for comets with impact velocity 50 km/s is D(sub m) less than 1600 m. Most of the atmospheric energy dissipation occurs in a fraction of a scale height, so large meteors appear to 'explode' or 'flare' at the end of their visible paths. This dissipation of energy in the atmosphere protects the earth from direct impact damage (e.g., craters), but it produces a blast wave that can do considerable damage. The area of destruction around the impact point in which the over-pressure in the blast wave exceeds 4 lb/sq in = 2.8 x 10(exp 5) dynes/cu cm, which is enough to knock over trees and destroy buildings, increases rapidly from zero for chondritic meteoroids less than 56 m in diameter (15 megatons) to about 200 sq km for those 80 m in diameter (48 megatons); the probable diameter of the tunguska impactor of 1908 is about 80 m. Crater formation and earthquakes are not significant in land impacts by stony asteroids less than about 200 m in diameter because of the air protection. A tsunami is probably the most devastating type of damage for asteroids 200 m to 1 km in diameter. An impact by an asteroid this size anywhere in the Atlantic would devastate coastal areas on both sides of the ocean. An asteroid a few kilometers across would produce a tsunami that would reach the foothills of the Appalachian Mountains in the upper half of the East Coast of the United States. Most of Florida is protected from a tsunami by the gradual slope of the ocean off its coast, which causes most of the tsunami energy to be reflected back into the Atlantic. The atmosphere plume produced by asteroids with diameters exceeding about 120 m cannot be contained by the atmosphere, so this bubble of high-temperature gas forms a new layer on top of the atmosphere. The dust entrapped in this hot gas is likely to have optical depths exceeding tau = 10 for asteroids with diameters exceeding about 0.5 to 1 km. The optical flux from asteroids 60 m or more in diameter is enough to ignite pine forests. However, the blast wave from an impacting asteroid goes beyond the radius in which the fire starts. The blast wave tends to blow out the fire, so it is likely that the impact will char the forest, as at Tunguska, but the impact will not produce a sustained fire. Because comets dissipate their energy much higher in the atmosphere than asteroids, they illuminate a much larger region and their blast wave is weaker. So they are much more effective in producing large fires. This suggests that the KT impactor was a comet rather than an asteroid.

Hills, J. G.↗

Fluid Creep Effects on Near-Wall Solute Transport for Non-Isothermal Ampoules

There is a growing practical and theoretical interest in developing accurate macroscopic modelling for flows arising in chemical or physical vapor transport (VT) crystal growth experiments, including those conducted in reduced gravity environments. Rosner was the first person to point out that previously neglected rarefield gas dynamics phenomena (Stefan and bouyancy-driven flows) become rather important sources of convection. In particular, the combination of rarefaction and strong gradients of temperature (and/or concentration) tangential to the side-walls of the ampoule induces convective flows known as thermal (and concentration) 'creep' respectively. His order-of-magnitude estimates revealed that thermal creep effects can be non-negligible even at normal gravitational levels. On the macroscopic level, the bulk fluid mechanics can be adequately described by the familiar macroscopic equations as long as the boundary conditions are modified to account for the integrated effect of kinetic boundary layers adjacent to solid boundaries. Motivated by the growing importance of these phenomena, we have embarked on a series of computational studies to elucidate these fundamental creep-induced effects for a rarefied gas in simple, two-dimensional confined geometries. However, unlike previous related studies, we resort to a microscopic description of the gas, mathematically expressed by the Boltzmann integro-differential equation. We employ the direct simulation Monte Carlo (DSMC) method of Bird, the theoretical foundations and several practical applications. In the case of thermally induced flows, the no-time counter method of Bird is used, as implemented for a hard-sphere gas. The scheme has been also extended to account for realistic molecular interaction models, an extension necessary if the diffusion physics underlying concentration creep are to be captured.

Papadopoulos, Dimitrios↗

Isolation of novel ribozymes that ligate AMP-activated RNA substrates

BACKGROUND: The protein enzymes RNA ligase and DNA ligase catalyze the ligation of nucleic acids via an adenosine-5'-5'-pyrophosphate 'capped' RNA or DNA intermediate. The activation of nucleic acid substrates by adenosine 5'-monophosphate (AMP) may be a vestige of 'RNA world' catalysis. AMP-activated ligation seems ideally suited for catalysis by ribozymes (RNA enzymes), because an RNA motif capable of tightly and specifically binding AMP has previously been isolated. RESULTS: We used in vitro selection and directed evolution to explore the ability of ribozymes to catalyze the template-directed ligation of AMP-activated RNAs. We subjected a pool of 10(15) RNA molecules, each consisting of long random sequences flanking a mutagenized adenosine triphosphate (ATP) aptamer, to ten rounds of in vitro selection, including three rounds involving mutagenic polymerase chain reaction. Selection was for the ligation of an oligonucleotide to the 5'-capped active pool RNA species. Many different ligase ribozymes were isolated; these ribozymes had rates of reaction up to 0.4 ligations per hour, corresponding to rate accelerations of approximately 5 x10(5) over the templated, but otherwise uncatalyzed, background reaction rate. Three characterized ribozymes catalyzed the formation of 3'-5'-phosphodiester bonds and were highly specific for activation by AMP at the ligation site. CONCLUSIONS: The existence of a new class of ligase ribozymes is consistent with the hypothesis that the unusual mechanism of the biological ligases resulted from a conservation of mechanism during an evolutionary replacement of a primordial ribozyme ligase by a more modern protein enzyme. The newly isolated ligase ribozymes may also provide a starting point for the isolation of ribozymes that catalyze the polymerization of AMP-activated oligonucleotides or mononucleotides, which might have been the prebiotic analogs of nucleoside triphosphates.

NASA Discipline Exobiology↗

Shear Strength and Release of Large Cryogenic Pellets from the Barrel of a Shattered Pellet Injector for Disruption Mitigation

Shattered pellet injection (SPI) has been chosen as the baseline disruption mitigation system on ITER due to its ability to rapidly inject material deep into the plasma to greatly increase the plasma density and radiate the thermal energy. SPI utilizes a mechanical punch or high-pressure gas to release and accelerate a pellet that has been cryogenically desublimated in the barrel of a pipe gun. Various material injection combinations could possibly be implemented during different phases of a disruption event to radiate plasma energy, reduce electromagnetic loads on machine components, avoid the formation of runaway electrons, or to dissipate runaway electrons that form. Each injection phase could possibly utilize combinations of deuterium, neon, or argon. In this paper we outline experimental measurements of pellet material shear strength at SPI operating temperatures to understand the force needed to release SPI pellets. Deuterium, neon, argon, and deuterium-neon mixture pellets with diameters of 8.5, 12.5, and 15.7 mm are formed at a range of relevant gas pressures and temperatures and dislodged from the cold zone with a slow-moving piston driven by a motor. The slow-moving piston is kept above the triple point temperature of the material while the pellet is forming, then cooled to below the triple point temperature before contacting the pellet to minimize any thermal conduction to the pellet. The piston incorporates a load cell to measure the force applied when the pellet breaks away from the cold zone in the barrel.The ability of the gas and punch methods to exceed the shear strength of the studied pellet materials for release has been analyzed. High pressure gas delivered by fast opening valves produce pressure shock to the pellet due to supersonic expansion of the propellant gas. Pressure (and therefore, force) oscillations are present due to transverse density propagation throughout the breech volume. Mechanical punches deliver an impact force through a high-kinetic energy impact. The effect of the mechanical shock on the pellet has been explored and will be presented in this paper. Scaling to larger ITER size SPI pellets will be described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Stability Analysis of Parameter Varying Genetic Toggle Switches Using Koopman Operators

The genetic toggle switch is a well known model in synthetic biology that represents the dynamic interactions between two genes that repress each other. The mathematical models for the genetic toggle switch that currently exist have been useful in describing circuit dynamics in rapidly dividing cells, assuming fixed or time-invariant kinetic rates. There is a growing interest in being able to model and extend synthetic biological function for growth conditions such as stationary phase or during nutrient starvation. As cells transition from one growth phase to another, kinetic rates become time-varying parameters. In this paper, we propose a novel class of parameter varying nonlinear models that can be used to describe the dynamics of genetic circuits, including the toggle switch, as they transition from different phases of growth. We show that there exists unique solutions for this class of systems, as well as for a class of systems that incorporates the microbial phenomena of quorum sensing. Further, we show that the domain of these systems, which is the positive orthant, is positively invariant. We also showcase a theoretical control strategy for these systems that would grant asymptotic monostability of a desired fixed point. We then take the general form of these systems and analyze their stability properties through the framework of time-varying Koopman operator theory. A necessary condition for asymptotic stability is also provided as well as a sufficient condition for instability. A Koopman control strategy for the system is also proposed, as well as an analogous discrete time-varying Koopman framework for applications with regularly sampled measurements.

Harrison, Jamiree (ORCID:0000000265351579)↗