Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesh Adaptation”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Multiscale Simulations of Magnetic Island Coalescence

We describe a new interactive parallel Adaptive Mesh Refinement (AMR) framework written in the Python programming language. This new framework, PyAMR, hides the details of parallel AMR data structures and algorithms (e.g., domain decomposition, grid partition, and inter-process communication), allowing the user to focus on the development of algorithms for advancing the solution of a systems of partial differential equations on a single uniform mesh. We demonstrate the use of PyAMR by simulating the pairwise coalescence of magnetic islands using the resistive Hall MHD equations. Techniques for coupling different physics models on different levels of the AMR grid hierarchy are discussed.

Dorelli, John C.↗
NASA NTRS2010

CME Initiation and Reconnection

Coronal mass ejections (CMEs) are the most massive explosions in the heliosphere, and the primary drivers of geoeffective space weather. This talk will be focused on fast CMEs, which travel at Alfvenic speeds as high as 2500 km/s. These ejections are associated with solar flares, prominence eruptions, and energetic particles accelerated near the Sun and in interplanetary space. CMEs require sufficient energy storage, in the form of magnetic stress, and rapid release of this energy. Although it is generally agreed that magnetic reconnection is the key to fast CME initiation, different models incorporate reconnection in different ways. One promising model --- the breakout scenario --- involves reconnection in two distinct yet interconnected locations: breakout reconnection ahead of the CME, and flare reconnect ion behind it. This model has been validated through 2D and 3D MHD simulations and favorable comparison with the observed properties of many fast CMEs. I will discuss what we have learned about the onset and evolution of breakout and flare reconnect ion from recent high-resolution 2D simulations of CME initiation with adaptive mesh refinement and numerical resistivity.

Karpen, Judy T.↗
NASA NTRS2010

Formation and Reconnection of Three-Dimensional Current Sheets in the Solar Corona

Current-sheet formation and magnetic reconnection are believed to be the basic physical processes responsible for much of the activity observed in astrophysical plasmas, such as the Sun s corona. We investigate these processes for a magnetic configuration consisting of a uniform background field and an embedded line dipole, a topology that is expected to be ubiquitous in the corona. This magnetic system is driven by a uniform horizontal flow applied at the line-tied photosphere. Although both the initial field and the driver are translationally symmetric, the resulting evolution is calculated using a fully three-dimensional magnetohydrodynamic (3D MHD) simulation with adaptive mesh refinement that resolves the current sheet and reconnection dynamics in detail. The advantage of our approach is that it allows us to apply directly the vast body of knowledge gained from the many studies of 2D reconnection to the fully 3D case. We find that a current sheet forms in close analogy to the classic Syrovatskii 2D mechanism, but the resulting evolution is different than expected. The current sheet is globally stable, showing no evidence for a disruption or a secondary instability even for aspect ratios as high as 80:1. The global evolution generally follows the standard Sweet- Parker 2D reconnection model except for an accelerated reconnection rate at a very thin current sheet, due to the tearing instability and the formation of magnetic islands. An interesting conclusion is that despite the formation of fully 3D structures at small scales, the system remains close to 2D at global scales. We discuss the implications of our results for observations of the solar corona. Subject Headings: Sun: corona Sun: magnetic fields Sun: reconnection

Edmondson, J. K.↗
NASA NTRS2011

CME Onset and Take-Off

For understanding and eventually predicting coronal mass ejections/eruptive flares, two critical questions must be answered: What is the mechanism for eruption onset, and what is the mechanism for the rapid acceleration? We address these questions in the context of the breakout model using 2.5D MHD simulations with adaptive mesh refinement (AMR). The AMR capability allowed us to achieve ultra-high numerical resolution and, thereby, determine the influence of the effective Lundquist number on the eruption. Our calculations show that, at least, for the breakout model, the onset of reconnection external to the highly sheared filament channel is the onset mechanism. Once this reconnection turns on, eruption is inevitable. However, as long as this is the only reconnection in the system, the eruption remains slow. We find that the eruption undergoes an abrupt "take-off" when the flare reconnection below the erupting plasmoid develops significant reconnection jets. We conclude that in fast CMEs, flare reconnection is the primary mechanism responsible for both flare heating and CME acceleration. We discuss the implications of these results for SDO observations and describe possible tests of the model.

