Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesh compression”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗
NASA NTRS1995

Methodology for simulation of unsteady flows to determine the time-dependent interference between stationary and moving boundaries

A new methodology is developed to simulate unsteady flows about prescribed and aerodynamically determined moving boundary problems. The method couples the fluid dynamics and rigid-body dynamics equations to capture the time-dependent interference between stationary and moving boundaries. The unsteady, compressible, inviscid (Euler) equations are solved on dynamic, unstructured grids by an explicit, finite-volume, upwind method. For efficiency, the grid adaptation is performed within a window around the moving object. The Eulerian equations of the rigid-body dynamics are solved by a Runge-Kutta method in a non-inertial frame of reference. The two-dimensional flow solver is validated by computing the flow past a sinusoidally-pitching airfoil and comparing these results with the experimental data. The overall methodology is used for two two-dimensional examples: the flow past an airfoil which is performing a three-degrees-of-freedom motion in a transonic freestream, and the free-fall of a store after separation from a wing-section. Then the unstructured mesh methodology is extended to three-dimensions to simulate unsteady flow past bodies in relative motion, where the trajectory is determined from the instantaneous aerodynamics. The flow solver and the adaptation scheme in three dimensions are validated by simulating the transonic, unsteady flow around a wing undergoing a forced, periodic, pitching motion, and comparing the results with the experimental data. To validate the trajectory code, the six-degrees-of-freedom motion of a store separating from a wing was computed using the experimentally determined force and moment fields, then comparing with an independently generated trajectory. Finally, the overall methodology was demonstrated by simulating the unsteady flowfield and the trajectory of a store dropped from a wing. The methodology, its computational cost notwithstanding, has proven to be accurate, automated, easy for dynamic gridding, and relatively efficient for the required man-hours.

Singh, Kamakhya Prasad↗
DOE OSTI2012

Collaborative Research: Continuous Dynamic Grid Adaptation in a Global Atmospheric Model - Application and Refinement

