Results for “Mesh Refinement” — Page 31

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 NTRSDate not supplied

Development of a Combined Cohesive and Virtual Crack-Closure Technique Approach to Represent R-Curves

Resistance curve (R-curve) effects due to fiber bridging, crack migration, and other blunting mechanisms are common in composite laminates. These mechanisms can dramatically increase the resistance to crack propagation but can be challenging to simulate. Delamination with R-curves can be analyzed using cohesive zone modeling (CZM) or the virtual crack closure technique (VCCT). Large fracture process zones can be simulated with CZM, but they require highly refined meshes. Coarser meshes can be used with VCCT, but this method is only applicable to small fracture process zones with R-curves defined as functions of position. Therefore, a technique with the computational efficiency of VCCT and the natural ability of cohesive elements to represent large fracture process zones is desirable. An approach is proposed that starts as a CZM in which the cohesive traction separation law (TSL) is separated into high-strength (HS) and low-strength (LS) components. The HS part, which is responsible for the mesh requirements of the CZM analysis, is replaced by VCCT. The combined method is evaluated by analyzing the response of double cantilevered beam specimens. Two specimen layup configurations, [0/90/90/0] 3s and [0/90/90/0] 9s , are evaluated. The R-curve response from each configuration is determined and applied to the constitutive properties of models built with CZM and with VCCT. The results and computational efficiency of the CZM, VCCT, and combined approaches are compared. The results are indicative that a combined cohesive/VCCT approach can enable progressive failure analyses to retain the computational efficiency of VCCT with the ability of the cohesive elements to capture R-curve effects.

Virtual Crack-Closure Technique↗

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

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

Investigation of Physical Mechanisms for Jet Noise Reduction by Plug Nozzle Porosity

This paper presents the results of a computational study on porous-plug nozzles for their eventual application in reducing takeoff noise levels for supersonic civil transport. The study explores a wide range of porous-plug geometries, encompassing variations in hole size, plug length and porosity. Computational fluid dynamics (CFD) analyses are employed in concert with an advanced automated mesh refinement (AMR) scheme, Sketch-to-Solution (S2S), to investigate porous plug nozzle design strategies that may diminish shock strength, broad-band shock noise (BBSN), flow separation, and other undesirable characteristics associated with high takeoff noise levels. A selected subset of the porous-plug nozzles that were analyzed computationally were manufactured at subscale using stereolithography 3D printing. These nozzles were tested experimentally, obtaining far-field noise spectral data and schlieren flow visualization pictures. Detailed comparisons between CFD predictions and experimental results were conducted. These comparisons focus on gaining a better understanding of the underlying physical mechanisms responsible for noise reduction with the porous plugs relative to corresponding solid plugs.

computational fluid dynamics↗

NASA NTRS2025

Physics-Based Modeling and Simulation of Self-Reacting Friction Stir Welding Using Computational Fluid Dynamics

A physics-based model was developed to simulate the behavior of material in a self-reacting friction stir welding (SR-FSW) process for the joining of metals. This steady-state model builds upon fundamental computational fluid dynamic (CFD) principles within Ansys Fluent to solve the discretized equations. The effective viscosity is calculated using a viscoplastic model using a Sheppard-Wright formulation of flow stress. Numerous advancements have been made in the incorporated physics including (1) temperature-dependent material properties; (2) locally adaptable flow and thermal boundary conditions; and (3) adapting material properties in nugget in response to microstructural changes. Simulation strategies to accelerate computation and improve numerical stability include adapting the mesh refinement and solver relaxation factors during simulation. The result is a highly robust and computationally efficient model capable of providing the material flow and temperature history across the domain. As material history determines the local microstructure and ultimately weld strength, an accurate and detailed physics-based model has the potential to accelerate SR-FSW process development. The model is highly adaptable to changes in process parameters, tool design, or alloy.

Process Modeling↗

NASA NTRSDate not supplied

Investigation of Physical Mechanisms for Jet Noise Reduction by Plug Nozzle Porosity

