Search NASA⌕ Search

SEARCH · Search NASA

Results for “Widely Used”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Introducing the GeoRePORT Resource Size Tool: Reporting on Geothermal Resource Size Estimations Using the Geothermal Resource Portfolio Optimization and Reporting Technique (GeoRePORT): Preprint

The Geothermal Resource Portfolio Optimization and Reporting Technique (GeoRePORT) was developed with funding from the U.S. Department of Energy Geothermal Technologies Office to assist in identifying and pursuing long-term investment strategies through the development of a resource reporting protocol. The assessment protocols used in GeoRePORT allow for comparison of project attributes across locations and geological settings to understand the feasibility of geothermal development. This work introduces the Resource Size Tool, a new feature within the GeoRePORT package that compiles two independent methods for estimating geothermal resource size in terms of energy capacity in MW. Energy production potential for twenty-three case studies was estimated with the Resource Size Tool in order to 1) generate a reasonable range of resource size estimates for a particular geothermal field; 2) illustrate the advantages and limitations of each methodology (such as data input requirements, estimate accuracy and precision, and the appropriate circumstances of use); and 3) test the ability of the resource size tool to provide useful and accurate information for geothermal stakeholders. The tool employs two methods widely used in the geothermal industry: (1) USGS Volumetric and (2) Power Density. Results from our case studies show general overlap between these two methods in terms of resource size estimates; however, they also reveal key differences between the two approaches that should be considered when using such estimates to drive development. First, the two methods rely on different input parameters and therefore one method may be more appropriate and/or accurate for a given project than the other. Second, the Power Density method was found to generate wider ranges of resource size predictions, more consistently aligning with actual power production of the field but with larger scales of error; whereas the USGS Volumetric method predicts narrower ranges but tends to overestimate when compared to current MW production. Future work will refine variables used in the methods with input data from other sections of GeoRePORT and modify uncertainty levels based on the particular datasets used for a given project.

geological↗
DOE OSTI2022

Meshfree Methods for Modeling Chemo-Mechanical Cathode Cracking in Li-Ion Batteries

Through repeated charging and discharging, the electrodes of a Li-ion battery experience cyclic loading causing swelling and contraction due to the movement of lithium, also known as intercalation. This research focuses specifically on the chemo-mechanical cracking in the cathode. Because cathode particles are comprised of many randomly-oriented grains, which have highly anisotropic material properties, the expansion and contraction is very non-uniform. As a result, stress concentrations tend to form between grains, which necessitates the modeling of crack propagation largely along grain boundaries. Chemo-mechanical models are simulated mainly by two methods: the cohesive zone model (CZM) and the continuous damage model (CDM). The CZM is more accurate at capturing the sharp discontinuities of a crack but is very computationally expensive and intractable for large-scale models as a result. Conversely, the CDM is easily computed but not well-suited to easily allow for discontinuous field variables, which are inherent across a crack. This study aims to improve the CDM's ability to capture discontinuous cracks. Current versions of the CDM use the finite element method (FEM), which is one of the most widely used approaches for spatial discretization. This research investigates the use of the Reproducing Kernel Particle Method (RKPM), a meshfree method, for spatial discretization and aims to achieve a chemo-mechanical crack propagation model that achieves both intermediate accuracy and computational time between a CDM and CZM. To gain a better understanding of how the meshfree model compares to the current FEM model in capturing cathode crack propagation, further investigation is needed.

cathode cracking↗
DOE OSTI2023

Meshfree Multiphysics Damage Modeling of Li-ion Battery Materials