This project had goals of advancing the performance capabilities of the numerical general circulation model EULAG and using it to produce a fully operational atmospheric global climate model (AGCM) that can employ either static or dynamic grid stretching for targeted phenomena. The resulting AGCM combined EULAG's advanced dynamics core with the "physics" of the NCAR Community Atmospheric Model (CAM). Effort discussed below shows how we improved model performance and tested both EULAG and the coupled CAM-EULAG in several ways to demonstrate the grid stretching and ability to simulate very well a wide range of scales, that is, multi-scale capability. We leveraged our effort through interaction with an international EULAG community that has collectively developed new features and applications of EULAG, which we exploited for our own work summarized here. Overall, the work contributed to over 40 peer-reviewed publications and over 70 conference/workshop/seminar presentations, many of them invited. 3a. EULAG Advances EULAG is a non-hydrostatic, parallel computational model for all-scale geophysical flows. EULAG's name derives from its two computational options: EULerian (flux form) or semi-LAGrangian (advective form). The model combines nonoscillatory forward-in-time (NFT) numerical algorithms with a robust elliptic Krylov solver. A signature feature of EULAG is that it is formulated in generalized time-dependent curvilinear coordinates. In particular, this enables grid adaptivity. In total, these features give EULAG novel advantages over many existing dynamical cores. For EULAG itself, numerical advances included refining boundary conditions and filters for optimizing model performance in polar regions. We also added flexibility to the model's underlying formulation, allowing it to work with the pseudo-compressible equation set of Durran in addition to EULAG's standard anelastic formulation. Work in collaboration with others also extended the demonstrated range of validity of soundproof models, showing that they are more broadly applicable than some had previously thought. Substantial testing of EULAG included application and extension of the Jablonowski-Williamson baroclinic wave test - an archetype of planetary weather - and further analysis of multi-scale interactions arising from collapse of temperature fronts in both the baroclinic wave test and simulations of the Held-Suarez idealized climate. These analyses revealed properties of atmospheric gravity waves not seen in previous work and further demonstrated the ability of EULAG to simulate realistic behavior over several orders of magnitude of length scales. Additional collaborative work enhanced capability for modeling atmospheric flows with adaptive moving meshes and demonstrated the ability of EULAG to move into petascale computing. 3b. CAM-EULAG Advances We have developed CAM-EULAG in collaboration with former project postdoc, now University of Cape Town Assistant Professor, Babatunde Abiodun. Initial study documented good model performance in aqua-planet simulations. In particular, we showed that the grid adaptivity (stretching) implemented in CAM-EULAG allows higher resolution in selected regions without causing anomalous behavior such as spurious wave reflection. We then used the stretched-grid version to analyze simulated extreme precipitation events in West Africa, comparing the precipitation and event environment with observed behavior. The model simulates fairly well the spatial scale and the interannual and intraseasonal variability of the extreme events, although its extreme precipitation intensity is weaker than observed. In addition, both observations and the simulations show possible forcing of extreme events by African easterly waves. 3c. Other Contributions Through our collaborations, we have made contributions to a wide range of outcomes. For research focused on terrestrial behavior, these have included (1) upwind schemes for gas dynamics, (2) a nonlinear perspective on the dynamics of the Madden-Julian Oscillation, (3) numerical realism of thermal convection, (4) extreme precipitation processes, (5) unstructured meshes, (6) porous-media analogies for simulating flow in urban environments, (7) Monge-Ampere enhancements of semi-Lagrangian methods and (8) further extensions of the generality of EULAG applications to geophysical problems. The generality of EULAG for multi-scale flow was also demonstrated through solar applications. Simulations of solar magnetohydrodynamics showed that EULAG produces the solar cycle fairly well. Although the model's solar cycle is significantly longer than observed, it shows remarkable agreement with observations in features of spatial evolution. These efforts in environments that are extreme from a geophysical perspective have helped to establish the robustness of EULAG as an all-scale code for geophysical applications.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2005

Navier-Stokes Computations With One-Equation Turbulence Model for Flows Along Concave Wall Surfaces

This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation turbulence model is an effective approach for turbulence modeling in the near solid wall surface region of flow over a concave wall.

Wang, Chi R.↗
NASA NTRS1985

Numerical studies of 2-dimensional flows

A formulation of the lambda scheme for the analysis of two dimensional inviscid, compressible, unsteady transonic flows is presented. The scheme uses generalized Riemann variables to determine the appropriate two point, one sided finite difference approximation for each derivative in the unsteady Euler equations. These finite differences are applied at the predictor and corrector levels with shock updating at each level. The weaker oblique shocks are captured, but strong near normal shocks are fitted into the flow using the Rankine-Hugoniot relations. This code is demonstrated with a numerical example of a duct flow problem with developing normal and oblique shock waves. The technique is implemented in a code which has been made efficient by streamlining to a minimal number of operations and by eliminating branch statements. The scheme is shown to provide an accurate analysis of the flow, including formation, motions, and interactions of shocks; the results obtained on a relatively coarse mesh are comparable to those obtained by other methods on much finer meshes.

Moretti, G.↗
NASA NTRS2001

Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computations

