Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model discrepancy”

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 1,441 records · Page 80

Solving inverse problems in stochastic models using deep neural networks and adversarial training

Inverse problems associated with stochastic models constitute a significant portion of scientific and engineering applications. In such cases the unknown quantities are distributions. The applicability of traditional methods is limited because of their demanding assumptions or prohibitive computational consumption; for example, maximum likelihood methods require closed-form density functions, and Markov Chain Monte Carlo needs a large number of simulations. We propose a new method that estimates the unknown distribution by matching the statistical properties between observed and simulated random processes. We leverage the expressive power of neural networks to approximate the unknown distribution and use a discriminative neural network for computing the statistical discrepancies between the observed and simulated random processes. Here we demonstrated numerically that the proposed methods can estimate both the model parameters and learn complicated unknown distributions.

42 ENGINEERING↗

Multipoint Observations of Energetic Particle Injections and Substorm Activity During a Conjunction Between Magnetospheric Multiscale (MMS) and Van Allen Probes

Here, this study examines multipoint observations during a conjunction between Magnetospheric Multiscale (MMS) and Van Allen Probes on 7 April 2016 in which a series of energetic particle injections occurred. With complementary data from Time History of Events and Macroscale Interactions during Substorms, Geotail, and Los Alamos National Laboratory spacecraft in geosynchronous orbit (16 spacecraft in total), we develop new insights on the nature of energetic particle injections associated with substorm activity. Despite this case involving only weak substorm activity (maximum AE <300 nT) during quiet geomagnetic conditions in steady, below-average solar wind, a complex series of at least six different electron injections was observed throughout the system. Intriguingly, only one corresponding ion injection was clearly observed. All ion and electron injections were observed at <600 keV only. MMS reveals detailed substructure within the largest electron injection. A relationship between injected electrons with energy <60 keV and enhanced whistler mode chorus wave activity is also established from Van Allen Probes and MMS. Drift mapping using a simplified magnetic field model provides estimates of the dispersionless injection boundary locations as a function of universal time, magnetic local time, and L shell. The analysis reveals that at least five electron injections, which were localized in magnetic local time, preceded a larger injection of both electrons and ions across nearly the entire nightside of the magnetosphere near geosynchronous orbit. The larger ion and electron injection did not penetrate to L < 6.6, but several of the smaller electron injections penetrated to L < 6.6. Due to the discrepancy between the number, penetration depth, and complexity of electron versus ion injections, this event presents challenges to the current conceptual models of energetic particle injections.

58 GEOSCIENCES↗

A redetermination of the succinonitrile-water phase diagram

The phase diagram of the water/succinonitrile system often used in modeling binary metallic alloys is determined experimentally using a laser scattering method. Triply distilled and zone-refined succinonitrile is placed in a jacketed water-cooled vessel with a magnetic mixer and a copper-constantan thermocouple, and the melting point is determined under vacuum and under Ar by observing the response of a photodiode placed opposite a laser source as the chamber is cooled and reheated. Then water is added, and the eutectic and monotectic temperatures and the entire phase diagram are determined under Ar. The results are presented in a table and graph and compared to those of Schreinemakers (1897): significant discrepancies are noted and attributed to impurities in the succinonitrile used by Schreinemakers. The importance of accurate phase diagrams for modeling binary alloys is stressed.

Smith, J. E., Jr.↗

Measurement and Modeling of Vertically Resolved Aerosol Optical Properties and Radiative Fluxes Over the ARM SGP Site

