Search NASA⌕ Search

SEARCH · Search NASA

Results for “reference cell model”

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 307 records · Page 17

Contamination Control for Ultra-Sensitive Life-Detection Missions

"This study report summarizes technological developments in science-required contamination control for cost-capped, life-detection missions in the Solar System. Technology advances focus on developing implementable strategies to restrict contamination of the spacecraft and instruments by the launch vehicle hardware. Secondary bake out operations during cruise further reduce molecular contaminants down to femtomolar levels in the sample path. The study validated a full-spacecraft, deployable barrier design as an effective strategy to isolate the spacecraft from the pre-launch processing and launch environment (e.g., fairing acoustic materials). It accommodates late mounting of RTGs and other required prelaunch activities. A new, high-fidelity physics, contamination-transport model for particles (including cells and their parts) and science-relevant molecules (e.g., possible biomolecules) was developed to validate the barrier concept. Importantly, this model takes into account the physics of extremely clean surfaces achieved by the best traditional contamination engineering. The ability to deal with very low levels of particulate and molecular contamination makes this model unique among standard contamination modeling approaches. Modeling also demonstrated that a high-temperature bake out during interplanetary cruise of the sample collector that comprises the largest surface area in the sample path for an Enceladus orbiter reference mission, was extremely effective. This study included the evaluation of other details necessary for an effective contamination control strategy for life detection missions including feasibility of achieving parts per trillion level cleanliness for semi-volatile hydrocarbons in nitrogen purge gas, identification of launch vehicle services required for the baseline contamination control, and defining additional contamination control measures to be taken during launch operations. Key results of the study are: • New spacecraft barrier design that accommodates RTGs is readily cleanable and repairable. Deployment of a 1/3 scale model was successfully demonstrated. • The barrier reduces particle contamination (likely biological) from fairing to spacecraft by 1e-2 to 1e-3. • On-cruise bake-out of critical surfaces significantly reduced molecular contamination (by as much as 1e-12). Subsequently, the probability of a surface contamination particle being transported to an instrument by an ice particle is less than 5.1e-5 (for microbes specifically, 4.39e-10). The report concludes that a full-spacecraft barrier designed specifically for life-detection missions that require stringent control of particulate and molecular contaminates is an effective means of mitigating risks of false positive and negative results for life signature investigations. Further, secondary cleaning steps for critical sections of the sample path can be highly effective at reducing molecular contaminants to low femtomolar levels."

contamination, life detection, barrier, modeling↗

Transcriptional networks underlying a primary ovarian insufficiency disorder in alligators naturally exposed to EDCs

Interactions between the endocrine system and environmental contaminants are responsible for impairing reproductive development and function. Despite the taxonomic diversity of affected species and attendant complexity inherent to natural systems, the underlying signaling pathways and cellular consequences are mostly studied in lab models. To resolve the genetic and endocrine pathways that mediate affected ovarian function in organisms exposed to endocrine disrupting contaminants in their natural environments, we assessed broad-scale transcriptional and steroidogenic responses to exogenous gonadotropin stimulation in juvenile alligators (Alligator missippiensis) originating from a lake with well-documented pollution (Lake Apopka, FL) and a nearby reference site (Lake Woodruff, FL). Here, we found that individuals from Lake Apopka are characterized by hyperandrogenism and display hyper-sensitive transcriptional responses to gonadotropin stimulation when compared to individuals from Lake Woodruff. Site-specific transcriptomic divergence appears to be driven by wholly distinct subsets of transcriptional regulators, indicating alterations to fundamental genetic pathways governing ovarian function. Consistent with broad-scale transcriptional differences, ovaries of Lake Apopka alligators displayed impediments to folliculogenesis, with larger germinal beds and decreased numbers of late-stage follicles. After resolving the ovarian transcriptome into clusters of co-expressed genes, most site-associated modules were correlated to ovarian follicule phenotypes across individuals. However, expression of two site-specific clusters were independent of ovarian cellular architecture and are hypothesized to represent alterations to cell-autonomous transcriptional programs. Collectively, our findings provide high resolution mapping of transcriptional patterns to specific reproductive function and advance our mechanistic understanding regarding impaired reproductive health in an established model of environmental endocrine disruption.

59 BASIC BIOLOGICAL SCIENCES↗

Exotic edge states of C 3 high-fold fermions in honeycomb lattices