The recently developed essentially fourth-order or higher low dissipative shock-capturing scheme of Yee, Sandham and Djomehri (1999) aimed at minimizing nu- merical dissipations for high speed compressible viscous flows containing shocks, shears and turbulence. To detect non smooth behavior and control the amount of numerical dissipation to be added, Yee et al. employed an artificial compression method (ACM) of Harten (1978) but utilize it in an entirely different context than Harten originally intended. The ACM sensor consists of two tuning parameters and is highly physical problem dependent. To minimize the tuning of parameters and physical problem dependence, new sensors with improved detection properties are proposed. The new sensors are derived from utilizing appropriate non-orthogonal wavelet basis functions and they can be used to completely switch to the extra numerical dissipation outside shock layers. The non-dissipative spatial base scheme of arbitrarily high order of accuracy can be maintained without compromising its stability at all parts of the domain where the solution is smooth. Two types of redundant non-orthogonal wavelet basis functions are considered. One is the B-spline wavelet (Mallat & Zhong 1992) used by Gerritsen and Olsson (1996) in an adaptive mesh refinement method, to determine regions where re nement should be done. The other is the modification of the multiresolution method of Harten (1995) by converting it to a new, redundant, non-orthogonal wavelet. The wavelet sensor is then obtained by computing the estimated Lipschitz exponent of a chosen physical quantity (or vector) to be sensed on a chosen wavelet basis function. Both wavelet sensors can be viewed as dual purpose adaptive methods leading to dynamic numerical dissipation control and improved grid adaptation indicators. Consequently, they are useful not only for shock-turbulence computations but also for computational aeroacoustics and numerical combustion. In addition, these sensors are scheme independent and can be stand alone options for numerical algorithm other than the Yee et al. scheme.

Sjoegreen, B.↗
DOE OSTI2024

Beryllium–tungsten graded density inner shells in double shell capsules for improved hydrodynamic stability

The outer surface of the high-Z inner shell in the double shell configuration of inertial confinement fusion experiments experiences Rayleigh–Taylor instability growth during the implosion process due to inverted density and pressure gradients between a highly compressed foam interstitial layer and the accelerating dense inner shell. Graded density layers have long been known to reduce instability growth rates. In this study, we employ high-fidelity radiation hydrodynamic simulations to demonstrate this improved stability when grading beryllium into tungsten. We first characterize the response to L-band preheat of these layers using a newly calibrated radiation drive. While graded layer capsules suffer reduced performance (here, measured as DD neutron yield from a CD foam fuel) in 1D simulations due to reduced kinetic energy coupling and reduced fuel compression, they suffer less of a performance drop when 2D instabilities are accounted for. With the improved stability of graded layers, we explore the performance of capsules with larger fuel radii and thinner shells as a preliminary study to find new designs in which graded layers produce the highest yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2018

Implementation of a Matrix Crack Spacing Parameter in a Continuum Damage Mechanics Finite Element Model

Continuum Damage Mechanics (CDM) based progressive damage and failure analysis (PDFA) methods have demonstrated success in a variety of finite element analysis (FEA) implementations. However, the technical maturity of CDM codes has not yet been proven for the full design space of composite materials in aerospace applications. CDM-based approaches represent the presence of damage by changing the local material stiffness definitions and without updating the original mesh or element integration schemes. Without discretely representing cracks and their paths through the mesh, damage in models with CDM-based materials is often distributed in a region of partially damaged elements ahead of stress concentrations. Having a series of discrete matrix cracks represented by a softened region may affect predictions of damage propagation and, thus, structural failure. This issue can be mitigated by restricting matrix damage development to discrete, fiber-aligned rows of elements; hence CDM-based matrix cracks can be implemented to be more representative of discrete matrix cracks. This paper evaluates the effect of restricting CDM matrix crack development to discrete, fiber-aligned rows where the spacing of these rows is controlled by a user-defined crack spacing parameter. Initially, the effect of incrementally increasing matrix crack spacing in a unidirectional center notch coupon is evaluated. Then, the lessons learned from the center notch specimen are applied to open-hole compression finite element models. Results are compared to test data, and the limitations, successes, and potential of the matrix crack spacing approach are discussed.

Hyder, Imran↗
DOE OSTI2023

The MOOSE fluid properties module