Antiochos, S. K.↗
NASA NTRS2005

NASA Tech Briefs, February 2005

Topics discussed include: Instrumentation for Sensitive Gas Measurements; Apparatus for Testing Flat Specimens of Thermal Insulation; Quadrupole Ion Mass Spectrometer for Masses of 2 to 50 Da; Miniature Laser Doppler Velocimeter for Measuring Wall Shear; Coherent Laser Instrument Would Measure Range and Velocity; Printed Microinductors for Flexible Substrates; Digital Receiver for Microwave Radiometry; Printed Antennas Made Reconfigurable by Use of MEMS Switches; Traffic-Light-Preemption Vehicle-Transponder Software Module; Intersection-Controller Software Module; Central-Monitor Software Module; Estimating Effects of Multipath Propagation on GPS Signals; Parallel Adaptive Mesh Refinement Library; Predicting Noise From Aircraft Turbine-Engine Combustors; Generating Animated Displays of Spacecraft Orbits; Diagnosis and Prognosis of Weapon Systems; Training Software in Artificial-Intelligence Computing Techniques; APGEN Version 5.0; Single-Command Approach and Instrument Placement by a Robot on a Target; Three-Dimensional Audio Client Library; Isogrid Membranes for Precise, Singly Curved Reflectors; Nickel-Tin Electrode Materials for Nonaqueous Li-Ion Cells; Photocatalytic Coats in Glass Drinking-Water Bottles; Fast Laser Shutters With Low Vibratory Disturbances; Series-Connected Buck Boost Regulators; Space Physics Data Facility Web Services; Split-Resonator, Integrated-Post Vibratory Microgyroscope; Blended Buffet-Load-Alleviation System for Fighter Airplane; Gifford-McMahon/Joule-Thomson Refrigerator Cools to 2.5 K; High-Temperature, High-Load-Capacity Radial Magnetic Bearing; Fabrication of Spherical Reflectors in Outer Space; Automated Rapid Prototyping of 3D Ceramic Parts; Tissue Engineering Using Transfected Growth-Factor Genes; Automation of Vapor-Diffusion Growth of Protein Crystals; Atom Skimmers and Atom Lasers Utilizing Them; Gears Based on Carbon Nanotubes; Patched Off-Axis Bending/Twisting Actuators for Thin Mirrors; and Improving Control in a Joule-Thomson Refrigerator.

Source record↗
NASA NTRS2004

Evaluating Performance of Components

Parallel Component Performance Benchmarks is a computer program developed to aid the evaluation of the Common Component Architecture (CCA) - a software architecture, based on a component model, that was conceived to foster high-performance computing, including parallel computing. More specifically, this program compares the performances (principally by measuring computing times) of componentized versus conventional versions of the Parallel Pyramid 2D Adaptive Mesh Refinement library - a software library that is used to generate computational meshes for solving physical problems and that is typical of software libraries in use at NASA s Jet Propulsion Laboratory.

Katz, Daniel↗
NASA NTRS2011

The Angular Momentum of Baryons and Dark Matter Halos Revisited

Recent theoretical studies have shown that galaxies at high redshift are fed by cold, dense gas filaments, suggesting angular momentum transport by gas differs from that by dark matter. Revisiting this issue using high-resolution cosmological hydrodynamics simulations with adaptive-mesh refinement (AMR), we find that at the time of accretion, gas and dark matter do carry a similar amount of specific angular momentum, but that it is systematically higher than that of the dark matter halo as a whole. At high redshift, freshly accreted gas rapidly streams into the central region of the halo, directly depositing this large amount of angular momentum within a sphere of radius r = 0.1R(sub vir). In contrast, dark matter particles pass through the central region unscathed, and a fraction of them ends up populating the outer regions of the halo (r/R(sub vir) > 0.1), redistributing angular momentum in the process. As a result, large-scale motions of the cosmic web have to be considered as the origin of gas angular momentum rather than its virialised dark matter halo host. This generic result holds for halos of all masses at all redshifts, as radiative cooling ensures that a significant fraction of baryons remain trapped at the centre of the halos. Despite this injection of angular momentum enriched gas, we predict an amount for stellar discs which is in fair agreement with observations at z=0. This arises because the total specific angular momentum of the baryons (gas and stars) remains close to that of dark matter halos. Indeed, our simulations indicate that any differential loss of angular momentum amplitude between the two components is minor even though dark matter halos continuously lose between half and two-thirds of their specific angular momentum modulus as they evolve. In light of our results, a substantial revision of the standard theory of disc formation seems to be required. We propose a new scenario where gas efficiently carries the angular momentum generated by large-scale structure motions deep inside dark matter halos, redistributing it only in the vicinity of the disc.

