Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic mechanical loading”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

An Axial-Torsional, Thermomechanical Fatigue Testing Technique

A technique for conducting strain-controlled, thermomechanical, axial-torsional fatigue tests on thin-walled tubular specimens was developed. Three waveforms of loading, namely, the axial strain waveform, the engineering shear strain waveform, and the temperature waveform were required in these tests. The phasing relationships between the mechanical strain waveforms and the temperature and axial strain waveforms were used to define a set of four axial-torsional, thermomechanical fatigue (AT-TMF) tests. Real-time test control (3 channels) and data acquisition (a minimum of 7 channels) were performed with a software program written in C language and executed on a personal computer. The AT-TMF testing technique was used to investigate the axial-torsional thermomechanical fatigue behavior of a cobalt-base superalloy, Haynes 188. The maximum and minimum temperatures selected for the AT-TMF tests were 760 and 316 C, respectively. Details of the testing system, calibration of the dynamic temperature profile of the thin-walled tubular specimen, thermal strain compensation technique, and test control and data acquisition schemes, are reported. The isothermal, axial, torsional, and in- and out-of-phase axial-torsional fatigue behaviors of Haynes 188 at 316 and 760 C were characterized in previous investigations. The cyclic deformation and fatigue behaviors of Haynes 188 in AT-TMF tests are compared to the previously reported isothermal axial-torsional behavior of this superalloy at the maximum and minimum temperatures.

Kalluri, Sreeramesh↗
NASA NTRS1996

Design and implementation of scalable tape archiver

In order to reduce costs, computer manufacturers try to use commodity parts as much as possible. Mainframes using proprietary processors are being replaced by high performance RISC microprocessor-based workstations, which are further being replaced by the commodity microprocessor used in personal computers. Highly reliable disks for mainframes are also being replaced by disk arrays, which are complexes of disk drives. In this paper we try to clarify the feasibility of a large scale tertiary storage system composed of 8-mm tape archivers utilizing robotics. In the near future, the 8-mm tape archiver will be widely used and become a commodity part, since recent rapid growth of multimedia applications requires much larger storage than disk drives can provide. We designed a scalable tape archiver which connects as many 8-mm tape archivers (element archivers) as possible. In the scalable archiver, robotics can exchange a cassette tape between two adjacent element archivers mechanically. Thus, we can build a large scalable archiver inexpensively. In addition, a sophisticated migration mechanism distributes frequently accessed tapes (hot tapes) evenly among all of the element archivers, which improves the throughput considerably. Even with the failures of some tape drives, the system dynamically redistributes hot tapes to the other element archivers which have live tape drives. Several kinds of specially tailored huge archivers are on the market, however, the 8-mm tape scalable archiver could replace them. To maintain high performance in spite of high access locality when a large number of archivers are attached to the scalable archiver, it is necessary to scatter frequently accessed cassettes among the element archivers and to use the tape drives efficiently. For this purpose, we introduce two cassette migration algorithms, foreground migration and background migration. Background migration transfers cassettes between element archivers to redistribute frequently accessed cassettes, thus balancing the load of each archiver. Background migration occurs the robotics are idle. Both migration algorithms are based on access frequency and space utility of each element archiver. To normalize these parameters according to the number of drives in each element archiver, it is possible to maintain high performance even if some tape drives fail. We found that the foreground migration is efficient at reducing access response time. Beside the foreground migration, the background migration makes it possible to track the transition of spatial access locality quickly.

Nemoto, Toshihiro↗
NASA NTRS1997

Aeroelastic Analysis of Counter Rotation Fans