In order to meet one of its goals - to relate observations of radiative fluxes and radiances to the atmospheric composition - the Department of Energy's Atmospheric Radiation Measurement (ARM) program has pursued measurements and modeling activities that attempt to determine how aerosols impact atmospheric radiative transfer, both directly and indirectly. However, significant discrepancies between aerosol properties measured in situ or remotely remain. One of the objectives of the Aerosol Intensive Operational Period (TOP) conducted by ARM in May 2003 at the ARM Southern Great Plains (SGP) site in north central Oklahoma was to examine and hopefully reduce these differences. The IOP involved airborne measurements from two airplanes over the heavily instrumented SGP site. We give an overview of airborne results obtained aboard the Center for Interdisciplinary Remotely-Piloted Aircraft Studies (CIRPAS) Twin Otter aircraft. The Twin Otter performed 16 research flights over the SGP site. The aircraft carried instrumentation to perform in-situ measurements of aerosol absorption, scattering, extinction and particle size. This included such novel techniques as the photoacoustic and cavity ring-down methods for in-situ absorption (675 nm) and extinction (675 and 1550 nm) and a new multiwavelength, filter-based absorption photometer (467, 530, 660 nm). A newly developed instrument measured cloud condensation nucleus concentration (CCN) concentrations at two supersaturation levels. Aerosol optical depth and extinction (354-2139 nm) were measured with the NASA Ames Airborne Tracking 14-channel sunphotometer. Furthermore, up-and downwelling solar (broadband and spectral) and infrared radiation were measured using seven individual radiometers. Three up-looking radiometers werer mounted on a newly developed stabilized platform, keeping the instruments level up to aircraft pitch and roll angles of approximately 10(exp 0). This resulted in unprecedented continuous vertical profiles of radiative fluxes, which we will compare to modeled fluxes using the aforementioned data as input.

Schmid, B.↗

Relating flow resistance to equivalent roughness

Describing flow resistance using the physical properties of an underlying surface is a recalcitrant problem in overland flow models. If discharge measurements are available, an equivalent roughness (e.g., Manning’s n) can be calibrated to represent the effects of surface properties within the domain with a single numerical value. Alternatively, the flow resistance can be estimated from discharge and velocity measured at a point, typically a runoff plot outlet. However, such experimental estimates are often inconsistent with the equivalent roughness determined from calibration to discharge, even if both derive from the same dataset. For example, if Manning’s equation is used to parameterize flow resistance, the Manning’s n obtained by calibrating a model to discharge differs from the value of n calculated from measured flow and velocity at the hillslope outlet. Here, this discrepancy is resolved by deriving a correction factor relating experimentally-determined flow resistance to the equivalent roughness. The derived correction factor is tested for four commonly-used resistance formulations using 129 rainfall simulator experiments. The correction factor is necessary to reproduce measured velocities, and yields minor improvements in discharge prediction. Plain Language Summary: Accurate runoff prediction is needed for land and water management in dryland regions, where sporadic and limited rainfall necessitate efficient water use and drought mitigation strategies. The skill of runoff models is known to be hindered by out ability to estimate flow resistance, which is the quantity that describes how energy is lost from flowing water to the underlying surface. Typically, models represent flow resistance with an equivalent roughness, e.g., Manning’s n, that is adjusted until the model can reproduce available discharge observations at watershed scale. However, the flow resistance measured in plot-scale experiments (1–10 m) often exceeds equivalent roughness coefficients by a factor of 10. This means that the direct use of plot-scale experimental data to parameterize runoff models could cause errors in discharge and runoff velocity predictions. Here, we resolve these differences by deriving an analytic correction factor that relates flow resistance to the equivalent roughness required for models to reproduce experimental velocity and discharge data. This correction factor is tested using rainfall simulator data from 129 experiments performed in the US Southwest covering a wide range of precipitation intensities, soil textures and vegetation types. Use of the correction factor substantially improves model prediction of flow velocity, which is needed for reproducing the timing of flood events and the estimation of erosion.

54 ENVIRONMENTAL SCIENCES↗

A regional simulation of marine boundary-layer clouds

