Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesh Adaptivity”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Ground Tests on the BOLT Geometry at Mach 7: Cross-Facility Comparison and Stability Analysis

Transition measurements have been obtained through two experimental campaigns conducted independently by the German Aerospace Center and the French Aerospace Lab with the French Alternative Energies and Atomic Energy Commission on subscale models of the BOLT-1 flight experiment geometry. This paper details a cross-facility comparison of measurements obtained at Mach 7, as well as subsequent computational analysis. Infrared (IR) thermography measurements obtained by both campaigns have facilitated a global comparison of the transition front across facilities at analagous conditions, which are found to be in good agreement. High-frequency surface pressure fluctuation data demonstrate significant amplification of instabilities with Mack-mode characteristics in the outboard regions of the acreage. These measurements are compared to stability analyses of varying fidelity. The computational methods employed to characterize the boundary-layer transition phenomena include the traditional line-marching implementation of the Parabolized Stability Equations (PSE), 2D eigenvalue analysis coupled with PSE, and Adaptive Mesh Refinement Wavepacket Tracking. While the 2D eigenvalue analysis is found to predict instabilities, which correlate in terms of frequency and acreage location to the experimental measurements, the predicted amplification for these instabilities is lower than would typically be expected for transition. Line-marching results for traveling crossflow produce the best match to the experimental transition front, with a consistent transition N factor of approximately 3-3.5.

Boundary Layer Transition↗
NASA NTRS2024

Computational Techniques to Generate Space Launch System Aerodynamic Databases

This document describes the reasoning and trade studies used to evaluate tools for constructing the aerodynamic lineload databases for the liftoff and transition phases of flight for the Space Launch System. Three computational fluid dynamics codes (USM3D, FUN3D, Kestrel) were investigated with various turbulence models, as well as detached eddy simulation variants for the launch vehicle in free air and in proximity to the tower. Decisions were made mostly based on results from brief developmental studies performed in response to specific, unforeseen challenges that were encountered in the analysis of a given configuration. The need to develop databases in a timely manner, as well as accurately capture the expected leeward-wake flowfield characteristics, led to the selection of the Kestrel flow solver with its delayed detached eddy simulation method, the Spalart-Allmaras turbulence model, and the adaptive mesh refinement capability in the off-body Cartesian grid region.

Karen A. Deere↗
NASA NTRS2024

FUN3D Manual: 14.1

This manual describes the installation and execution of FUN3D version 14.1, including optional dependent packages. FUN3D is a suite of computational fluid dynamics simulation and design tools that uses mixed-element unstructured grids in a large number of formats, including structured multiblock and overset grid systems. A discretely-exact adjoint solver may be used for formal design optimization, error estimation, and mesh adaptation. FUN3D also offers a reacting, real-gas capability and provides GPU acceleration of many common simulation options.

William K Anderson↗
NASA NTRSDate not supplied

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept during Supersonic Retropropulsion

Simulations of unsteady supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle were performed using the OVERFLOW Computational Fluid Dynamics (CFD) solver. High-fidelity flow solver techniques, including Detached Eddy Simulation (DES) turbulence modeling and Adaptive Mesh Refinement (AMR), were employed to obtain improved realism in CFD predictions. Simulation conditions and geometry configurations were designed to match specific runs in the Descent System Study (DSS) wind tunnel testing (WTT) campaign. The accuracy of each simulation is assessed by direct comparison to experimental data. Comparisons of computational predictions of the SRP flowfield and bow shock shape to experimental schlieren imaging show reasonable prediction of mean shock shape, with approximately 10% similarity in shock standoff distance for selected conditions, as well as similarity in local, time-varying fluctuations of the shock-plume interaction. Comparisons of discrete measurements of surface pressure coefficient (Cp) indicate CFD accuracy within approximately 10% of the experiment across the majority of the model heatshield, with larger variations at some of the heatshield edge locations with stronger flow unsteadiness. Simulated unsteadiness of these chaotic flows, which were highly dynamic and multi-modal, was shown to be within 20-40% of experimentally-measured pressure standard deviation (SD) for the majority of the sampled locations.

Supersonic Retropropulsion↗
NASA NTRSDate not supplied

Enhanced Simulation Techniques in Predicting Sonic Boom Loudness Using CFD