Aeroelastic problems in turbomachinery and propfans can be static or dynamic in nature. The analysis of static aeroelastic problems is involved primarily with determination: (a) of the shape of the blades and the steady aerodynamic loads on the blades (which are inter-dependent), (b) of the resultant steady stresses and (c) of the static instability (divergence) margin, if applicable. In this project, we were concerned exclusively with dynamic aeroelastic behavior. The analysis of dynamic aeroelastic problems is involved with the determination: (a) of the unsteady aerodynamic loads on blades and the dynamic motion of the blades (which are again inter-dependent), (b) of the resultant dynamic stresses and their effect on fatigue life and (c) of the dynamic instability (flutter), if applicable. There are two primary dynamic aeroelastic phenomena of interest to designers of turbomachinery and propfans: flutter and forced response. Flutter generally refers to the occurrence of rapidly growing self-excited oscillations leading to catastrophic failure of the blade. When certain nonlinear phenomena are present, flutter response may lead to a potentially dangerous limit cycle oscillation rather than an immediate catastrophic failure. Forced response generally refers to the steady-state oscillations that occur as a consequence of excitations external to the rotor in question. These excitations typically result from the presence of upstream obstructions, inflow distortions, downstream obstructions, or mechanical sources such as tip-casing contact or shaft and gear meshing. Significant forced response leads to blade fatigue, and at design conditions, generally contributes to a degradation of blade life. At other operating conditions, forced response may lead to catastrophic failure due to severe blade fatigue in a short duration of time.

Keith, Theo G., Jr.↗
NASA NTRS2010

Modeling Geometry and Progressive Failure of Material Interfaces in Plain Weave Composites

A procedure combining a geometrically nonlinear, explicit-dynamics contact analysis, computer aided design techniques, and elasticity-based mesh adjustment is proposed to efficiently generate realistic finite element models for meso-mechanical analysis of progressive failure in textile composites. In the procedure, the geometry of fiber tows is obtained by imposing a fictitious expansion on the tows. Meshes resulting from the procedure are conformal with the computed tow-tow and tow-matrix interfaces but are incongruent at the interfaces. The mesh interfaces are treated as cohesive contact surfaces not only to resolve the incongruence but also to simulate progressive failure. The method is employed to simulate debonding at the material interfaces in a ceramic-matrix plain weave composite with matrix porosity and in a polymeric matrix plain weave composite without matrix porosity, both subject to uniaxial cyclic loading. The numerical results indicate progression of the interfacial damage during every loading and reverse loading event in a constant strain amplitude cyclic process. However, the composites show different patterns of damage advancement.

Hsu, Su-Yuen↗
NASA NTRS2004

Glow Discharge Plasma Demonstrated for Separation Control in the Low-Pressure Turbine

Flow separation in the low-pressure turbine (LPT) is a major barrier that limits further improvements of aerodynamic designs of turbine airfoils. The separation is responsible for performance degradation, and it prevents the design of highly loaded airfoils. The separation can be delayed, reduced, or eliminated completely if flow control techniques are used. Successful flow control technology will enable breakthrough improvements in gas turbine performance and design. The focus of this research project was the development and experimental demonstration of active separation control using glow discharge plasma (GDP) actuators in flow conditions simulating the LPT. The separation delay was shown to be successful, laying the foundation for further development of the technologies to practical application in the LPT. In a fluid mechanics context, the term "flow control" means a technology by which a very small input results in a very large effect on the flow. In this project, the interest is to eliminate or delay flow separation on LPT airfoils by using an active flow control approach, in which disturbances are dynamically inserted into the flow, they interact with the flow, and they delay separation. The disturbances can be inserted using a localized, externally powered, actuating device, examples are acoustic, pneumatic, or mechanical devices that generate vibrations, flow oscillations, or pulses. A variety of flow control devices have been demonstrated in recent years in the context of the external aerodynamics of aircraft wings and airframes, where the incoming flow is quiescent or of a very low turbulence level. However, the flow conditions in the LPT are significantly different because there are high levels of disturbances in the incoming flow that are characterized by high free-stream turbulence intensity. In addition, the Reynolds number, which characterizes the viscous forces in the flow and is related to the flow speed, is very low in the LPT passages.

Ashpis, David e.↗
DOE OSTI2024

Hydrogen effects on the deformation and slip localization in a single crystal austenitic stainless steel