The fluid properties module in MOOSE~\cite{lindsay2022moose} serves a variety of fluid simulation applications based on MOOSE, including the MOOSE Navier Stokes module~\cite{moose_ns}, Pronghorn~\cite{pgh}, the MOOSE thermal hydraulics module, SAM~\cite{sam}, RELAP-7~\cite{relap7}, Sockeye~\cite{sockeye}, Pronghorn-subchannel~\cite{subchannel} and the MOOSE porous flow module~\cite{porous}. These applications are used to solve coarse mesh multi-dimensional thermal-hydraulics~\cite{pgh}, 1D systems analysis~\cite{sam,relap7} in nuclear reactor analysis, heat pipe modeling~\cite{sockeye} and porous flow simulations~\cite{porous} for underground gas storage and water seepage. The use of consistent fluid properties across fluid flow applications facilitates coupled flow simulations~\cite{anl_sam_pgh,osti_1889653}. Each application has historically driven the implementation of several fluid properties, which were later extended to be compatible with other applications. The unique diversity of applications of the module, due to its presence in MOOSE, has driven its expansion to new fluids, such as advanced nuclear reactor coolants and, more recently, arbitrary functions or tables-based property definitions, as detailed in section~\ref{content}, as well as numerous thermophysical properties and variable sets, as detailed in subsection~\ref{sec:prop}. The need for different discretizations of flow equations based on the compressibility of the fluid has motivated support for both primitive (pressure- and temperature-based) and conservative (internal energy- and specific volume-based) flow variables in the module. Thermodynamic relations are used to convert between these two formulations, as needed. The module serves a dual purpose of providing fluid properties for direct use in flow simulations and facilitating the implementation of user-specific fluid properties. Contributions of new properties for existing fluids or new fluids are strongly encouraged.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

CFD and X-ray Analysis of Gaseous Direct Injection from an Outward Opening Injector

Using natural gas in an internal combustion engine (ICE) is emerging as a promising way to improve thermal efficiency and reduce exhaust emissions. In the development of such engine platforms, computational fluid dynamics (CFD) plays a fundamental role in the optimization of geometries and operating parameters. One of the most relevant issues in the simulation of direct injection (DI) gaseous processes is the accurate prediction of the gas jet evolution. The simulation of the injection process for a gaseous fuel does not require complex modeling, nevertheless properly describing high-pressure gas jets remains a challenging task. At the exit of the nozzle, the injected gas is under-expanded, the flow becomes supersonic and shocks occur due to compressibility effects. These phenomena lead to challenging computational requirements resulting from high grid resolution and low computational time-steps. In this paper, the simulation of gaseous injection from an outward opening injector was carried out with the commercial CFD software CONVERGE. Argon was injected into a quiescent chamber filled with Nitrogen. Two free-jet cases with injection pressure of 15 and 10 bar, and ambient pressure of 1 and 2 bar respectively, were tested. For the 15 bar injection pressure and 1 bar ambient pressure, a jet impinging on a barrier case was performed as well. The numerical results were validated against experimental mass measurements using X-ray radiography. The results, presented in terms of 2D maps and transverse sections of equivalent path length of argon, showed good agreement between modeling and experiments in terms of jet penetration, shape and gas mass distribution. These emphasized the significant impact of mesh size on the results for the simulation of this type of injector. Due to the small geometrical dimensions a fine mesh is required in order to deliver an accurate description of the mixture formation. Numerical results were discussed with a particular emphasis on the application to engine simulations, where a trade-off between accuracy and computational time is a necessity.

Bartolucci, Lorenzo↗
DOE OSTI2019

A new micromechanics based full field numerical framework to simulate the effects of dynamic recrystallization on the formability of HCP metals

