Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear embedding”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Spatio-Temporal Clustering of Seismicity Enabled by Off-Fault Plasticity

While significant progress has been made in understanding earthquake source processes in linear elastic domains, the effect of more realistic rheologies including plasticity is poorly understood. Here, in this work, we simulate the sequence of earthquake and aseismic slip of a 2D antiplane rate-and-state fault embedded in a full-space elastic-plastic bulk. We show that off-fault plasticity may lead to partial ruptures as well as temporal clustering of seismic events. Furthermore, the interaction of fault slip and off-fault plasticity results in pockets of slip deficit. While the energy dissipated through plastic deformation remains a small fraction of the total energy budget, its impact on the source characteristics is disproportionally large through the redistribution of stresses and viscous relaxation. Our results suggest a new mechanism of dynamic heterogeneity in earthquake physics that may have important implications on earthquake size distribution and energy budget.

58 GEOSCIENCES↗
DOE OSTI2026

Topology-Dependent Performance of Free-Space Photonic Quantum Networks Under Noise

Photonic quantum communication enables secure and high-fidelity information transfer beyond classical limits, with direct relevance to emerging quantum networks operating in free-space environments. While physical-layer models of depolarizing noise, Gamma–Gamma turbulence statistics, entanglement swapping, and decoy-state QKD security bounds are individually well established, prior work typically treats these components in isolation or under fixed network assumptions. In this work, we develop a unified topology-aware analytical framework that simultaneously integrates free-space optical link budgets, turbulence-induced visibility degradation, depolarizing qubit noise, multi-hop entanglement cascade dynamics, teleportation fidelity thresholds, CHSH nonlocality certification, and asymptotic decoy-state secret key rate bounds across star, mesh, and ring graph structures. Rather than introducing new physical channel models, we demonstrate that identical physical links exhibit fundamentally different end-to-end performance once embedded within different network topologies. Mesh architectures minimize visibility cascade through hop-count reduction but incur quadratic hardware scaling. Star topologies minimize link count but concentrate noise and synchronization overhead at the hub. Ring configurations offer linear hardware scaling with multiplicative fidelity degradation. The results establish topology as a first-order design parameter in near-term free-space quantum networks operating without full quantum repeater infrastructures. While motivated by distributed multi-agent architectures, the framework applies broadly to terrestrial, airborne, and satellite-assisted photonic quantum communication systems.

QKD↗
NASA NTRS1980

Analytic and approximate boundary-fitted coordinate systems for fluid flow simulation

A method for constructing a direct functional relation between the computational and physical domains such that boundaries in the computational domain map into boundaries in the physical domain is presented. This method which is called the 'two-boundary technique' is based on generating functions that represent two disconnected boundaries in the physical domain and an independent function which connects one boundary with another to generate grids. The suggested connecting functions are simple linear and cubic parametric polynomials with parameters based on position and derivatives at the disconnected boundaries. The boundaries can be defined by analytical, approximate, and tension spline functions; the grid density in specified regions is controlled by embedding contraction functions in either the boundary representation or the connecting function.

Smith, R. E.↗
NASA NTRS2018

Electrical Properties of Unidirectional Carbon Fiber Reinforced Epoxy-Graphene Composites

This paper aims at establishing connections between electrical resistivity and stresses in carbon fiber reinforced epoxy-graphene composites. To avoid the interconnection of multiple parameters, we examined composite specimens consisting of a single strand or tow (12,000 carbon fiber filaments) embedded in an epoxy matrix containing up to 1.0 wt% of graphene (TG679) particles. The specimens show resistance increase with the tensile stress. The characteristic of the dependence varies with the content of graphene particles, tending to be linear as the percentage of graphene content increases. We did not observe any significant changes in the ultimate strength of the composite with an increase of the graphene particles concentration.

Daniel Wentzel↗
NASA NTRS2019

Evaluation of System Reliability and Heavy Ion Angular Effects

Microsemi (Microchip) RTG4 embedded triple modular redundant (TMR) phase-locked-loop (PLL) SEU data is presented. SEU data analysis includes: 1) Evaluation of heavy-ion beam angular effects (rectangular parallel pipe (RPP) or no RPP), 2) Importance of finding linear energy transfer (LET) onset (L0), 3) Comparison of prediction rate techniques.

Space radiation environments↗
DOE OSTI2020

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗
DOE OSTI2023

Star Block Copolymers at Homopolymer Interfaces: Conformation and Compatibilization