Hydrogen is known to embrittle austenitic stainless steels, which are widely used in high-pressure hydrogen storage and delivery systems, but the mechanisms that lead to such material degradation are still being elucidated. The current work investigates the deformation behavior of single crystal austenitic stainless steel 316L through combined uniaxial tensile testing, characterization and atomistic simulations. Thermally precharged hydrogen is shown to increase the critical resolved shear stress (CRSS) without previously reported deviations from Schmid’s law. Molecular dynamics simulations further expose the statistical nature of the hydrogen and vacancy contributions to the CRSS in the presence of alloying. Slip distribution quantification over large in-plane distances (> 1 mm), achieved via atomic force microscopy (AFM), highlights the role of hydrogen increasing the degree of slip localization in both single and multiple slip configurations. The most active slip bands accumulate significantly more deformation in hydrogen precharged specimens, with potential implications for damage nucleation. For $\langle$110$\rangle$ tensile loading, slip localization further enhances the activity of secondary slip, increases the density of geometrically necessary dislocations and leads to a distinct lattice rotation behavior compared to hydrogen-free specimens, as evidenced by electron backscatter diffraction (EBSD) maps. Finally, the results of this study provide a more comprehensive picture of the deformation aspect of hydrogen embrittlement in austenitic stainless steels.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Dynamic Nonlinear Elastic Stability of Helicopter Rotor Blades in Hover and in Forward Flight

Equations for large coupled flap-lag motion of hingeless elastic helicopter blades are consistently derived. Only torsionally-rigid blades excited by quasi-steady aerodynamic loads are considered. The nonlinear equations of motion in the time and space variables are reduced to a system of coupled nonlinear ordinary differential equations with periodic coefficients, using Galerkin's method for the space variables. The nonlinearities present in the equations are those arising from the inclusion of moderately large deflections in the inertia and aerodynamic loading terms. The resulting system of nonlinear equations has been solved, using an asymptotic expansion procedure in multiple time scales. The stability boundaries, amplitudes of nonlinear response, and conditions for existence of limit cycles are obtained analytically. Thus, the different roles played by the forcing function, parametric excitation, and nonlinear coupling in affecting the solution can be easily identified, and the basic physical mechanism of coupled flap-lag response becomes clear. The effect of forward flight is obtained with the requirement of trimmed flight at fixed values of the thrust coefficient.

Friedmann, P.↗
DOE OSTI2020

Real-time damage characterization for GFRCs using high-speed synchrotron X-ray phase contrast imaging

Here, we report the application of high-speed synchrotron X-ray phase contrast imaging (PCI) in real-time damage characterization for glass fiber reinforced composites (GFRCs) subjected to dynamic loading. Dynamic single-edge notched bending (DSENB) experiments on pre-notched S-2 GFRCs were performed on a modified Kolsky compression bar. During loading, the synchrotron X-ray beam penetrated through the composite specimen from the side to detect damage evolution inside the material. Entire dynamic events were recorded by a high-speed camera as image sequences.0° and 90° unidirectional and cross-ply composites were investigated. An optical imaging technique was also employed to capture similar dynamic events in comparison with the radiographic imaging. It is demonstrated that high-speed X-ray PCI had sufficient phase contrast to characterize a crack initiation at a 20-μm spatial resolution within 920 ns and track the crack geometry during propagation, thereby providing reliable data to quantify the dynamic damage resistance of GFRCs. Furthermore, being capable of recognizing microscopic damage-related features at a sub-10-μm resolution, high-speed X-ray PCI provided fundamental material failure mechanisms to reveal the essential of macroscale structural failure of composites. It can also track the damage evolution inside and between individual plies of laminated composites. However, current high-speed X-ray PCI technique only supports in-situ observation and the high timing and spatial resolutions are limited within a field of view of ~2.5 mm in square, preventing its application in the three-dimensional and larger-area damage detection for GFRC structures.

42 ENGINEERING↗
DOE OSTI2025

