Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesh Adaptation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27

Betelgeuse as a Merger of a Massive Star with a Companion

We investigate the merger between a 16M ⊙ star, on its way to becoming a red supergiant (RSG), and a 4M ⊙ main-sequence companion. Our study employs three-dimensional hydrodynamic simulations using the state-of-the-art adaptive mesh refinement code Octo-Tiger. The initially corotating binary undergoes interaction and mass transfer, resulting in the accumulation of mass around the companion and its subsequent loss through the second Lagrangian point (L2). The companion eventually plunges into the envelope of the primary, leading to its spin-up and subsequent merger with the helium core. We examine the internal structural properties of the post-merger star, as well as the merger environment and the outflow driven by the merger. Our findings reveal the ejection of approximately ∼0.6 M ⊙ of material in an asymmetric and somewhat bipolar outflow. We import the post-merger stellar structure into the MESA stellar evolution code to model its long-term nuclear evolution. In certain cases, the post-merger star exhibits persistent rapid equatorial surface rotation as it evolves in the H–R diagram toward the observed location of Betelgeuse. These cases demonstrate surface rotation velocities of a similar magnitude to those observed in Betelgeuse, along with a chemical composition resembling that of Betelgeuse. In other cases, efficient rotationally induced mixing leads to slower surface rotation. This pioneering study aims to model stellar mergers across critical timescales, encompassing dynamical, thermal, and nuclear evolutionary stages.

79 ASTRONOMY AND ASTROPHYSICS↗

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗

Time-dependent finite-difference simulation of unsteady interactive flows

The solution of the time-dependent, Reynolds-averaged, Navier-Stokes equations for unsteady, interacting flows by finite-difference algorithms is discussed. Specific examples include (1) unsteady transonic flow over a thick biconvex airfoil, (2) determination of buffet boundaries for a transonic lifting airfoil, (3) the simulation of aileron buzz and (4) dynamic stall. Algorithms considered include explicit methods, mixed (or hybrid) methods, and fully implicit methods. Consideration of time scales for computational stability, computational accuracy, and physical accuracy and the use of time-dependent adaptive meshing to realize computational efficiency are also discussed.

Deiwert, G. S.↗

Adaptive grid methods for partial differential equations

A number of techniques for constructing adaptive mesh generators for use in solving partial differential equations are reviewed in this paper. Techniques reviewed include methods based on steady grid generation schemes and those which are explicitly designed to determine grid speeds in a time-dependent or space-marching problem. Results for candidate methods are included and suggestions for areas of future research are suggested.

Anderson, D. A.↗

An expanding vortex site for the r-process in rotating stellar collapse

The astrophysical site of r-process nucleosynthesis is investigated theoretically using two-dimensional expanding-vortex stellar-collapse simulations based on the Eulerian adaptive-mesh MHD code of Symbalisty (1984). The results are presented graphically, and it is found that the classical r-process can be explained as the result of the collapse of the highly rotating iron core of a 15-solar-mass star, in which the angular velocity along the rotation axis reaches a maximum of 1000 rad/s, corresponding to an angular momentum of 3.4 x 10 to the 48th erg/s for the inner 1.5 solar mass. The ejected jets are shown to yield about 0.0004 solar mass per supernova, sufficient to explain the observed abundances of r-process products.

Symbalisty, E. M. D.↗

Refined numerical solution of the transonic flow past a wedge

A numerical procedure combining the ideas of solving a modified difference equation and of adaptive mesh refinement is introduced. The numerical solution on a fixed grid is improved by using better approximations of the truncation error computed from local subdomain grid refinements. This technique is used to obtain refined solutions of steady, inviscid, transonic flow past a wedge. The effects of truncation error on the pressure distribution, wave drag, sonic line, and shock position are investigated. By comparing the pressure drag on the wedge and wave drag due to the shocks, a supersonic-to-supersonic shock originating from the wedge shoulder is confirmed.

Liang, S.-M.↗

An adaptive finite element procedure for compressible high speed flows

A practical finite element-based solution procedure for high-speed inviscid compressible flow problems is described. The method provides time-accurate solutions to the coompressible Euler equations, and is computationally more efficient than the one-step Taylor-Galerkin approach and better suited for implementation on the modern generation of vector computers. The method is coupled to an adaptive mesh refinement process that enables steady state solutions of improved quality to be obtained.

Loehner, R.↗

Vectorizable algorithms for adaptive schemes for rapid analysis of SSME flows

An initial study into vectorizable algorithms for use in adaptive schemes for various types of boundary value problems is described. The focus is on two key aspects of adaptive computational methods which are crucial in the use of such methods (for complex flow simulations such as those in the Space Shuttle Main Engine): the adaptive scheme itself and the applicability of element-by-element matrix computations in a vectorizable format for rapid calculations in adaptive mesh procedures.

Oden, J. Tinsley↗

Recent developments in FEM-CFD

The current status of CFD with regard to unstructured grids employing finite element methods and Eulerian frames is reviewed. Algorithms suitable for the computation of large three-dimensional problems involving flow past arbitrary geometries are developed. Adaptive mesh refinement strategy is reviewed, and domain splitting or local time-stepping are briefly addressed. The development of search algorithms of optimal order, variable time-stepping Jacobi smoothers for elliptic problems, and transport concepts for hyperbolics to help achieve good performance for unstructured multigrid processes is discussed. As examples, transient supersonic flow in a channel, regular shock reflection of a wall, viscous flow past a protruberance, potential flow past a cylinder, and Burgers equation are considered.

Loehner, R.↗

Computational solid mechanics (finite elements and boundary elements) - Present status and future directions