Star block copolymers (s-BCPs), comprised of multiple linear diblock copolymers joined at a central point, are shown to segregate to the interface between two immiscible homopolymers that are identical to the blocks of the s-BCPs. The s-BCPs undergo a configurational transition at the interface, with different blocks of copolymers being embedded in their respective homopolymers, thereby bridging the interface and promoting adhesion. A series of 4-arm s-BCPs were synthesized with hydrogenated or deuterated polystyrene (PS/dPS) as the core block and poly(2-vinylpyridine) (P2VP) as the corona block, which was directly placed at the interface between the two homopolymers. Neutron reflectivity (NR) was used to determine the concentration profiles of the PS homopolymer, s-BCP core blocks, and P2VP total segments under equilibrium. The investigation varies the molecular weight (MW) and the total number of s-BCPs at the interface. Self-consistent-field theory (SCFT) was also employed to calculate the concentration profiles of the components at the interface, which were in excellent agreement with experimental results. The NR showed that the interfacial width between the homopolymers increased with the increasing number of s-BCPs at the interface up to a saturation limit. Beyond this limit, additional s-BCPs were released into the corona-miscible phase as unimolecular micelles. For a comparable interlayer thickness of s-BCPs at the interface, lower MW s-BCPs generated a broader interface. SCFT analysis suggested that, at the same packing density, the arms of the low MW s-BCPs align more parallel to the interface, while the arms of high MW s-BCPs adopt a more normal orientation, like their linear BCP counterparts. Furthermore, it was also observed that the core blocks, constrained by the junction points, were oriented more parallel and closer to the interface than the corona blocks. The phase behavior of the polymer blends revealed that s-BCP additives can efficiently reduce the domain size, with the low MW yielding smaller domain sizes due to the greater reduction in the interfacial energy and the high MW arresting phase separation due to their higher binding energy and a jamming of the interfacial assemblies. Asymmetric double cantilever beam (ADCB) tests demonstrated that s-BCPs promoted adhesion more efficiently than their linear BCP counterparts due to stronger binding energy per molecule, suggesting a more efficient compatibilizer for polymer upcycling. In conclusion, the results from these studies provide fundamental insights into the assembly of s-BCPs at homopolymer interfaces, the reduction of domain size, and promotion of adhesion, providing a strategy for the use of s-BCPs as stealth surfactants and universal compatibilizers.

interface↗
DOE OSTI2021

Materials representation and transfer learning for multi-property prediction

The adoption of machine learning in materials science has rapidly transformed materials property prediction. Hurdles limiting full capitalization of recent advancements in machine learning include the limited development of methods to learn the underlying interactions of multiple elements as well as the relationships among multiple properties to facilitate property prediction in new composition spaces. To address these issues, we introduce the Hierarchical Correlation Learning for Multi-property Prediction (H-CLMP) framework that seamlessly integrates: (i) prediction using only a material's composition, (ii) learning and exploitation of correlations among target properties in multi-target regression, and (iii) leveraging training data from tangential domains via generative transfer learning. The model is demonstrated for prediction of spectral optical absorption of complex metal oxides spanning 69 three-cation metal oxide composition spaces. H-CLMP accurately predicts non-linear composition-property relationships in composition spaces for which no training data are available, which broadens the purview of machine learning to the discovery of materials with exceptional properties. This achievement results from the principled integration of latent embedding learning, property correlation learning, generative transfer learning, and attention models. The best performance is obtained using H-CLMP with transfer learning [H-CLMP(T)] wherein a generative adversarial network is trained on computational density of states data and deployed in the target domain to augment prediction of optical absorption from composition. H-CLMP(T) aggregates multiple knowledge sources with a framework that is well suited for multi-target regression across the physical sciences.

36 MATERIALS SCIENCE↗
NASA NTRS2004

New Effective Material Couple--Oxide Ceramic and Carbon Nanotube-- Developed for Aerospace Microsystem and Micromachine Technologies