A generalization of the graphene honeycomb model to the case where each site in the honeycomb lattice contains a n -fold degenerate set of eigenstates of the C 3 symmetry has been recently proposed to describe several systems, including triangulene crystals and photonic lattices. These generalized honeycomb models are defined by ( n a , n b ) , the number of C 3 eigenstates in the a and b sites of the unit cell, resulting in n a + n b bands. Thus, the (1,1) case gives the coventional honeycomb model that describes the two low-energy bands in graphene. Generalizations, such as (2,1), (2,2), and (3,3) display several nontrivial features, such as coexisting graphenelike Dirac cones with flat bands, both at zero and finite energy, as well as robust degeneracy points where a flat band and a parabolic band meet at the Γ point. Here we explore the edge states of this class of crystals, using as reference triangulene crystals, and we find several types of edge states absent in the conventional (1,1) honeycomb case, associated to the nontrivial features of the two-dimensional bands of the high-fold case. First, we find dispersive edge states associated to the finite-energy flat bands, that occur both at the armchair and zigzag termination. Second, in the case of noncentrosymmetric triangulene crystals that lead to a S = 1 Dirac band, we have a bonding-antibonding pair of dispersive edge states, localized in the same edge so that their energy splitting is reduced as their localization increases, opposite to the conventional behavior of pairs of states localized in opposite edges. Third, for the (3,3) case, that hosts a gap separating a pair of flat conduction and valence bands, we find nondispersive edge states with E = 0 in all edge terminations. Published by the American Physical Society 2024

Madail, L. (ORCID:0000000208457748)↗

Pressure-Volume-Temperature (PVT) Gauging of an Isothermal Cryogenic Propellant Tank Pressurized with Gaseous Helium

Results are presented for pressure-volume-temperature (PVT) gauging of a liquid oxygen/liquid nitrogen tank pressurized with gaseous helium that was supplied by a high-pressure cryogenic tank simulating a cold helium supply bottle on a spacecraft. The fluid inside the test tank was kept isothermal by frequent operation of a liquid circulation pump and spray system, and the propellant tank was suspended from load cells to obtain a high-accuracy reference standard for the gauging measurements. Liquid quantity gauging errors of less than 2 percent of the tank volume were obtained when quasi-steady-state conditions existed in the propellant and helium supply tanks. Accurate gauging required careful attention to, and corrections for, second-order effects of helium solubility in the liquid propellant plus differences in the propellant/helium composition and temperature in the various plumbing lines attached to the tanks. On the basis of results from a helium solubility test, a model was developed to predict the amount of helium dissolved in the liquid as a function of cumulative pump operation time. Use of this model allowed correction of the basic PVT gauging calculations and attainment of the reported gauging accuracy. This helium solubility model is system specific, but it may be adaptable to other hardware systems.

VanDresar, Neil T.↗

Charged-cell periodic DFT simulations via an impurity model based on density embedding: Application to the ionization potential of liquid water

Calculations of charged systems in periodic boundary conditions (PBC) are problematic because there are spurious interactions between the charges in different periodic images that can affect the physical picture. In addition, the intuitive limit of Coulomb interactions decaying to zero as the interacting charges are placed at infinite separation no longer applies, and for example total energies become undefined. Leveraging subsystem density functional theory (also known as density embedding) we define an impurity model that embeds a finite neutral or charged subsystem within an extended (infinite) surrounding subsystem. The combination of the impurity model and a consistent choice of the Coulomb reference provides us with an algorithm for evaluating the ionization potential (IP) in extended systems. Additionally, we demonstrate our approach in a pilot calculation of the IP of liquid water, based on a configuration from a prior ab initio molecular dynamics (AIMD) simulation of liquid water (Genova et al., J. Chem. Phys. 2016, 144, 234105). The calculations with the impurity model capture the broadening on the ionization energies introduced by the interactions between the water molecules. Furthermore, the calculated average IP value (10.5 eV) compare favorably to experiments (9.9-10.06 eV) and very recent simulations based on the GW approximation (10.55 eV), while at the same time outperforming density embedding calculations carried out with a naïve handling of the electrostatic interactions (about 7 eV).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spatially Mapping the CO 2 Alkaline Sorbent Diffuse Microenvironment Using Operando Raman Spectroscopy