This paper presents the results of a computational study on porous-plug nozzles for their eventual application in reducing takeoff noise levels for supersonic civil transport. The study explores a wide range of porous-plug geometries, encompassing variations in hole size, plug length and porosity. Computational fluid dynamics (CFD) analyses are employed in concert with an advanced automated mesh refinement (AMR) scheme, Sketch-to-Solution (S2S), to investigate porous plug nozzle design strategies that may diminish shock strength, broad-band shock noise (BBSN), flow separation, and other undesirable characteristics associated with high takeoff noise levels. A selected subset of the porous-plug nozzles that were analyzed computationally were manufactured at subscale using stereolithography 3D printing. These nozzles were tested experimentally, obtaining far-field noise spectral data and schlieren flow visualization pictures. Detailed comparisons between CFD predictions and experimental results were conducted. These comparisons focus on gaining a better understanding of the underlying physical mechanisms responsible for noise reduction with the porous plugs relative to corresponding solid plugs.

Nozzles↗

DOE OSTI2020

Progressive collapse analysis of stiffened panel structure with consideration of actual welding distortion and residual stresses

To date, inclusion of welding distortion and residual stresses in determining the strength of marine structures is based on empirical formulations due to welding complexity. In this study, a more feasible numerical model was developed from the perspective of inherent strain theory. The accuracy of this simplified numerical model was confirmed by transient thermo-elastic–plastic finite element analysis using a novel dynamic mesh refining and iterative substructure method which is 50 times faster than conventional finite element method. Welding distortion by the two numerical approaches was validated with experimental measurement. The collapse analyses indicated that the inherent strain model predicted a compressive strength and load-displacement curve with high accuracy, as validated by the transient thermo-mechanical model. The global stiffness and the ultimate strength of the welded structure would be overestimated by about 22% if no welding effects were considered. The proposed computation procedure for determining the ultimate strength of welded structures is accurate, simple, and practical and enables improvement of ultimate strength through optimization of welding joint design to reduce the welding residual stress and distortion.

Huang, Hui↗

DOE OSTI2022

MFIX-Exa: Performance prediction of multiphase energy conversion devices

MFIX-Exa targets the high-fidelity CFD-DEM model in which particles are unresolved by the fluid grid (typically using a fluid mesh approximately twice the particle diameter) but discrete particle dynamics including collisions are resolved with a simple linear-spring dashpot model. The fluid field is solved with a modern low Mach number formation of a cell-centered, nodal-pressure approximate projection method using a Godunov scheme . Physics capabilities include: complex geometries via embedded boundaries (EBs), open and closed system ideal gas equations of state, species transport and heterogeneous chemi- cal reactions. MFIX-Exa is built on the AMReX software framework, the ECP Block-Structured Adaptive Mesh Refinement (AMR) Co-Design Center, which allows the code to be portable and performant.Further integration into the ECP ecosystem includes HYPRE linear solvers and Ascent for in situ visualization. Current project efforts focus on scaling up realistic simulations to the KPP of the challenge problem: NETs 50kW chemical looping reactor, shakedown and scaling on JSLE and OLCF TDS systems and reaching for any remaining performance improvements.

Musser, Jordan↗

DOE OSTI2022

Novel Solver Algorithms for Nearly Singular Linear Systems Arising in Combustion Modelling

Direct Numerical Simulations of realistic combustion devices are extremely challenging due to the wide separation of scales in the simulation, for example an internal combustion (IC) engine chamber, and the flame thickness of a high-pressure flame. The PeleLMeX solver uses adaptive mesh refinement (AMR) to evolve multi-species reacting flows in the low Mach number limit at the Exascale and relies on an embedded boundary (EB) approach to represent complex geometries. In that framework, the EB geometries often give rise to very small cut-cells along the boundary, which translate into extreme ill-conditioning of the pressure-projection, with eigenvalues that span 15-16 orders of magnitude. In this talk, we focus on the case of a typical IC piston bowl geometry for which we present on a novel approach towards solving these nearly singular linear systems with ILU-based, C-AMG smoothers on massively parallel architectures. In particular, we use scaling and equilibration algorithms to handle the non-normality of the upper triangular factors. This enables us to approximate the highly sequential triangular solve algorithm, embedded in the AMG smoothing-solve phase, with Jacobi iterations. This approximation can be written as a convergent Neumann series whose terms are composed of highly parallel sparse matrix vector multiplications. The result is an algorithm that substantially decreases setup and solve time, compared to state-of-the-art, for these challenging linear systems.