This work presents a new full-field, efficient and mesh-free numerical framework to model microstructure evolution, dynamic recrystallization (DRX) and formability in hexagonal closed-packed (HCP) metals such as magnesium alloys at warm temperatures. A rate tangent-fast Fourier transform-based elasto-viscoplastic crystal plasticity constitutive model for HCP metals (RTCP-FFT-HCP) is coupled with a probabilistic cellular automata (CA) approach to model DRX (CA-DRX). Furthermore, this new model is coupled with the Marciniak-Kuczynski (M-K) approach to model formability of magnesium alloys at elevated temperatures. The RTCP-FFT-HCP model computes macro stress-strain, twinning volume fraction, micromechanical fields, texture evolution and local dislocation density. Nucleation of new grains and their subsequent growth is modeled using the cellular automata approach with probabilistic state switching rule. First, the proposed RTCP-FFT-HCP model is validated by comparing the predicted stress-strain responses and texture evolution under uniaxial tension and compression with experimental measurements for AZ31 sheet alloy at room temperature. The coupled CA-RTCPFFT-HCP model is further validated by comparing the predicted stress-strain responses and texture evolution in uniaxial compression with experimental measurements at 100 °C, 200 °C and 300 °C for the AZ31 sheet alloy. Next, the forming limit diagrams (FLDs) with and without including the effects of DRX are simulated at 100 °C, 200 °C and 300 °C for AZ31 sheet alloy. The predicted FLDs with DRX show very good agreement with the experimental measurements and clearly demonstrate the need for an accurate DRX model. The study reveals that the DRX strongly affects the deformed grain structure, grain size and texture evolution and also highlights the importance accounting for DRX during FLD simulations at high temperatures.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Simulation of Complex Cracking in Plain Weave C/SiC Composite under Biaxial Loading

Finite element analysis is performed on a mesh, based on computed geometry of a plain weave C/SiC composite with assumed internal stacking, to reveal the pattern of internal damage due to biaxial normal cyclic loading. The simulation encompasses intertow matrix cracking, matrix cracking inside the tows, and separation at the tow-intertow matrix and tow-tow interfaces. All these dissipative behaviors are represented by traction-separation cohesive laws. Not aimed at quantitatively predicting the overall stress-strain relation, the simulation, however, does not take the actual process of fiber debonding into account. The fiber tows are represented by a simple rule-of-mixture model where the reinforcing phase is a hypothetical one-dimensional material. Numerical results indicate that for the plain weave C/SiC composite, 1) matrix-crack initiation sites are primarily determined by large intertow matrix voids and interlayer tow-tow contacts, 2) the pattern of internal damage strongly depends on the loading path and initial stress, 3) compressive loading inflicts virtually no damage evolution. KEY WORDS: ceramic matrix composite, plain weave, cohesive model, brittle failure, smeared crack model, progressive damage, meso-mechanical analysis, finite element.

Cheng, Ron-Bin↗
DOE OSTI2019

Pronghorn

Pronghorn is designed as a multi-dimensional coarse-mesh conjugate heat transfer simulator (tightly coupled fluid and solid phase heat transfer) for nuclear reactor core applications. Pronghorn’s design philosophy of is that: 1. The tool will be agnostic of reactor concept, including VHTRs, SFRs, FHRs, and LWRs. 2. Must be higher fidelity than traditional nuclear systems/safety analysis (0D-1D). 3. Avoid lower-length scale phenomena, such as boundary layer theory, i.e., rely on closure relations. 4. Apply mixture theory (homogenization) wherever possible. 5. Conserve mass and energy (solve all physics on single mesh). 6. Must be easily coupled to both MOOSE-based and external physics applications. 7. Does not require a super computer. Thus, Pronghorn’s algorithmic philosophy is to resolve large scale flow features and model lower-length scale phenomena, such as boundary layers. The coarse mesh approach for lower-length scale modeling is to develop a multi-component multi-phase homogenization of solid and liquid phases into a mixture formulation inside a finite element. For example, a pebble bed core is represented as a two-phase mixture of coolant and fuel pebbles, where the fuel phase is stationary. Thus, the solid phase exerts a drag force on the fluid phase, composed of form and viscous drag, and both phases exchange internal energy through conduction and convective heat transfer. For Pronghorn, the volume fraction of fluid can vary in space, but not in time. Pronghorn has been validated against the VHTR PBMR400 benchmark problems, used for rod ejection simulations have been conducted for thermal-fluids/neutronics for both pebble-bed and prismatic gas-cooled reactors and simple thermal-neutronics coupled LWR benchmark problems. Pronghorn flow models cover a wide range of transient and steady flow conditions, including, porous fluid flow, compressible flow, and nearly incompressible, thermally expansive flow. The physics can be solved in three-dimensional Cartesian coordinates and two-dimensional cylindrical (r-z) coordinates.

