SEARCH · Engineering Papers

Results for “finite scale equations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1978

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗

NASA NTRS1991

Development of turbulence models for shear flows by a double expansion technique

Turbulence models are developed by supplementing the renormalization group (RNG) approach of Yakhot and Orszag with scale expansions for the Reynolds stress and production of dissipation terms. The additional expansion parameter (eta) is the ratio of the turbulent to mean strain time scale. While low-order expansions appear to provide an adequate description of the Reynolds stress, no finite truncation of the expansion for the production of dissipation term in powers of eta suffices - terms of all orders must be retained. Based on these ideas, a new two-equation model and Reynolds stress transport model are developed for turbulent shear flows. The models are tested for homogeneous shear flow and flow over a backward facing step. Comparisons between the model predictions and experimental data are excellent.

Yakhot, V.↗

NASA NTRS1992

Development of turbulence models for shear flows by a double expansion technique

Yakhot, V.↗

DOE OSTI2026

Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity

Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1985

A note on the finite differencing of the linearized primitive equations' lower boundary condition

This note examines the accuracy of finite difference solutions of the midlatitude primitive equations and the quasi-geostrophic equation. First order accurate forward differencing of the equations' lower boundary condition is shown to poorly simulate the radiating wave response to midlatitude heating. Forward differencing always exaggerates the magnitude of the radiating response. For a realistic heating height scale and for a reasonable mesh size this exaggeration is on the order of 50 percent. Central differencing of the lower boundary condition gives an error of only about 3 percent.

Jacqmin, D.↗

NASA NTRS1994

Nonlinear hydrodynamics of cosmological sheets. 1: Numerical techniques and tests

We present the numerical techniques and tests used to construct and validate a computer code designed to study the multidimensional nonlinear hydrodynamics of large-scale sheet structures in the universe, especially the fragmentation of such structures under various instabilities. This code is composed of two codes, the hydrodynamical code ZEUS-2D and a particle-mesh code. The ZEUS-2D code solves the hydrodynamical equations in two dimensions using explicit Eulerian finite-difference techniques, with modifications made to incorporate the expansion of the universe and the gas cooling due to Compton scattering, bremsstrahlung, and hydrogen and helium cooling. The particle-mesh code solves the equation of motion for the collisionless dark matter. The code uses two-dimensional Cartesian coordinates with a nonuniform grid in one direction to provide high resolution for the sheet structures. A series of one-dimensional and two-dimensional linear perturbation tests are presented which are designed to test the hydro solver and the Poisson solver with and without the expansion of the universe. We also present a radiative shock wave test which is designed to ensure the code's capability to handle radiative cooling properly. And finally a series of one-dimensional Zel'dovich pancake tests used to test the dark matter code and the hydro solver in the nonlinear regime are discussed and compared with the results of Bond et al. (1984) and Shapiro & Struck-Marcell (1985). Overall, the code is shown to produce accurate and stable results, which provide us a powerful tool to further our studies.

Anninos, Wenbo Y.↗

NASA NTRS2001

Thermal Analysis of a Carbon Fiber Rope Barrier for Use in the Reusable Solid Rocket Motor Nozzle Joint-2

This study provides for development and verification of analysis methods used to assess performance of a carbon fiber rope (CFR) thermal barrier system that is currently being qualified for use in Reusable Solid Rocket Motor (RSRM) nozzle joint-2. Modeled geometry for flow calculations considers the joint to be vented with the porous CFR barriers placed in the "open' assembly gap. Model development is based on a 1-D volume filling approach where flow resistances (assembly gap and CFRs) are defined by serially connected internal flow and the porous media "Darcy" relationships. Combustion gas flow rates are computed using the volume filling code by assuming a lumped distribution total joint fill volume on a per linear circumferential inch basis. Gas compressibility, friction and heat transfer are included in the modeling. Gas-to-wall heat transfer is simulated by concurrent solution of the compressible flow equations and a large thermal 2-D finite element (FE) conduction grid. The derived numerical technique loosely couples the FE conduction matrix with the compressible gas flow equations, Free constants that appear in the governing equations are calibrated by parametric model comparison to hot fire subscale test results. The calibrated model is then used to make full-scale motor predictions using RSRM aft dome environments. Model results indicate that CFR thermal barrier systems will provide a thermally benign and controlled pressurization environment for the RSRM nozzle joint-2 primary seal activation.