This paper outlines advancements in predicting sonic boom loudness within the Launch, Ascent, and Vehicle Aerodynamics (LAVA) computational framework. Traditionally, a two step process consisting of a steady state computational fluid dynamics problem for near-field analysis and a far-field propagation solver for calculation of loudness metrics has been used. Improvements to this process made in this work include utilizing a high-order space marching method for mid-field computations, developing a novel output-based mesh adaptation method targeting error in near-field pressure sig-natures, and developing a robust scripting system using curvilinear grids to increase robustness and simplify the process of running large databases of simulation cases. These advancements are detailed and applied to the simulation of the X-59, presenting comparative cost and timing analyses between the prior two step workflow and the current three step procedure. We achieve increased accuracy and robustness for loudness predictions with at least a50%computational cost reduction.

CST↗
NASA NTRSDate not supplied

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept during Supersonic Retropropulsion

Simulations of unsteady supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle were performed using the OVERFLOW Computational Fluid Dynamics (CFD) solver. High-fidelity flow solver techniques, including Detached Eddy Simulation (DES) turbulence modeling and Adaptive Mesh Refinement (AMR), were employed to obtain improved realism in CFD predictions. Simulation conditions and geometry configurations were designed to match specific runs in the Descent System Study (DSS) wind tunnel testing (WTT) campaign. The accuracy of each simulation is assessed by direct comparison to experimental data. Comparisons of computational predictions of the SRP flowfield and bow shock shape to experimental schlieren imaging show reasonable prediction of mean shock shape, with approximately 10% similarity in shock standoff distance for selected conditions, as well as similarity in local, time-varying fluctuations of the shock-plume interaction. Comparisons of discrete measurements of surface pressure coefficient (Cp) indicate CFD accuracy within approximately 10% of the experiment across the majority of the model heatshield, with larger variations at some of the heatshield edge locations with stronger flow unsteadiness. Simulated unsteadiness of these chaotic flows, which were highly dynamic and multi-modal, was shown to be within 20-40% of experimentally-measured pressure standard deviation (SD) for the majority of the sampled locations.

Supersonic Retropropulsion↗
NASA NTRSDate not supplied

HLPW-5: Overview and Workshop Summary

The Fifth AIAA CFD High-Lift Prediction Workshop was held with the goal of assessing the numerical prediction capability of current-generation computational fluid dynamics (CFD) technology for swept, medium/high-aspect-ratio wings in high-lift configurations. A key aspect of this endeavor was the use of Technology Focus Groups (TFG), which included both mesh generation and flow-solver experts working together to accelerate advancements for their particular CFD methodology, by addressing key questions of importance prior to the workshop. The high-lift version of the NASA Common Research Model (CRM-HL) configuration was the focus of this workshop, and was used for three unique test cases. Wind-tunnel data were available for comparison for one of the test cases. Altogether, 365 datasets of CFD results were submitted by 47 teams, with 41 teams contributing to the multiple configurations of Case 1, 40 to Case 2, and 18 to Case 3. This paper provides a high-level summary of the results and conclusions from the workshop. As concluded from past workshops, fixed-grid Reynolds-averaged Navier-Stokes methods continued to be inaccurate and inconsistent for high-lift flows near maximum lift. However, application of mesh-adaptation technology helped to achieve improved consistency. Scale-resolving methods appeared most promising for predicting high-lift flow physics, particularly at maximum lift. Best practices for these methods were refined over the course of the workshop and new challenges were identified.

Adam M Clark↗
NASA NTRSDate not supplied

The Effect of Stochastic Acceleration on the Dynamics of Solar Energetic Particles in the Heliosphere

Solar Energetic Particles (SEPs) are crucial in space weather phenomena, impacting space-based technologies and human activities in space. Understanding the dynamics of SEPs, including their acceleration and transport mechanisms, is essential for improving predictive models and mitigating adverse effects. This study investigates the impact of stochastic acceleration (SA) on the dynamics of SEPs within the heliosphere. By integrating stochastic acceleration into a SEP transport model, we assess its impact relative to diffusive shock acceleration (DSA). The investigation utilizes the Adaptive Mesh Particle Simulator (AMPS) within the Space Weather Modeling Framework (SWMF) to simulate SEP dynamics, incorporating solar wind parameters and turbulence from the Alfven Wave Solar Model (AWSoM). To do this, we solve the focused transport equation to model the propagation of SEPs along magnetic field lines, extending these lines to 5 AU to incorporate the potential effects of pitch angle scattering beyond 1 AU. This aspect could notably impact the decay phase of a SEP event. Modeling the transport of SEPs is coupled with the simulation of solar wind dynamics, the interplanetary magnetic field, and the characteristics of Alfven wave turbulence. This presentation covers the modeling approach used in this research, the characterization of the impact of the stochastic acceleration, and the role of pitch angle scattering at various heliocentric distances on the dynamics of the SEP events' decay phase.

SEP↗
NASA NTRSDate not supplied

Sonic Boom Ground Noise Minimization Via the Adjoint Method