When designing a chemical process, the local balance of transport and kinetics, collectively referred to as the diffuse microenvironment, plays a critical role in performance but is difficult to directly observe. This work demonstrates a method of two-dimensional spatial chemical mapping of the diffuse microenvironment in the context of alkali metal hydroxide direct air capture of carbon dioxide using a custom operando gas-absorption flow cell along with confocal Raman spectroscopy. Notably, we observe the concentration boundary layer near the gas–liquid interface and elucidate the interplay of carbonate and bicarbonate ions within it while inferring local hydroxide depletion through continuum modeling. These first of their kind observations provide a technique to compare the performance of direct air capture solvents based on diffuse microenvironment dynamics while also providing metrics important for air contactor design such as boundary layer thickness. Overall, this work showcases a new experimental platform to study interfacial diffuse microenvironments in and outside of the field of direct air capture of carbon dioxide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Decoupling electrode kinetics to elucidate reaction mechanisms in alkaline water electrolysis

Alkaline water electrolysis (AWE) presents key advantages, including reduced material costs, enhanced operational stability, and compatibility with non-precious metal catalysts, positioning it as a scalable route for hydrogen production. In this study, we introduce a minimally invasive single-cell configuration incorporating a reference electrode via diaphragm extension to form an internal ion channel. This setup, combined with an interfaced potentiostat and auxiliary electrometer, enables real-time, independent monitoring of anode and cathode behavior, offering high-resolution electrochemical diagnostics. While it is well established that the hydrogen evolution reaction (HER) exhibits sluggish kinetics in alkaline media, our study reveals that this limitation persists even in practical AWE systems where nickel-based substrates are used as electrodes. This observation is supported by both experimental data and voltage breakdown modeling. Arrhenius-type analysis reveals that localized electric fields induced by catalysts shift the reaction kinetics from classical Butler–Volmer behavior toward a Marcus-like regime, where interfacial molecular dynamics and bimolecular charge transfer dominate. We propose a semi-empirical model and a surficial reaction mechanism to describe these dynamics. This work underscores the critical need for cathode innovation and provides a rational framework for designing advanced catalysts and electrode architectures to optimize AWE performance.

08 HYDROGEN↗

Langmuir Probe Measurements in an Inductively Coupled GEC Reference Cell Plasma

Measurements of electron number density, electron temperature, and electron energy distribution function (EEDF) using a compensated Langmuir probe have been performed on an inductively (transformer ) coupled Gaseous Electronics Conference (GEC) reference cell plasma. The plasma source is operated with CH4, CF4, or their mixtures with argon. The effect of independently driving the electrode supporting the wafer on the probe data is studied. In particular, we find that the plasma structure depends on the phase in addition to the magnitude of the power coupled to the electrode relative to that of the transformer coil. The Langmuir probe is translated in a plane parallel to the electrode to investigate the spatial structure of the plasma. The probe data is also compared with fluid model predictions.

Ji, J. S.↗

Numerical methods for instability mitigation in the modeling of laser wakefield accelerators in a Lorentz-boosted frame

Modeling of laser-plasma wakefield accelerators in an optimal frame of reference has been shown to produce orders of magnitude speed-up of calculations from first principles. Obtaining these speedups required mitigation of a high-frequency instability that otherwise limits effectiveness. In this paper, methods are presented which mitigated the observed instability, including an electromagnetic solver with tunable coefficients, its extension to accommodate Perfectly Matched Layers and Friedman's damping algorithms, as well as an efficient large bandwidth digital filter. It is observed that choosing the frame of the wake as the frame of reference allows for higher levels of filtering or damping than is possible in other frames for the same accuracy. Detailed testing also revealed the existence of a singular time step at which the instability level is minimized, independently of numerical dispersion. A combination of the techniques presented in this paper prove to be very efficient at controlling the instability, allowing for efficient direct modeling of 10 GeV class laser plasma accelerator stages. The methods developed in this paper may have broader application, to other Lorentz-boosted simulations and Particle-In-Cell simulations in general.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of homogenized and enhanced diffusion solutions of model PWR problems

