Engineering Papers⌕ Search

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.

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

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Simulation and Modeling of the Elliptic Streamline Flow

Direct numerical simulation (DNS) was used to study homogeneous turbulence with a mean flow having elliptic streamlines. This flow combines the effects of rotation and strain on the turbulence. There are many important reasons for studying the elliptic streamline flow. This flow contains the effects of both rotation and strain and is therefore similar to the mean flow in a vortex strained in the plane perpendicular to its axis. Such flows provide insight into the fundamental vortical interactions within turbulence. A strained vortex also occurs in airplane wakes, in which each wingtip vortex induces a strain field on the other. The strain field can affect the stability of these vortices and thereby their turbulent structure downstream. The ability to understand and predict the turbulent structure of the vortices is important to the wake hazard problem which is of major concern for the safety of commercial aircraft.

Blaisdell, Gregory A.↗
NASA NTRS1994

Radiative Heat Transfer in Finite Cylindrical Enclosures with Nonhomogeneous Participating Media

Results of a numerical solution for radiative heat transfer in homogeneous and nonhomogeneous participating media are presented. The geometry of interest is a finite axisymmetric cylindrical enclosure. The integral formulation for radiative transport is solved by the YIX method. A three-dimensional solution scheme is applied to two-dimensional axisymmetric geometry to simplify kernel calculations and to avoid difficulties associated with treating boundary conditions. As part of the effort to improve modeling capabilities for turbulent jet diffusion flames, predicted distributions for flame temperature and soot volume fraction are used to calculate radiative heat transfer from soot particles in such flames. It is shown that the nonhomogeneity of radiative property has very significant effects. The peak value of the divergence of radiative heat flux could be underestimated by 2 factor of 7 if a mean homogeneous radiative property is used. Since recent studies have shown that scattering by soot agglomerates is significant in flames, the effect of magnitude of scattering is also investigated and found to be nonnegligible.

Hsu, Pei-Feng↗
NASA NTRS2010

Comprehensive T-Matrix Reference Database: A 2007-2009 Update

The T-matrix method is among the most versatile, efficient, and widely used theoretical techniques for the numerically exact computation of electromagnetic scattering by homogeneous and composite particles, clusters of particles, discrete random media, and particles in the vicinity of an interface separating two half-spaces with different refractive indices. This paper presents an update to the comprehensive database of T-matrix publications compiled by us previously and includes the publications that appeared since 2007. It also lists several earlier publications not included in the original database.

Mishchenko, Michael I.↗
NASA NTRS1990

An experimental study of the effects of rapid rotation on turbulence

Experiments, large eddy simulations and direct numerical simulations all show that rapid rotation has a striking effect on homogeneous nearly isotropic turbulence. The cascade process is effectively inhibited by rotation, and thus dissipation is greatly reduced. Some attempts have been made to incorporate this effect in turbulence models. Numerical simulation further showed the somewhat surprising result that anisotropic homogeneous turbulence subjected to rotation tended towards an isotropic state; however, the residual anisotropy was not zero. Reynolds (1989) performed a detailed analysis using Rapid Distortion Theory (RDT) and showed that a reduction in the anisotropy is indeed expected and if the anisotropy is produced by irrotational strain then the anisotropy tensor b(sub ij) is asymptotically driven to half its initial value. Our objective is to extend the work of the experiments mentioned above to lower turbulent Rossby numbers and to confirm some of the results obtained by Reynolds (1989) for anisotropic turbulence.

Veeravalli, S. V.↗
NASA NTRS2019

Orbits Through Polytropes

We describe how orbital tunnels could be used to transport payloads through the Earth. If you use a brachistochrone for the tunnel, the body forces in the tunnel become overwhelmingly large for small angular distances traveled. Projectiles move along an orbital tunnel faster than they would along abrachistochrone connecting the same points but the body force components cancel. We describe how parabolic Keplerian orbits outside the object merge onto quasi-Keplerian orbits inside the object. We use models of the interior of the Earth with three values of the polytropic index (n) to calculate interior or bits that travel between surface points. The n¼3 results are also scaled to the Sun. Numerical integrations of the equations describing polytropes were used to generate the initial models. Numerical integration of the equations of motion are then used to calculate the angular distance you can travel along the surface and the traversal time as a function of the parabolic periaps is distance for each model. Trajectories through objects of low central condensation show a focusing effect that decreases as the central condensation increases. Analytic solutions for the trajectories in a homogeneous sphere are derived and compared to the numeric results. The results can be scaled to other planets, stars, or even globular clusters.

Gjerlov, Amalia↗
DOE OSTI2025

Hydrodynamic Coupling to a Homogenized Radiation Transport Method based on Young Measures

