SEARCH · Engineering Papers

Results for “Numerical homogenization”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS2023

Universal Two-Phase Convection Heat Transfer Correlations for Cryogenic Pipe Chilldown

This paper presents a set of universal two-phase convection heat transfer correlations for modeling boiling heat transfer during pipe chilldown fit over the widest available range of cryogenic fluids and thermodynamic conditions. The correlations improve upon prior correlations that were developed separately for liquid nitrogen (LN2) and liquid hydrogen (LH2) pipe quenching datasets. The new correlations include equations to calculate the single-phase vapor heat transfer, film boiling heat transfer, transition boiling heat transfer, nucleate boiling heat transfer, single-phase liquid heat transfer, bulk vapor temperature during high quality film boiling, Leidenfrost temperature, critical heat flux, critical heat flux temperature, and the onset of nucleate boiling temperature. The correlations were validated against LH2, LN2, liquid methane, liquid oxygen, and liquid argon pipe quenching datasets. The eight datasets cover the following parameter ranges: pipe lengths of 0.1 to 6.5m; outer pipe diameters of 12.7 to 25.4 mm; pipe wall thicknesses from 0.51 to 1.64 mm; flow directions of upward, downward, and horizontal; gravity levels of 1g and 0g±0.01g. A numerical model assumes homogeneous mixing between the vapor and liquid and implicitly integrates the coupled energy equations for the pipe and fluid, as well as continuity for the fluid. The model was used to estimate the equilibrium quality and fluid mass accumulation along the pipes during chilldown from estimates of the fluid-pipe heat transfer extracted from temperature measurements. The correlations can be implemented in lumped parameter codes such as SINDA/FLUINT and the Generalized Fluid System Simulation Program (GFSSP) to improve the accuracy of chilldown time and chilldown boiloff mass predictions. Such predictions are useful for designing ground or in-space cryogenic liquid transfer systems.

Universal Correlations↗

NASA NTRS1979

Numerical study of an explosion in a non-homogeneous medium with and without magnetic fields

A version of the two step Lax-Wendroff difference method with second order accuracy is used to seek solutions of the unsteady magnetohydrodynamic (MHD) equations for the study of an explosion in a non-homogeneous medium with and without magnetic fields. The explosion is initiated by introducing a finite amount of energy in a small volume of gas which leads to an instantaneous increase of gas temperature. Two types of magnetic field configurations (i.e. open and closed) are considered to illustrate the dependence and differences of magnetic and a non-magnetic flow motion. Numerical results show the strong dependence of induced MHD flow field on the magnetic field configuration and strength. The development and decay of the fast and slow MHD shock waves, as well as the ordinary gasdynamic shock waves, are also represented well by the present computer simulations. In conclusion, we have demonstrated that the numerical scheme we have used is a reliable one for studying 2-dimensional MHD problems in non-homogeneous medium.

Han, S. M.↗

NASA NTRS2011

Broken Ergodicity in MHD Turbulence in a Spherical Domain

Broken ergodicity (BE) occurs in Fourier method numerical simulations of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence. Although naive statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that low-wave-number coefficients have non-zero mean values that can be very large compared to the associated standard deviation. In other words, large-scale coherent structure (i.e., broken ergodicity) in homogeneous MHD turbulence can spontaneously grow out of random initial conditions. Eigenanalysis of the modal covariance matrices in the probability density functions of ideal statistical theory leads to a theoretical explanation of observed BE in homogeneous MHD turbulence. Since dissipation is minimal at the largest scales, BE is also relevant for resistive magnetofluids, as evidenced in numerical simulations. Here, we move beyond model magnetofluids confined by periodic boxes to examine BE in rotating magnetofluids in spherical domains using spherical harmonic expansions along with suitable boundary conditions. We present theoretical results for 3-D and 2-D spherical models and also present computational results from dynamical simulations of 2-D MHD turbulence on a rotating spherical surface. MHD turbulence on a 2-D sphere is affected by Coriolus forces, while MHD turbulence on a 2-D plane is not, so that 2-D spherical models are a useful (and simpler) intermediate stage on the path to understanding the much more complex 3-D spherical case.

Shebalin, John V.↗

NASA NTRS1993

The nonlinear Galerkin method: A multi-scale method applied to the simulation of homogeneous turbulent flows

Using results of Direct Numerical Simulation (DNS) in the case of two-dimensional homogeneous isotropic flows, the behavior of the small and large scales of Kolmogorov like flows at moderate Reynolds numbers are first analyzed in detail. Several estimates on the time variations of the small eddies and the nonlinear interaction terms were derived; those terms play the role of the Reynolds stress tensor in the case of LES. Since the time step of a numerical scheme is determined as a function of the energy-containing eddies of the flow, the variations of the small scales and of the nonlinear interaction terms over one iteration can become negligible by comparison with the accuracy of the computation. Based on this remark, a multilevel scheme which treats differently the small and the large eddies was proposed. Using mathematical developments, estimates of all the parameters involved in the algorithm, which then becomes a completely self-adaptive procedure were derived. Finally, realistic simulations of (Kolmorov like) flows over several eddy-turnover times were performed. The results are analyzed in detail and a parametric study of the nonlinear Galerkin method is performed.