Maximizing long-term biohydrogen production with Clostridium thermocellum for high solids conversion of lignocellulosic biomass

Biological hydrogen production from lignocellulosic biomass sustainably couples organic waste reduction with renewable energy generation. Efficient conversion is challenged by the structural complexity of lignocellulose and resulting recalcitrance to enzymatic degradation. Clostridium thermocellum natively breaks down biomass with highly effective hemi-/cellulases systems (i.e., cellulosomes) and generates hydrogen in anaerobic cultivation, creating a compelling platform for lignocellulosic biohydrogen production. Achieving commercially viable production rates requires balancing high biomass loading and throughput against uniform mixing conditions required for enzyme dispersion, pH and temperature control, and efficient hydrogen and metabolite removal in continuous operation. To address these barriers to process intensification, we implemented novel reactor and process designs for high-solids lignocellulosic biomass fermentations using the C. thermocellum KJC19-9 strain, genetically engineered for co-utilization of cellulose and hemicellulose sugars (i.e., xylose). Via computational fluid dynamics (CFD) modeling and experimental validation, we achieved a >50% improvement in biohydrogen production with an improved anchor-type impeller morphology, coupled to a threefold reduction in agitation rate. To further reduce rheological constraints and accumulation of toxic metabolites, we then transitioned the process to sequencing fed-batch operation. The resulting process generated 24.87 L H 2 L −1 from 160 g L −1 of deacetylated and mechanically refined (DMR)-pretreated corn stover biomass over 16 days while solubilizing >95% of influent cellulose and hemicellulose, setting a new performance benchmark for continuous production of biohydrogen from lignocellulose.

08 HYDROGEN↗
DOE OSTI2022

304L Can Crush Validation Studies

Accurate prediction of ductile behavior of structural alloys up to and including failure is essential in component or system failure assessment, which is necessary for nuclear weapons alteration and life extensions programs of Sandia National Laboratories. Modeling such behavior requires computational capabilities to robustly capture strong nonlinearities (geometric and material), rate- dependent and temperature-dependent properties, and ductile failure mechanisms. This study's objective is to validate numerical simulations of a high-deformation crush of a stainless steel can. The process consists of identifying a suitable can geometry and loading conditions, conducting the laboratory testing, developing a high-quality Sierra/SM simulation, and then drawing comparisons between model and measurement to assess the fitness of the simulation in regards to material model (plasticity), finite element model construction, and failure model. Following previous material model calibration, a J 2 plasticity model with a microstructural BCJ failure model is employed to model the test specimen made of 304L stainless steel. Simulated results are verified and validated through mesh and mass-scaling convergence studies, parameter sensitivity studies, and a comparison to experimental data. The converged mesh and degree of mass-scaling are the mesh discretization with 140,372 elements, and a mass scaling with a target time increment of 1.0e-6 seconds and time step scale factor of 0.5, respectively. Results from the coupled thermal-mechanical explicit dynamic analysis are comparable to the experimental data. Simulated global force vs displacement (F/D) response predicts key points such as yield, ultimate, and kinks of the experimental F/D response. Furthermore, the final deformed shape of the can and field data predicted from the analysis are similar to that of the deformed can, as measured by 3D optical CMM scans and DIC data from the experiment.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Large-scale wind turbine structures

The purpose of this presentation is to show how structural technology was applied in the design of modern wind turbines, which were recently brought to an advanced stage of development as sources of renewable power. Wind turbine structures present many difficult problems because they are relatively slender and flexible; subject to vibration and aeroelastic instabilities; acted upon by loads which are often nondeterministic; operated continuously with little maintenance in all weather; and dominated by life-cycle cost considerations. Progress in horizontal-axis wind turbines (HAWT) development was paced by progress in the understanding of structural loads, modeling of structural dynamic response, and designing of innovative structural response. During the past 15 years a series of large HAWTs was developed. This has culminated in the recent completion of the world's largest operating wind turbine, the 3.2 MW Mod-5B power plane installed on the island of Oahu, Hawaii. Some of the applications of structures technology to wind turbine will be illustrated by referring to the Mod-5B design. First, a video overview will be presented to provide familiarization with the Mod-5B project and the important components of the wind turbine system. Next, the structural requirements for large-scale wind turbines will be discussed, emphasizing the difficult fatigue-life requirements. Finally, the procedures used to design the structure will be presented, including the use of the fracture mechanics approach for determining allowable fatigue stresses.