A regional version of a two-layer model is used to simulate marine boundary-layer clouds over an area (20 deg - 42 deg and 117 - 145 deg W) of the eastern North Pacific. The large-scale conditions requried for the model are provided by European Centre for Medium-Range Weather Forecasts analyses and the Comprehensive Ocean-Atmosphere Data Set. The simulated cloud-top height and fractional cloudiness are compared with satellite data. The model predicts a realistic pattern of cloud-top height, although the heights appear to be overestimated in the southwest and underestimated along the eastern border of the domain. The overall simulated cloudiness is qualitatively comparable to the satellite derived, but the model predicts excessive cloudiness over cold ocean surfaces and insufficient cloudiness over an area centered at 30 deg N, 130 deg W. The possible causes for these discrepancies are discussed. Budget studies show that horizontal advection contributes substantially to the cloud-top height and cloud-layer moisture budgets. Thus, one may need to consider the effects of horizontal advection when validating steady-state solutions of a one-dimensional model. The simulated cloudiness is shown to be sensitive to large-scale subsidence, moisture above the clouds, and drizzle. An increase in the large-scale divergence by 2 x 10(exp -6)/sec decreases average cloudiness for the area by 30%. An increase in the moisture at 850 mb by 2 g/kg increases the cloudiness by 20%. The suppression of drizzle in the model increases the cloudiness by 34% and significantly changes its pattern.

Wang, Shouping↗

Comparison of ozone formation attribution techniques in the northeastern United States

The Integrated Source Apportionment Method (ISAM) has been revised in the Community Multiscale Air Quality (CMAQ) model. This work updates ISAM to maximize its flexibility, particularly for ozone (O 3 ) modeling, by providing multiple attribution options, including products inheriting attribution fully from nitrogen oxide reactants, fully from volatile organic compound (VOC) reactants, equally from all reactants, or dynamically from NO x or VOC reactants based on the indicator gross production ratio of hydrogen peroxide (H 2 O 2 ) to nitric acid (HNO 3 ). The updated ISAM has been incorporated into the most recent publicly accessible versions of CMAQ (v5.3.2 and beyond). This study's primary objective is to document these ISAM updates and demonstrate their impacts on source apportionment results for O 3 and its precursors. Additionally, the ISAM results are compared with the Ozone Source Apportionment Technology (OSAT) in the Comprehensive Air-quality Model with Extensions (CAMx) and the brute-force method (BF). All comparisons are performed for a 4 km horizontal grid resolution application over the northeastern US for a selected 2 d summer case study (9 and 10 August 2018). General similarities among ISAM, OSAT, and BF results add credibility to the new ISAM algorithms. However, some discrepancies in magnitude or relative proportions among tracked sources illustrate the distinct features of each approach, while others may be related to differences in model formulation of chemical and physical processes. Despite these differences, OSAT and ISAM still provide useful apportionment data by identifying the geographical and temporal contributions of O 3 and its precursors. Both OSAT and ISAM attribute the majority of O 3 and NO x contributions to boundary, mobile, and biogenic sources, whereas the top three contributors to VOCs are found to be biogenic, boundary, and area sources.

54 ENVIRONMENTAL SCIENCES↗

CLUST - EVAP Monte Carlo Simulation Applications for Determining Effective Energy Deposition in Silicon by High Energy Protons

The CLUST-EVAP is a Monte Carlo simulation of the interaction of high energy (25 - 400 MeV) protons with silicon nuclei. The initial nuclear cascade stage is modeled using the CLUST model developed by Indiana University over 30 years ago. The second stage, in which the excited nucleus evaporates particles in random directions, is modeled according to the evaporation algorithm provided by H. H. K. Tang of IBM. Using the CLUST-EVAP code to model fragment produ6tion and the Vavilov-Landau theory to model fluctuations in direct ionization in thin silicon layers, we have predicted energy deposition in silicon components for various geometrical configurations. We have compared actual measurements with model predictions for geometry's such as single, thin silicon particle detectors, telescopic particle detectors flown in space to measure the environment, and thin sensitive volumes of modern micro-electronic components. We have recently compared the model predictions with actual measurements made by the DOSTEL spectrometer flown in the Shuttle payload bay on STS-84. The model faithfully reproduces the features and aids in interpretation of flight results of this instrument. We have also applied the CLUST-EVAP model to determine energy deposition in the thin sensitive volumes of modern micro-electronic components. We have accessed the ability of high energy (200 MeV) protons to induce latch-up in certain devices that are known to latch up in heavy ion environments. However, some devices are not nearly as susceptible to proton induced latch-up as expected according to their measured heavy ion latch-up cross sections. The discrepancy is believed to be caused by the limited range of the proton-silicon interaction fragments. The CLUST-EV AP model was used to determine a distribution of these fragments and their range and this is compared to knowledge of the ranges required based on the known device structure. This information is especially useful in accessing the risk to on-orbit perfonnance in a heavy ion environment based on testing performed with only protons.