combustion modelling↗

DOE OSTI2022

Simulations of Low Mach Number Reactive Flows Coupled with Electric Fields

Applying electric fields to flames can provide an effective, non-invasive control of the combustion process. In this work, an algorithm recently introduced to handle the interactions between reactive flows and an electric field is extended to handle multiple dimensions and adaptive mesh refinement (AMR), and applied to a series of laminar flames. Experimental data allows to validate the ability of the algorithm to capture the mechanisms driving the flame response to electric fields.

electric fields↗

DOE OSTI2022

Simulations of Low Mach Number Reactive Flows Coupled with Electric Fields: Preprint

Applying electric fields to flames can provide an effective, non-invasive control of the combustion process. In this work, an algorithm recently introduced to handle the interactions between reactive flows and an electric field is extended to handle multiple dimensions and adaptive mesh refinement (AMR), and applied to a series of laminar flames. Experimental data allows to validate the ability of the algorithm to capture the mechanisms driving the flame response to electric fields.

electric field↗

DOE OSTI2022

DNS of Hydrodynamic Instabilities of Laminar H2/O2/N2 Flames at Elevated Pressure

The atypical combustion properties of hydrogen have long provided an interest for hydrogen-enriched combustion, emphasized in recent years by increasingly stringent emission regulations. Significant work has been focused on the small hydrocarbon(C1-C3)/hydrogen blends, but much less is known about the effect of hydrogen on heavier hydrocarbon combustion, of interest for the aviation and automotive industry. In this work, we perform direct numerical simulations (DNS) of spherically expanding laminar n-dodecane/H2 flames in a constant-volume vessel at elevated pressure and a range of hydrogen seeding levels. A low Mach number model is used to represent the gas and flame dynamics allowing for a temporally varying, spatially homogeneous pressure field, while also enabling numerical evolution of the system numerically at time step sizes governed by the advective CFL limit rather than acoustic processes. We use an adaptive mesh refinement (AMR) approach to tackle the large separation of scales associated with high pressure premixed flame surfaces propagating in a laboratory-scale closed domain. We incorporate detailed chemistry and transport models for the reacting flow and focus the analysis on the effects of the pressure rise on the flame propagation characteristics and morphology as the hydrogen content is increased. The presence of fast diffusive H2 and the pressure rise results in dramatic thinning of the flame and the potential onset of hydrodyanmic and theromdiffusive instabilities. We assess the potential role of these instabilities on the ability to accurately measure laminar burning speed experimentally based on the mean flame surface propagation speed.

high pressure↗

DOE OSTI2022

(Towards) DNS of a Laboratory Lean CH4/H2 Low-Swirl Flame Impinging on an Inclined Wall

Due to downsizing trends, flame-wall interaction (FWI) is increasingly prominent in gas turbines (GTs). FWI has direct consequences on flame stabilization and pollutant emissions, but it is not well understood in turbulent flows representative of GTs. We present results from a direct numerical simulation (DNS) of a turbulent CH4/H2 model GT low-swirl laboratory-scale flame interacting with an inclined wall. The results from the laboratory flame include simultaneous measurements of velocity using stereo particle imaging velocimetry and OHxCH2O planar laser induced fluorescence. The adaptive-mesh refinement solver PeleLMeX is used, with 24-species, 105-reaction reduced Aramco chemical kinetics mechanism. The premixed fuel-air mixture consists of hydrogen-enriched methane with 70% hydrogen volume fraction and 0.4 equivalence ratio. The inflow is prescribed to match experimental measurements at the burner exit. Karlovitz and turbulent Reynolds numbers are 300 and 400, respectively. The simulation and experimental results show excellent agreement. The flame features a bowl-shape stabilization, with a corrugated, continuous flame front at the leading edge, followed by fragmented reaction zones downstream. A large diffuse cloud of CH2O is formed downstream of the quenching point. The simulation results indicate that the cloud of CH2O is the result of incomplete methane combustion, with CH2O "leaking" from the locally quenched reaction zones.The DNS provides fine-grain resolution of turbulence-flame-wall interaction that cannot be captured with experimental measurements. With access to the entire solution vector at each cell of the computational domain, the local quenching.