Spera, David A.↗
DOE OSTI2022

Modeling Distributed Energy Resources for Analyzing Distribution Systems with High Renewable Penetration: Preprint

Increasing levels of inverter-based distributed energy resources (IBDERs) impact the legacy overcurrent distribution protection systems. The fault current injections of IBDERs are limited to between 1-2 p.u. of the rated current. The combination of the limited fault current and varying system load makes it increasingly difficult to reliably set the overcurrent protection. To address this challenge, researchers are proposing adaptive overcurrent protection mechanisms that can adapt to changing system conditions. Accurate modeling of IBDERs is required to develop reliable protection schemes. This paper presents a photovoltaics-based IBDER model in OpenDSS for adaptive overcurrent protection. The paper presents the validation of the proposed IBDER model against a detailed model in electromagnetic transient programs. Protection analysis of the Electric Power Research Institute J1 feeder with the proposed IBDER model is also presented.

adaptive↗
DOE OSTI2026

A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs

In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.

creep-fatigue↗
DOE OSTI2024

A recyclable self-healing composite with advanced sensing property

Polymer-based composites frequently encounter damage, often lurking beneath the surface and proving challenges to their early detection and repair. While material-based sensors show promise for encoding self-sensing properties within these composites, their in situ healing and reprocessability remain significant challenges. Therefore, the overarching goal of this study is the creation of a reprocessable polymeric composite encoded with self-healing attributes and the ability to autonomously sense damage. At the core of this innovation are vitrimers, a polymeric material characterized by a covalently adaptive dynamic network responsive to external factors such as heat. They combine thermoset-like resilience with thermoplastic-like flowability on demand under external stimuli. We nanoengineer a polyester-based vitrimeric polymer by incorporating piezoresistive carbon nanotubes (CNTs) as reinforcing elements that not only enhance its mechanical strength but also create a percolation network within the composite, thereby enabling piezoresistive self-sensing properties, all the while preserving the intrinsic self-healing capabilities offered by the vitrimeric matrix. The fabrication process of the composite involves a solvent-free in situ polymerization method that combines epoxy and anhydride-containing monomers with ~ 0.1 wt.% of CNTs. Once it was established that the introduction of CNTs into the polymeric matrix did not compromise the mechanical properties of the composite, their strain-sensing properties were characterized by applying cyclic loading while measuring their electrical resistance. Strikingly, CNT-enhanced vitrimer composite consistently retains its mechanical and sensing properties through repeated cycles of reshaping and reprocessing, underscoring its potential as a robust distributed strain sensor. This polyester-based vitrimeric composite is also easily recyclable without harsh chemical treatments. Preliminary findings from this study conclusively demonstrate that the bulk composite boasts both self-sensing capabilities and in situ detect healing properties, charting a promising course towards the development of a mechanically resilient multifunctional composite that seamlessly integrates selfsensing and healing capabilities.

Rohewal, Sargun Singh↗
NASA NTRS2006

Multiscale Analysis of Delamination of Carbon Fiber-Epoxy Laminates with Carbon Nanotubes

