Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Volume”

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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

NEAMS Technical Area Support in MOOSE

The MOOSE framework is a foundational capability used by the NEAMS program to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and other new capabilities. The work transcribed in this report was conducted in direct support of the simulation tools and has already been deployed. The capabilities outlined in this report include enabling selective polynomial basis refinement, implementing a custom convergence system, building a scalable preconditioner for saddle-point problems, and much more.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Increased accuracy of multiphysics simulations through flexible execution, transient algorithms, and modular physics

The MOOSE framework is a foundational capability used by the NEAMS program to create over 15 different simulation tools for advanced nuclear reactors. Due to MOOSE’s broad use, improvements to the framework in support of modeling and simulation goals are critical to the program. Such improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and new capabilities. The work described in this report was conducted in direct support of the simulation tools and has already been deployed. The capabilities were implemented in the same order as they are covered in this report: multiple time integrators in the same input file, initial design of framework Components, an input file Application block, extension of NetGen to 3D geometries in MOOSE, and deployment of executors in the multi-system paradigm. These five additions are fundamental capabilities that will be leveraged by many NEAMS applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

NEAMS Technical Area Support in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) framework is a foundational capability used by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and other new capabilities. The work described in this report was conducted in direct support of the simulation tools and has already been deployed. The capabilities outlined in this report include implementing hash table matrix assembly for efficient sparsity pattern construction for contact in BISON, developing re-step testing infrastructure for ensuring the viability of overlapping domain coupling between SAM and Pronghorn, allowing unique preconditioners for single-input multi-system solves, supporting multi-system in MOOSE’s workhorse executioners, and many more smaller feature enhancements and bug fixes.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Benchmarking numerical methods for impact and cratering applications

Large scale computational models are important for studying impact cratering events that are prevalent both on Earth and, more broadly, in this solar system. To address these problems, models must reliably account for both large length scales (e.g., kilometers) and relatively long time scales (hundreds of seconds). This work benchmarks two such approaches, a more traditional hydrodynamics approach and a finite-discrete element method (FDEM), for impact cratering applications. Both 2D and 3D results are discussed for two different impact velocities, 5 km/s and 20 km/s, striking normal to the target and, for 3D simulations, 45° from vertical. In addition, comparisons to previously published data are presented. Finally, differences in how these methods model damage are discussed. Ultimately, both approaches show successful modeling of several different impact scenarios.

58 GEOSCIENCES↗
DOE OSTI2020

Numerical Considerations in the Modeling of a High Explosive Cylinder Experiment Using an ALE Continuum Mechanics Code

A series of experiments involving the detonation of PBX 9501 encased in a copper cylinder are modeled with the objective of evaluating a proposed set of phenomenological parameters for the Wescott–Stewart–Davis reactive burn model. The numerical analysis is conducted using the Los Alamos continuum mechanics code FLAG. Numerical considerations pertaining to various aspects of modeling the experiments using FLAG are discussed. It is shown that use of the proposed set of phenomenological parameters results in predictions of free-surface velocity that match empirically measured velocities reasonably well.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Localized Exponential Time Differencing Method for Shallow Water Equations: Algorithms and Numerical Study

Here, we explore the performance of the exponential time differencing (ETD) method applied to the rotating shallow water equations. Comparing with explicit time stepping of the same order accuracy in time, the ETD algorithms could reduce the computational time in many cases by allowing the use of large time step sizes while still maintaining numerical stability. To accelerate the ETD simulations, we propose a localized approach that synthesizes the ETD method and overlapping domain decomposition. By dividing the original problem into many subdomain problems of smaller sizes and solving them locally, the proposed approach could speed up the calculation of matrix exponential vector products. Several standard test cases for shallow water equations of one or multiple layers are considered. The results show great potential of the localized ETD method for high-performance computing because each subdomain problem can be naturally solved in parallel at every time step.

58 GEOSCIENCES↗
NASA NTRS1972

Three-dimensional plasma diffusion in a very strong magnetic field.

The thermal equilibrium coefficient of spatial diffusion transverse to a strong uniform dc magnetic field is calculated for a fully ionized plasma. The particles are assumed to move transverse to the field only as a consequence of the E x B drift, but to move freely parallel to it. The particles interact only electrostatically. The calculation is done at large, but fixed and finite, plasma volume. It is shown that, as magnetic field density approaches infinity, the leading term of the coefficient of transverse spatial diffusion falls off in a prescribed manner, but contains a multiplicative factor which goes to zero as the plasma volume becomes infinite. The method of calculation fails for unbounded plasmas.

Montgomery, D.↗
NASA NTRS1979

TF34 engine compression system computer study