Kimm, Taysun↗
NASA NTRS2011

Consequences of the Breakout Model for Particle Acceleration in CMEs and Flares

The largest and most efficient particle accelerators in the solar system are the giant events consisting of a fast coronal mass ejection (CME) and an intense X-class solar flare. Both flares and CMEs can produce l0(exp 32) ergs or more in nonthermal particles. Two general processes are believed to be responsible: particle acceleration at the strong shock ahead of the CME, and reconnection-driven acceleration in the flare current sheet. Although shock acceleration is relatively well understood, the mechanism by which flare reconnection produces nonthermal particles is still an issue of great debate. We address the question of CME/flare particle acceleration in the context of the breakout model using 2.5D MHD simulations with adaptive mesh refinement (AMR). The AMR capability allows us to achieve ultra-high numerical resolution and, thereby, determine the detailed structure and dynamics of the flare reconnection region. Furthermore, we employ newly developed numerical analysis tools for identifying and characterizing magnetic nulls, so that we can quantify accurately the number and location of magnetic islands during reconnection. Our calculations show that flare reconnection is dominated by the formation of magnetic islands. In agreement with many other studies, we find that the number of islands scales with the effective Lundquist number. This result supports the recent work by Drake and co-workers that postulates particle acceleration by magnetic islands. On the other hand, our calculations also show that the flare reconnection region is populated by numerous shocks and other indicators of strong turbulence, which can also accelerate particles. We discuss the implications of our calculations for the flare particle acceleration mechanism and for observational tests of the models.

Antiochos, S. K.↗
NASA NTRS2012

Numerical Predictions of Sonic Boom Signatures for a Straight Line Segmented Leading Edge Model

A sonic boom wind tunnel test was conducted on a straight-line segmented leading edge (SLSLE) model in the NASA Langley 4- by 4- Foot Unitary Plan Wind Tunnel (UPWT). The purpose of the test was to determine whether accurate sonic boom measurements could be obtained while continuously moving the SLSLE model past a conical pressure probe. Sonic boom signatures were also obtained using the conventional move-pause data acquisition method for comparison. The continuous data acquisition approach allows for accurate signatures approximately 15 times faster than a move-pause technique. These successful results provide an incentive for future testing with greatly increased efficiency using the continuous model translation technique with the single probe to measure sonic boom signatures. Two widely used NASA codes, USM3D (Navier-Stokes) and CART3D-AERO (Euler, adjoint-based adaptive mesh), were used to compute off-body sonic boom pressure signatures of the SLSLE model at several different altitudes below the model at Mach 2.0. The computed pressure signatures compared well with wind tunnel data. The effect of the different altitude for signature extraction was evaluated by extrapolating the near field signatures to the ground and comparing pressure signatures and sonic boom loudness levels.

Elmiligui, Alaa A.↗
NASA NTRS2012

Adaptive Immersed Boundary Simulations for the Launch Environment

A high-fidelity computational fluid dynamics simulation of a next generation heavy lift space vehicle during launch is presented. The purpose of the simulation is to evaluate the acoustic overpressures during ignition to permit re-design of the launch site to safely handle heavy lift vehicles. The simulation is performed using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) code, an immersed boundary block-structured Cartesian adaptive mesh refinement based solver. A verification and validation study of LAVA in the launch environment context is also performed, comparing to flight data and previous simulations of a Space Shuttle launch

Barad, Michael F.↗
NASA NTRS2012

The Mechanisms for the Onset and Explosive Eruption of Coronal Mass Ejections and Eruptive Flares

