2014 U.S. DRIVE Highlight - Local ALE Method or Moving Parts in Engine Simulations
Abstract Not Provided
2014 U.S. DRIVE Highlight Local ALE Method or Movi↗
SEARCH · Search NASA
Results for “ALE methods”
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.
A 3D finite element ALE method using an approximate Riemann solution
Arbitrary Lagrangian–Eulerian finite volume methods that solve a multidimensional Riemann-like problem at the cell center in a staggered grid hydrodynamic (SGH) arrangement have been proposed. This research proposes a new 3D finite element arbitrary Lagrangian–Eulerian SGH method that incorporates a multidimensional Riemann-like problem. Here, two different Riemann jump relations are investigated. A new limiting method that greatly improves the accuracy of the SGH method on isentropic flows is investigated. A remap method that improves upon a well-known mesh relaxation and remapping technique in order to ensure total energy conservation during the remap is also presented. Numerical details and test problem results are presented.
97 MATHEMATICS AND COMPUTING↗
Second-invariant-preserving Remap of the 2D deviatoric stress tensor in ALE methods
Not Available
Klíma, Matěj (ORCID:0000000324967614)↗
Towards an IMEX Monolithic ALE Method with Integrated UQ for Multiphysics Shock-hydro
A number of critical science, engineering and advanced technology applications require predictive analysis of complex shock-hydrodynamics of fluid/solid materials with possible electromagnetic interaction. The physical mechanisms include wave-phenomena, material-transport, diffusion, chemical reactions, and electromagnetics. The highly nonlinear multiple-time and length-scale response of these mechanisms include discontinuities formed from shocks, contact surfaces, and complex tabular equations-of-state (EOS). Current dominant computational solution strategies use ad-hoc combinations of operator- splitting, semi-implicit, and explicit time-integration methods and decoupled nonlinear-solvers. While these approaches have enabled progress in forward simulation, the inherited mathematical structure has not provided stability, accuracy and efficiency to resolve all the dynamical time-scales of interest, nor has it enabled integrated fast sensitivity analysis and uncertainty quantification (UQ). This draft report describes initial progress towards developing a new multiphysics shock-hydro capability that is intended to be more robust, mathematically well-structured and can readily be combined with advanced higher-order implicit/explicit (IMEX) time integration and efficient adjoint-enhanced uncertainty quantification (UQ) techniques.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
A coupled ALE-AMR method for shock hydrodynamics
We present a numerical method combining adaptive mesh refinement (AMR) with arbitrary Lagrangian-Eulerian (ALE) mesh motion for the simulation of shock hydrodynamics on unstructured grids. The primary goal of the coupled method is to use AMR to reduce numerical error in ALE simulations at reduced computational expense relative to uniform fine mesh calculations, in the same manner that AMR has been used in Eulerian simulations. We also identify deficiencies with ALE methods that AMR is able to mitigate, and discuss the unique coupling challenges. The coupled method is demonstrated using three-dimensional unstructured meshes of up to O(10 7 ) tetrahedral cells. Convergence of ALE-AMR solutions towards both uniform fine mesh ALE results and analytic solutions is demonstrated. Speed-ups of 5-10× for a given level of error are observed relative to uniform fine mesh calculations.
42 ENGINEERING↗
Challenges in atomic layer etching of gallium nitride using surface oxidation and ligand-exchange
Two atomic layer etching (ALE) methods were studied for crystalline GaN, based on oxidation, fluorination, and ligand exchange. Etching was performed on unintentionally doped GaN grown by hydride vapor phase epitaxy. For the first step, the GaN surfaces were oxidized using either water vapor or remote O 2 -plasma exposure to produce a thin oxide layer. Removal of the surface oxide was addressed using alternating exposures of hydrogen fluoride (HF) and trimethylgallium (TMG) via fluorination and ligand exchange, respectively. Several HF and TMG super cycles were implemented to remove the surface oxide. Each ALE process was monitored in situ using multiwavelength ellipsometry. X-ray photoelectron spectroscopy was employed for the characterization of surface composition and impurity states. Additionally, the thermal and plasma-enhanced ALE methods were performed on patterned wafers and transmission electron microscopy (TEM) was used to measure the surface change. The x-ray photoelectron spectroscopy measurements indicated that F and O impurities remained on etched surfaces for both ALE processes. Ellipsometry indicated a slight reduction in thickness. TEM indicated a removal rate that was less than predicted. In conclusion, we suggest that the etch rates were reduced due to the ordered structure of the oxide formed on crystalline GaN surfaces.
36 MATERIALS SCIENCE↗
Three-dimensional local ALE-FEM method for fluid flow in domains containing moving boundaries/objects interfaces
A three-dimensional finite element method for the numerical simulations of fluid flow in domains containing moving rigid objects or boundaries is developed. The method falls into the general category of Arbitrary Lagrangian Eulerian methods; it is based on a fixed mesh that is locally adapted in the immediate vicinity of the moving interfaces and reverts to its original shape once the moving interfaces go past the elements. The moving interfaces are defined by separate sets of marker points so that the global mesh is independent of interface movement and the possibility of mesh entanglement is eliminated. The results is a fully robust formulation capable of calculating on domains of complex geometry with moving boundaries or devises that can also have a complex geometry without danger of the mesh becoming unsuitable due to its continuous deformation thus eliminating the need for repeated re-meshing and interpolation. Moreover, the boundary conditions on the interfaces are imposed exactly. This work is intended to support the internal combustion engines simulator KIVA developed at Los Alamos National Laboratories. The model's capabilities are illustrated through application to incompressible flows in different geometrical settings that show the robustness and flexibility of the technique to perform simulations involving moving boundaries in a three-dimensional domain.
16 TIDAL AND WAVE POWER↗
Dynamic mesh adaptation for front evolution using discontinuous Galerkin based weighted condition number relaxation
A new mesh smoothing method designed to cluster cells near a dynamically evolving interface is presented. The method is based on weighted condition number mesh relaxation with the weight function computed from a level set representation of the interface. The weight function is expressed as a Taylor series based discontinuous Galerkin projection, which makes the computation of the derivatives of the weight function needed during the condition number optimization process a trivial matter. For cases when a level set is not available, a fast method for generating a low-order level set from discrete cell-centered fields, such as a volume fraction or index function, is provided. Results show that the low-order level set works equally well as the actual level set for mesh smoothing. Meshes generated for a number of interface geometries are presented, including cases with multiple level sets. Lastly, dynamic cases with moving interfaces show the new method is capable of maintaining a desired resolution near the interface with an acceptable number of relaxation iterations per time step, which demonstrates the method's potential to be used as a mesh relaxer for arbitrary Lagrangian Eulerian (ALE) methods.
97 MATHEMATICS AND COMPUTING↗
Adaptive reconnection-based arbitrary Lagrangian Eulerian method
We present a new adaptive Arbitrary Lagrangian Eulerian (ALE) method. This method is based on the reconnection-based ALE (ReALE) methodology of Refs. [35], [34] and [6]. The main elements in a standard ReALE method are: an explicit Lagrangian phase on an arbitrary polygonal (in 2D) mesh in which the solution and positions of grid nodes are updated; a rezoning phase in which a new grid is defined by changing the connectivity (using Voronoi tessellation) but not the number of cells; and a remapping phase in which the Lagrangian solution is transferred onto the new grid. Furthermore, in the standard ReALE method, the rezoned mesh is smoothed by using one or several steps toward centroidal Voronoi tessellation, but it is not adapted to the solution in any way.
97 MATHEMATICS AND COMPUTING↗
Dynamic Mesh Adaptation for Front Evolution Using Discontinuous Galerkin Based Weighted Condition Number Mesh Relaxation
A new mesh smoothing method designed to cluster mesh cells near a dynamically evolving interface is presented. The method is based on weighted condition number mesh relaxation with the weight function being computed from a level set representation of the interface. The weight function is expressed as a Taylor series based discontinuous Galerkin projection, which makes the computation of the derivatives of the weight function needed during the condition number optimization process a trivial matter. For cases when a level set is not available, a fast method for generating a low-order level set from discrete cell-centered elds, such as a volume fraction or index function, is provided. Results show that the low-order level set works equally well for the weight function as the actual level set. Meshes generated for a number of interface geometries are presented, including cases with multiple level sets. Dynamic cases for moving interfaces are presented to demonstrate the method's potential usefulness to arbitrary Lagrangian Eulerian (ALE) methods.
97 MATHEMATICS AND COMPUTING↗
Arbitrary Lagrangian–Eulerian methods for modeling high-speed compressible multimaterial flows
This paper reviews recent developments in Arbitrary Lagrangian Eulerian (ALE) methods for modeling high speed compressible multimaterial flows in complex geometry on general polygonal meshes. We only consider the indirect ALE approach which consists of three key stages: a Lagrangian stage, in which the solution and the computational mesh are updated; a rezoning stage, in which the nodes of the computational mesh are moved to improve grid quality; and a remapping stage, in which the Lagrangian solution is transferred to the rezoned mesh.
97 MATHEMATICS AND COMPUTING↗
Compatible, energy conserving, bounds preserving remap of hydrodynamic fields for an extended ALE scheme
From the very origins of numerical hydrodynamics in the Lagrangian work of von Neumann and Richtmyer [83], the issue of total energy conservation as well as entropy production has been problematic. Because of well known problems with mesh deformation, Lagrangian schemes have evolved into Arbitrary Lagrangian–Eulerian (ALE) methods [39] that combine the best properties of Lagrangian and Eulerian methods. Energy issues have persisted for this class of methods. We believe that fundamental issues of energy conservation and entropy production in ALE require further examination. The context of the paper is an ALE scheme that is extended in the sense that it permits cyclic or periodic remap of data between grids of the same or differing connectivity. The principal design goals for a remap method then consist of total energy conservation, bounded internal energy, and compatibility of kinetic energy and momentum. We also have secondary objectives of limiting velocity and stress in a non-directional manner, keeping primitive variables monotone, and providing a higher than second order reconstruction of remapped variables. Particularly, the new contributions fall into three categories associated with: energy conservation and entropy production, reconstruction and bounds preservation of scalar and tensor fields, and conservative remap of nonlinear fields. Our paper presents a derivation of the methods, details of implementation, and numerical results for a number of test problems. The methods requires volume integration of polynomial functions in polytopal cells with planar facets, and the requisite expressions are derived for arbitrary order.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions
A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Constrained optimization framework for interface-aware sub-scale dynamics models for voids closure in Lagrangian hydrodynamics
In hydrocodes, voids are used to represent vacuum and model free boundaries between vacuum and real materials. We give a systematic description of a new treatment of void closure in the framework of the multimaterial arbitrary Lagrangian–Eulerian (ALE) methods. This includes a new formulation of the interface-aware sub-scale-dynamics (IA-SSD) closure model for multimaterial cells with voids, which is used in the Lagrangian stage of our indirect ALE scheme. The results of the comprehensive testing of the new model are presented for one- and two-dimensional multimaterial calculations in the presence of voids. Finally, we also present a sneak peek of a realistic shaped charge calculation in the presence of voids and solids.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
A massively-parallel, unstructured overset method to simulate moving bodies in turbulent flows
An unstructured overset method capable of performing direct numerical simulation (DNS) and large eddy simulation (LES) of many (O(10 5 )) moving bodies, utilizing many computational cores (O(10 5 )), in turbulent, incompressible fluid flow is presented. Unstructured meshes are attached to bodies and placed within a fixed background domain. Body meshes are allowed to arbitrarily overlap and move throughout the domain. Within each mesh a high resolution, unstructured, non-dissipative finite volume method is used to solve for the flow field. Boundary conditions for each mesh are provided by interpolation from flow solutions on overlapping meshes. When many unstructured meshes of different resolution overlap, care is required in the connection between the different flow solutions. An interpolant is created which seeks to preserve volume conservation of flow quantities between meshes regardless of overlapping mesh differences. An implicit fractional step method is used for time advancement, requiring the calculation of a predicted fluid velocity and corrector pressure field. For the predictor step, the resulting interpolation is directly introduced into the implicit equations for the predicted flow field. For the corrected pressure field, the continuity between meshes is weakly enforced using a penalty formulation. The pressure formulation is symmetric, positive-definite and non-singular resulting in a formulation which is readily solvable using traditional iterative matrix inversion techniques. An Arbitrary Euler-Lagrangian (ALE) method coupled to a 6 degrees of freedom rigid body equation system (6-DOF) is used for body motion. For rotation, a quaternion representation is used to solve Euler's equations of rigid body motion. A linear spring damper model, which uses geometry information readily available from the overset assembly process, is used for collisions. Validation of the method for canonical flow fields is presented including assessment of order of accuracy and kinetic energy conservation properties. Particle-resolved direct numerical simulation (PR-DNS) of single particles in various flow fields are presented for validation. PR-DNS results of 50,000 spherical particles freely moving within turbulent channel flow are shown as a demonstration of the method at full scale. LES results of a marine propeller under crashback conditions are shown to demonstrate the ability to simulate highly unsteady turbulent flows over complex, moving geometries.
97 MATHEMATICS AND COMPUTING↗
Adaptive laboratory evolution for increased temperature tolerance of the diatom Nitzschia inconspicua
Abstract Outdoor microalgal cultivation for the production of valuable biofuels and bioproducts typically requires high insolation and strains with high thermal (>37°C) tolerance. While some strains are naturally thermotolerant, other strains of interest require improved performance at elevated temperatures to enhance industrial viability. In this study, adaptive laboratory evolution (ALE) was performed for over 300 days using consecutive 0.5°C temperature increases in a constant temperature incubator to attain greater thermal tolerance in the industrially relevant diatom Nitzschia inconspicua str. Hildebrandi. The adapted strain was able to grow at a constant temperature of 37.5°C; whereas this constant temperature was lethal to the parental control, which had an upper‐temperature boundary of 35.5°C before adaptive evolution. Several high‐temperature clonal isolates were obtained from the evolved population following ALE, and increased temperature tolerance was observed in the clonal, parent, and non‐clonal adapted cultures. This ALE method demonstrates the development of enhanced industrial algal strains without the production of genetically modified organisms (GMOs).
59 BASIC BIOLOGICAL SCIENCES↗
A CFD Study Assisted with Experimental Confirmation for Liquid Shape Control of Electromagnetically Levitated Bulk Materials
Two types of optimized electromagnetic levitators were designed to achieve the free surface shape control of metallic melts as the combined results of the Lorentz force, sample gravity, and surface tension. The levitation behavior of bulk melts was investigated using computational fluid dynamics (CFD) modeling coupled with high-frequency electromagnetic field analysis and the arbitrary Lagrangian–Eulerian (ALE) method. The difference in the oscillation behavior between the solid and molten samples was explained by the damping of the electromagnetic induction and liquid viscosity. The motion of the mass center, melt shape, flow pattern, and Lorentz force in the levitated melt were determined within a wide excitation current range. With the increase in the applied current, the melt’s centroid position rose sharply in the area with low current but displayed a slow increase in the area with high current. Meanwhile, the stable shape of the bulk melt showed the typical transition from a long taper through a short taper and then into a rhombus. The internal flow pattern transformed from a simple double-loop structure to a complex configuration with three or four loops. The dependence of the deformation on the Bond number was analyzed in the two types of levitators. In addition, the stable shape and swing process of the bulk Al melt within the two types of electromagnetic levitation (EML) systems were quantitatively studied under protective inert gas conditions. The melt contour could be well described by the 10th Legendre polynomial function with a deviation of less than 0.5 pct.
36 MATERIALS SCIENCE↗
Establishing mesh topology in multi-material cells: enabling technology for robust and accurate multi-material simulations
Real world problems are typically multi-material, combining materials such as gases, liquids and solids that have very different properties. The material interfaces may be fixed in time or can be a part of the solution, as in fluid-structure interactions or air-water dynamics, and therefore move and change shape. In such problems the computational mesh may be non-conformal to interfaces due to complexity of these interfaces, presence of small fractions of materials, or because the mesh does not move with the flow, as in the arbitrary Lagrangian–Eulerian (ALE) methods. In order to solve problems of interest on such meshes, interface reconstruction methods are usually used to recover an approximation of material regions within the cells. For a cell intersecting multiple material regions, these approximations of contained subregions can be considered as single-material subcells in a local mesh that we call a minimesh. In this paper, we discuss some of the requirements that discretization methods have on topological information in the resulting hierarchical meshes and present an approach that allows incorporating the buildup of sufficiently detailed topology into the nested dissections based PLIC-type reconstruction algorithms (e.g. Volume-of-Fluid, Moment-of-Fluid) in an efficient and robust manner. Specifically, we describe the X-MOF interface reconstruction algorithm in 2D, which extends the Moment-Of-Fluid (MOF) method to include the topology of minimeshes created inside of multi-material cells and parent-child relations between corresponding mesh entities on different hierarchy levels. X-MOF retains the property of being local to a cell and not requiring external communication, which makes it suitable for massively parallel applications. Here, we demonstrate some scaling results for the X-MOF implementation in Tangram, a modern interface reconstruction framework for exascale computing.
97 MATHEMATICS AND COMPUTING↗