Through repeated charging and discharging cycles, the electrodes of a Li-ion battery experience significant swelling and contraction due to the movement of lithium, also known as intercalation. This research focuses specifically on the chemo-mechanical cracking in the cathode. Because cathode particles are comprised of many randomly-oriented grains, which have highly anisotropic material properties, the expansion and contraction is very non-uniform. As a result, stress concentrations tend to form between grains, which necessitates the modeling of crack propagation largely along grain boundaries and material interfaces. Chemo-mechanical damage models are generally simulated by one of two methods: the cohesive zone model (CZM) and the continuous damage model (CDM). The CZM is more accurate at capturing the sharp discontinuities of a crack but is very computationally expensive and intractable for large-scale models as a result. Conversely, the CDM is easily computed but not well-suited to easily allow for discontinuous field variables, which are inherent across a crack. This study aims to improve the CDM's ability to capture discontinuous cracks. Current versions of the CDM use the finite element method (FEM), which is one of the most widely used approaches for spatial discretization. This research investigates the use of the Reproducing Kernel Particle Method (RKPM), a meshfree method, for spatial discretization and aims to achieve a coupled chemo-mechanical crack propagation model that enhances accuracy while maintaining high computational efficiency. This model features a fully coupled, iterative electrochemistry solution, which informs the meshfree damage model of impending crack formation. To gain a better understanding of how the meshfree model compares to the current FEM model in capturing cathode crack propagation, further investigation is needed.

damage modeling↗
DOE OSTI2025

Tree-level carbon stock estimations across diverse species using multi-source remote sensing integration

Forests are critical carbon sinks, and remote sensing has been increasingly widely used for forest monitoring and biomass estimations. However, species-specific tree-level studies remain limited. In this study, we demonstrated the feasibility of integrating UAV-based LiDAR with high-resolution optical satellite imagery (0.5 m) to estimate biomass for individual trees across different species. The proposed method accurately estimated biomass for 53 trees (R² = 0.82, rRMSE = 0.44), with species-specific datasets, showing an average 25.2% increase in R² and a 14.8% reduction in rRMSE. A novel vegetation index combining forest structure parameters with vegetation indices (VIs) was developed using high-resolution multispectral satellite data (3 m) to explore its relationship with individual tree biomass. Combining forest structural parameters with VIs further improved estimation accuracy, achieving an R²of 0.89 and an rRMSE of 0.34. Species-specific datasets show an 11.6% increase in R²compared to methods without VIs, and a 22.2% improvement over methods using only VIs. SHapley Additive exPlanations (SHAP) analysis shows that the volume feature played a key role in model performance and remained stable throughout the training process. Altogether, the proposed approach enhances individual tree biomass and carbon sink estimations, showing great potential for large-scale precise forest carbon monitoring using multi-source remote sensing data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Synthesis of PHI-type zeolite harmotome using barium as an inorganic structure-directing agent inspired from natural mineral compositions

Zeolites are nanoporous aluminosilicates widely used as commercial adsorbents, heterogeneous catalysts, and ion-exchange materials due to their unique porosity, acidity, and thermal stability. The vast majority of synthetically-realized zeolite structures are prepared under hydrothermal conditions in alkaline media, and often in the presence of an organic structure-directing agent that facilitates the crystallization of diverse porous networks. Due to economic and environmental disadvantages of organic-based syntheses, it is often desirable to produce zeolites in organic-free media using alkali metal ions as the most commonly employed inorganic structure-directing agents (ISDAs). Here, in this study, we were inspired by the elemental compositions of natural zeolite minerals to examine the use of an alkaline earth metal ion as the ISDA to prepare zeolites with the PHI framework type. Two natural minerals exhibiting this crystal structure are phillipsite and harmotome. Both structures contain sodium and potassium as extra-framework cations, but only harmotome contains barium cations. Here we examined how the combination of alkali and alkaline earth ISDAs promote the formation of PHI-type zeolites. Our findings reveal that each isostructure is prepared under different synthesis conditions where simply adding barium to a growth mixture does not always lead to harmotome formation. Harmotome is formed in more dilute and alkaline growth media; however, time-resolved analyses of zeolite crystallization reveal that harmotome can be synthesized at low temperatures (ca. 65 °C) with rates that are 5-times faster than syntheses of phillipsite. This seemingly indicates that barium reduces the energetic barrier for PHI-type zeolite nucleation in ways that are not fully understood. Moreover, harmotome crystals exhibit an average particle size of ca. 60 nm, which is nearly 40-fold smaller than phillipsite crystals prepared under similar conditions. Collectively, these findings identify routes to control crystallization of PHI-type zeolites and may provide a broader strategy for using the composition of natural zeolite minerals to optimize the preparation of synthetic analogues.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Using active matter to introduce spatial heterogeneity to the susceptible infected recovered model of epidemic spreading