Debussche, A.↗

NASA NTRS1987

Albedos of homogeneous semi-infinite canopies - Comparison of two-stream analytic and numerical solutions

The albedo of plant canopies is treated as a problem in radiative transfer. Albedos calcualted from an iterative multistream numerical model are compared with those calculated with an analytic two-stream solution. With the assumption of a randomly homogeneous distribution of leaf positions and orientations and isotropic scattering by individual leaves, the single-scattering albedo of the canopy can be found analytically. This single-scattering solution is incorporated into the two-stream solution and used to benchmark the multistream numerical model in the single-scattering limit. Relative errors so established in the multistream model are O(0.3 percent) or less. The two-stream model is also found to be remarkably accurate, with the error in multiply scattered radiation O(5 percent) or less, corresponding to absolute errors in visible albedo of less than 0.001 and near-infrared albedo of less than or equal to 0.01. Thus the two-stream model should be adequate for many purposes, such as climate modeling, provided the assumptions of homogeneous canopy and isotropic scattering are not too unrealistic.

Dickinson, Robert E.↗

NASA NTRS1987

Study of one-dimensional spectral dynamic equations of the Reynolds stresses in homogeneous anisotropic turbulence: Application to split-spectrum modeling

The CTR numerical data base generated by direct simulation of homogeneous anisotropic turbulence was used to calculate all of the terms in the spectral balance equations for the turbulent Reynolds stresses. The aim in not only to test the main closure assumptions used in the split-spectrum models, but also to try to devise improved hypotheses deduced from the statistical information. Numerical simulations of turbulent flows provide a large amount of data, a thought provoking wealth of information. The main advantage of this type of comparison is that a great variety of flows can be considered, and this is necessary to test closure hypotheses. Moreover various initial conditions can be introduced in the calculation, even if they are not experimentally feasible. All the terms in the spectral equations can be calculated. The limited Reynolds numbers of the simulations and the statistical noise caused by a small sample, particularly at the large scales, causes some difficulty in the interpretation of the results, but the method of approach proved to be a powerful tool for testing and improving spectral closures.

Schiestel, R.↗

NASA NTRS1989

Cirrus crystal nucleation by homogeneous freezing of solution droplets

A numerical model consisting of a system of differential equations is used to study cirrus crystal nucleation in a rising parcel containing a distribution of cloud condensation nuclei. The evolution of the particle population and the thermodynamic variables in the parcel are examined. The results suggest that, if homogeneous freezing is not considered, liquid water should be detected below -40 C. If homogeneous freezing is considered, the rapid growth of ice crystals and vapor depletion prevent water saturation from being reached. It is shown that the likelihood of a droplet being frozen is increased by lower temperatures, larger droplet diameter, or lower solution density.

Heymsfield, Andrew J.↗

NASA NTRS1980

Stress singularities at the vertex of a cylindrically anisotropic wedge

The plane elasticity problem for a cylindrically anisotropic solid is formulated. The form of the solution for an infinite wedge shaped domain with various homogeneous boundary conditions is derived and the nature of the stress singularity at the vertex of the wedge is studied. The characteristic equations giving the stress singularity and the angular distribution of the stresses around the vertex of the wedge are obtained for three standard homogeneous boundary conditions. The numerical examples show that the singular behavior of the stresses around the vertex of an anisotropic wedge may be significantly different from that of the isotropic material. Some of the results which may be of practical importance are that for a half plane the stress state at r = 0 may be singular and for a crack the power of stress singularity may be greater or less than 1/2.

Delale, F.↗

NASA NTRS1985

Electromagnetic electron temperature anisotropy instabilities

This paper considers electromagnetic Vlasov instabilities driven by electron temperature anisotropies in a homogeneous, nonrelativistic magnetized plasma. Numerical solutions of the full linear dispersion equation for bi-Maxwellian distribution functions and instabilities propagating parallel to the magnetic field are presented. Parametric dependences of the maximum growth rates of the electron fire hose and whistler anisotropy instabilities are given.

Gary, S. P.↗

NASA NTRS1986

Turbulent solutions of the Navier-Stokes equations

Analytical and numerical approaches to the mechanics of turbulent flow are examined in a general introduction and illustrated with graphs. Topics discussed include the averaged and unaveraged basic equations, numerical solutions and methods, homogeneous fluctuations and turbulence with no mean flow, uniformly sheared fluctuations and turbulence, inhomogeneous fluctuations and turbulence in a developing shear layer, and steady-state homogeneous turbulence with a spatially periodic body force.