A multidisciplinary design optimization methodology to directly minimize the ground-level noise generated by the sonic boom of a high altitude supersonic body is presented. A Cartesian Euler flow solver is coupled with an atmospheric propagation tool to create a ground-noise analysis capability for supersonic bodies. Adjoint formulations for both the flow solver and propagation tool are also coupled to compute sensitivities of shape variations in the body in a highly efficient manner. A gradient-based optimizer is then introduced to forge a valuable design capability. Output-based mesh adaptation that is driven directly by ground-level noise is employed to increase accuracy and provide error estimation. The design method is demonstrated first on a simple axisymmetric body with few design variables to evaluate the efficacy of the optimization scheme. Guided by the results of this initial case, the problem is then repeated with somewhat different design variables to further demonstrate the capabilities of the design method. Finally, the method is applied to a real-world problem by optimizing control surface deflection settings of a low-boom aircraft to minimize ground noise while maintaining trimmed, level flight.

ARMD↗
DOE OSTI2013

Adaptive Implicit Non-Equilibrium Radiation Diffusion

We describe methods for accurate and efficient long term time integra- tion of non-equilibrium radiation diffusion systems: implicit time integration for effi- cient long term time integration of stiff multiphysics systems, local control theory based step size control to minimize the required global number of time steps while control- ling accuracy, dynamic 3D adaptive mesh refinement (AMR) to minimize memory and computational costs, Jacobian Free Newton-Krylov methods on AMR grids for efficient nonlinear solution, and optimal multilevel preconditioner components that provide level independent solver convergence.

Philip, Bobby↗
DOE OSTI2014

Parameter Analysis of the VPIN (Volume synchronized Probability of Informed Trading) Metric

VPIN (Volume synchronized Probability of Informed trading) is a leading indicator of liquidity-induced volatility. It is best known for having produced a signal more than hours before the Flash Crash of 2010. On that day, the market saw the biggest one-day point decline in the Dow Jones Industrial Average, which culminated to the market value of $1 trillion disappearing, but only to recover those losses twenty minutes later (Lauricella 2010). The computation of VPIN requires the user to set up a handful of free parameters. The values of these parameters significantly affect the effectiveness of VPIN as measured by the false positive rate (FPR). An earlier publication reported that a brute-force search of simple parameter combinations yielded a number of parameter combinations with FPR of 7%. This work is a systematic attempt to find an optimal parameter set using an optimization package, NOMAD (Nonlinear Optimization by Mesh Adaptive Direct Search) by Audet, le digabel, and tribes (2009) and le digabel (2011). We have implemented a number of techniques to reduce the computation time with NOMAD. Tests show that we can reduce the FPR to only 2%. To better understand the parameter choices, we have conducted a series of sensitivity analysis via uncertainty quantification on the parameter spaces using UQTK (Uncertainty Quantification Toolkit). Results have shown dominance of 2 parameters in the computation of FPR. Using the outputs from NOMAD optimization and sensitivity analysis, We recommend A range of values for each of the free parameters that perform well on a large set of futures trading records.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Akuna - Integrated Toolsets Supporting Advanced Subsurface Flow and Transport Simulations for Environmental Management

A next generation open source subsurface simulator and user environment for environmental management is being developed through a collaborative effort across Department of Energy National Laboratories. The flow and transport simulator, Amanzi, will be capable of modeling complex subsurface environments and processes using both unstructured and adaptive meshes at very fine spatial resolutions that require supercomputing-scale resources. The user environment, Akuna, provides users with a range of tools to manage environmental and simulator data sets, create models, manage and share simulation data, and visualize results. Underlying the user interface are core toolsets that provide algorithms for sensitivity analysis, parameter estimation, and uncertainty quantification. Akuna is open-source, cross platform software that is initially being demonstrated on the Hanford BC Cribs remediation site. In this paper, we describe the emerging capabilities of Akuna and illustrate how these are being applied to the BC Cribs site.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Progress in Fast, Accurate Multi-scale Climate Simulations

We present a survey of physical and computational techniques that have the potential to con- tribute to the next generation of high-fidelity, multi-scale climate simulations. Examples of the climate science problems that can be investigated with more depth include the capture of remote forcings of localized hydrological extreme events, an accurate representation of cloud features over a range of spatial and temporal scales, and parallel, large ensembles of simulations to more effectively explore model sensitivities and uncertainties. Numerical techniques, such as adaptive mesh refinement, implicit time integration, and separate treatment of fast physical time scales are enabling improved accuracy and fidelity in simulation of dynamics and allow more complete representations of climate features at the global scale. At the same time, part- nerships with computer science teams have focused on taking advantage of evolving computer architectures, such as many-core processors and GPUs, so that these approaches which were previously considered prohibitively costly have become both more efficient and scalable. In combination, progress in these three critical areas is poised to transform climate modeling in the coming decades.