ONeill, Pat M.↗

How does uncertainty of soil organic carbon stock affect the calculation of carbon budgets and soil carbon credits for croplands in the U.S. Midwest?

Cropland carbon budget depicts the amount of carbon flowing in and out of agroecosystems and the changes in carbon stocks of soil and living biomass during the same period. Soil carbon credit is the additional change in soil carbon stock under certain farming practices compared with the business-as-usual practices. Accurately calculating cropland carbon budget and soil carbon credit is critical to assessing climate change mitigation potential in agroecosystems. The calculation of cropland carbon budget and soil carbon credit is sensitive to local soil and climatic conditions, especially initial soil organic carbon (SOC) stock, which is determined by both SOC concentration (SOC%) and bulk density (Bulk_Density). SOC stock data are either from soil sampling or gridded public survey data. In agroecosystem models, SOC stock data are a key model input for quantifying cropland carbon budget and soil carbon credit. However, various types and degrees of uncertainties exist in SOC stock datasets, which propagate to the quantification of SOC stock change. In particular, a large discrepancy is found in two widely used SOC stock datasets — Rapid Carbon Assessment dataset (RaCA) and Gridded Soil Survey Geographic Database (gSSURGO) — in the U.S. Midwest, with a relative difference (quantified using Normalized Root Mean Square Error, NRMSE) of 48.0% for 0–30 cm SOC stock between the two datasets. It remains largely unclear how uncertainty in SOC stocks affects the calculation of cropland carbon budget and soil carbon credit. To address this question, we used a well-validated process-based agroecosystem model, ecosys, to assess the impacts of SOC stock uncertainty on carbon budget and soil carbon credit calculation in the U.S. Midwestern corn-soybean rotation systems. Our results reveal the following findings: (1) A sizable discrepancy exists in simulated cropland carbon budget between using gSSURGO and using RaCA for their SOC% and Bulk_Density as model inputs, with a Pearson correlation coefficient (r) of only 0.4 for simulated change of SOC stock (ΔSOC) using these two different soil datasets. (2) Simulated cropland carbon budget components were more sensitive to initial SOC% than to Bulk_Density. For example, the upper and lower quartiles of multi-year averaged ΔSOC were –29.8 and 4.8 gC/m 2 /year for the selected counties respectively, with an uncertainty of 13.7 and 0.7 gC/m 2 /year induced by uncertainties in initial SOC% and Bulk_Density, respectively. (3) Both simulated ΔSOC and its uncertainty were negatively correlated with initial SOC%, whereas ΔSOC was negatively correlated with air temperature, and ΔSOC uncertainty was positively correlated with air temperature. (4) The uncertainty of calculated soil carbon credits was much smaller compared with the uncertainty of calculated absolute carbon budgets assuming the same SOC stock uncertainty level in the inputs. Specifically, in our assessment comparing planting cover crops vs no cover crop, the uncertainty of calculated soil carbon credits induced by initial SOC% uncertainty was less than 4% (relative to the quantified value of the soil carbon credits) for 90% of the cases. Our analysis highlights that high accuracy measurement of SOC% as inputs is needed for the calculation of cropland carbon budgets; however, soil carbon credit quantification is much less sensitive to the initial SOC% inputs, and the current publicly available soil datasets (e.g., gSSURGO) are largely suitable for the calculation of soil