The prime driving force for using microsystem and micromachine technologies in transport vehicles, such as spacecraft, aircraft, and automobiles, is to reduce the weight, power consumption, and volume of components and systems to lower costs and increase affordability and reliability. However, a number of specific issues need to be addressed with respect to using microsystems and micromachines in aerospace applications--such as the lack of understanding of material characteristics; methods for producing and testing the materials in small batches; the limited proven durability and lifetime of current microcomponents, packaging, and interconnections; a cultural change with respect to system designs; and the use of embedded software, which will require new product assurance guidelines. In regards to material characteristics, there are significant adhesion, friction, and wear issues in using microdevices. Because these issues are directly related to surface phenomena, they cannot be scaled down linearly and they become increasingly important as the devices become smaller. When microsystems have contacting surfaces in relative motion, the adhesion and friction affect performance, energy consumption, wear damage, maintenance, lifetime and catastrophic failure, and reliability. Ceramics, for the most part, do not have inherently good friction and wear properties. For example, coefficients of friction in excess of 0.7 have been reported for ceramics and ceramic composite materials. Under Alternate Fuels Foundation Technologies funding, two-phase oxide ceramics developed for superior high-temperature wear resistance in NASA's High Operating Temperature Propulsion Components (HOTPC) project and new two-layered carbon nanotube (CNT) coatings (CNT topcoat/iron bondcoat/quartz substrate) developed in NASA's Revolutionary Aeropropulsion Concepts (RAC) project have been chosen as a materials couple for aerospace applications, including micromachines, in the nanotechnology lubrication task because of their potential for superior friction and wearf properties in air and in an ultrahigh vacuum, spacelike environment. At the NASA Glenn Research Center, two-phase oxide ceramic eutectics, Al2O3/ZrO2(Y2O3), were directionally solidified using the laser-float-zone process, and carbon nanotubes were synthesized within a high-temperature tube furnace at 800 C. Physical vapor deposition was used to coat all quartz substrates with 5-nm-thick iron as catalyst and bondcoat, which formed iron islands resembling droplets and serving as catalyst particles on the quartz. A series of scanning electron micrographs showing multiwalled carbon nanotubes directionally grown as aligned "nanograss" on quartz is presented. Unidirectional sliding friction eperiments were conducted at Glenn with the two-layered CNT coatings in contact with the two-phase Al2O3/ZrO2(Y2O3) eutectics in air and in ultrachigh vacuum. The main criteria for judging the performance of the materials couple for solid lubrication and antistick applications in a space environment were the coefficient of friction and the wear resistance (reciprocal of wear rate), which had to be less than 0.2 and greater than 10(exp 5) N(raised dot)/cubic millimetes, respectively, in ultrahigh vacuum. In air, the coefficient of friction for the CNT coatings in contact with Al2O3/ZrO2 (Y2O3) eutectics was 0.04, one-fourth of that for quartz. In an ultrahigh vacuum, the coefficient of friction for CNT coatings in contact with Al2O3/ZrO2 (Y2O3) was one-third of that for quartz. The two-phase Al2O3/ZrO2 (Y2O3) eutectic coupled with the two-layered CNT coating met the coefficient of friction and wear resistance criteria both in air and in an ultrahigh vacuum, spacelike environment. This material's couple can dramatically improve the stiction (or adhesion), friction, and wear resistance of the contacting surfaces, which are major issues for microdevices and micromachines.

Miyoshi, Kazuhisa↗
DOE OSTI2021

Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction

Stochastic density functional theory (sDFT) is becoming a valuable tool for studying ground-state properties of extended materials. The computational complexity of describing the Kohn-Sham orbitals is replaced by introducing a set of random (stochastic) orbitals leading to linear and often sub-linear scaling of certain ground-state observables at the account of introducing a statistical error. Schemes to reduce the noise are essential, for example, for determining the structure using the forces obtained from sDFT. Recently, we have introduced two embedding schemes to mitigate the statistical fluctuations in the electron density and resultant forces on the nuclei. Both techniques were based on fragmenting the system either in real space or slicing the occupied space into energy windows, allowing for a significant reduction in the statistical fluctuations. For chemical accuracy, further reduction of the noise is required, which could be achieved by increasing the number of stochastic orbitals. However, the convergence is relatively slow as the statistical error scales as 1 / N χ according to the central limit theorem, where N χ is the number of random orbitals. In this paper, we combined the embedding schemes mentioned above and introduced a new approach that builds on overlapped fragments and energy windows. The new approach significantly lowers the noise for ground-state properties, such as the electron density, total energy, and forces on the nuclei, as demonstrated for a G-center in bulk silicon.

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

Computation of multi-dimensional viscous supersonic flow

A method has been developed for two- and three-dimensional computations of viscous supersonic jet flows interacting with an external flow. The approach employs a reduced form of the Navier-Stokes equations which allows solution as an initial-boundary value problem in space, using an efficient noniterative forward marching algorithm. Numerical instability associated with forward marching algorithms for flows with embedded subsonic regions is avoided by approximation of the reduced form of the Navier-Stokes equations in the subsonic regions of the boundary layers. Supersonic and subsonic portions of the flow field are simultaneously calculated by a consistently split linearized block implicit computational algorithm. The results of computations for a series of test cases associated with supersonic jet flow is presented and compared with other calculations for axisymmetric cases. Demonstration calculations indicate that the computational technique has great promise as a tool for calculating a wide range of supersonic flow problems including jet flow. Finally, a User's Manual is presented for the computer code used to perform the calculations.