Computational finite-element and boundary-element methods are reviewed, and their application to the mechanics of solids is discussed. Stability conditions for general FEMs are considered in addition to the use of least-order, stable, invariant, or hybrid/mixed isoparametric elements as alternatives to the displacement-based isoparametric elements. The use of symbolic manipulation, adaptive mesh refinement, transient dynamic response, and boundary-element methods for linear elaslticity and finite-strain problems of inelastic materials are also discussed.

Atluri, S. N.↗

An adaptive finite element method for high speed flows

The solution of the equations of compressible high speed flow, on unstructured triangular grids in 2D and tetrahedral grids in 3D, is considered. Solution methods based upon both Taylor-Galerkin and Runge-Kutta time-stepping techniques are presented and the incorporation of the ideas of flux corrected transport (FCT) is discussed. These methods are combined with an adaptive mesh regeneration procedure and are employed in the solution of several examples, consisting of Euler flows in both 2D and 3D and Navier-Stokes flows in 2D.

Peraire, J.↗

Convection dominated problems

The paper surveys the last ten years of activity of the INME Swansea, dealing with problems of convection dominated flow. The basic explicit/implicit characteristic Galerkin process and its application to adaptive mesh refinement used in the solution of realistic problems is focused on.

Peraire, J.↗

Finite element methods of analysis for high speed viscous flows

Over the past three years a finite element based procedure for the solution of high speed viscous compressible flows was developed. The approach followed was to compute steady state solutions via a false transient, using an explicit time stepping scheme, and to attempt to improve the solution quality by incorporating adaptive mesh procedures. The main thrust of the work was to continue on the extension of the approach to the solution of some realistic compressible viscous flows. When flows at high Reynolds number are investigated, is soon becomes apparent that explicit techniwues have to be supplemented if they are to deal effectively with the large variations in element size and aspect ratio which characterize the computational grids necessary for adequate resolutin of the primary flow features. For this reason, the Taylor-Galerkin solution algorithm was rewritten in an explicit/implicit form. Solutions were computed for the problems of a flow past a flat plate, M-3, Re-100; shock/boundary layer interaction, M-2, Re-296000; flow over a compression corner, M-11.68, Re246000; and unifrom flow past a circular cylinder, M-6.34, Re-39770. A summary of the results is included and demonstrates the numerical performance of the scheme.

Source record↗

Adaptive remeshing for compressible flow computations

The present, quality-enhancing adaptive-mesh procedure for two-dimensional Euler equation steady state solutions is implemented by means of linear triangular elements and an explicit time-stepping scheme, in conjunction with a finite element solution algorithm. The meshes thus generated typically take the form of stretched elements in the vicinity of one-dimensional flow features; a considerable variation in element size may thereby emerge which allows the desired high-quality solutions to be obtained with commensurately high computational efficiency.

Peraire, J.↗

An error estimation procedure for plate bending elements

Procedures for identifying and eliminating errors inherent in individual finite elements and those due to the discretization of the continuum are presented. The elemental errors are identified through the use of an element formulation procedure based on physically interpretable strain gradient interpolation functions. The use of physically interpretable notation allows these errors to be eliminated using rational arguments. The discretization errors are identified by comparing the finite-element solution with a smoothed superconvergent solution. The errors thus identified are used to guide an adaptive mesh refinement procedure which produces improved results.

Dow, John O.↗

Numerical investigation of chemically reacting flows in ramjet dump combustors

The time-dependent Navier-Stokes equations, including second-order turbulence model, are numerically integrated by using four-stage Runge-Kutta scheme to predict the steady-state supersonic flow structures in ramjet dump combustors. The formulation is derived for reacting flows with finite-rate chemistry. In the present study, it is firstly attempted to assess the accuracy of existing high-order turbulence model in supersonic flows. The comparison shows reasonable agreement between calculated and measured data in terms of velocity distributions. It is indicated that a modified constant C-mu for calculating turbulent eddy viscosity is needed in the supersonic flow regime and the adaptive meshing is preferred to capture the recirculation zone. In the reacting flow calculation, the results from a test case of hydrogen and air combustion at premixed conditon show that the rearward facing step is able to increase flow residence time and stabilize the flame in supersonic flows.

Hsieh, Kwang-Chung↗

Finite element analysis in fluids; Proceedings of the Seventh International Conference on Finite Element Methods in Flow Problems, University of Alabama, Huntsville, Apr. 3-7, 1989

Recent advances in computational fluid dynamics are examined in reviews and reports, with an emphasis on finite-element methods. Sections are devoted to adaptive meshes, atmospheric dynamics, combustion, compressible flows, control-volume finite elements, crystal growth, domain decomposition, EM-field problems, FDM/FEM, and fluid-structure interactions. Consideration is given to free-boundary problems with heat transfer, free surface flow, geophysical flow problems, heat and mass transfer, high-speed flow, incompressible flow, inverse design methods, MHD problems, the mathematics of finite elements, and mesh generation. Also discussed are mixed finite elements, multigrid methods, non-Newtonian fluids, numerical dissipation, parallel vector processing, reservoir simulation, seepage, shallow-water problems, spectral methods, supercomputer architectures, three-dimensional problems, and turbulent flows.

Chung, T. J.↗

Numerical simulations in combustion

This paper reviews numerical simulations in reacting flows in general and combustion phenomena in particular. It is shown that use of implicit schemes and/or adaptive mesh strategies can improve convergence, stability, and accuracy of the solution. Difficulties increase as turbulence and multidimensions are considered, particularly when finite-rate chemistry governs the given combustion problem. Particular attention is given to the areas of solid-propellant combustion dynamics, turbulent diffusion flames, and spray droplet vaporization.

Chung, T. J.↗