The stability of the fan and the compressor components was examined individually using linearized and time dependent, one dimensional stability analysis techniques. The stability of the fan core integrated compression system was investigated using a two dimensional compression system model. The analytical equations on which this model was based satisfied the mass, axial momentum, radial momentum, and energy conservation equations for flow through a finite control volume. The results gave an accurate simulation of the flow through the compression system. The speed lines of the components were reproduced; the points of instability were accurately predicted; the locations where the instability was initiated in the fan and the core were indicated; and the variation of the bypass ratio during flow throttling was calculated. The validity of the analytical techniques was then established by comparing these results with test data and with results obtained from the steady state cycle deck.

Hosny, W. M.↗
NASA NTRS1986

Determination of U, V, and W from single station Doppler radar radial velocities

The ST/MST (stratosphere troposphere/mesosphere stratosphere troposphere) clear air Doppler radar, or wind profiler, is an important tool in observational meteorology because of its capability to remote observe dynamic parameters of the atmosphere. There are difficulties in transforming the observed radial velocities into meteorological wind components. How this problem has been treated in the past is reviewed, and some of the analysis is recast to a form more suited to the high diagnostic abilities of a number of fixed beam configurations with reference to a linear wind field. The results, in conjunction with other works which treats problems such as the effects of finite sample volumes in the presence of nonhomogeneous atmospheric reflectivity, have implications important to the design of both individual MST/ST radars and MST/ST radar networks. The key parameters to uncoupling terms in the scaling equations are w sub x and w sub y. Whenever the stratiform condition, which states that these two parameters are negligible, is satisfied, a five beam ST radar may determine unbiased values of u, v, and w for sample volumes directly above the radar. The divergence and partial deformation of the flow may also be determined. Three beam systems can determine w and w sub z, but are unable to obtain u and v wind components uncontaminated by vertical sheer terms, even when the stratiform condition is satisfied.

Clark, W. L.↗
NASA NTRS1989

Control-volume based Navier-Stokes equation solver valid at all flow velocities

A control-volume based finite difference method to solve the Reynolds averaged Navier-Stokes equations is presented. A pressure correction equation valid at all flow velocities and a pressure staggered grid layout are used in the method. Example problems presented herein include: a developing laminar channel flow, developing laminar pipe flow, a lid-driven square cavity flow, a laminar flow through a 90-degree bent channel, a laminar polar cavity flow, and a turbulent supersonic flow over a compression ramp. A k-epsilon turbulence model supplemented with a near-wall turbulence model was used to solve the turbulent flow. It is shown that the method yields accurate computational results even when highly skewed, unequally spaced, curved grids are used. It is also shown that the method is strongly convergent for high Reynolds number flows.

Kim, S.-W.↗
NASA NTRS1989

Numerical computation of shock wave-turbulent boundary layer interaction in transonic flow over an axisymmetric curved hill

A control-volume based finite difference computation of a turbulent transonic flow over an axisymmetric curved hill is presented. The numerical method is based on the SIMPLE algorithm, and hence the conservation of mass equation is replaced by a pressure correction equation for compressible flows. The turbulence is described by a k-epsilon turbulence model supplemented by a near-wall turbulence model. In the method, the dissipation rate in the region very close to the wall is obtained from an algebraic equation and that for the rest of the flow domain is obtained by solving a partial differential equation for the dissipation rate. The other flow equations are integrated up to the wall. It is shown that the present turbulence model yields the correct location of the compression shock. The other computational results are also in good agreement with experimental data.

Kim, S.-W.↗
NASA NTRS1989

Calculation of reattaching shear layers in divergent channel with a multiple-time-scale turbulence model

Numerical calculations of turbulent reattaching shear layers in a divergent channel are presented. The turbulence is described by a multiple-time-scale turbulence model. The turbulent flow equations are solved by a control-volume based finite difference method. The computational results are compared with those obtained using k-epsilon turbulence models and algebraic Reynolds stress turbulence models. It is shown that the multiple-time-scale turbulence model yields significantly improved computational results than the other turbulence models in the region where the turbulence is in a strongly inequilibrium state.

Kim, S.-W.↗
NASA NTRS1989

Numerical investigation of an internal layer in turbulent flow over a curved hill