A multi-scale analysis is presented to parametrically describe the Mode I delamination of a carbon fiber/epoxy laminate. In the midplane of the laminate, carbon nanotubes are included for the purposes of selectively enhancing the fracture toughness of the laminate. To analyze carbon fiber epoxy carbon nanotube laminate, the multi-scale methodology presented here links a series of parameterizations taken at various length scales ranging from the atomistic through the micromechanical to the structural level. At the atomistic scale molecular dynamics simulations are performed in conjunction with an equivalent continuum approach to develop constitutive properties for representative volume elements of the molecular structure of components of the laminate. The molecular-level constitutive results are then used in the Mori-Tanaka micromechanics to develop bulk properties for the epoxy-carbon nanotube matrix system. In order to demonstrate a possible application of this multi-scale methodology, a double cantilever beam specimen is modeled. An existing analysis is employed which uses discrete springs to model the fiber bridging affect during delamination propagation. In the absence of empirical data or a damage mechanics model describing the effect of CNTs on fracture toughness, several tractions laws are postulated, linking CNT volume fraction to fiber bridging in a DCB specimen. Results from this demonstration are presented in terms of DCB specimen load-displacement responses.

Riddick, Jaret C.↗
DOE OSTI2025

In situ deformation of antigorite-olivine two-phase mixtures: Implications for dynamics and seismic anisotropy in the mantle wedge

Water released from hydrous minerals in subducting slabs reacts with the overlying plate, resulting in widespread serpentinization in the mantle wedge. Deformation of serpentinized peridotites has been invoked to explain forearc seismic anisotropy, yet studies of the mechanical properties and deformation behaviors of serpentine-bearing multiphase aggregates remain limited. Here we deformed olivine-antigorite mixtures containing 70, 50, and 20 vol.% of antigorite at 2.5 – 7.6 GPa, 673 K and strain rates of ∼10 –5 –10 –4 s –1 . Elasto-viscoplastic self-consistent simulations, constrained by synchrotron X-ray diffraction (XRD) data, were used to estimate lattice strain, stress–strain partitioning, crystallographic preferred orientations (CPO) development, and aggregate strength. Selected run products were also analyzed by electron backscatter diffraction for comparison with the CPO results obtained from XRD experiments. We found olivine transitions from A- or B-type to C-type when antigorite fraction drops to 20 vol.%, coinciding with a microstructural change from interconnected weak layers to a load-bearing framework (LBF). An additional run on a sample Atg50/Ol50 with preexisting microstructures suggested the formation of LBF was promoted by these microstructures, although the preexisting antigorite CPO has been overprinted at 20.8 % strain and could be erased completely by subsequent deformation in nature. Estimated viscosity of the two-phase mixtures suggests that low-degree serpentinization (≤20 %) in the mantle wedge may increase the strength of olivine-rich peridotite and hinder slab-mantle decoupling, whereas high-degree serpentinization (≥50–70 %) weakens the peridotite and favors decoupling if sufficient viscosity contrast (>10) develops. Seismic anisotropy shows a nonlinear dependence on antigorite fraction: antigorite CPO governs the anisotropy of the mixtures with ≥50 vol.% antigorite, whereas olivine CPO dominates at low fractions (∼20 vol.%). The presence of pre-existing microstructures reduces seismic anisotropy of the deformed mixtures, however the persistence of pre-existing CPO in actively subducting slabs remains uncertain, making their significance over geological timescales questionable.

Crystallographic preferred orientation↗
DOE OSTI2023

In-Situ X-Ray Imaging High Strain Rate Compression of Laminate Al-Graphene Composite and Mechanical Property Characterization

Owing to the opaque nature of the laminated structures, traditional high-speed optical camera cannot be used to detect the dynamic process of sub-surface deformation. In this article, we report a study of using high speed X-ray imaging to study the high strain rate deformation in laminated Al structures. We used a Kolsky bar apparatus to apply dynamic compression and a high-speed synchrotron X-ray phase contrast imaging (PCI) setup to conduct the in situ X-ray imaging study. The in situ X-ray imaging captures the shock wave propagation in the laminated structures. After shock compression, we characterized the microstructures by using transmission electron microscopy (TEM), which demonstrates an increase of dislocation density. In conclusion, the micro-pillar compression tests show that the yield strength at 0.2% offset of laminated Al-graphene composite has a significant increase of 67%, from 30 to 50 MPa, compared to laminate Al after shock loading.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