Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D equilibrium”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Progress in computing nozzle/plume flow fields

The long-term goal is to develop the capability to predict chemically-reacting, multi-stream nozzle and plume flow fields. Two basic Navier-Stokes solvers, including the widely used F-3D code, are upgraded to include several upwind difference schemes and portable chemistry packages. Current computational capabilities for solving equilibrium single-stream and multi-stream, frozen gas, and finite rate chemistry problems are described. A variety of complex nozzle and plume flows were computed. Solutions presented include axisymmetric plume flow for ideal and equilibrium air, 3-D NASP nozzle/afterbody flow, and an internal nozzle calculation comparing various finite-rate chemistry packages.

Ruffin, Stephen M.↗
NASA NTRS1993

A global three-dimensional model of the stratospheric sulfuric acid layer

A 3D model which encompasses SO2 production from OCS, followed by its oxidation to gaseous H2SO4, the condensation-evaporation equilibrium of gaseous and particulate H2SO4, and finally particle condensation and rainout, is presently used to study processes maintaining the nonvolcanically-perturbed stratosphere's sulfuric acid layer. A comparison of the results thus obtained with remotely sensed stratospheric aerosol extinction data shows the model to simulate the general behavior of stratospheric aerosol extinction.

Golombek, Amram↗
NASA NTRS2008

The CME Flare Arcade and the Width of the CME in the Outer Corona