54 ENVIRONMENTAL SCIENCES↗

SDYN-GANs: Adversarial learning methods for multistep generative models for general order stochastic dynamics

We introduce adversarial learning methods for data-driven generative modeling of dynamics of nth-order stochastic systems. Our approach builds on Generative Adversarial Networks (GANs) with generative model classes based on stable m-step stochastic numerical integrators. From observations of trajectory samples, we introduce methods for learning long-time predictors and stable representations of the dynamics. Our approaches use discriminators based on Maximum Mean Discrepancy (MMD), training protocols using both conditional and marginal distributions, and methods for learning dynamic responses over different time-scales. We show how our approaches can be used for modeling physical systems to learn force-laws, damping coefficients, and noise-related parameters. Our adversarial learning approaches provide methods for obtaining stable generative models for dynamic tasks including long-time prediction and developing simulations for stochastic systems.

• Artificial intelligence (AI) / machine learning ↗

Advanced Diagnostic System on Earth Observing One

In this infusion experiment, the Livingstone 2 (L2) model-based diagnosis engine, developed by the Computational Sciences division at NASA Ames Research Center, has been uploaded to the Earth Observing One (EO-1) satellite. L2 is integrated with the Autonomous Sciencecraft Experiment (ASE) which provides an on-board planning capability and a software bridge to the spacecraft's 1773 data bus. Using a model of the spacecraft subsystems, L2 predicts nominal state transitions initiated by control commands, monitors the spacecraft sensors, and, in the case of failure, isolates the fault based on the discrepant observations. Fault detection and isolation is done by determining a set of component modes, including most likely failures, which satisfy the current observations. All mode transitions and diagnoses are telemetered to the ground for analysis. The initial L2 model is scoped to EO-1's imaging instruments and solid state recorder. Diagnostic scenarios for EO-1's nominal imaging timeline are demonstrated by injecting simulated faults on-board the spacecraft. The solid state recorder stores the science images and also hosts: the experiment software. The main objective of the experiment is to mature the L2 technology to Technology Readiness Level (TRL) 7. Experiment results are presented, as well as a discussion of the challenging technical issues encountered. Future extensions may explore coordination with the planner, and model-based ground operations.

Hayden, Sandra C.↗

Global-scale multidecadal variability missing in state-of-the-art climate models

Reliability of future global warming projections depends on how well climate models reproduce the observed climate change over the twentieth century. In this regard, deviations of the model-simulated climate change from observations, such as a recent “pause” in global warming, have received considerable attention. Such decadal mismatches between model-simulated and observed climate trends are common throughout the twentieth century, and their causes are still poorly understood. Here we show that the discrepancies between the observed and simulated climate variability on decadal and longer timescale have a coherent structure suggestive of a pronounced Global Multidecadal Oscillation. Surface temperature anomalies associated with this variability originate in the North Atlantic and spread out to the Pacific and Southern oceans and Antarctica, with Arctic following suit in about 25–35 years. While climate models exhibit various levels of decadal climate variability and some regional similarities to observations, none of the model simulations considered match the observed signal in terms of its magnitude, spatial patterns and their sequential time development. These results highlight a substantial degree of uncertainty in our interpretation of the observed climate change using current generation of climate models.

54 ENVIRONMENTAL SCIENCES↗

SWIFT Code Assessment for Two Similar Transonic Compressors