climate modeling↗
DOE OSTI2013

miniGMG

miniGMG is a compact geometric multigrid (MG benchmark designed to proxy the performance characteristics of the solves found in adaptive mesh refinement multigrid (AMR MG) applications. It solves the equation a*alpha*u - b div beta grad u = f on a large cubical domain using cell-centered values and periodic boundaries. a and b are scalar constants, alpha and beta are space varying constants, f is the right hand side, and u is the solution. The cubical domain is divided into cubical subdomains which are distributed across the supercomputer. The righthand side f is generated from a continuous function for u. Thus, one can solve this linear system of equations and verify the error properties of the discrete solution u by comparing to a sampling of the continuous u. miniGMG implements both v- and f- cycles, as well as two constant coeffficient discretications of the equation above (one 27pt 2nd order, and one 13pt 4th order). In both cases, the (embedded) v-cycles are truncated when a subdomain reaces 4^3, 2^3, or 1^3 cells (a u-cycle). The multigrd solver then swtiches to one of seven coarse grid (bottom) solves. relaxation, BiCGStab, CABiCGStab, TelescopingCABiCGStab, CG, or CACG. The code exploits two forms of communication-avoiding. The first is exemplified by the referenced SC'12 paper in which DRAM data movment is avoided when smoothing. The second, described in the IPDPS'14 paper avoids collectives in a bottom solve (the CA bottom solves). There are a number of implementations of the operators for both CPUs and NVIDIA GPUs. The CPUs have been optimized to exploit a wavefront approach to communication-avoiding. The GPU versions include a number of optimized and tunable implementations. The code has been demonstrated to scale to 46K nodes on the BGQ Mira (750K cores), 9K nodes on Edison (111K cores), and 14K GPUs on Titan (37M cuda cores). It thus acts as a scalable testbed for a wide range of computer science and applied math research. Nominally, one invokes the benchmark as ./run b x y z alpha beta gamma Each subdomain is (2^b)^3 cells Each MPI process receives a x-by-y-by-z collection of subdomains. There are alpha*beta*gamma total MPI processes For correctness (a cubical domain), *alpha==y*beta==z*gamma.

Williams, Sameul↗
DOE OSTI2014

Snow Micro-Structure Model

PIKA is a MOOSE-based application for modeling micro-structure evolution of seasonal snow. The model will be useful for environmental, atmospheric, and climate scientists. Possible applications include application to energy balance models, ice sheet modeling, and avalanche forecasting. The model implements physics from published, peer-reviewed articles. The main purpose is to foster university and laboratory collaboration to build a larger multi-scale snow model using MOOSE. The main feature of the code is that it is implemented using the MOOSE framework, thus making features such as multiphysics coupling, adaptive mesh refinement, and parallel scalability native to the application. PIKA implements three equations: the phase-field equation for tracking the evolution of the ice-air interface within seasonal snow at the grain-scale; the heat equation for computing the temperature of both the ice and air within the snow; and the mass transport equation for monitoring the diffusion of water vapor in the pore space of the snow.

Micah Johnson, Andrew Slaughter↗
DOE OSTI2014

miniAMR

miniAMR is a "miniapp" which will be added to the Mantevo project. It provides a platform to experiment with different aspects of Adaptive Mesh Refinement (AMR). The calculations are performed on random data and the refinement is driven by simple objects moving through the mesh, so no real problem is solved. By having a small code that runs quickly, we can use to study different aspects of AMR and easily port the code to new testbed machines

Source record↗
DOE OSTI2013

Modeling Solar Wind Flow with the Multi-Scale Fluid-Kinetic Simulation Suite

Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS) is a package of numerical codes capable of performing adaptive mesh refinement simulations of complex plasma flows in the presence of discontinuities and charge exchange between ions and neutral atoms. The flow of the ionized component is described with the ideal MHD equations, while the transport of atoms is governed either by the Boltzmann equation or multiple Euler gas dynamics equations. We have enhanced the code with additional physical treatments for the transport of turbulence and acceleration of pickup ions in the interplanetary space and at the termination shock. In this article, we present the results of our numerical simulation of the solar wind (SW) interaction with the local interstellar medium (LISM) in different time-dependent and stationary formulations. Numerical results are compared with the Ulysses, Voyager, and OMNI observations. Finally, the SW boundary conditions are derived from in-situ spacecraft measurements and remote observations.

79 ASTRONOMY AND ASTROPHYSICS↗