Resolving radiation transport fields subject to opacity profiles with strong, oscillatory line structure while potentially falling under intermediate optical depth conditions presents a numerical challenge in radiation transport modeling. The Young measure-based homogenization technique formulated by Haut et al. (2017) was investigated as a candidate method for resolving radiation fields under such conditions more accurately. The method was compared against frequently-utilized mean opacity methods as the Rosseland and Planck formulations. In this work, all methods were tested through radiation slab calculations separately comprised of aluminum, copper, and krypton, each for different thermodynamic conditions. Following these offline radiation slab calculations, demonstrations shifted towards the SCEPTRE radiation transport code and, subsequently, the multiphysics ALEGRA code for approximately-coupled radiation-material simulations. Throughout all the simulations shown in this study, for a fixed computational cost, the homogenized method was observed to be more accurate than any of the solutions determined through traditional mean opacity approaches.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Geometry and dynamics of passive scalar structures in compressible turbulent mixing

We present a structure-based numerical analysis of passive scalar mixing in decaying homogeneous isotropic turbulence (DHIT) and shock-turbulence interaction canonical configurations. The analysis focuses on the temporal evolution of ensembles of passive scalar structures, initialized as spheres of different sizes relative to the Taylor microscale. An algorithm is introduced to track the evolution of each individual structure and the interactions with other structures in the ensemble, relating changes in the surface geometry and the underlying physical processes (turbulent transport, scalar dissipation, and shock compression). The tracking algorithm is applied to datasets from shock-capturing direct numerical simulations of DHIT, with Taylor microscale Reynolds number Re λ =40 and turbulence Mach number M t =0.2 ⁠, and STI cases in which the turbulence is processed by a shock wave at Mach numbers M = 1.5 and 3.0. Temporal surface convolution increases for initially larger structures, resulting in a higher probability of locally hyperbolic geometries where breakup into smaller structures occurs. Shock-induced deformation of the structures amplifies breakup processes, enhancing mixing, particularly for larger structures. Mixing enhancement by the shock is manifested as an amplification of the surface-averaged scalar gradient, which increases for initially larger structures. The alignment between the scalar gradient and the most extensional strain-rate eigendirection on the scalar isosurfaces also increases across the shock. Larger magnitudes of the scalar gradient and its alignment with the most compressive strain-rate eigendirection correlate with flatter surface regions. Shock-induced structure compression increases the area coverage of flat regions, where the amplification of scalar gradient is localized.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Compressibility Effects on the Growth and Structure of Homogeneous Turbulent Shear Flow

Compressibility effects within decaying isotropic turbulence and homogeneous turbulent shear flow have been studied using direct numerical simulation. The objective of this work is to increase our understanding of compressible turbulence and to aid the development of turbulence models for compressible flows. The numerical simulations of compressible isotropic turbulence show that compressibility effects are highly dependent on the initial conditions. The shear flow simulations, on the other hand, show that measures of compressibility evolve to become independent of their initial values and are parameterized by the root mean square Mach number. The growth rate of the turbulence in compressible homogeneous shear flow is reduced compared to that in the incompressible case. The reduced growth rate is the result of an increase in the dissipation rate and energy transfer to internal energy by the pressure-dilatation correlation. Examination of the structure of compressible homogeneous shear flow reveals the presence of eddy shocklets, which are important for the increased dissipation rate of compressible turbulence.

Blaisdell, G. A.↗
NASA NTRSDate not supplied

Numerical Analysis of Cavitation-Induced Two-Phase Fluid Transients in Cryogenic Feedlines Using Finite Volume Procedure