One goal of the NASA Fundamental Aeronautics Program is the assessment of computational fluid dynamic (CFD) codes used for the design and analysis of many aerospace systems. This paper describes the assessment of the SWIFT turbomachinery analysis code for two similar transonic compressors, NASA rotor 37 and stage 35. The two rotors have identical blade profiles on the front, transonic half of the blade but rotor 37 has more camber aft of the shock. Thus the two rotors have the same shock structure and choking flow but rotor 37 produces a higher pressure ratio. The two compressors and experimental data are described here briefly. Rotor 37 was also used for test cases organized by ASME, IGTI, and AGARD in 1994-1998. Most of the participating codes over predicted pressure and temperature ratios, and failed to predict certain features of the downstream flowfield. Since then the AUSM+ upwind scheme and the k- turbulence model have been added to SWIFT. In this work the new capabilities were assessed for the two compressors. Comparisons were made with overall performance maps and spanwise profiles of several aerodynamic parameters. The results for rotor 37 were in much better agreement with the experimental data than the original blind test case results although there were still some discrepancies. The results for stage 35 were in very good agreement with the data. The results for rotor 37 were very sensitive to turbulence model parameters but the results for stage 35 were not. Comparison of the rotor solutions showed that the main difference between the two rotors was not blade camber as expected, but shock/boundary layer interaction on the casing.

Chima, Rodrick V.↗

The role of radical-radical chain-propagating pathways in the phenyl + propargyl reaction

Well-skipping radical-radical reactions can provide a chain-propagating pathway for formation of polycyclic radicals implicated in soot inception. Here we use controlled pyrolysis in a microreactor to isolate and examine the role of well-skipping channels in the phenyl (C 6 H 5 ) + propargyl (C 3 H 3 ) radical-radical reaction at temperatures of 800-1600 K and pressures near 25 Torr. The temperature and concentration dependence of the closed-shell (C 9 H 8 ) and radical (C 9 H 7 ) products are observed using electron-ionization mass spectrometry. The flow in the reactor is simulated using a boundary layer model employing a chemical mechanism based on recent rate coefficient calculations. Further, comparison between simulation and experiment shows reasonable agreement, within a factor of 3, while suggesting possible improvements to the model. In contrast, eliminating the well-skipping reactions from the chemistry mechanism causes a much larger discrepancy between simulation and experiment in the temperature dependence of the radical concentration, revealing that the well-skipping pathways, especially to form indenyl radical, are significant at temperatures of 1200 K and higher. While most C 9 H 7 forms by well-skipping at 25 Torr, an additional simulation indicates that the wellskipping channels only contribute around 3% of the C 9 H x yield at atmospheric pressure, thus indicating a negligible role of the well-skipping pathways at atmospheric and higher pressures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Reaction of Criegee Intermediate CH 2 OO with Water Dimer: Primary Products and Atmospheric Impact

The rapid reaction of the smallest Criegee intermediate, CH 2 OO, with water dimers is the dominant removal mechanism for CH 2 OO in the Earth's atmosphere, but its products are not well understood. This reaction was recently suggested as a significant source of the most abundant tropospheric organic acid, formic acid (HCOOH), which is consistently underpredicted by atmospheric models. Furthermore, using time-resolved measurements of reaction kinetics by UV absorption and product analysis by photoionization mass spectrometry, we show that the primary products of this reaction are formaldehyde and hydroxymethyl hydroperoxide (HMHP), with direct HCOOH yields of less than 10%. Incorporating our results into a global chemistry-transport model further reduces HCOOH levels by 10–90%, relative to previous modeling assumptions, which indicates that the reaction CH 2 OO + water dimer by itself cannot resolve the discrepancy between the measured and predicted HCOOH levels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Endothermic self-interacting dark matter in Milky Way-like dark matter haloes

ABSTRACT Self-interacting dark matter (SIDM) offers the potential to mitigate some of the discrepancies between simulated cold dark matter (CDM) and observed galactic properties. We introduce a physically motivated SIDM model to understand the effects of self interactions on the properties of Milky Way and dwarf galaxy sized haloes. This model consists of dark matter with a nearly degenerate excited state, which allows for both elastic and inelastic scattering. In particular, the model includes a significant probability for particles to up-scatter from the ground state to the excited state. We simulate a suite of zoom-in Milky Way-sized N-body haloes with six models with different scattering cross sections to study the effects of up-scattering in SIDM models. We find that the up-scattering reaction greatly increases the central densities of the main halo through the loss of kinetic energy. However, the physical model still results in significant coring due to the presence of elastic scattering and down-scattering. These effects are not as apparent in the subhalo population compared to the main halo, but the number of subhaloes is reduced compared to CDM.