Martineau, RichardC.↗
DOE OSTI2012

Progress towards a high-gain and robust target design for heavy ion fusion

Recently [E. Henestroza et al., Phys. Plasmas 18, 032702 (2011)], a new inertial-fusion target configuration, the X-target, using one-sided axial illumination has been explored. This class of target uses annular and solid-profile heavy ion beams to compress and ignite deuterium-tritium (DT) fuel that fills the interior of metal cases that have side-view cross sections in the shape of an 'X.' X-targets using all-DT-filled metal cases imploded by three annular ion beams resulted in fuel densities of {approx}50 g/cm{sup 3} at peak compression, and fusion gains of {approx}50, comparable to heavy ion driven hohlraum targets [D. A. Callahan-Miller and M. Tabak, Phys. Plasmas 7, 2083 (2000)]. This paper discusses updated X-target configurations that incorporate inside the case a propellant (plastic) and a pusher (aluminum) surrounding the DT fuel. The updated configurations are capable of assembling higher fuel areal densities {approx}2 g/cm{sup 2} using two annular beams to implode the target to peak DT densities {approx}100 g/cm{sup 3}, followed by a fast-ignition solid ion beam which heats the high-density fuel to thermonuclear temperatures in {approx}200 ps to start the burn propagation, obtaining gains of {approx}300. These targets have been modeled using the radiation-hydrodynamics code HYDRA [M. M. Marinak et al., Phys. Plasmas 8, 2275 (2001)] in two- and three- dimensions to study the properties of the implosion as well as the ignition and burn propagation phases. At typical Eulerian mesh resolutions of a few microns, the aluminum-DT interface shows negligible Rayleigh-Taylor (RT) and Richtmyer-Meshkov instability growth; also, the shear flow of the DT fuel as it slides along the metal X-target walls, which drives the RT and Kelvin Helmholtz instabilities, does not have a major effect on the burning rate. An analytic estimate of the RT instability process at the Al-DT interface shows that the aluminum spikes generated during the pusher deceleration phase would not reach the ignition zone in time to affect the burning process. Also, preliminary HYDRA calculations, using a higher resolution mesh to study the shear flow of the DT fuel along the X-target walls, indicate that metal-mixed fuel produced near the walls would not be transferred to the DT ignition zone (at maximum {rho}R) located at the vertex of the X-target.

36 MATERIALS SCIENCE↗
DOE OSTI2014

The matryoshka run: A Eulerian refinement strategy to study the statistics of turbulence in virialized cosmic structures

We study the statistical properties of turbulence driven by structure formation in a massive merging galaxy cluster at redshift z = 0. The development of turbulence is ensured as the largest eddy turnover time is much shorter than the Hubble time independent of mass and redshift. We achieve a large dynamic range of spatial scales through a novel Eulerian refinement strategy where the cluster volume is refined with progressively finer uniform nested grids during gravitational collapse. This provides an unprecedented resolution of 7.3 h {sup –1} kpc across the virial volume. The probability density functions of various velocity-derived quantities exhibit the features characteristic of fully developed compressible turbulence observed in dedicated periodic-box simulations. Shocks generate only 60% of the total vorticity within the R {sub vir}/3 region and 40% beyond that. We compute second- and third-order longitudinal and transverse structure functions for both solenoidal and compressional components in the cluster core, virial region, and beyond. The structure functions exhibit a well-defined inertial range of turbulent cascade. The injection scale is comparable to the virial radius but increases toward the outskirts. Within R {sub vir}/3, the spectral slope of the solenoidal component is close to Kolmogorov's, but for the compressional component is substantially steeper and close to Burgers's; the flow is mostly solenoidal and statistically rigorously, which is consistent with fully developed homogeneous and isotropic turbulence. Small-scale anisotropy appears due to numerical artifact. Toward the virial region, the flow becomes increasingly compressional, the structure functions become flatter, and modest genuine anisotropy appear particularly close to the injection scale. In comparison, mesh adaptivity based on Lagrangian refinement and the same finest resolution leads to a lack of turbulent power on a small scales, an excess thereof on large scales, and unreliable density-weighted structure functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Computational Fluid Dynamics Simulation of Lignocellulosic Biomass Transport in a Compression-Screw Feeder