Abstract The widely used susceptible-infected-recovered (S-I-R) epidemic model assumes a uniform, well-mixed population, and incorporation of spatial heterogeneities remains a major challenge. Understanding failures of the mixing assumption is important for designing effective disease mitigation approaches. We combine a run-and-tumble self-propelled active matter system with an S-I-R model to capture the effects of spatial disorder. Working in the motility-induced phase separation regime both with and without quenched disorder, we find two epidemic regimes. For low transmissibility, quenched disorder lowers the frequency of epidemics and increases their average duration. For high transmissibility, the epidemic spreads as a front and the epidemic curves are less sensitive to quenched disorder; however, within this regime it is possible for quenched disorder to enhance the contagion by creating regions of higher particle densities. We discuss how this system could be realized using artificial swimmers with mobile optical traps operated on a feedback loop.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Phase retrieval for refraction-enhanced x-ray radiography using a deep neural network

X-ray refraction-enhanced radiography (RER) or phase contrast imaging is widely used to study internal discontinuities within materials. The resulting radiograph captures both the decrease in intensity caused by material absorption along the x-ray path, as well as the phase shift, which is highly sensitive to gradients in density. A significant challenge lies in effectively analyzing the radiographs to decouple the intensity and phase information and accurately ascertain the density profile. Conventional algorithms often yield ambiguous and unrealistic results due to difficulties in including physical constraints and other relevant information. We have developed an algorithm that uses a deep neural network to address these issues and applied it to extract the detailed density profile from an experimental RER. To generalize the applicability of our algorithm, we have developed a technique that quantitatively evaluates the complexity of the phase retrieval process based on the characteristics of the sample and the configuration of the experiment. Accordingly, this evaluation aids in the selection of the neural network architecture for each specific case. Beyond RER, the model has potential applications for other diagnostics where phase retrieval analysis is required.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A miniature multi-anvil apparatus using diamond as anvils—MDAC: Multi-axis diamond anvil cell

The diamond anvil cell (DAC) has been widely used in high-pressure research. Despite significant progress over the past five decades, the opposed anvil geometry in the DAC inevitably leads to a disk-shaped sample configuration at high pressure. This intrinsic limitation is largely responsible for the large pressure and temperature gradients in the DAC, which often compromise precise experiments and their characterizations. Here, we designed and fabricated a multi-axis diamond anvil cell (MDAC) by adopting the concept of a multi-anvil apparatus but using single crystal diamonds as the anvil material. Preliminary data show that the MDAC can generate extreme pressure conditions above 100 GPa. The advantages of the MDAC over a traditional opposed anvil DAC include thicker, voluminous samples, quasi-hydrostatic, or designed deviatoric stress conditions, and multidirectional access windows for optical applications and x-ray probes. In this article, we present the design and performance of a prototype MDAC, as well as the application prospects in high-pressure research.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Sub-millisecond keyhole pore detection in laser powder bed fusion using sound and light sensors and machine learning

Laser powder bed fusion is a mainstream additive manufacturing technology widely used to manufacture complex parts in prominent sectors, including aerospace, biomedical, and automotive industries. However, during the printing process, the presence of an unstable vapor depression can lead to a type of defect called keyhole porosity, which is detrimental to the part quality. In this study, we developed an effective approach to locally detect the generation of keyhole pores during the printing process by leveraging machine learning and a suite of optical and acoustic sensors. Simultaneous synchrotron x-ray imaging allows the direct visualization of pore generation events inside the sample, offering high-fidelity ground truth. A neural network model adopting SqueezeNet architecture using single-sensor data was developed to evaluate the fidelity of each sensor for capturing keyhole pore generation events. Our comparative study shows that the near infrared images gave the highest prediction accuracy, followed by 100 kHz and 20 kHz microphones, and the photodiode sensitive to processing laser wavelength had the lowest accuracy. Using a single sensor, over 90% prediction accuracy can be achieved with a temporal resolution as short as 0.1 ms. A data fusion scheme was also developed with features extracted using SqueezeNet neural network architecture and classification using different machine learning algorithms. Our work demonstrates the correlation between the characteristic optical and acoustic emissions and the keyhole oscillation behavior, and thereby provides strong physics support for the machine learning approach.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Improving the National Solar Radiation Database (NSRDB) Using a Physics-Based Direct Normal Irradiance (DNI) Model