We have investigated the onset and acceleration of coronal mass ejections (CMEs) and eruptive flares. To isolate the eruption physics, our study uses the breakout model, which is insensitive to the energy buildup process leading to the eruption. We performed 2.5D simulations with adaptive mesh refinement that achieved the highest overall spatial resolution to date in a CME/eruptive flare simulation. The ultra-high resolution allows us to separate clearly the timing of the various phases of the eruption. Using new computational tools, we have determined the number and evolution of all X- and O-type nulls in the system, thereby tracking both the progress and the products of reconnection throughout the computational domain. Our results show definitively that CME onset is due to the start of fast reconnection at the breakout current sheet. Once this reconnection begins, eruption is inevitable; if this is the only reconnection in the system, however, the eruption will be slow. The explosive CME acceleration is triggered by fast reconnection at the flare current sheet. Our results indicate that the explosive eruption is caused by a resistive instability, not an ideal process. Moreover, both breakout and flare reconnections begin first as a form of weak tearing characterized by a slowly evolving plasmoids, but eventually transition to a fast form with well-defined Alfvenic reconnection jets and rapid flux transfer. This transition to fast reconnection is required for both CME onset and explosive acceleration. We discuss the key implications of our results for CME/flare observations and for theories of magnetic reconnection.

Karpen, Judith T.↗
NASA NTRS2012

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in Flash

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

Dubey, A.↗
NASA NTRS2013

ISSM: Ice Sheet System Model

In order to have the capability to use satellite data from its own missions to inform future sea-level rise projections, JPL needed a full-fledged ice-sheet/iceshelf flow model, capable of modeling the mass balance of Antarctica and Greenland into the near future. ISSM was developed with such a goal in mind, as a massively parallelized, multi-purpose finite-element framework dedicated to ice-sheet modeling. ISSM features unstructured meshes (Tria in 2D, and Penta in 3D) along with corresponding finite elements for both types of meshes. Each finite element can carry out diagnostic, prognostic, transient, thermal 3D, surface, and bed slope simulations. Anisotropic meshing enables adaptation of meshes to a certain metric, and the 2D Shelfy-Stream, 3D Blatter/Pattyn, and 3D Full-Stokes formulations capture the bulk of the ice-flow physics. These elements can be coupled together, based on the Arlequin method, so that on a large scale model such as Antarctica, each type of finite element is used in the most efficient manner. For each finite element referenced above, ISSM implements an adjoint. This adjoint can be used to carry out model inversions of unknown model parameters, typically ice rheology and basal drag at the ice/bedrock interface, using a metric such as the observed InSAR surface velocity. This data assimilation capability is crucial to allow spinning up of ice flow models using available satellite data. ISSM relies on the PETSc library for its vectors, matrices, and solvers. This allows ISSM to run efficiently on any parallel platform, whether shared or distrib- ISSM: Ice Sheet System Model NASA's Jet Propulsion Laboratory, Pasadena, California uted. It can run on the largest clusters, and is fully scalable. This allows ISSM to tackle models the size of continents. ISSM is embedded into MATLAB and Python, both open scientific platforms. This improves its outreach within the science community. It is entirely written in C/C++, which gives it flexibility in its design, and the power/speed that C/C++ allows. ISSM is svn (subversion) hosted, on a JPL repository, to facilitate its development and maintenance. ISSM can also model propagation of rifts using contact mechanics and mesh splitting, and can interface to the Dakota software. To carry out sensitivity analysis, mesh partitioning algorithms are available, based on the Scotch, Chaco, and Metis partitioners that ensure equal area mesh partitions can be done, which are then usable for sampling and local reliability methods.

Larour, Eric↗
NASA NTRS2014

Numerical Simulations For the F-16XL Aircraft Configuration