Clayton, J. Louie↗

NASA NTRS2002

Thermal Analysis of a Carbon Fiber Rope Barrier for Use in the Reusable Solid Rocket Motor Nozzle Joint-2

This study provides development and verification of analysis methods used to assess performance of a carbon fiber rope (CFR) thermal barrier system that is currently being qualified for use in Reusable Solid Rocket Motor (RSRM) nozzle joint-2. Modeled geometry for flow calculations considers the joint to be vented with the porous CFR barriers placed in the 'open' assembly gap. Model development is based on a 1-D volume filling approach where flow resistances (assembly gap and CFRs) are defined by serially connected internal flow and the porous media 'Darcy' relationships. Combustion gas flow rates are computed using the volume filling code by assuming a lumped distribution total joint fill volume on a per linear circumferential inch basis. Gas compressibility, friction and heat transfer are included in the modeling. Gas-to-wall heat transfer is simulated by concurrent solution of the compressible flow equations and a large thermal 2-D finite element (FE) conduction grid. The derived numerical technique loosely couples the FE conduction matrix with the compressible gas flow equations. Free constants that appear in the governing equations are calibrated by parametric model comparison to hot fire subscale test results. The calibrated model is then used to make full-scale motor predictions using RSRM aft dome environments. Model results indicate that CFR thermal barrier systems will provide a thermally benign and controlled pressurization environment for the RSRM nozzle joint-2 primary seal activation.

Clayton, J. Louie↗

NASA NTRS2007

An Engineering Solution for Solving Mesh Size Effects in the Simulation of Delamination with Cohesive Zone Models

This paper presents a methodology to determine the parameters to be used in the constitutive equations of Cohesive Zone Models employed in the simulation of delamination in composite materials by means of decohesion finite elements. A closed-form expression is developed to define the stiffness of the cohesive layer. A novel procedure that allows the use of coarser meshes of decohesion elements in large-scale computations is also proposed. The procedure ensures that the energy dissipated by the fracture process is computed correctly. It is shown that coarse-meshed models defined using the approach proposed here yield the same results as the models with finer meshes normally used for the simulation of fracture processes.

FROM↗

NASA NTRS1996

Comparative Study of Advanced Turbulence Models for Turbomachinery

A computational study has been undertaken to study the performance of advanced phenomenological turbulence models coded in a modular form to describe incompressible turbulent flow behavior in two dimensional/axisymmetric and three dimensional complex geometry. The models include a variety of two equation models (single and multi-scale k-epsilon models with different near wall treatments) and second moment algebraic and full Reynolds stress closure models. These models were systematically assessed to evaluate their performance in complex flows with rotation, curvature and separation. The models are coded as self contained modules that can be interfaced with a number of flow solvers. These modules are stand alone satellite programs that come with their own formulation, finite-volume discretization scheme, solver and boundary condition implementation. They will take as input (from any generic Navier-Stokes solver) the velocity field, grid (structured H-type grid) and computational domain specification (boundary conditions), and will deliver, depending on the model used, turbulent viscosity, or the components of the Reynolds stress tensor. There are separate 2D/axisymmetric and/or 3D decks for each module considered. The modules are tested using Rocketdyn's proprietary code REACT. The code utilizes an efficient solution procedure to solve Navier-Stokes equations in a non-orthogonal body-fitted coordinate system. The differential equations are discretized over a finite-volume grid using a non-staggered variable arrangement and an efficient solution procedure based on the SIMPLE algorithm for the velocity-pressure coupling is used. The modules developed have been interfaced and tested using finite-volume, pressure-correction CFD solvers which are widely used in the CFD community. Other solvers can also be used to test these modules since they are independently structured with their own discretization scheme and solver methodology. Many of these modules have been independently tested by Professor C.P. Chen and his group at the University of Alabama at Huntsville (UAH) by interfacing them with own flow solver (MAST).

Hadid, Ali H.↗

DOE OSTI2023

Impact of Aerodynamic Modeling Assumptions on Flutter Speeds of Vertical-Axis Wind Turbines