flame-wall interactions↗

DOE OSTI2023

Examination of a Methane/Diesel RCCI Engine Using Pele: Preprint

Multi-fuel, advanced injection strategies have become increasingly promising as a strategy to mitigate the emissions generated from internal combustion engines. By carefully controlling the combustion phasing in-cylinder, these new multi-pulse, multi-fuel injection strategies are able to burn in the low-temperature combustion regime where both NOx and soot are not readily produced, reducing the need for extensive exhaust gas recirculation systems. In this study, we examine a reactivity-controlled compression ignition (RCCI) strategy that uses an early pre-filled methane-air mixture with low turbulence background as the low-reactivity fuel and a direct injection of four discrete dodecane jets as a surrogate for the high-reactivity diesel fuel. We use the Pele software suite, a highly optimized, exascale-ready, adaptive mesh refinement codebase to perform high-resolution numerical simulations of a scaled down, single cylinder from the RCCI engine. Here, we resolve the ignition kernels down to micrometer scales and present several statistical quantities evaluating the development of the flow and detailing the onset of ignition and subsequent flame development. Particular attention is paid to the conditions surrounding the onset of the first ignition kernels and discussing what led to the development of those conditions.

CFD↗

DOE OSTI2023

Simulating Catalysis with Realistic Pellet Geometries Using Mesoflow: A Case Study of Catalytic Propane Dehydrogenation

We present a case study of catalytic propane dehydrogenation with our open-source multiphysics solver, Mesoflow. The solver was developed to simulate reactive flow coupled to heterogeneous catalytic reactions and deactivation in the context of complex, mesoscale geometry. The method leverages cartesian block-structured adaptive mesh refinement to capture realistic catalyst microstructural features acquired directly from X-ray computed tomography data. A kinetic model for propane dehydrogenation and catalyst deactivation was developed based on temporal analysis of products (TAP) reactor experiments. The TAP reactor experiments allow for precise characterization of intrinsic kinetic reaction steps which are implemented into Mesoflow simulations to model the spatial and temporal evolution of reactants, products, and catalyst active sites. The short-term and long-term deactivation behavior is studied by using XCT data collected from fresh and aged catalyst pellets, which exhibit different microstructural features. This study employs time-splitting algorithms to connect disparate reaction and flow timescales, enabling the simulations to achieve realistic deactivation timescales on the order of minutes while the flow time-scales for small particles (100 microns) are several milliseconds. We also introduce a flexible automated python script that writes the necessary files to construct a Mesoflow simulation from user-created chemical mechanisms. We will also introduce a few new features that are added to Mesoflow such as higher order schemes, implicit chemistry integrators and the ability to run on AMD and NVIDIA graphics-processing-units.

AMReX↗

DOE OSTI2024

Influence of Real Gas Effects on Chemical Kinetics in Oxycombustion in Supercritical Carbon Dioxide

Oxycombustion in supercritical carbon dioxide is an integral part of the Allam Cycle, a technology that enables carbon-neutral use of fossil-fuels and carbon-negative use of biofuels. We simulate oxycombustion in a realistic combustor geometry for two sets of fuel conditions: pure methane, and a 40 percent methane/60 percent carbon dioxide blend, both at 343.15 K. The fuel jet mixes with a preheated swirler of 20 percent oxygen and 80 percent carbon dioxide at 1005.35 K, with a 100 percent carbon dioxide coflow at 783.15 K. The entire system operates at a pressure of 300 bar, putting the entire system above the critical temperature and pressure of carbon dioxide. Simulations are performed using PeleC, a compressible block-structured adaptive-mesh refinement (AMR) reacting flow code. Direct comparison of results using thermodynamically self-consistent implementations of the ideal gas equation of state and the Soave-Redlich-Kwong (SRK) equation of state allow quantification of the real gas impacts on combustion kinetics, flow temperatures, and pollutant formation.

carbon dioxide↗