Model problem comparisons in slab geometry are made between two forms of homogenized diffusion theory and enhanced diffusion theory. The pin-cell discontinuity factors for homogenized diffusion calculations are derived from homogenized variational nodal P1 response matrices and from standard finite differencing. Enhanced diffusion theory consists of applying quasi-reflected interface conditions to reduce variational nodal Pn response matrices to one degree of freedom per interface, without homogenization within the cell. As expected both homogenized diffusion methods preserve reaction rates exactly if the discontinuity factors are derived from the P 11 reference solutions. If no reference lattice solution is available, discontinuity factors may be approximated from single cells with reflected boundary conditions; the computational effort is then comparable to calculating the enhanced diffusion response matrices. In this situation enhanced diffusion theory gives the most accurate results and finite difference discontinuity factors the least accurate. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Active Reconditioning of Retired Lithium-ion Battery Packs from Electric Vehicles for Second Life Applications

Utilizing the remaining capacity in retired lithium-ion (Li-ion) batteries from electric vehicles (EVs) for second-life applications has shown economic and environmental benefits. However, achieving homogeneity among the capacities of cells before their second life is critical to exploit the benefits. This article proposes a new active reconditioning approach with the potential to make short-term reconditioning of batteries before second life feasible. A control objective map determines each cell's state-of-charge (SOC) operating window based on its capacity relative to other cells. The SOC reference translates into distinct differential currents through the cells, wherein the higher-capacity cells undergo more frequent and deep charge and discharge cycles than their lower-capacity counterparts. The proposed solution achieves capacity homogeneity within the battery pack with low reconditioning time and minimal fade in the overall pack capacity. The feasibility of the reconditioning approach under varying load and environmental conditions is assessed through simulations, encompassing factors such as the number of cycles per day, depth-of-discharge, battery pack temperature, and cell resting time at different SOCs. Furthermore, the simulation model employs a battery pack with sixteen series-connected 75 Ah Kokam lithium nickel manganese cobalt oxide (NMC) cells with a 3.6% initial capacity imbalance. A reconditioning time of 1.3 months is achieved with a final capacity imbalance of 0.1% and an overall capacity fade of 0.005%, thereby confirming the viability of the reconditioning process. Moreover, experimental validation using eight retired battery cells from a Nissan Leaf demonstrates a substantial decrease in the capacity imbalance of cells from 9.4% to 2.15% within 78 days, effectively affirming the efficacy of the proposed reconditioning scheme.

25 ENERGY STORAGE↗

Spindown of Pulsars Interacting with Companion Winds: The Impact of Magnetospheric Compression

Pulsars in binary systems with strong companion winds can have the magnetopause separating their magnetosphere from the wind located well within their light cylinder. This bow-like enclosure effectively creates a waveguide that confines the pulsar's electromagnetic fields and can significantly alter its spindown. In this paper, we study the spindown of compressed pulsar magnetospheres in such systems. We parameterize the confinement as the ratio between the equatorial position of the magnetopause (or standoff distance) R m and the pulsar's light cylinder R LC . Using particle-in-cell simulations, we quantify the pulsar spindown for a range of compressions, R m /R LC = 1/3–1, and inclination angles, χ = 0°...90°, between magnetic and rotation axes. Our strongly confined models (R m /R LC = 1/3) show two distinct limits. For χ = 0°, the spindown of a compressed pulsar magnetosphere is enhanced by approximately a factor of three compared to an isolated pulsar due to the increased number of open magnetic field lines. Conversely, for χ = 90°, the compressed pulsar spins down at less than 40% of the rate of an isolated reference pulsar due to the mismatch between the pulsar wind stripe wavelength and the waveguide size. We apply our analysis to the 2.77 s oblique rotator (χ = 60°) in the double-pulsar system PSR J0737-3039. With the numerically derived spindown estimate, we constrain its surface magnetic field to B* ≈ (7.3 ± 0.2) × 10 11 G. We discuss the time modulation of its period derivative, the effects of compression on its braking index, and implications for the radio eclipse in PSR J0737-3039.

79 ASTRONOMY AND ASTROPHYSICS↗

Cell biology experiments conducted in space

A review of cell biology experiments conducted during the first two decades of space flight is provided. References are tabulated for work done with six types of living test system: isolated viruses, bacteriophage-host, bacteria, yeasts and filamentous fungi, protozoans, and small groups of cells (such as hamster cell tissue and fertilized frog eggs). The general results of studies involving the survival of cells in space, the effect of space flight on growing cultures, the biological effects of multicharged high-energy particles, and the effects of space flight on the genetic apparatus of microorganisms are summarized. It is concluded that cell systems remain sufficiently stable during space flight to permit experimentation with models requiring a fixed cell line during the space shuttle era.