Theodorsen’s unsteady aerodynamic theory has been used extensively in flutter speed prediction of wind turbine blades. In this study, three key assumptions of Theodorsen theory (thin airfoil, flat wake, and small angle of attack) have been revisited, and all three assumptions have been addressed in combination to obtain new lift and moment equations that are subsequently applied in flutter calculations for full-scale vertical-axis wind turbine (VAWT) rotors, not just of an individual airfoil or blade. Furthermore, edgewise aerodynamics terms are added to the lift and moment equations to include their effects on flutter speeds. The newly obtained equations were implemented in the OWENS (Offshore Wind ENergy Simulation) toolkit, which is an FEM (Finite Element Method)-based toolkit for aeroelastic analysis of VAWTs. The effect of modifying each of these assumptions has been studied for the flutter RPM prediction of three primary modes of flutter of VAWTs: propeller, butterfly, and tower modes. For the land-based case, the most change was observed for the flutter RPM of the propeller mode, with a maximum increase of 2.62% for the two-bladed UTD 5 MW VAWT case. For the floating offshore case, the primary flutter modes (tower and platform pitch) were not significantly affected.

Engineering↗

NASA NTRS1991

Numerical analysis of ignition transients in solid rocket motors

A model to analyze the unstready multidimensional turbulent flow in the head-end, star slot section of a solid rocket motor during the ignition transient is developed. The present paper examines the fluid dynamic aspects of the starting transient. A finite-difference solution of the unsteady, compressible, full Navier-Stokes equations, together with a two-equation k-epsilon turbulence model, is obtained utilizing the MacCormack explicit, predictor-corrector technique. Computed results for the flow in the slot are compared with experimental data obtained in cold flow tests of a scaled model of the Space Shuttle SRM head-end, star grain section. The agreement between the calculated flowfield and existing experimental data is very good.

Ciucci, A.↗

NASA NTRS1997

Large-Eddy Simulation of a High Reynolds Number Flow Around a Cylinder Including Aeroacoustic Predictions

The dynamic subgrid-scale model is employed in large-eddy simulations of flow over a cylinder at a Reynolds number, based on the diameter of the cylinder, of 90,000. The Centric SPECTRUM(trademark) finite element solver is used for the analysis. The far field sound pressure is calculated from Lighthill-Curle's equation using the computed fluctuating pressure at the surface of the cylinder. The sound pressure level at a location 35 diameters away from the cylinder and at an angle of 90 deg with respect to the wake's downstream axis was found to have a peak value of approximately 110 db. Slightly smaller peak values were predicted at the 60 deg and 120 deg locations. A grid refinement study suggests that the dynamic model demands mesh refinement beyond that used here.

Spyropoulos, Evangelos T.↗

NASA NTRS2008

Direct-View Multi-Point Two-Component Interferometric Rayleigh Scattering Velocimeter

This paper describes an instantaneous velocity measurement system based on the Doppler shift of elastically scattered laser light from gas molecules (Rayleigh scattering) relative to an incident laser. The system uses a pulsed laser as the light source, direct-viewing optics to collect the scattered light, an interferometer to analyze spectrally the scattered light mixed with the incident laser light, and a CCD camera to capture the resulting interferogram. The system is capable of simultaneous, spatially (approximately 0.2 mm(exp 3)) and temporally (approximately 40 ns) resolved, multiple point measurements of two orthogonal components of flow velocity in the presence of background scattered light, acoustic noise and vibrations, and flow particulates. Measurements in a large-scale axi-symmetric Mach 1.6 H2-air combustion-heated jet running at a flow sensible enthalpy specific to Mach 5.5 hypersonic flight are performed to demonstrate the technique. The measurements are compared with CFD calculations using a finite-volume discretization of the Favre-averaged Navier-Stokes equations (VULCAN code).

Bivolaru, Daniel↗

DOE OSTI2020

Matrix-free preconditioning for high-order H (curl) discretizations