Moore, Sterling, & Suess (2007, ApJ, 668, 1221) present evidence that (1) a CME is typically a magnetic bubble, a low-beta gplasmoid with legs h having roughly the 3D shape of a light bulb, and (2) in the outer corona the CME plasmoid is in lateral pressure equilibrium with the ambient magnetic field. They present three CMEs observed by SOHO/LASCO, each from a very different source located near the limb. One of these CMEs came from a compact ejective eruption from a small part of a sunspot active region, another came from a large quiet-region filament eruption, and the third CME, an extremely large and fast one, was produced in tandem with an X20 flare arcade that was centered on a huge delta sunspot. Each of these CMEs had more or less the classic lightbulb silhouette and attained a constant heliocentric angular width in the outer corona. This indicates that the CME plasmoid attained lateral magnetic pressure balance with the ambient radial magnetic field in the outer corona. This lateral pressure balance, together with the standard scenario for CME production by the eruption of a sheared-core magnetic arcade, yields the following simple estimate of the strength B(sub Flare) of the magnetic field in the flare arcade produced together with the CME: B(sub Flare) 1.4(theta CME/theta Flare)sup 2 G, where theta (sub CME) is the heliocentric angular width of the CME plasmoid in the outer corona and theta (sub Flare) is the heliocentric angular width of the full-grown flare arcade. Conversely, theta (sub CME) approximately equal to (R(sub Sun)sup -1(phi(sub Flare)/1.4)sup 1/2 radians, where Flare is the magnetic flux covered by the full-grown flare arcade. In addition to presenting the three CMEs of Moore, Sterling, & Suess (2007) and their agreement with this relation between CME and Flare, we present a further empirical test of this relation. For CMEs that erupt from active regions, the co-produced flare arcade seldom if ever covers the entire active region: if AR is the total magnetic flux of the active region, Flare . AR, and we predict that CME. (R(sub Sun))sup -1(theta AR/1.4)sup 1/2 radians. For a random sample of 31 CMEs that erupted from active regions within 30 of the limb, for each CME we have measured CME from LASCO/C3 and have measured AR from a SOHO/MDI magnetogram of the source active region when it was within 30 of disk center. We find that each CME obeys the above predicted inequality, none having width greater than half of the upper bound given by theta(sub AR). Thus, an active region's magnetic flux content, together with its location on the solar disk, largely determines whether the active region can possibly produce a CME that is wide enough to intercept the Earth.

Moore, Ron↗
DOE OSTI2025

Novel Electrochemical Technique: Serial Open Circuit Potentiometry to Examine the Separability of Lanthanides in Molten Salts

The accumulation of lanthanide fission products in molten salts reduces the efficiency of pyroprocessing used nuclear fuel. Lanthanide fission products are notoriously difficult to separate due to their similar physical and electrical properties. We developed a novel electrochemical technique called serial open circuit potentiometry (SOCP) which allows us to delineate between metal deposition reactions with similar reduction potentials to great precision. SOCP has several advantages over traditional open circuit potentiometry including its ability to create 3D plots that describe all metal deposition reactions in a system. This information can be used to precisely determine equilibrium potentials of metal deposition reactions and quantify their error, determine nucleation potentials, and determine the nucleation and growth mechanism. In this work, we demonstrated the effectiveness of SOCP by studying molten salt mixtures containing CeCl3, GdCl3, and both CeCl3 and GdCl3 in LiCl-KCl eutectic.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Axisymmetric simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

A benchmark exercise for the modeling of vertical displacement events (VDEs) is presented and applied to the 3D nonlinear magnetohydrodynamic codes M3D-C 1 , JOREK, and NIMROD. The simulations are based on a vertically unstable NSTX equilibrium enclosed by an axisymmetric resistive wall with a rectangular cross section. A linear dependence of the linear VDE growth rates on the resistivity of the wall is recovered for sufficiently large wall conductivity and small temperatures in the open field line region. We report the benchmark results show good agreement between the VDE growth rates obtained from linear NIMROD and M3D-C 1 simulations and from the linear phase of axisymmetric nonlinear JOREK, NIMROD, and M3D-C 1 simulations. Axisymmetric nonlinear simulations of a full VDE performed with the three codes are compared, and an excellent agreement is found regarding the plasma location and plasma currents, as well as eddy and halo currents in the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Electromechanical coupling of the solar atmosphere; Proceedings of the OSL Workshop, Capri, Italy, May 27-31, 1991

The present conference discusses the role of magnetic flux tubes as communication channels, flux tube sizes and their temporal evolution, magnetic field line topology in the solar active regions, weak solar magnetic fields, explosive events and magnetic reconnection in the solar atmosphere, and 3D kinematic reconnection of plasmoids with nulls. Also discussed are coronal heating mechanisms, coronal heating through a lack of MHD equilibrium, Alfven waves in current-carrying inhomogeneous plasmas, hydrostatic models of X-ray coronal loops, MHD turbulence in an expanding atmosphere, and hot mass transport in the solar active prominence.

Spicer, Daniel S.↗
NASA NTRS2005

Estimating Aeroheating of a 3D Body Using a 2D Flow Solver

A method for rapidly estimating the aeroheating, shear stress, and other properties of hypersonic flow about a three-dimensional (3D) blunt body has been devised. First, the geometry of the body is specified in Cartesian coordinates. The surface of the body is then described by its derivatives, coordinates, and principal curvatures. Next, previously relatively simple equations are used to find, for each desired combination of angle of attack and meridional angle, a scaling factor and the shape of an equivalent axisymmetric body. These factors and equivalent shapes are entered as inputs into a previously developed computer program that solves the two-dimensional (2D) equations of flow in a non-equilibrium viscous shock layer (VSL) about an axisymmetric body. The coordinates in the output of the VSL code are transformed back to the Cartesian coordinates of the 3D body, so that computed flow quantities can be registered with locations in the 3D flow field of interest. In tests in which the 3D bodies were elliptic paraboloids, the estimates obtained by use of this method were found to agree well with solutions of 3D, finite-rate-chemistry, thin-VSL equations for a catalytic body.

Scott, Carl D.↗
DOE OSTI2026

Hydride and Seek: Comparing Crystallographic Hydride Placement Techniques with an Open-Shell Cobalt Complex

Locating hydrides is crucial in organometallic chemistry but difficult to do accurately using X-ray diffraction. Electron diffraction has been proposed as a way to overcome this problem but has not been systematically compared to neutron diffraction and to quantum crystallography (Hirshfeld atom refinement, HAR) to test this hypothesis. Here, we present a comparative analysis of methods for a terminal cobalt hydride complex by comparing a single-crystal neutron diffraction reference structure to results from single-crystal X-ray diffraction with and without Hirshfeld atom refinement (HAR, NoSpherA2), density functional theory (DFT), and electron diffraction (3D-ED/MicroED) refined under kinematical and dynamical formalisms. Conventional X-ray diffraction gives lower precision than neutron diffraction as expected. Despite expected improvements, HAR gives systematic deviation from the neutron benchmark. Interestingly, optimized DFT equilibrium geometries are closer to the neutron value than the value from HAR. On the other hand, electron diffraction with a high-quality data set coupled with dynamical refinement localizes the hydride in difference maps and gives excellent agreement with the neutron data. Dynamical refinement is crucial, as kinematical refinement does not allow assignment of a hydride peak. This cross-modal comparison defines the conditions under which 3D-ED/MicroED delivers high-precision metal–hydride distances for this open-shell cobalt hydride.

anions↗
DOE OSTI2023

Grain boundary effects in high-temperature liquid-metal dealloying: a multi-phase field study

Abstract A multi-phase field model is employed to study the microstructural evolution of an alloy undergoing liquid dealloying, specifically considering the role of grain boundaries. A semi-implicit time-stepping algorithm using spectral methods is implemented, which enables simulating large 2D and 3D domains over long time scales while still maintaining a realistic interfacial thickness. Simulations reveal a mechanism of coupled grain–boundary migration to maintain equilibrium contact angles with the topologically complex solid–liquid interface, which locally accelerates diffusion-coupled growth of a liquid channel into the precursor. This mechanism asymmetrically disrupts the ligament connectivity of the dealloyed structure in qualitative agreement with published experimental observations. The grain boundary migration-assisted corrosion channels form even for precursors with small amounts of the dissolving alloy species, below the parting limit . The activation of this grain boundary dealloying mechanism depends strongly on grain boundary mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2024

FLARE: field line analysis and reconstruction for 3D boundary plasma modeling

The FLARE code is a magnetic mesh generator that is integrated within a suite of tools for the analysis of the magnetic geometry in toroidal fusion devices. A magnetic mesh is constructed from field line segments and permits fast reconstruction of field lines in 3D boundary plasma codes such as EMC3-EIRENE. Both intrinsically non-axisymmetric configurations (stellarators) and those with symmetry breaking perturbations of an axisymmetric equilibrium (tokamaks) are supported. The code itself is written in Modern Fortran with MPI support for parallel computing, and it incorporates object-oriented programming for the definition of the magnetic field and the material surface geometry. Extended derived types for a number of different magnetohydrodynamic equilibrium and plasma response models are implemented. The core element of FLARE is a field line tracer with adaptive step-size control, and this is integrated into tools for the construction of Poincaré maps and invariant manifolds of X-points. A collection of high-level procedures that generate output files for visualization is build on top of that. The analysis modules are build with Python frontends that facilitate customization of tasks and/or scripting of parameter scans.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Molecular mechanisms of thickness-dependent water desalination in polyamide reverse-osmosis membranes

Here, using non-equilibrium molecular simulations, we systematically elaborate the relationships among synthesis (membrane's thickness and morphology), atomic-scale transport mechanism, and separation performance (permeability and selectivity) for 3D-printed polyamide (PA) membranes. Results indicated that water diffusion, swelling ratio, water flux, and water permeance proportionally decrease with increasing membrane thickness (4.0–32.5 nm). PA membranes with different thicknesses can achieve almost 100% salt rejection over the simulation time. Importantly, water permeability exponentially decreases with increasing thickness, and 15 nm is identified as the critical membrane thickness for efficient water transport. The discontinuous water-available space spreads all over PA membranes with thicknesses greater than 5 nm, allowing water molecules to jump by way of the temporary open-and-close pores. However, the connected water-useable space exists in PA membranes with thicknesses below 5 nm, offering the continuous channels to dominate water transport. More significantly, pore distribution is more homogeneous as the thickness increases. The applied high pressures can lead to membrane compaction during reverse osmosis and the thicker membranes show a lower compression ratio. In short, these investigations provide molecular insights for effectively designing and manufacturing PA membranes for water desalination and treatment at the molecular level.

42 ENGINEERING↗
DOE OSTI2020

RELAP5-3D Modeling of High Temperature Test Facility (HTTF) Test PG-26

The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled nuclear reactors with prismatic fuel and reflector blocks. Several tests have been completed, and more are still planned to at the HTTF, including depressurized conduction cooldown (DCC) and pressurized conduction cooldown (PCC) transients. This report analyses test PG-26, a progression of the Double Ended Inlet-Outlet Crossover Duct Break transient that is referred to as a DCC. PG-26 has been performed at the HTTF between May 30 and June 30, 2019. Core initial conditions (i.e., before the DCC started) have been met using low power (<100 kW) and two of ten available electric heaters. The DCC transient was initiated during the 50 th hour of the test. The break valves were opened, and hot helium from the core and cold helium from the reactor cavity simulation tank (RCST) started mixing. The gases flowed in a countercurrent fashion, where the top half of the hot duct contained hot helium that flowed in one direction and cold helium that flowed in the other direction in the bottom half of the duct. After the pressure and density reached equilibrium, the event entered a diffusion mode. The onset of a reverse natural circulation was not observed during the DCC period of the test. Version 4.4.2ie of the RELAP5-3D computer code has been used to model the HTTF PG-26 test, and results have been compared to available high-quality measured data. The model used in this study is the quality-controlled HTTF RELAP5-3D model (HTTF base 2018-04-19 QA), originally developed by P. Bayless. The report includes RELAP5-3D results of the “base calculations” as well as some sensitivities to important uncertain model inputs, such as primary helium mass flow rate, core ceramic thermal properties, as well as heat evacuation and loop friction models. Using the base RELAP5-3D model predicts a countercurrent helium flow in the hot duct observed at the beginning of the DCC, but instead of going into a molecular diffusion mode, the model predicts the onset of natural convection. Increasing friction in the core and hot duct prevents the natural convection from happening in some of the simulations. Although some temperatures are well predicted (and even overpredicted), the general tendency is to underpredict the ceramic and helium temperatures and heat removal rates during the DCC, resulting in many of the assessment findings being in minimal or insufficient agreement with the data. It is worth noting that the described discrepancies between measured data and RELAP5-3D predictions are not RELAP5-3D code limitations. More so, they reflect limitations in boundary condition and thermal property knowledge. While the RELAP5-3D calculations of the test provide some insights into what happens during the transient, and point to missing or potentially uncertain data to which the experimenters can direct their attention, the principal conclusion is that the PG-26 test data are insufficient for a system code assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Optically induced changes in the band structure of the Weyl charge-density-wave compound (TaSe 4 ) 2 I

Collective modes are responsible for the emergence of novel quantum phases in topological materials. In the quasi-one dimensional (1D) Weyl semimetal (TaSe 4 ) 2 I , a charge density wave (CDW) opens band gaps at the Weyl points, thus turning the system into an axionic insulator. Melting the CDW would restore the Weyl phase, but 1D fluctuations extend the gapped regime far above the 3D transition temperature (T CDW = 263 K), thus preventing the investigation of this topological phase transition with conventional spectroscopic methods. Here we use a non-equilibrium approach: we perturb the CDW phase by photoexcitation, and we monitor the dynamical evolution of the band structure by time- and angle-resolved photoelectron spectroscopy. We find that, upon optical excitation, electrons populate the linearly dispersing states at the Fermi level (E F ), and fill the CDW gap. The dynamics of both the charge carrier population and the band gap renormalization (BGR) show a fast component with a characteristic time scale of a few hundreds femtoseconds. However, the BGR also exhibits a second slow component on the µs time scale. The combination of an ultrafast response and of persistent changes in the spectral weight at E$_\mathrm{F}$, and the resulting sensitivity of the linearly dispersing states to optical excitations, may explain the high performances of (TaSe 4 ) 2 I as a material for broadband infrared photodetectors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-Temperature Gas-Cooled Pebble-Bed Reactors Running In And Transient Modeling Capabilities Demonstration

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor (HTGR) applications. The first part of this effort focuses on the running-in verification of Griffin's multiphysics capabilities, specifically for simulating the evolution of Pebble Bed reactor cores from startup to equilibrium. In the absence of validation data, code-to-code comparisons are conducted with Kugelpy, showing good agreement for key quantities like maximum power density and fresh core k-eigenvalue predictions. However, discrepancies in equilibrium core predictions suggest potential issues with cross sections, underscoring the need for further refinement and evaluation. The HTTF system analysis code benchmark involves RELAP5-3D, SAM, and GAMMA+ to assess their predictive capabilities for HTTF behavior under both normal operation and pressurized conduction cooldown (PCC) transient conditions. While there is good agreement in predicting major parameters such as coolant temperature, solid temperature, and flow distribution, discrepancies in transient behavior highlight differences in modeling approaches, nodalizations, and heat transfer models. The HTTF lower plenum CFD benchmark employs nekRS to simulate flow mixing phenomena, successfully capturing relevant flow physics and demonstrating mesh independence in complex geometries. Preliminary results suggest a relatively uniform temperature field but significant unsteadiness in the flow, requiring time-averaging analyses. The GPBR200 system analysis code benchmark uses SAM's core channel and porous media models, incorporating an RCCS loop for decay heat removal. During steady-state and transient conditions, including protected de-pressurized and pressurized loss of forced cooling (DLOFC and PLOFC), both models show good agreement in predicting temperature profiles and key parameters. Notably, while the core channel model underpredicts convective heat transfer effects, both models maintain temperatures well below the TRISO fuel safety limit. These benchmarking efforts collectively enhance the predictive capabilities of the tools used in HTGR design and safety analysis, guiding developments to improve their accuracy and applicability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRSDate not supplied

3D Material Response of the MSL Heatshield Using NuSil-Coated PICA

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield that used NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performance for protecting probes during planetary entry. Data from the Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. A space-grade silicone-based coating commercially known as NuSil CV-1144-0 [3] was applied to the entire MSL heatshield, including the MEDLI plugs, to mitigate the spread of dust from PICA. Modeling the thermal response of PICA-NuSil (PICA-N) system is still an open challenge. Ground testing of PICA-N models exhibited surface temperature jumps of the order of 150 K due to oxide scale formation and sub-sequent NuSil burn-off. It is therefore critical to include a validated model for the material response of the coating in engineering codes. A test campaign has been conducted at the NASA’s Langley HyMETS [4] facility to screen the response of PICA-N and gather detailed data on its behavior [5]. A first model of PICA-N thermal response has been developed using the Hy-METS experiments [6]. The objective of this work is to analyze the material response of the latest PICA-N model compared to the engineering model used to simulate the entry of MSL. The environment and material response around the MSL aeroshell during Mars atmospheric entry is simulated using a collection of tools. The Direct Simulation Monte Carlo SPARTA code [7] is used in the rarefied regime, the Data Parallel Line Relaxation (DPLR) code [8] is used in the continuum regime and radiative heating conditions are provided by the Nonequilibrium air radiation (NEQAIR) code [9] to estimate the environmental conditions. The thermal response inside the material is computed using the Porous material Analysis Toolbox based on Open-FOAM (PATO) [10,11,12]. Thermodynamic and chemistry properties are estimated using the Mutation++ library [13]. The approach implemented in PATO as a first cut PICA-N thermal response model is outlined in Figure 1. While the recession is less than the coating thickness, the Surface mass and energy balance Boundary Condition (SBC) uses the NuSil B’ tables. Once the recession removes the coating, the usual PICA B’ tables are used for the SBC. The B’ tables are computed using an equilibrium solver implemented in Mutation++, given the temperature, pressure, blowing rate, composition of the pyrolysis and environment gases, and the condensed species at the surface. Preliminary results of the 3D material response of the MSL heat-shield at the peak heating (80 sec after Entry Interface) are shown in Figure 2. Current NASA’s mission to Mars, Mars 2020, used the spare heatshield of MSL for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the PATO high-fidelity material response capability was benchmarked against flight data from MEDLI. This effort represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes. This bench-marking is awaiting the final release of the MEDLI-2 data.

Aerospace↗
NASA NTRS2023

Reconnection in a Pinch

A recently published analysis of current sheets has updated the classic Harris 1D static solution by considering multiple classes of charged particle trajectories in a generalized and dynamic current sheet. It uses a 1D PIC simulation to describe dynamic pinching and bifurcation of the current sheet. These 1D results strongly suggest that properties of the inflowing plasma, including the plasma beta, have an important effect on the equilibrium thickness of the pinched current sheet. Since 1D studies cannot describe magnetic reconnection, the time appears right to carry such 1D studies over to 2D or 3D simulations to explore current sheet thickness effects on reconnection. The Magnetospheric Multiscale Mission (MMS), with its well-resolved multipoint measurements of collisionless plasma and fields, has found that collisionless reconnection is accompanied by non-adiabatic motions of electrons that only occur in magnetic structures with a narrow scale comparable to electron inertial lengths (d e ). The recent 1D studies suggest that a plasma pinch to such scales may only occur for inflowing magnetized plasmas with relatively low plasma beta. We conclude that a parametric exploration of simulated and observed reconnection inflow conditions, particularly plasma beta, should shed light on the enablement of reconnection in collisionless plasmas.

reconnection↗
DOE OSTI2023

A Sustainable Multi-Dimensional Printable Material

Polymeric materials are usually tailored for specific functionality. A single polymer exhibiting multiple simultaneous functionalities often requires intricate molecular architecture, which is difficult to manufacture at scale because of its complex synthesis routes. In this work, a facile, partly renewable composition—prepared via reactive melt processing—that induces tunable functionalities such as 3D printability, shape recovery, and self-healing while exhibiting satisfactory mechanical properties is reported. The system with a hydrogen-bonded 3D network consists of thermally reversible nano-scale agglomerates of sustainable, rigid phenolic oligomers and crystallizable flexible polymer. Local molecular mobility and temperature-dependent relaxation and recovery of the non-equilibrium networked states enable exploiting these simultaneous functionalities. Transitions involving solidification and structure stabilization at ambient temperature spanning several hours after preheating only at 70 °C directly contrast typical thermoplastic or thermoplastic elastomer behaviors. Results from this study can inform the design of future rheology modifiers and materials for soft robotics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2010

3D Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability: Instability of a Static Column - 1

We have investigated the development of current-driven (CD) kink instability through three-dimensional relativistic MHD simulations. A static force-free equilibrium helical magnetic configuration is considered in order to study the influence of the initial configuration on the linear and nonlinear evolution of the instability. We found that the initial configuration is strongly distorted but not disrupted by the kink instability. The instability develops as predicted by linear theory. In the non-linear regime the kink amplitude continues to increase up to the terminal simulation time, albeit at different rates, for all but one simulation. The growth rate and nonlinear evolution of the CD kink instability depends moderately on the density profile and strongly on the magnetic pitch profile. The growth rate of the kink mode is reduced in the linear regime by an increase in the magnetic pitch with radius and the non-linear regime is reached at a later time than for constant helical pitch. On the other hand, the growth rate of the kink mode is increased in the linear regime by a decrease in the magnetic pitch with radius and reaches the non-linear regime sooner than the case with constant magnetic pitch. Kink amplitude growth in the non-linear regime for decreasing magnetic pitch leads to a slender helically twisted column wrapped by magnetic field. On the other hand, kink amplitude growth in the non-linear regime nearly ceases for increasing magnetic pitch.

Mizuno, Yosuke↗