Astronomy & Astrophysics↗

First detection of chromospheric magnetic field changes during an X1-flare

Stepwise changes of the photospheric magnetic field, which often becomes more horizontal, have been observed during many flares. Previous interpretations include coronal loops that contract, and it has been speculated that such jerks could be responsible for sunquakes. Here we report the detection of stepwise chromospheric line-of-sight magnetic field (B{sub LOS}) changes obtained through spectropolarimetry of Ca ii 8542 Å with DST/IBIS during the X1-flare SOL20140329T17:48. These changes are stronger (<640 Mx cm{sup −2}) and appear in larger areas than their photospheric counterparts (<320 Mx cm{sup −2}). The absolute value of B{sub LOS} more often decreases than increases. Photospheric changes are predominantly located near a polarity inversion line, and chromospheric changes near the footpoints of loops. The locations of changes are near, but not exactly co-spatial to hard X-ray emission and neither to enhanced continuum emission nor to a small sunquake. Enhanced chromospheric and coronal emission is observed in nearly all locations that exhibit changes of B{sub LOS}, but the emission also occurs in many locations without any B{sub LOS} changes. Photospheric and chromospheric changes of B{sub LOS} show differences in timing, sign, and size and seem independent of each other. A simple model of contracting loops yields changes of the opposite sign to those observed. An explanation for this discrepancy could be increasing loop sizes or loops that untwist in a certain direction during the flare. It is yet unclear which processes are responsible for the observed changes and their timing, size, and location, especially considering the incoherence between the photosphere and the chromosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Laboratory Study of the Torus Instability Threshold in Solar-relevant, Line-tied Magnetic Flux Ropes

Coronal mass ejections (CMEs) occur when long-lived magnetic flux ropes (MFRs) anchored to the solar surface destabilize and erupt away from the Sun. This destabilization is often described in terms of an ideal magnetohydrodynamic instability called the torus instability. It occurs when the external magnetic field decreases sufficiently fast such that its decay index, ${n}_{}=-z\,\partial (\mathrm{ln}{B}_{})/\partial z$, is larger than a critical value, $n\gt {n}_{\mathrm{cr}}^{}$, where ${n}_{\mathrm{cr}}^{}=1.5$ for a full, large aspect ratio torus. However, when this is applied to solar MFRs, a range of conflicting values for ${n}_{\mathrm{cr}}^{}$ is found in the literature. To investigate this discrepancy, we have conducted laboratory experiments on arched, line-tied flux ropes and applied a theoretical model of the torus instability. Our model describes an MFR as a partial torus with foot points anchored in a conducting surface and numerically calculates various magnetic forces on it. This calculation yields better predictions of ${n}_{\mathrm{cr}}^{}$ that take into account the specific parameters of the MFR. Here, we describe a systematic methodology to properly translate laboratory results to their solar counterparts, provided that the MFRs have a sufficiently small edge safety factor or, equivalently, a large enough twist. After this translation, our model predicts that ${n}_{\mathrm{cr}}^{}$ in solar conditions falls near ${n}_{\mathrm{cr}}^{\mathrm{solar}}\sim 0.9$ and within a larger range of ${n}_{\mathrm{cr}}^{\mathrm{solar}}\sim (0.7,1.2)$, depending on the parameters. The methodology of translating laboratory MFRs to their solar counterparts enables quantitative investigations of CME initiation through laboratory experiments. These experiments allow for new physics insights that are required for better predictions of space weather events but are difficult to obtain otherwise.

79 ASTRONOMY AND ASTROPHYSICS↗