The greater arithmetic intensity of high-order finite element discretizations makes them attractive for implementation on next-generation hardware, but assembly of high-order finite element operators as matrices is prohibitively expensive. As a result, the development of general algebraic solvers for such operators has been an open research challenge. Fast matrix-free application of high-order operators has received significant attention in the literature in the context of Poisson-type problems, but preconditioners and solvers for inverting more general operators are not very well-developed. In this paper, we consider the problem of preconditioning a definite Maxwell operator at high polynomial order without assembling a matrix. We show that given efficient preconditioners for high-order H 1 finite element problems on the same mesh, efficient H(curl) preconditioners can be constructed in an auxiliary space framework. We demonstrate the resulting preconditioners in a practical setting with tensor-product basis functions on an unstructured mesh of quadrilaterals. Overall, our approach uses a sparsified H 1 solver constructed on a low-order mesh of the nodal points of the underlying high-order space, and we show that the resulting H(curl) preconditioner is effective at very high polynomial orders for two-dimensional model problems with complicated geometry, varying piecewise constant coefficients, and curved elements. The resulting preconditioner scales with nearly optimal O(p d+1 ) floating point operation count and optimal O(p d ) memory transfer requirements, outperforming existing Maxwell preconditioners in the high-order regime.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1987

Experiments on the applicability of MAE techniques for predicting sound diffraction by irregular terrains

The sound field diffracted by a single smooth hill of finite impedance is studied both analytically, within the context of the theory of Matched Asymptotic Expansions (MAE), and experimentally, under laboratory scale modeling conditions. Special attention is given to the sound field on the diffracting surface and throughout the transition region between the illuminated and the shadow zones. The MAE theory yields integral equations that are amenable to numerical computations. Experimental results are obtained with a spark source producing a pulse of 42 microsec duration and about 130 Pa at 1 m. The insertion loss of the hill is inferred from measurements of the acoustic signals at two locations in the field, with subsequent Fourier analysis on an IBM PC/AT. In general, experimental results support the predictions of the MAE theory, and provide a basis for the analysis of more complicated geometries.

Berthelot, Yves H.↗

NASA NTRSDate not supplied

Understanding the Role of CFD, Including RANS, DES, LES and DNS, in Hypersonic Propulsion Design and Analysis

Computational fluid dynamics (CFD) has matured to the point where it is becoming ubiquitous in the design and analysis of aerospace vehicles and their propulsion systems. The current state-of-the-art of CFD in hypersonic propulsion design and analysis as well as a historical perspective of its development over the last 50 years is presented. The current state-of-the-art for the design of finite rate chemically reacting turbulent hypersonic flow is shown to currently be Reynolds averaged simulations (RAS) where the Reynolds averaged Navier-Stokes (RAS) equations are solved to obtain steady state solutions. An overview of the physical modeling techniques employed by RAS for the simulation of hypersonic chemically reacting flow are presented. However, scale resolving simulation techniques, such as detached eddy simulation (DES), large eddy simulations (LES), hybrid RAS/LES and direct numerical simulations (DNS), are becoming becoming powerful tools for point design analysis. One example of such a scale resolving technique, hybrid RAS/LES, will be briefly described as well as how it can be used to help address model form uncertainties in RAS.

Hypersonic↗

DOE OSTI2020

Thermal conductivity of two-dimensional disordered fibrous materials defined by interfiber thermal contact conductance and intrinsic conductivity of fibers

A general framework for the theoretical analysis and numerical calculations of the effective thermal conductivity of two-dimensional homogeneous and isotropic disordered fibrous materials is developed in this work based on the model of random contacts between soft-core spherocylinders. The analysis accounts for the interfiber contact conductance and intrinsic conductivity of fibers and is performed in a wide space of governing parameters that includes the fiber aspect ratio, Biot number calculated for a single thermal contact between fibers, and material density ranging from values corresponding to the percolation threshold up to those characteristic of dense fiber networks. For dense networks, exact theoretical equations for the thermal conductivity of materials composed of spherocylinders with an arbitrary aspect ratio and Biot number are derived. The effect of the intrinsic conductivity of fibers on the thermal transport in fibrous materials is found to depend on the density of contacts and can be significant in sufficiently dense fiber networks even if the Biot number for a single thermal contact is small. Semiempirical corrections to the theoretical equations are derived for small and moderate fiber densities. The power law exponent describing the approximate dependence of the conductivity on the density parameter is found to vary from values close to 1 up to values exceeding 2 when evaluated within different finite ranges of the density parameter. This finding explains the variability of scaling laws for thermal conductivity of fibrous materials suggested in the literature based on numerical simulations performed in different regions of the space of material parameters.

36 MATERIALS SCIENCE↗