The National Solar Radiation Database (NSRDB) is a widely used resource providing satellite-derived solar data across the United States and globally. While the NSRDB employs a physical model for computing global horizontal irradiance (GHI), its current method for estimating cloudy-sky direct normal irradiance (DNI) relies on surface observations and empirical models. Recently, a novel physics-based approach, the Fast All-Sky Radiation Model for Solar applications with DNI (FARMS-DNI), was developed to enhance the DNI forecasting. FARMS-DNI incorporates both direct and scattered solar radiation within the circumsolar region, resulting in improved day-ahead DNI predictions when integrated into the Weather Research and Forecasting model with Solar extensions (WRF-Solar). This study integrates FARMS-DNI into the NSRDB algorithm to generate high-resolution DNI data from satellite resources. Our findings reveal that FARMS-DNI effectively mitigates the substantial DNI overestimation present in the conventional NSRDB across surface sites, particularly in conditions categorized as cloudy overcast. Consequently, this innovative model substantially enhances the overall accuracy of the NSRDB.

Xie, Yu↗
NASA NTRS1993

Bone geometry, structure and mineral distribution using Dual energy X ray Absorptiometry (DXA)

Dual energy x-ray absorptiometry (DXA) is currently the most widely used method of analyzing regional and whole body changes in bone mineral content (BMC) and areal (g/sq cm) bone mineral density (BMD). However, BMC and BMD do not provide direct measures of long bone geometry, structure, or strength nor do regional measurements detect localized changes in other regions of the same bone. The capabilities of DXA can be enhanced significantly by special processing of pixel BMC data which yields cross-sectional geometric and structural information. We have extended this method of analysis in order to develop non-uniform structural beam models of long bones.

Whalen, Robert↗
NASA NTRS1997

Evaluation of a High-Accuracy MacCormack-Type Scheme Using Benchmark Problems

Due to their inherent dissipation and stability, the MacCormack scheme and its variants have been widely used in the computation of unsteady flow and acoustic problems. However, these schemes require many points per wavelength in order to propagate waves with a reasonable amount of accuracy. In this work, the linear wave propagation characteristics of MacCormack-type schemes are shown by solving several of the CAA Benchmark Problems.

Hixon, R.↗
NASA NTRS2014

Validation of Land Surface Temperature Products Derived from the Visible Infrared Imaging Radiometer Suite (VIIRS) Using Ground-based and Heritage Satellite Measurements

Thermal infrared satellite observations of the Earth's surface are widely used to retrieve Land Surface Temperature (LST) and monitor LST changes around the world. Since January 2012, the Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar-Orbiting Partnership (S-NPP) has provided daily observations of LST with a spatial resolution of 750 m at nadir. Comparison of the standard VIIRS LST product with the equivalent daily standard product from the Moderate Resolution Imaging Spectroradiometer (MODIS) collection-5 and with ground-based measurements over vegetated and inland water surfaces showed good agreement. Analysis indicated the accuracy and precision of the VIIRS product over these cover types was 0.2 K and 2.0 K respectively provided the analyses included appropriate compensation for any spatial heterogeneity in LST within the validation site. However, comparisons between in situ LST and the VIIRS and MODIS LST.

Land surface temperature↗
NASA NTRS2023

Regional Analysis of the 2015–16 Lower Mekong River Basin Drought Using NASA Satellite Observations

Study region Lower Mekong River Basin (LMRB) Study focus Satellite remote sensing products are widely used for monitoring droughts. Using NASA satellite sensors of precipitation (Global Measurement Mission, GPM), soil moisture (Soil Moisture Active and Passive, SMAP), and terrestrial water storage (Gravity Recovery and Climate Experiment, GRACE), this study evaluates the historical drought in the LMRB during 2015–16. SMAP soil moisture was validated against in-situ soil moisture, and GPM precipitation and SMAP soil moisture were cross-validated with streamflow observations. The spatiotemporal dynamics of soil moisture were also examined in different ranges of catchment areas. In performing the analysis, we used lagged correlations between hydrological variables and the indices of the Standardized Precipitation Index (SPI) and Standardized Streamflow Index (SSI). New hydrological insights for the regions Spatio-temporal patterns of drought in 2015–16 were examined from the entire basin to small watersheds. A mismatch occurs when using GRACE data to study droughts in small watersheds (many of the small watersheds would be a fraction of the few 100 km2 spatial resolutions of GRACE pixel). In smaller watersheds, hydrological drought (SSI) was closely defined with SMAP soil moisture downscaled to 1 km rather than the meteorological drought index (SPI). By leveraging satellite-based observations across a range of spatial scales, this study highlights the utility of Earth observations in informing water resources and land management decisions at the regional scale.