The development of an internal layer in a turbulent boundary layer flow over a curved hill is investigated numerically. The turbulence field of the boundary layer flow over the curved hill is compared with that of a turbulent flow over a symmetric airfoil (which has the same geometry as the curved hill except that the leading and trailing edge plates were removed) to study the influence of the strongly curved surface on the turbulence field. The turbulent flow equations are solved by a control-volume based finite difference method. The turbulence is described by a multiple-time-scale turbulence model supplemented with a near-wall turbulence model. Computational results for the mean flow field (pressure distributions on the walls, wall shearing stresses and mean velocity profiles), the turbulence structure (Reynolds stress and turbulent kinetic energy profiles), and the integral parameters (displacement and momentum thicknesses) compared favorably with the measured data. Computational results show that the internal layer is a strong turbulence field which is developed beneath the external boundary layer and is located very close to the wall. Development of the internal layer was more obviously observed in the Reynolds stress profiles and in the turbulent kinetic energy profiles than in the mean velocity profiles. In this regard, the internal layers is significantly different from wall-bounded simple shear layers in which the mean velocity profile characterizes the boundary layer most distinguishably. Development of such an internal layer, characterized by an intense turbulence field, is attributed to the enormous mean flow strain rate caused by the streamline curvature and the strong pressure gradient. In the turbulent flow over the curved hill, the internal layer begin to form near the forward corner of the hill, merges with the external boundary layer, and develops into a new fully turbulent boundary layer as the fluid flows in the downstream direction. For the flow over the symmetric airfoil, the boundary layer began to form from almost the same location as that of the curved hill, grew in its strength, and formed a fully turbulent boundary layer from mid-part of the airfoil and in the downstream region. Computational results also show that the detailed turbulence structure in the region very close to the wall of the curved hill is almost the same as that of the airfoil in most of the curved regions except near the leading edge. Thus the internal layer of the curved hill and the boundary layer of the airfoil were also almost the same. Development of the wall shearing stress and separation of the boundary layer at the rear end of the curved hill mostly depends on the internal layer and is only slightly influenced by the external boundary layer flow.

Kim, S-W.↗
NASA NTRS1989

A numerical analysis of Stefan problems for generalized multi-dimensional phase-change structures using the enthalpy transforming model

An enthalpy transforming scheme is proposed to convert the energy equation into a nonlinear equation with the enthalpy, E, being the single dependent variable. The existing control-volume finite-difference approach is modified so it can be applied to the numerical performance of Stefan problems. The model is tested by applying it to a three-dimensional freezing problem. The numerical results are in agreement with those existing in the literature. The model and its algorithm are further applied to a three-dimensional moving heat source problem showing that the methodology is capable of handling complicated phase-change problems with fixed grids.

Cao, Yiding↗
NASA NTRS1990

Calculation of reattaching shear layers in divergent channel with a multiple-time-scale turbulence model

Numerical calculations of turbulent reattaching shear layers in a divergent channel are presented. The turbulence is described by a multiple-time-scale turbulence model. The turbulent flow equations are solved by a control-volume based finite difference method. The computational results are compared with those obtained using k-epsilon turbulence models and algebraic Reynolds stress turbulence models. It is shown that the multiple-time-scale turbulence model yields significantly improved computational results than the other turbulence models in the region where the turbulence is in a strongly inequilibrium state.

Kim, S.-W.↗
NASA NTRS1990

Thermal buoyancy and Marangoni convection in a two fluid layered system - A numerical study

Thermal buoyancy and surface tension driven convection is numerically investigated in a system with two, immiscible fluids. The geometry investigated, has an open cavity configuration with the lighter fluid situated on top of the heavier fluid, forming a stable layered system. The upper fluid meniscus and the interface are assumed to be flat and undeformable in the calculations. The governing equations and boundary/interface conditions are solved by a control volume based finite difference scheme for two pairs of immiscible fluids; the water-hexadecane system and a so called generic system. The steady state calculations show, that dramatically different flows are predicted when the interfacial tension effects are included or excluded from the system model. These differences are particularly appreciable in surface tension dominated flows, that are typical of microgravity situations. Complex flow patterns, with induced secondary flows are noticed in both the fluids. In general, the overall system heat transfer is found to increase with increases in the thermal buoyancy and surface tension effects, but the behavior of the system flow and thermal fields is not easy to characterize, when different combinations of these forces are considered.

Ramachandran, N.↗
NASA NTRS1990

Small domes on Venus - Characteristics and origin

The areal and size-frequency distribution, abundance, individual characteristics, and geologic and regional associations of small domes in the Venera data set are examined. It is noted that, because of their numbers and widespread occurence, regardless of origin, these structures will be important in the local and global geographic interpretation of the surface. It is concluded that small domes of 20 km in diameter occur in numbers of order 10,000 on the northern quarter of the surface of Venus, and that small domes are interpreted to be predominantly low shield volcanos and to represent multiple centralized effusive eruptions of discrete volumes over finite periods. It is noted that small domes occur throughout the northern latitudes, but attain maximum regional concentrations in two areas, one located north of Thetis Regio in Akkruva Colles, and another located on the northeast flanks of Beta Regio in Guinevere Planitia. The presence of two global concentrations implies that magma production and eruption were enhanced within these two areas.

Aubele, Jayne C.↗