Buggeln, R. C.↗
NASA NTRS1993

The embedded young stars in the Taurus-Auriga molecular cloud. II - Models for scattered light images

We describe NIR imaging observations of embedded young stars in the Taurus-Auriga molecular cloud. We find a large range in J-K and H-K colors for these class I sources. The bluest objects have colors similar to the reddest T Tauri stars in the cloud; redder objects lie slightly above the reddening line for standard ISM dust and have apparent K extinctions of up to 5 mag. Most of these sources also show extended NIR emission on scales of 10-20 arcsec which corresponds to linear sizes of 1500-3000 AU. The NIR colors and nebular morphologies for this sample and the magnitude of linear polarization in several sources suggest scattered light produces most of the NIR emission in these objects. We present modeling results that suggest mass infall rates that agree with predictions for cold clouds and are generally consistent with rates estimated from radiative equilibrium models. For reasonable dust grain parameters, the range of colors and extinctions require flattened density distributions with polar cavities evacuated by bipolar outflows. These results support the idea that infall and outflow occur simultaneously in deeply embedded bipolar outflow sources. The data also indicate fairly large centrifugal radii and large inclinations to the rotational axis for a typical source.

Kenyon, Scott J.↗
DOE OSTI2022

Exact quantum scars in the chiral nonlinear Luttinger liquid

While the chiral linear Luttinger liquid is integrable via bosonization, its nonlinear counterpart does not admit for an analytic solution. In this work, we find a subextensive number of exact eigenstates for a large family of density-density interaction terms. These states are embedded in a continuum of strongly correlated excited states. The real-space entanglement entropy of some exact states scales logarithmically with system size while that of others has volume-law scaling. We introduce momentum-space entanglement as an unambiguous differentiator between these exact states and the remaining excited states. With regard to momentum space, the exact states behave as bona fide quantum many-body scars: they exhibit identically zero momentum-space entanglement, while typical eigenstates behave thermally. We corroborate this finding by a level statistics analysis. Here we detail the general formalism for systematically finding all interaction terms and associated exact states, and present a number of infinite exact state sequences extending to arbitrarily high energies. Unlike many previous examples of quantum many-body scars, the exact states uncovered here do not lie at equidistant energies and do not follow from a special operator algebra. Instead, they are uniquely enabled by the interplay of Fermi statistics and chirality.

36 MATERIALS SCIENCE↗
NASA NTRS1987

A computational code for two-dimensional unsteady magnetohydrodynamics by the method of characteristics

A computational code for solving two-dimensional, time-dependent MHD equations by the method of characteristics is presented. Its capabilities are demonstrated by solving two very different problems for which analytical solutions exist: linearized, standing MHD wave motions in a magnetized cylindrical plasma, and nonlinear self-similar expansion of a magnetized plasma ball. The nonlinear development of standing MHD wave solutions in a cylindrical plasma is also studied. The method can be naturally embedded in the computational architecture of massively parallel processors.

Lou, Y. Q.↗
NASA NTRS2009

Aerodynamic and Aeroelastic Characteristics of a Tension Cone Inflatable Aerodynamic Decelerator

The supersonic aerodynamic and aeroelastic characteristics of a tension cone inflatable aerodynamic decelerator were investigated by wind tunnel testing. Two sets of tests were conducted: one using rigid models and another using textile models. Tests using rigid models were conducted over a Mach number range from 1.65 to 4.5 at angles of attack from -12 to 20 degrees. The axial, normal, and pitching moment coefficients were found to be insensitive to Mach number over the tested range. The axial force coefficient was nearly constant (C(sub A) = 1.45 +/- 0.05) with respect to angle of attack. Both the normal and pitching moment coefficients were nearly linear with respect to angle of attack. The pitching moment coefficient showed the model to be statically stable about the reference point. Schlieren images and video showed a detached bow shock with no evidence of large regions of separated flow and/or embedded shocks at all Mach numbers investigated. Qualitatively similar static aerodynamic coefficient and flow visualization results were obtained using textile models at a Mach number of 2.5. Using inflatable textile models the torus pressure required to maintain the model in the fully-inflated configuration was determined. This pressure was found to be sensitive to details in the structural configuration of the inflatable models. Additional tests included surface pressure measurements on rigid models and deployment and inflation tests with inflatable models.

Clark, Ian G.↗
NASA NTRS2012