Mekong↗
NASA NTRS2009

Performance of MODIS Thermal Emissive Bands On-orbit Calibration Algorithms

Two nearly identical copies of the Moderate Resolution Imaging Spectroradiometer (MODIS) are currently operated on-board the Terra and Aqua spacecrafts, launched in December 1999 and May 2002, respectively. Together, they have produced an unprecedented amount of science data products, which are widely used for the studies of changes in the Earth's system of land, oceans, and atmosphere. MODIS is a cross-track scanning radiometer, which uses a two-sided scan mirror and collects data continuously over a wide scan angle range (+/-55 degree relative to the instrument nadir) each scan of 1.47 seconds. It has 36 spectral bands with wavelengths ranging from visible (VIS) to long-wave infrared (LWIR). MODIS bands 1-19 and 26 are the reflective solar bands (RSB) and bands 20-25 and 27-36 are the thermal emissive bands (TEB). MODIS was developed and designed with improvements made over its heritage sensors (such as AVHRR and Landsat) and, in particular, with more stringent calibration requirements. Because of this, MODIS was built with a set of state-of-art on-board calibrators (OBC), which include a solar diffuser (SD), a solar diffuser stability monitor (SDSM), a blackbody (BB), a spectroradiometric calibration assembly (SRCA), and a space view (SV) port. With the exception of view angle differences, MODIS OBC measurements and the Earth View (EV) observations are made via the same optical path. MODIS TEB have a total of 160 individual TEB detectors (10 per band), which are located on two cold focal plane assemblies (CFPA). For nominal on-orbit operation, the CFPA temperature is controlled at 83K via a passive radiative cooler. For the TEB, the calibration requirements at specified typical scene radiances are less than or equal to 1% with an exception for the fire detection (low gain) band. MODIS TEB on-orbit calibration is performed on a scan-by-scan basis using a quadratic calibration algorithm, and data collected from sensor responses to the onboard BB and SV. The BB serves as the thermal calibration source and the SV provides measurements for the sensor's background and offsets. MODIS on-board BB is a v-grooved plate with its temperature measured using 12 platinum resistive thermistors (PRT) uniformly embedded in the BB substrate. All the BB thermistors were characterized pre-launch with reference to the NIST temperature standards. Unlike typical BB operations in many heritage sensors, which have no temperature control capability, the MODIS on-board BB can be operated at any temperatures between instrument ambient (about 270K) and 315K and can also be varied continuously within this range. This feature has significantly enhanced the MODIS' capability of tracking and updating the TEB nonlinear calibration coefficients over its entire mission. Following a brief description of MODIS TEB on-orbit calibration methodologies and its onboard BB operational activities, this paper provides a comprehensive performance assessment of MODIS TEB quadratic calibration algorithm. It examines the scan-by-scan, orbit-by-orbit, daily, and seasonal variations of detector responses and associated impact due changes in the CFPA and instrument temperatures. Specifically, this paper will analyze the contribution by each individual thermal emissive source term (BB, scan cavity, and scan mirror), the impact on the Level 1 B data product quality due to pre-launch and on-orbit calibration uncertainties. A comparison of Terra and Aqua TEB on-orbit performance, lessons learned, and suggestions for future improvements will also be made.

Xiong, Xiaoxiong↗
DOE OSTI2025

Electrostatically Enhanced Buried Interface Binding of Self‐Assembled Monolayers for Efficient And Stable Inverted Perovskite Solar Cells