Deissler, R. G.↗

NASA NTRS1987

Helicity fluctuations in incompressible turbulent flows

Results from direct numerical simulations of several homogeneous flows and fully developed turbulent channel flow indicate that the probability distribution function (pdf) of relative helicity density exhibits at most a 20 percent deviation from a flat distribution. Isotropic flows exhibit a slight helical nature but the presence of mean strain in homogeneous turbulence suppresses helical behavior. All the homogeneous turbulent flows studied show no correlation between relative helicity density and the dissipation of turbulent kinetic energy. The channel flow simulations indicate that, except for low-dissipation regions near the outer edge of the buffer layer, there is no tendency for the flow to be helical. The strong peaks in the relative helicity density pdf and the association of these peaks with regions of low dissipation found in previous simulations by Pelz et al.(1985) are not observed.

Rogers, Michael M.↗

NASA NTRS1993

Parameterization of subgrid-scale stress by the velocity gradient tensor

The objective of this work is to construct and evaluate subgrid-scale models that depend on both the strain rate and the vorticity. This will be accomplished by first assuming that the subgrid-scale stress is a function of the strain and rotation rate tensors. Extensions of the Caley-Hamilton theorem can then be used to write the assumed functional dependence explicitly in the form of a tensor polynomial involving products of the strain and rotation rates. Finally, use of this explicit expression as a subgrid-scale model will be evaluated using direct numerical simulation data for homogeneous, isotropic turbulence.

Lund, Thomas S.↗

NASA NTRS1994

The stabilizing effect of compressibility in turbulent shear flow

Direct numerical simulation of turbulent homogeneous shear flow is performed in order to clarify compressibility effects on the turbulence growth in the flow. The two Mach numbers relevant to homogeneous shear flow are the turbulent Mach number M(t) and the gradient Mach number M(g). Two series of simulations are performed where the initial values of M(g) and M(t) are increased separately. The growth rate of turbulent kinetic energy is observed to decrease in both series of simulations. This 'stabilizing' effect of compressibility on the turbulent energy growth rate is observed to be substantially larger in the DNS series where the initial value of M(g) is changed. A systematic companion of the different DNS cues shows that the compressibility effect of reduced turbulent energy growth rate is primarily due to the reduced level of turbulence production and not due to explicit dilatational effects. The reduced turbulence production is not a mean density effect since the mean density remains constant in compressible homogeneous shear flow. The stabilizing effect of compressibility on the turbulence growth is observed to increase with the gradient Mach number M(g) in the homogeneous shear flow DNS. Estimates of M(g) for the mixing and the boundary layer are obtained. These estimates show that the parameter M(g) becomes much larger in the high-speed mixing layer relative to the high-speed boundary layer even though the mean flow Mach numbers are the same in the two flows. Therefore, the inhibition of turbulent energy production and consequent 'stabilizing' effect of compressibility on the turbulence (over and above that due to the mean density variation) is expected to be larger in the mixing layer relative to the boundary layer in agreement with experimental observations.

Sarkar, S.↗

NASA NTRS2003

Large-Eddy Simulation of Conductive Flows at Low Magnetic Reynolds Number

In this paper we study the LES method with dynamic procedure in the context of conductive flows subject to an applied external magnetic field at low magnetic Reynolds number R(sub m). These kind of flows are encountered in many industrial applications. For example, in the steel industry, applied magnetic fields can be used to damp turbulence in the casting process. In nuclear fusion devices (Tokamaks), liquid-lithium flows are used as coolant blankets and interact with the surrounding magnetic field that drives and confines the fusion plasma. Also, in experimental facilities investigating the dynamo effect, the flow consists of liquid-sodium for which the Prandtl number and, as a consequence, the magnetic Reynolds number is low. Our attention is focused here on the case of homogeneous (initially isotropic) decaying turbulence. The numerical simulations performed mimic the thought experiment described in Moffatt in which an initially homogeneous isotropic conductive flow is suddenly subjected to an applied magnetic field and freely decays without any forcing. Note that this flow was first studied numerically by Schumann. It is well known that in that case, extra damping of turbulence occurs due to the Joule effect and that the flow tends to become progressively independent of the coordinate along the direction of the magnetic field. Our comparison of filtered direct numerical simulation (DNS) predictions and LES predictions show that the dynamic Smagorinsky model enables one to capture successfully the flow with LES, and that it automatically incorporates the effect of the magnetic field on the turbulence. Our paper is organized as follows. In the next section we summarize the LES approach in the case of MHD turbulence at low R(sub m) and recall the definition of the dynamic Smagorinsky model. In Sec. 3 we describe the parameters of the numerical experiments performed and the code used. Section 4 is devoted to the comparison of filtered DNS results and LES results. Conclusions are presented in Sec. 5.