There is continued interest in lignocellulosic biomass feedstocks for producing power, heat, and transportation liquid fuels in order to increase energy diversification and security. However, there are a few problems associated with the biomass production that must be resolved before biorefineries are considered economical. One of these problematic operations is the process of feeding biomass to reactor systems. Compression-screw feeders have been widely used as a way to fulfill this purpose. However, due to variability of feedstock properties and its complex challenges, production plants operate at less than 50% of their design capacity due to process upsets along with mechanical wear and failure. The objective of this study is to numerically investigate the screw-feeder operation and performance for various biomass physical properties and under different operating conditions. Therefore, the findings from this study will reduce the number of physical experiments needed to improve the plant design.

09 BIOMASS FUELS↗
NASA NTRS1993

Unstructured multigrid through agglomeration

In this work the compressible Euler equations are solved using finite volume techniques on unstructured grids. The spatial discretization employs a central difference approximation augmented by dissipative terms. Temporal discretization is done using a multistage Runge-Kutta scheme. A multigrid technique is used to accelerate convergence to steady state. The coarse grids are derived directly from the given fine grid through agglomeration of the control volumes. This agglomeration is accomplished by using a greedy-type algorithm and is done in such a way that the load, which is proportional to the number of edges, goes down by nearly a factor of 4 when moving from a fine to a coarse grid. The agglomeration algorithm has been implemented and the grids have been tested in a multigrid code. An area-weighted restriction is applied when moving from fine to coarse grids while a trivial injection is used for prolongation. Across a range of geometries and flows, it is shown that the agglomeration multigrid scheme compares very favorably with an unstructured multigrid algorithm that makes use of independent coarse meshes, both in terms of convergence and elapsed times.

Venkatakrishnan, V.↗
DOE OSTI2022

Comparison of three state-of-the-art crystal plasticity based deformation twinning models for magnesium alloys

In magnesium alloys, deformation twinning and its interactions with dislocation slip are responsible for a sigmoidal shape stress–strain behavior and an asymmetrical tension–compression yield strength in magnesium alloys. The sensitivity of twinning to the underlying microstructure renders the crystal plasticity method the most commonly adopted modeling approach for magnesium-twinning. This paper compares three state-of-the-art crystal plasticity-based twinning models from the literature, namely the elastic-viscoplastic self-consistent twinning-detwinning (EVPSC-TDT) model, crystal plasticity finite element model based on enhanced predominate twin reorientation approach (CPFE-ePTR), and the crystal plasticity finite element model based on “discrete twinning” approach (CPFE-DT). A polycrystalline microstructure is simulated with all three methods to compare the resulting stress–strain curves and lattice strains to those from the experimentally measured data. All three methods showed the capability of capturing the experimental results with different levels of accuracy. Additionally, the EVPSC-TDT method avoids solving the finite element matrices and showed the highest computational efficiency. The CPFE-ePTR model shows a higher accuracy in capturing the lattice strain. The CPFE-DT relies on high-resolution finite element mesh and is much slower than the other two methods, but it captured the local deformation concentration and stress reversal phenomena near the twin band, which was not possible with the other two methods. Based on the comparison, guidance for the selection of the appropriate model based on the specific modeling target is provided in this paper.

36 MATERIALS SCIENCE↗