Inverted p‐i‐n structure perovskite solar cells (PSCs) have outperformed traditional n‐i‐p PSCs in recent years. A key advancement is the use of self‐assembled monolayers (SAMs) as hole transport layers. One class of widely used SAMs is carbazole‐based phosphonic acids. However, it is found that these SAMs lack strong binding with transparent conducting oxides (TCO) and perovskite. The weak binding strength results in suboptimal interfacial adhesion of the buried interface, which limits the device's stability. Here, interfacial binding is enhanced by increasing the dipole moment that creates a strong interfacial electric field that enhances electrostatic interactions at the TCO/perovskite interface, while incorporating tailored functional groups in SAMs to improve chemical anchoring to TCO and binding to perovskite. Specifically, the donor‐acceptor SAM molecule 4‐(7‐(4‐(bis(4‐methoxyphenyl)amino)‐2,5‐difluorophenyl)benzo[c][1,2,5]thiadiazol‐4‐yl)benzoic acid (PAFTB) is employed, which features an enhanced dipole moment along with electron‐donating and electron‐withdrawing functional groups to optimize interfacial interactions. Compared to extensively used [2‐(9H‐carbazol‐9‐yl)ethyl]phosphonic acid (2PACz), PAFTB enhances total interfacial adhesion by 2.8 times, thereby improving the thermal stability of the layer. Using this approach, PSCs are demonstrated with a certified quasi‐steady‐state power conversion efficiency of 24.9% and maintain 80% of the initial efficiency after 900 h of maximum power point tracking at 85 °C.

Interfacial adhesion↗
DOE OSTI2021

Direct comparison of gamma, electron beam and X-ray irradiation doses on characteristics of low-density polyethylene, polypropylene homopolymer, polyolefin elastomer and chlorobutyl rubber medical device polymers

There is a growing need for increased efficiency in the sterilization of single use medical devices and other products that contain polymer components. Gamma radiation is widely used for devices suited for radiation sterilization; however, safety, throughput and cobalt-60 source availability are challenging the prospect of relying on gamma radiation to meet the anticipated needs of the industry. Use of electron beam (e-beam) and X-rays as alternatives to gamma for radiation sterilization have been hampered in part by a concern that these modalities may adversely affect polymer integrity and performance relative to the gamma method, for which the industry has had much more experience. To address this concern, the effects of sterilization-relevant doses of ebeam, X-ray and gamma radiation were directly compared using common medical device polymers found in two prototypical commercial devices currently sterilized using cobalt-60 gamma irradiation. The Becton, Dickinson and Company (BD) Vacutainer™ Plus tube contains low-density polyethylene and chlorobutyl rubber components, while the BD Vacutainer™ Push Button Blood Collection Set contains polypropylene homopolymer and polyolefin elastomer components. Injection-molded samples prepared from the polymers used in these products were exposed to target doses of 15, 35, 50 and 80 kGy using gamma, e-beam and X-ray radiation. Changes in coloration, tensile properties and hardness were measured for each condition, and the effects of e-beam and X-ray irradiation compared with the effects of gamma irradiation on these properties. Both e-beam and X-ray appear as viable alternatives to gamma irradiation for sterilization of the polymers tested.

sterilization, gamma, X-ray, electron beam, radiat↗
DOE OSTI2022

Stable and unstable supersonic stagnation of an axisymmetric rotating magnetized plasma

The Naval Research Laboratory ‘Mag Noh problem’, described in this paper, is a self-similar magnetized implosion flow, which contains a fast magnetohydrodynamic (MHD) outward propagating shock of constant velocity. We generalize the classic Noh ( OSTI Tech. Rep. 577058, 1983) problem to include azimuthal and axial magnetic fields as well as rotation. Our family of ideal MHD solutions is five parametric, each solution having its own self-similarity index, gas gamma, magnetization, the ratio of axial to the azimuthal field and rotation. While the classic Noh problem must have a supersonic implosion velocity to create a shock, our solutions have an interesting three-parametric special case with zero initial velocity in which magnetic tension, instead of implosion flow, creates the shock at $t=0+$ . Our self-similar solutions are indeed realized when we solve the initial value MHD problem with the finite volume MHD code Athena. We numerically investigated the stability of these solutions and found both stable and unstable regions in parameter space. Stable solutions can be used to test the accuracy of numerical codes. Unstable solutions have also been widely used to test how codes reproduce linear growth, transition to turbulence and the practically important effects of mixing. Now we offer a family of unstable solutions featuring all three elements relevant to magnetically driven implosions: convergent flow, magnetic field and a shock wave.

Mechanics↗