This article presents a numerical investigation of the cavitation induced fluid transient in cryogenic feedline due to sudden closure of the valve A one dimensional numerical model based on finite volume procedure consisting of nodes and branches using Generalized Fluid System Simulation Program ( has been employed for the prediction of pressure variations The unsteady friction formulation accurately determines wall shear stress, and a homogeneous mixture model is used for two phase flow modeling Numerical results obtained with the present formulation are compared with experimental results and Method of Characteristic ( based results It is observed that low residual in flow variable calculation improves the solution providing numerically stable results GFSSP is able to reproduce the fluid transient peak due to vapor collapse Subsequently, it can retain the shape of pressure wave for two phase fluid transients, but it is dampening faster compared to experiments Further investigation is required with set of experiments to validate the present numerical approach for different experimental data.

Fluid Transient↗
NASA NTRS1989

Solitary waves in mantle plumes

Numerical calculations in two dimensions show that solitary wave disturbances can propagate along thermal plumes in a homogeneous, viscous fluid with a thermally activated rheology similar to the mantle. Comparison of the numerical results with analytic results from an idealized two-fluid model indicates that the two-fluid model can be applied to solitary waves in the thermally activated plume. A two-fluid model of solitary waves on cylindrical conduits is accordingly applied to mantle plumes to estimate the propagation speeds, time durations, and pulse lengths of solitary waves in mantle plumes as a function of background mantle viscosity, plume flux, plume density deficit, plume viscosity, and the volume of material transported by the solitary wave. It is inferred that mantle plumes may be eposodic via solitary waves and that these disturbances might account for observed 0(10 m.y.) variations in the durations and spacings of episodes of enhanced hotspot volcanic activity. Solitary waves in mantle plumes could be generated by interactions among plumes and between plumes and the large-scale, time-dependent mantle circulation.

Schubert, G.↗
DOE OSTI2021

Experimental validation of shock propagation through a foam with engineered macro-pores

The engineered macro-pore foam provides a new way to study thermonuclear burn physics by utilizing capsules containing deuterated (D) foam and filling tritium (T) gas in the engineered macro-pores. The implosion of a thermonuclear capsule filled with an engineered macro-pore foam will be complex due to the interaction of a shock wave with the engineered macro-pores. It is our goal to quantify how substantially complex foam structures affect the shape of shock and bulk shock speed. A cylinder-shape shock tube experiment has been designed and performed at the Omega Laser Facility. In order to examine how a foam structure will affect shock propagation, we performed several tests varying (1) engineered macro-pore size, (2) average foam density, and (3) with/without neopentane (C 5 H 12 ) gas. X-ray radiographic data indicate that shock speed through engineered macro-pore foams depends strongly on average foam density and less on pore size. In this work, experimental shock propagation data helped guide two numerical simulation approaches: (1) a 2D simulation with homogenizing foams rather than explicitly simulating engineered macro-pores and (2) a 2D toroidal-pore approximation adopting a toroidal-tube geometry to model engineered macro-pores.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Electrostatically driven pattern formation in mixed charged–neutral multicomponent elastic membranes

Multicomponent crystalline and amorphous elastic shells exhibit heterogeneous surface patterns that provide distinctive functionalities in cellular environments. Such patterning typically arises from the competition between short-range attractive and long-range repulsive interactions in membranes. Here, we demonstrate that the intrinsic competition between electrostatic repulsion and elastic deformation is sufficient to drive spontaneous surface patterning in elastic shells, requiring no additional attractive interactions. Using numerical simulations, we demonstrate pattern formation in mechanically homogeneous membranes with heterogeneous surface charge composition across different topologies, including spheres, discs, and flat periodic membranes. We also examine patterns in crystalline and amorphous shells of coassembled charged and neutral components with different bending rigidities. At low charge fraction, discrete charged surface domains form. At intermediate charge fraction, the competition between electrostatics and elasticity leads to elongated domains (rods) of the charged component, which results in lamellar patterns at nearly equal fraction of the charged and neutral components. At high charge fraction, nanodomains of the neutral component form. Amorphous shells exhibit similar progressions but with disordered structures rather than ordered lamellar patterns. These pattern morphologies are observed in both the closed shells and flat membranes. As salt concentration increases, all patterns coarsen due to the screening of electrostatic interactions.

charge heterogeneity↗
DOE OSTI2026

Radial oscillations of viscous neutron stars: Zero diffusion case

In this article, the spectrum of radial oscillations of neutron stars is systematically studied within two frameworks of viscous relativistic hydrodynamics: the relativistic Navier-Stokes and Israel-Stewart theories. A correspondence is established between the discrete stellar eigenmodes and the continuous dispersion relation of perturbations around a homogeneous fluid, providing a basis for interpreting our numerical results. We analyze the Newtonian limit and assess the impact of relativistic corrections, such as the gravitational redshifting of microscopic relaxation timescales. We show that bulk viscosity can significantly affect the behavior of both hydrodynamic and nonhydrodynamic fundamental modes, and that, depending on the magnitude of the viscous effects, it is the nonhydrodynamic mode that becomes unstable beyond the turning point in a sequence of equilibrium configurations. These results provide a useful step toward systematic studies of neutron star quasinormal modes in the presence of viscosity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Shear Assisted Processing and Extrusion (ShAPE) of Aluminum Alloy 7075, 2024, and Al-12.4TM

The most common aluminum alloys utilized in the aerospace industry are 7075 and 2024 due to their high strength-to-weight ratio compared to advanced high strength steels and other aluminum alloys. Despite excellent performance, these aluminum alloys have seen limited use outside of the aerospace industry due in part to high cost. If high-performance aluminum extrusions could be made more cost effectively by eliminating energy intensive process steps typical of conventional extrusion, then numerous opportunities exist for more widespread adoption. A key reason for the high cost of 7075 and 2024 extrusions (25-75% higher than 6061) is their slow extrusion speed. 7075 and 2024 are limited to 2 m/min and 3.5 m/min respectively, in contrast to 6061 which can be extruded at 20–80 m/min. In addition to slow speed, aluminum alloys require numerous thermal treatments throughout the extrusion process including homogenization and pre-heating prior to extrusion, and solution heat treating and artificial aging after extrusion. Each of these steps contribute to the total energy consumed during manufacturing of extruded components. This project investigates the use of ShAPE to improve extrusion speed and reduce, or even eliminate, the typical thermal treatments for high strength aluminum alloys, all while improving material performance. The overarching goal of this project was to demonstrate that Shear Assisted Processing and Extrusion (ShAPE) can manufacture high-performance aluminum alloy tubing with lower manufacturing energy and improved mechanical properties compared to conventional extrusion. Unlike conventional extrusion where the billet is rammed against a stationary die using a strictly linear motion, the ShAPE process superimposes a rotational shear force by spinning the die while the billet is plunged. Compared to conventional linear extrusion, the ShAPE process imparts significantly more strain into the feedstock material, which enables the formation of novel microstructures. These microstructures manifest an array of property and process improvements for extrusion of high-performance aluminum alloys. The following accomplishments were achieved for this project: Extrusion of 7075 at 12.2 meters/min compared to 2 meters/min for conventional extrusion; Elimination of 7075 billet homogenization (430 °C for 20 hours) which is required prior to conventional extrusion; Elimination of 7075 billet pre-heating (400 °C for 1 hour) in a separate furnace which is required prior to conventional extrusion; Achieved 7075-T6with yield strength = 595 MPa, ultimate tensile strength = 531MPa, and elongation = 17.4% for extrusions made from unhomogenized billets. Exceeds the ASTM and ASM standard, and typical industry values; Achieved 7075-T5 (i.e., no solution heat treatment) with yield strength = 588 MPa, ultimate tensile strength = 535 MPa, and elongation = 14.8% for extrusions made from homogenized billets; Extrusion of 2024 at 7.4 meters/min compared to 3.5 meters/min for conventional extrusion; Achieved 2024-T8510 yield strength = 522 MPa, ultimate strength = 510MPa, and elongation = 7.1% for extrusions made from wrought billets. Exceeds the ASTM and ASM standard, and typical industry values; Extrusion of Al-12.4TM high-performance aluminum powder directly into tubing, in a single step, which eliminates process steps typical of powder metallurgy extrusion.

36 MATERIALS SCIENCE↗
NASA NTRS1973

On the relationship between acoustic energy density flux near the jet axis and far field acoustic intensity

By measurement and analysis, the relationship between the distribution of the outflow of acoustic energy over the jet boundary and the far-field intensity is considered. The physical quantity used is the gradient of the pressure evaluated on a geometrical plane at the smallest possible radial distance from the jet axis, but outside the vortical region, in the area where the homogeneous wave equation is reasonably well satisfied. The numerical and experimental procedures involved have been checked out by using a known source. Results indicate that the acoustic power output per unit length of the jet, in the region from which the sound emanates, peaks at approximately 9 diameters downstream. The acoustic emission for a jet Strouhal number of about 0.3 exceeds the emission for all other Strouhal numbers nearly everywhere along the measurement plane. However, the far-field peak intensity distribution obtained from the contribution of each station was found to depend on the spatial extent of the region where sound emanates from the jet, which, in turn, depends more on the far-field angle than on the Strouhal number.

Maestrello, L.↗
NASA NTRS1972

Diffraction of a plane pulse by a three-dimensional corner.

The conical solutions for the incidence of a plane pulse on a three-dimensional corner are presented. The corner is represented by a trihedron with one edge perpendicular to the other two. Both the boundary condition of the first kind and that of the second kind are considered. Outside the characteristic sphere of the vertex of the corner, the solution is represented by the well-known conical solutions in two variables. Inside the characteristic sphere, the problem involves three conical variables. By the separation of variables, the problem is reduced to that of an eigenvalue problem with an irregular boundary which is in turn reduced to a system of homogeneous algebraic equations. The eigenvalues are then determined numerically. By the superposition of the conical solutions for plane pulses, the solution for the incidence of a plane wave is obtained. Numerical examples simulating the incidence of a sonic boom on the corner of a structure are presented.

Ting, L.↗
NASA NTRS1990

Reduced and simplified chemical kinetics for air dissociation using Computational Singular Perturbation

The Computational Singular Perturbation CSP methods is employed (1) in the modeling of a homogeneous isothermal reacting system and (2) in the numerical simulation of the chemical reactions in a hypersonic flowfield. Reduced and simplified mechanisms are constructed. The solutions obtained on the basis of these approximate mechanisms are shown to be in very good agreement with the exact solution based on the full mechanism. Physically meaningful approximations are derived. It is demonstrated that the deduction of these approximations from CSP is independent of the complexity of the problem and requires no intuition or experience in chemical kinetics.

Goussis, D. A.↗