Numerical simulations of flow around the F-16XL are presented as a contribution to the Cranked Arrow Wing Aerodynamic Project International II (CAWAPI-II). The NASA Tetrahedral Unstructured Software System (TetrUSS) is used to perform numerical simulations. This CFD suite, developed and maintained by NASA Langley Research Center, includes an unstructured grid generation program called VGRID, a postprocessor named POSTGRID, and the flow solver USM3D. The CRISP CFD package is utilized to provide error estimates and grid adaption for verification of USM3D results. A subsonic high angle-of-attack case flight condition (FC) 25 is computed and analyzed. Three turbulence models are used in the calculations: the one-equation Spalart-Allmaras (SA), the two-equation shear stress transport (SST) and the kε turbulence models. Computational results, and surface static pressure profiles are presented and compared with flight data. Solution verification is performed using formal grid refinement studies, the solution of Error Transport Equations, and adaptive mesh refinement. The current study shows that the USM3D solver coupled with CRISP CFD can be used in an engineering environment in predicting vortex-flow physics on a complex configuration at flight Reynolds numbers.

Elmiligui, Alaa A.↗
NASA NTRS2014

Parallel Adaptive High-Order CFD Simulations Characterizing Cavity Acoustics for the Complete SOFIA Aircraft

This paper presents one-of-a-kind MPI-parallel computational fluid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft of a Boeing 747SP. These simulations focus on how the unsteady flow field inside and over the cavity interferes with the optical path and mounting of the telescope. A temporally fourth-order Runge-Kutta, and spatially fifth-order WENO-5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh refinement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32,000 cores and 4 billion cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregularities caused by the highly complex geometry. Limits to scaling beyond 32K cores are identified, and targeted code optimizations are discussed.

Acoustics↗
NASA NTRS2015

Los Alamos RAGE Simulations of the HAIV Mission Concept

The mitigation of potentially hazardous objects (PHOs) can be accomplished by a variety of methods including kinetic impactors, gravity tractors and several nuclear explosion options. Depending on the available lead time prior to Earth impact, non- nuclear options can be very effective at altering a PHOs orbit. However if the warning time is short nuclear options are generally deemed most effective at mitigating the hazard. The NIAC mission concept for a nuclear mission has been presented at several meetings, including the last PDC (2013).We use the adaptive mesh hydrocode RAGE to perform detailed simulations of this Hypervelocity Asteroid Intercept Vehicle (HAIV) mission concept. We use the RAGE code to simulate the crater formation by the kinetic impactor as well as the explosion and energy coupling from the follower nuclear explosive device (NED) timed to detonate below the original surface to enhance the energy coupling. The RAGE code has been well validated for a wide variety of applications. A parametric study will be shown of the energy and momentum transfer to the target 100 m diameter object: 1) the HAIV mission as planned; 2) a surface explosion and 3) a subsurface (contained) explosion; both 2) and 3) use the same source energy as 1).Preliminary RAGE simulations show that the kinetic impactor will carve out a surface crater on the object and the subsequent NED explosion at the bottom of the crater transfers energy and momentum to the target effectively moving it off its Earth crossing orbit. Figure 1 shows the initial (simplified) RAGE 2D setup geometry for this study. Figure 2 shows the crater created by the kinetic impactor and Figure 3 shows the time sequence of the energy transfer to the target by the NED.

Asteroid Mitigation↗
NASA NTRS2015

Simulation of Atmospheric-Entry Capsules in the Subsonic Regime

The accuracy of Computational Fluid Dynamics predictions of subsonic capsule aerodynamics is examined by comparison against recent NASA wind-tunnel data at high-Reynolds-number flight conditions. Several aspects of numerical and physical modeling are considered, including inviscid numerical scheme, mesh adaptation, rough-wall modeling, rotation and curvature corrections for eddy-viscosity models, and Detached-Eddy Simulations of the unsteady wake. All of these are considered in isolation against relevant data where possible. The results indicate that an improved predictive capability is developed by considering physics-based approaches and validating the results against flight-relevant experimental data.

Simulation↗
NASA NTRS2016

Parallel Adjective High-Order CFD Simulations Characterizing SOFIA Cavity Acoustics

This paper presents large-scale MPI-parallel computational uid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft fuselage of a Boeing 747SP. These simulations focus on how the unsteady ow eld inside and over the cavity interferes with the optical path and mounting structure of the telescope. A temporally fourth-order accurate Runge-Kutta, and spatially fth-order accurate WENO- 5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh re nement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32k CPU cores and 4 billion compu- tational cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregular numerical cost associated with blocks con- taining boundaries. Limits to scaling beyond 32k cores are identi ed, and targeted code optimizations are discussed.

CFD↗