Knaepen, B.↗

NASA NTRS1992

Homogeneous buoyancy-generated turbulence

Using a theoretical analysis of fundamental equations and a numerical simulation of the flow field, the statistically homogeneous motion that is generated by buoyancy forces after the creation of homogeneous random fluctuations in the density of infinite fluid at an initial instant is examined. It is shown that analytical results together with numerical results provide a comprehensive description of the 'birth, life, and death' of buoyancy-generated turbulence. Results of numerical simulations yielded the mean-square density mean-square velocity fluctuations and the associated spectra as functions of time for various initial conditions, and the time required for the mean-square density fluctuation to fall to a specified small value was estimated.

Batchelor, G. K.↗

NASA NTRS2017

Scale Dependence of Multiplier Distributions for Particle Concentration, Enstrophy, and Dissipation in the Inertial Range of Homogeneous Turbulence

Turbulent flows preferentially concentrate inertial particles depending on their stopping time or Stokes number, which can lead to significant spatial variations in the particle concentration. Cascade models are one way to describe this process in statistical terms. Here, we use a direct numerical simulation (DNS) dataset of homogeneous, isotropic turbulence to determine probability distribution functions (PDFs) for cascade multipliers, which determine the ratio by which a property is partitioned into subvolumes as an eddy is envisioned to decay into smaller eddies. We present a technique for correcting effects of small particle numbers in the statistics. We determine multiplier PDFs for particle number, flow dissipation, and enstrophy, all of which are shown to be scale dependent. However, the particle multiplier PDFs collapse when scaled with an appropriately defined local Stokes number. As anticipated from earlier works, dissipation and enstrophy multiplier PDFs reach an asymptote for sufficiently small spatial scales. From the DNS measurements, we derive a cascade model that is used it to make predictions for the radial distribution function (RDF) for arbitrarily high Reynolds numbers, Re, finding good agreement with the asymptotic, infinite Re inertial range theory of Zaichik and Alipchenkov [New J. Phys. 11, 103018 (2009)]. We discuss implications of these results for the statistical modeling of the turbulent clustering process in the inertial range for high Reynolds numbers inaccessible to numerical simulations.

turbulent flows↗

NASA NTRS1992

Homogeneous quantum electrodynamic turbulence

The electromagnetic field equations and Dirac equations for oppositely charged wave functions are numerically time-integrated using a spatial Fourier method. The numerical approach used, a spectral transform technique, is based on a continuum representation of physical space. The coupled classical field equations contain a dimensionless parameter which sets the strength of the nonlinear interaction (as the parameter increases, interaction volume decreases). For a parameter value of unity, highly nonlinear behavior in the time-evolution of an individual wave function, analogous to ideal fluid turbulence, is observed. In the truncated Fourier representation which is numerically implemented here, the quantum turbulence is homogeneous but anisotropic and manifests itself in the nonlinear evolution of equilibrium modal spatial spectra for the probability density of each particle and also for the electromagnetic energy density. The results show that nonlinearly interacting fermionic wave functions quickly approach a multi-mode, dynamic equilibrium state, and that this state can be determined by numerical means.

Shebalin, John V.↗

NASA NTRS1992

Search for subgrid scale parameterization by projection pursuit regression

The dependence of subgrid-scale stresses on variables of the resolved field is studied using direct numerical simulations of isotropic turbulence, homogeneous shear flow, and channel flow. The projection pursuit algorithm, a promising new regression tool for high-dimensional data, is used to systematically search through a large collection of resolved variables, such as components of the strain rate, vorticity, velocity gradients at neighboring grid points, etc. For the case of isotropic turbulence, the search algorithm recovers the linear dependence on the rate of strain (which is necessary to transfer energy to subgrid scales) but is unable to determine any other more complex relationship. For shear flows, however, new systematic relations beyond eddy viscosity are found. For the homogeneous shear flow, the results suggest that products of the mean rotation rate tensor with both the fluctuating strain rate and fluctuating rotation rate tensors are important quantities in parameterizing the subgrid-scale stresses. A model incorporating these terms is proposed. When evaluated with direct numerical simulation data, this model significantly increases the correlation between the modeled and exact stresses, as compared with the Smagorinsky model. In the case of channel flow, the stresses are found to correlate with products of the fluctuating strain and rotation rate tensors. The mean rates of rotation or strain do not appear to be important in this case, and the model determined for homogeneous shear flow does not perform well when tested with channel flow data. Many questions remain about the physical mechanisms underlying these findings, about possible Reynolds number dependence, and, given the low level of correlations, about their impact on modeling. Nevertheless, demonstration of the existence of causal relations between sgs stresses and large-scale characteristics of turbulent shear flows, in addition to those necessary for energy transfer, provides important insight into the relation between scales in turbulent flows.

Meneveau, C.↗