Taylor, G. R.↗

Investigation of electrolyte measurement in diluted whole blood using spectroscopic and chemometric methods

The feasibility of using near-infrared (NIR) spectroscopy in combination with partial least-squares (PLS) regression was explored to measure electrolyte concentration in whole blood samples. Spectra were collected from diluted blood samples containing randomized, clinically relevant concentrations of Na+, K+, and Ca2+. Sodium was also studied in lysed blood. Reference measurements were made from the same samples using a standard clinical chemistry instrument. Partial least squares (PLS) was used to develop calibration models for each ion with acceptable results (Na+, R2 = 0.86, CVSEP = 9.5 mmol/L; K+, R2 = 0.54, CVSEP = 1.4 mmol/L; Ca2+, R2 = 0.56, CVSEP = 0.18 mmol/L). Slightly improved results were obtained using a narrower wavelength region (470-925 nm) where hemoglobin, but not water, absorbed indicating that ionic interaction with hemoglobin is as effective as water in causing measurable spectral variation. Good models were also achieved for sodium in lysed blood, illustrating that cell swelling, which is correlated with sodium concentration, is not required for calibration model development.

Non-NASA Center↗

Design and application of a kinetic model of lipid metabolism in Saccharomyces cerevisiae

Lipid biosynthesis plays a vital role in living cells and has been increasingly engineered to overproduce various lipid-based chemicals. However, owing to the tightly constrained and interconnected nature of lipid biosynthesis, both understanding and engineering of lipid metabolism remain challenging, even with the help of mathematical models. Here we report the development of a kinetic metabolic model of lipid metabolism in Saccharomyces cerevisiae that integrates fatty acid biosynthesis, glycerophospholipid metabolism, sphingolipid metabolism, storage lipids, lumped sterol synthesis, and the synthesis and transport of relevant target-chemicals, such as fatty acids and fatty alcohols. The model was trained on lipidomic data of a reference S. cerevisiae strain, single knockout mutants, and lipid overproduction strains reported in literature. The model was used to design mutants for fatty alcohol overproduction and the lipidomic analysis of the resultant mutant strains coupled with model-guided hypothesis led to discovery of a futile cycle in the triacylglycerol biosynthesis pathway. In addition, the model was used to explain successful and unsuccessful mutant designs in metabolic engineering literature. Thus, this kinetic model of lipid metabolism can not only enable the discovery of new phenomenon in lipid metabolism but also the engineering of mutant strains for overproduction of lipids.

59 BASIC BIOLOGICAL SCIENCES↗

Damage to cellular DNA from particulate radiations, the efficacy of its processing and the radiosensitivity of mammalian cells. Emphasis on DNA double strand breaks and chromatin breaks