NASA Tech Briefs, September 2012

Topics covered include: Beat-to-Beat Blood Pressure Monitor; Measurement Techniques for Clock Jitter; Lightweight, Miniature Inertial Measurement System; Optical Density Analysis of X-Rays Utilizing Calibration Tooling to Estimate Thickness of Parts; Fuel Cell/Electrochemical Cell Voltage Monitor; Anomaly Detection Techniques with Real Test Data from a Spinning Turbine Engine-Like Rotor; Measuring Air Leaks into the Vacuum Space of Large Liquid Hydrogen Tanks; Antenna Calibration and Measurement Equipment; Glass Solder Approach for Robust, Low-Loss, Fiber-to-Waveguide Coupling; Lightweight Metal Matrix Composite Segmented for Manufacturing High-Precision Mirrors; Plasma Treatment to Remove Carbon from Indium UV Filters; Telerobotics Workstation (TRWS) for Deep Space Habitats; Single-Pole Double-Throw MMIC Switches for a Microwave Radiometer; On Shaft Data Acquisition System (OSDAS); ASIC Readout Circuit Architecture for Large Geiger Photodiode Arrays; Flexible Architecture for FPGAs in Embedded Systems; Polyurea-Based Aerogel Monoliths and Composites; Resin-Impregnated Carbon Ablator: A New Ablative Material for Hyperbolic Entry Speeds; Self-Cleaning Particulate Prefilter Media; Modular, Rapid Propellant Loading System/Cryogenic Testbed; Compact, Low-Force, Low-Noise Linear Actuator; Loop Heat Pipe with Thermal Control Valve as a Variable Thermal Link; Process for Measuring Over-Center Distances; Hands-Free Transcranial Color Doppler Probe; Improving Balance Function Using Low Levels of Electrical Stimulation of the Balance Organs; Developing Physiologic Models for Emergency Medical Procedures Under Microgravity; PMA-Linked Fluorescence for Rapid Detection of Viable Bacterial Endospores; Portable Intravenous Fluid Production Device for Ground Use; Adaptation of a Filter Assembly to Assess Microbial Bioburden of Pressurant Within a Propulsion System; Multiplexed Force and Deflection Sensing Shell Membranes for Robotic Manipulators; Whispering Gallery Mode Optomechanical Resonator; Vision-Aided Autonomous Landing and Ingress of Micro Aerial Vehicles; Self-Sealing Wet Chemistry Cell for Field Analysis; General MACOS Interface for Modeling and Analysis for Controlled Optical Systems; Mars Technology Rover with Arm-Mounted Percussive Coring Tool, Microimager, and Sample-Handling Encapsulation Containerization Subsystem; Fault-Tolerant, Real-Time, Multi-Core Computer System; Water Detection Based on Object Reflections; SATPLOT for Analysis of SECCHI Heliospheric Imager Data; Plug-in Plan Tool v3.0.3.1; Frequency Correction for MIRO Chirp Transformation Spectroscopy Spectrum; Nonlinear Estimation Approach to Real-Time Georegistration from Aerial Images; Optimal Force Control of Vibro-Impact Systems for Autonomous Drilling Applications; Low-Cost Telemetry System for Small/Micro Satellites; Operator Interface and Control Software for the Reconfigurable Surface System Tri-ATHLETE; and Algorithms for Determining Physical Responses of Structures Under Load.

Source record↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
NASA NTRS1993

Three-dimensional magnetohydrodynamics of the emerging magnetic flux in the solar atmosphere

The nonlinear evolution of an emerging magnetic flux tube or sheet in the solar atmosphere is studied through 3D MHD simulations. In the initial state, a horizontal magnetic flux sheet or tube is assumed to be embedded at the bottom of MHD two isothermal gas layers, which approximate the solar photosphere/chromosphere and the corona. The magnetic flux sheet or tube is unstable against the undular mode of the magnetic buoyancy instability. The magnetic loop rises due to the linear and then later nonlinear instabilities caused by the buoyancy enhanced by precipitating the gas along magnetic field lines. We find by 3D simulation that during the ascendance of loops the bundle of flux tubes or even the flux sheet develops into dense gas filaments pinched between magnetic loops. The interchange modes help produce a fine fiber flux structure perpendicular to the magnetic field direction in the linear stage, while the undular modes determine the overall buoyant loop structure. The expansion of such a bundle of magnetic loops follows the self-similar behavior observed in 2D cases studied earlier. Our study finds the threshold flux for arch filament system (AFS) formation to be about 0.3 x 10 exp 20 Mx.

Matsumoto, R.↗