For several years, it has been evident that cellular radiation biology is in a necessary period of consolidation and transition (Lett 1987, 1990; Lett et al. 1986, 1987). Both changes are moving apace, and have been stimulated by studies with heavy charged particles. From the standpoint of radiation chemistry, there is now a consensus of opinion that the DNA hydration shell must be distinguished from bulk water in the cell nucleus and treated as an integral part of DNA (chromatin) (Lett 1987). Concomitantly, sentiment is strengthening for the abandonment of the classical notions of "direct" and "indirect" action (Fielden and O'Neill 1991; O'Neill 1991; O'Neill et al. 1991; Schulte-Frohlinde and Bothe 1991 and references therein). A layer of water molecules outside, or in the outer edge of, the DNA (chromatin) hydration shell influences cellular radiosensitivity in ways not fully understood. Charge and energy transfer processes facilitated by, or involving, DNA hydration must be considered in rigorous theories of radiation action on cells. The induction and processing of double stand breaks (DSBs) in DNA (chromatin) seem to be the predominant determinants of the radiotoxicity of normally radioresistant mammalian cells, the survival curves of which reflect the patterns of damage induced and the damage present after processing ceases, and can be modelled in formal terms by the use of reaction (enzyme) kinetics. Incongruities such as sublethal damage are neither scientifically sound nor relevant to cellular radiation biology (Calkins 1991; Lett 1990; Lett et al. 1987a). Increases in linear energy transfer (LET infinity) up to 100-200 keV micron-1 cause increases in the extents of neighboring chemical and physical damage in DNA denoted by the general term DSB. Those changes are accompanied by decreasing abilities of cells normally radioresistant to sparsely ionizing radiations to process DSBs in DNA and chromatin and to recover from radiation exposure, so they make significant contributions to the relative biological effectiveness (RBE) of a given radiation.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗

Orbiter Capability for Providing Water to the International Space Station according to the Most Probable Flight Attitudes

Water to be generated by, delivered to, and processed by the International Space Station (ISS) is a critical Environmental Control and Life Support (ECLS) element, especially for the early ISS missions. A significant portion of the water required by the ISS shall be provided by the Shuttle Transportation System (STS) Orbiter. The balance of water generated by the Orbiter Fuel Cells (FC), minus that water consumed by the Orbiter Flash Evaporator System (FES) and crew, is available for transfer to the ISS. During later missions, crew respired and perspired water, as well as effluent water from the Orbiter LiOH canisters, will be collected as condensate and available for transfer to the ISS. Orbiter radiator performance provides the most variance in determining the amount of net Orbiter water available for transfer to the ISS. As radiator performance decreases, the dependence upon the FES (and FC water) increases for rejecting Orbiter waste heat. Generally, radiator performance decreases as the ISS assembly size increases (especially as solar arrays are added), and also as beta angle increases. ISS solar array deployment necessitates the use of models with articulating solar arrays (for Earth local-vertical attitudes), as array position dramatically affects Orbiter radiator performance. Recent developments in the relaxation of beta angle limitations have also increased the complexity and difficulty of providing water to the ISS. Other factors that may hinder the ability to transfer water are the number of empty Contingency Water Containers (CWCs) available, duration of open-hatch time, crew activity timeline, and full CWC storage capability. A parametric study has been accomplished that provides a quick-reference table for determining expected water generation rates for ISS missions 2A.2 through 7A.1. An hourly Orbiter water generation rate is reported according to a matrix that consists of: (1) (six) significant changes in ISS assembly configuration; (2) (four) beta angles (0 deg. , +37 deg., +53 deg. , and +75 deg.); (3) the (three) most representative ISS attitudes (XPOP-O, XPOP-180 and +XVV); (4) (four) Orbiter radiator configurations (both stowed, starboard deployed, port deployed, and both deployed) and (5) the (two) conditions (radiator inlet temperatures and fuel cell power) most consistent with sleep and wake periods. Those permutations of higher probability of occurrence than others have been identified. Another parametric study has been accomplished that provides a quick-reference table for determining expected water generation rates for ISS assembly complete missions. An hourly Orbiter water generation rate is reported according to a matrix that consists of: (1) (seven) beta angles (-75 deg., -60 deg., -30 deg., 0 deg., +30 deg., +60 deg., and +75 deg.); (2) the (nine) PYR angles that define the corners of the envelope; (3) (four) Orbiter radiator configurations (both stowed, starboard deployed, port deployed, and both deployed) and (4) the (two) conditions (radiator inlet temperatures and fuel cell power) most consistent with sleep and wake periods.

Dunaway, Brian↗

Locally adaptive artificial viscosity strategies for Lagrangian hydrodynamics

To accurately model inviscid flow with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity is introduced to convert kinetic energy into internal energy, thereby providing a mechanism to generate the required entropy increase across shocks. In this paper, we propose a new method for constructing an adaptive, artificial viscosity in the context of one-dimensional, staggered-grid Lagrangian hydrodynamics. Our adaptive, artificial viscosity is defined in terms of two parameters that depend on density locally in a neighborhood around a shock, and hence, vary cell-by-cell. Our methodology is based on building a reference set of pre-computed optimal, globally constant artificial viscosity coefficients for a family of isolated shock test problems. For arbitrary flows, the evaluation of the unknown coefficients is automated by first estimating shock intensity locally, and second computing the corresponding adaptive parameter value via interpolation over the data from the pre-computed reference set of optimal, globally constant coefficient values. To illustrate the performance of our new approach, we compare our results against two existing methods. The first method is a limiter-based approach, which relies on estimating velocity gradients of the flow, and the second method is utilized in several commercial codes. Finally, we demonstrate that our new adaptive methodology produces more accurate results for a variety of tests with propagating shock waves, as well as for the aforementioned family of isolated shock problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