Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct numerical simulation”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Streaming Instability with Multiple Dust Species – I. Favourable Conditions for the Linear Growth

A recent study suggests that the streaming instability, one of the leading mechanisms for driving the formation of planetesimals, may not be as efficient as previously thought. Under some disc conditions, the growth time-scale of the instability can be longer than the disc lifetime when multiple dust species are considered. To further explore this finding, we use both linear analysis and direct numerical simulations with gas fluid and dust particles tomutually validate and study the unstablemodes of the instability in more detail.We extend the previously studied parameter space by one order of magnitude in both the range of the dust-size distribution [Ts, min, Ts, max] and the total solid-to-gas mass ratio ε and introduce a third dimension with the slope q of the size distribution. We find that the fast-growth regime and the slow-growth regime are distinctly separated in the ε–Ts, max space, while this boundary is not appreciably sensitive to q or Ts, min.With a wide range of dust sizes present in the disc (e.g. Ts, min  10−3), the growth rate in the slow-growth regime decreases as more dust species are considered.With a narrow range of dust sizes (e.g. Ts, max/Ts, min = 5), on the other hand, the growth rate in most of the ε–Ts, max space is converged with increasing dust species, but the fast and the slow growth regimes remain clearly separated. Moreover, it is not necessary that the largest dust species dominate the growth of the unstable modes, and the smaller dust species can affect the growth rate in a complicated way. In any case, we find that the fast-growth regime is bounded by ε  1 or Ts, max  1, which may represent the favourable conditions for planetesimal formation.

Zhaohuan Zhu↗
NASA NTRS2022

Interaction of a Tunnel-like Acoustic Disturbance Field with a Blunt Cone Boundary Layer at Mach 8

The existing measurements of laminar-to-turbulent transition over circular cones in conventional (i.e., “noisy”) hypersonic wind tunnels have established that the transition location moves downstream when the nose radius is increased from zero. However, this initially downstream movement slows down and ultimately reverses beyond a critical value of the nose radius, and may be related to external forcing in the form of freestream disturbances and/or surface roughness. To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, hypersonic boundary-layer receptivity to broadband freestream acoustic disturbances from the nozzle wall of a digital conventional wind tunnel is investigated by both direct numerical simulations (DNS) and modal and nonmodal stability analysis. A Mach 8 flow over a 7 deg half-angle cone with a nose radius of 𝑅 𝑛 = 5.2mm and freestream Reynolds number of 12.2 × 10 6 m -1 is considered. The results show that the broadband tunnel noise in the free stream of a convectional hypersonic wind tunnel (i.e., outside of the nozzle-wall turbulent boundary layer) can be well represented by an acoustic model with an ansatz of slow acoustic waves. With successful calibration of the model parameters against the precursor tunnel DNS, such an acoustic ansatz can successfully reproduce both the frequency-wavenumber spectra and the temporal evolution of the broadband tunnel noise radiated from the nozzle wall. Additionally, the DNS of the Mach 8 blunt cone with tunnel-like acoustic input above the bow shock showed that the spectra of wall-pressure and heat-transfer fluctuations recovers the signature of the axisymmetric waves predicted by the nonmodal analysis. Furthermore, the azimuthal wavenumber and frequency spectrum of the temperature fluctuations as a function of the wall-normal distance show higher amplitudes for three-dimensional waves above the boundary-layer edge. The numerical schlieren contours show the inclined structures commonly observed in blunt cone experiments, demonstrating that they correspond to three-dimensional structures due to freestream disturbances in the presence of an entropy layer.

Boundary layer transition↗
NASA NTRS2022

Numerical Study of the Turbulent Prandtl Number in Supersonic Plane-Channel Flow – the Effect of Thermal Boundary Conditions

Direct numerical simulation is used to investigate the turbulent Prandtl number Prt above cold (isothermal) and hot (adiabatic) walls in a family of supersonic channel flows. A range of Reynolds numbers and mean temperature/density variations is considered. The value of Prt away from the wall approaches 0.85 above both the isothermal and adiabatic walls, with the width of the Prt = 0.85 region increasing with wall-variable-based Reynolds number Reτw (the latter strongly affected by the thermal boundary condition). The variable, near-wall Prt profiles from both the present and previous, passive-scalar simulations collapse as a function of the semilocal y∗ wall scaling proposed by Huang et al., with the y∗-dependence weakly parameterized by Reτw.

Turbulence↗
NASA NTRSDate not supplied

Computational Study of Low-Speed Boundary-Layer Laminar-Turbulent Transition over a 30º Swept Plate in a Channel

Reliable boundary-layer laminar-turbulent transition predictions are critical for the design of aircraft that use natural or controlled laminar flow to achieve low drag and better fuel efficiency. During cruise flight through low turbulence, the three-dimensional boundary layer over swept surfaces can transition prematurely due to stationary crossflow instabilities induced by surface roughness. Parabolized stability equations analyses have been successfully used to explain the different phases of the transition process under well controlled perturbations such as steps or discrete roughness elements. However, methodologies are still lacking to model the boundary layer receptivity to distributed random roughness, and the following non-linear instability stages leading to transition. For engineering applications, several semi-empirical methods based on linear amplification factors or local helicity have found growing adoption, but their reliability is limited by their level of simplification and the sparsity of data they are based on. To address this gap, a new study of the effect of distributed random surface roughness on crossflow transition was started, using a 30◦ swept plate as a test vehicle. The plate, which is mounted in the 2 Foot × 3 Foot Low-Speed Boundary-Layer Channel at the NASA Langley Research Center, has been tested with micro-scale roughness of varied geometry near the leading edge. In order to analyze the flow in the test section and prepare Direct Numerical Simulations with resolved roughness, Reynolds-Averaged Navier-Stokes (RANS) simulations of the flow over the plate installed in the wind tunnel test section were performed using the Langtry-Menter transition model with stationary crossflow extension. This paper presents the results of these simulations, with comparison to the experiment.

ARMD↗
NASA NTRSDate not supplied

DeepONet-Assisted Optimization of Surface Topography for Transition Delay in A Mach 4.5 Boundary Layer

We use deep learning, an ensemble variationaltechnique (EnVar), and direct numerical simulations(DNS) to design an optimal topography for a two-dimensional roughness element that delays the on-set of laminar-turbulent transition in a Mach 4.5 flat-plate boundary layer. Deep operator networks (Deep-ONets), which have the known ability to learn com-plex nonlinear operators within dynamical systems,are used for machine learning. For the baseline config-uration of a smooth flat plate, the second-mode wavesat the DNS inflow cause a quick nonlinear breakdownof the high-speed boundary layer within the computa-tional domain. Results reported in the present studyvalidate the ability of DeepONets to model the tran-sition delay via a given topography of the roughnesselement. The computing cost to optimize the rough-ness element for minimal skin-friction drag is substan-tially lowered by the DeepONets-based reduced-ordermodel. In comparison to the baseline method of EnVaroptimization based on DNS alone, the DeepONets-based EnVar optimizer is able to delay transition pastthe outflow boundary of the computational domainwhile utilizing almost 5–6 times fewer DNS.

Machine Learning↗
NASA NTRSDate not supplied

DeepONet-Assisted Optimization of Surface Topography for Transition Delay in a Mach 4.5 Boundary Layer

We use deep learning, an ensemble variational technique (EnVar), and direct numerical simulations(DNS) to design an optimal topography for a two-dimensional roughness element that delays the on-set of laminar-turbulent transition in a Mach 4.5 flat-plate boundary layer. Deep operator networks (DeepONets), which have the known ability to learn complex nonlinear operators within dynamical systems, are used for machine learning. For the baseline configuration of a smooth flat plate, the second-mode waves at the DNS inflow cause a quick nonlinear breakdown of the high-speed boundary layer within the computational domain. Results reported in the present study validate the ability of DeepONets to model the transition delay via a given topography of the roughness element. The computing cost to optimize the rough-ness element for minimal skin-friction drag is substantially lowered by the DeepONets-based reduced-order model. In comparison to the baseline method of EnVar optimization based on DNS alone, the DeepONets-based EnVar optimizer is able to delay transition past the outflow boundary of the computational domain while utilizing almost 5–6 times fewer DNS.

Machine Learning↗
NASA NTRSDate not supplied

Hypersonic Boundary-Layer Transition over a Blunt Circular Cone in a Mach 8 Digital Wind Tunnel

To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, three-dimensional direct numerical simulations (DNS) were conducted to investigate the interaction of a tunnel-like acoustic disturbance field with a Mach 8 laminar boundary layer over a 7 deg half-angle blunt cone with a nose radius of 𝑅 𝑛 = 5.2 mm. The flow conditions and the axisymmetric cone geometry matched those measured in the Sandia Hypersonic Wind Tunnel at Mach 8 (Sandia HWT-8), and the incident freestream acoustic disturbances in the DNS were generated by acoustic radiation from the nozzle-wall turbulent boundary layer in the same tunnel. Consistent with the predictions of the nonmodal instability analysis, the study shows the emergence of a second spectral peak for moderately oblique waves between 150 kHz and 200 kHz in the aft part of the cone. The spectral energy of temperature and streamwise velocity fluctuations is initially concentrated within the entropy layer. These fluctiations are eventually “swallowed” by the boundary layer as the entropy-layer thickness 𝛿 𝑆 becomes similar to the boundary-layer thickness 𝛿h𝑡 at downstream locations. The numerical schlieren contours show the same inclined structures commonly observed in blunt cone experiments, and these structures are found to be three dimensional and advect with a constant speed approximately equal to the mean flow speed near the edge of the boundary layer.

Hypersonic Flow↗
NASA NTRSDate not supplied

Hypersonic Boundary-Layer Transition over a Blunt Circular Cone in a Mach 8 Digital Wind Tunnel

To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, three-dimensional direct numerical simulations (DNS) were conducted to investigate the interaction of a tunnel-like acoustic disturbance field with a Mach 8 laminar boundary layer over a 7 deg half-angle blunt cone with a nose radius of 𝑅 𝑛 = 5.2 mm. The flow conditions and the axisymmetric cone geometry matched those measured in the Sandia Hypersonic Wind Tunnel at Mach 8 (Sandia HWT-8), and the incident freestream acoustic disturbances in the DNS were generated by acoustic radiation from the nozzle-wall turbulent boundary layer in the same tunnel. Consistent with the predictions of the nonmodal instability analysis, the study shows the emergence of a second spectral peak for moderately oblique waves between 150 kHz and 200 kHz in the aft part of the cone. The spectral energy of temperature and streamwise velocity fluctuations is initially concentrated within the entropy layer. These fluctiations are eventually “swallowed” by the boundary layer as the entropy-layer thickness 𝛿 𝑆 becomes similar to the boundary-layer thickness 𝛿h𝑡 at downstream locations. The numerical schlieren contours show the same inclined structures commonly observed in blunt cone experiments, and these structures are found to be three dimensional and advect with a constant speed approximately equal to the mean flow speed near the edge of the boundary layer.

Hypersonic Flow↗
NASA NTRSDate not supplied

New Capabilities and Improvements to the High-Order Glenn Flux Reconstruction Code

The Glenn Flux Reconstruction (GFR) code is a computational fluid dynamics (CFD) code under development at NASA Glenn Research Center. GFR is based on the high-order flux reconstruction (FR) method and provides a large-eddy simulation (LES) capability that is both accurate and efficient for complex aeropropulsion flows. Three significant new capabilities have been added to the code that improve its performance and functionality. First, a variety of explicit Runge-Kutta methods, including some with adaptive time stepping, were added to GFR with two methods offering a 33% improvement in time-to-solution. Second, GFR can now utilize fully unstructured, mixed-element meshes to more easily facilitate the grid generation process for complex geometries. Finally, a rotating reference frame capability has been added to GFR for solving rotating turbomachinery problems. A selection of results demonstrating these new capabilities are presented in this work. The Taylor-Green vortex problem is used to verify the new unstructured capability by showing similar accuracy and resolution for all element types. LES of the Turbulent Heat Flux Phase III (THX3) experiment with comparison to another high-order LES code and a popular Reynolds-averaged Navier-Stokes (RANS) code demonstrates the accuracy of the code for complex aeropropulsion flows. Finally, LES of a spacecraft cabin ventilation fan shows the ability of GFR to efficiently establish a fan performance map and identify operating points for further analysis at high orders of accuracy.

High-Order Methods↗
NASA NTRS1990

A Numerical Method for Direct Simulation of Turbulence in Complex Geometries

The ultimate goal of this work is to study the flow inside a channel with riblets on one of the two walls. The method has been tested for two dimensional flows in the presence of bodies with a geometrical singularity and for three dimensional flows inside domains described by Cartesian coordinates. The results have been compared with previous numerical simulations and with experimental results. The cases considered are: (1) the growth of Orr-Sommerfield waves in plane Poiseulle flow; (2) the flow over a backward facing step; (3) the flow past a wedge; and (4) the flow inside a narrow channel. Finally, the case of a channel with two large riblets on a wall has been simulated. In this case a limited number of grid points is sufficient, and in spite of the slow convergence for the pressure solver, one is able to obtain solutions with a reasonable amount of computer time. At present solutions with very fine grids in all three directions can not be obtained, due to lack of a fast pressure solver for general curvilinear coordinates.

P Orlandi↗
NASA NTRS1986

Spacecraft charging in the auroral plasma: Progress toward understanding the physical effects involved

The work is presented in four parts. First, main differences between the plasma environments in geostationary orbit and low polar orbit with regard to high-voltage charging situations are reviewed. Next, results are presented from a calculation of secondary-electron escape currents from negatively-charged spacecraft surfaces having various orientations relative to the local magnetic-field direction. It is shown that for finite ranges of combinations of electric and magnetic field directions, secondary-electron escape is completely suppressed and therefore cannot help to discharge the spacecraft. In such circumstances, secondary electrons may travel distances many times their gyroradii before reimpacting, and this may produce greatly increased secondary-electron surface currents. Thirdly, a simple rough estimate of the required conditions for high-voltage auroral-zone charging is developed. The results suggest that for any given spacecraft, surface potentials are likely to depend more strongly on the ratio of ambient flux of high-energy electrons to that of all ions, than on any other environmental parameter. Finally, preliminary results are presented of numerical simulation work directed toward testing this hypothesis. Numerical instabilities encountered in doing this simulation work probably are closely related to physical sensitivities inherent in the physics of the ion wake behind the spacecraft, and especially to beam-like constituents of the ion population in the wake.

Laframboise, J. G.↗
NASA NTRS1999

Direct-Numerical and Large-Eddy Simulations of a Non-Equilibrium Turbulent Kolmogorov Flow

A non-equilibrium form of turbulent Kolmogorov flow is set up by making an instantaneous change in the amplitude of the spatially-periodic forcing. It is found that the response of the flow to this instantaneous change becomes more dramatic as the wavenumber of the forcing is increased, and, at the same time, that the faithfulness with which the large-eddy-simulation results agree with the direct-numerical results decreases.

Woodruff, S. L.↗
NASA NTRS1993

One-equation near-wall turbulence modeling with the aid of direct simulation data

The length scales appearing in the relations for the eddy viscosity and dissipation rate in one-equation models were evaluated from direct numerical (DNS) simulation data for developed channel and boundary-layer flow at two Reynolds numbers each. To prepare the ground for the evaluation, the distribution of the most relevant mean-flow and turbulence quantities is presented and discussed, also with respect to Reynolds-number influence and to differences between channel and boundary-layer flow. An alternative model is tested as near wall component of a two-layer model by application to developed-channel, boundary-layer and backward-facing-step flows.

Rodi, W.↗
NASA NTRS1981

Numerical simulation of a compressible homogeneous, turbulent shear flow

A direct, low Reynolds number, numerical simulation was performed on a homogeneous turbulent shear flow. The full compressible Navier-Stokes equations were used in a simulation on the ILLIAC IV computer with a 64,000 mesh. The flow fields generated by the code are used as an experimental data base, to examine the behavior of the Reynols stresses in this simple, compressible flow. The variation of the structure of the stresses and their dynamic equations as the character of the flow changed is emphasized. The structure of the tress tensor is more heavily dependent on the shear number and less on the fluctuating Mach number. The pressure-strain correlation tensor in the dynamic uations is directly calculated in this simulation. These correlations are decomposed into several parts, as contrasted with the traditional incompressible decomposition into two parts. The performance of existing models for the conventional terms is examined, and a model is proposed for the 'mean fluctuating' part.

Feiereisen, W. J.↗
NASA NTRS2011

Geomagnetic Secular Variation Prediction with Thermal Heterogeneous Boundary Conditions

It has long been conjectured that thermal heterogeneity at the core-mantle boundary (CMB) affects the geodynamo substantially. The observed two pairs of steady and strong magnetic flux lobes near the Polar Regions and the low secular variation in the Pacific over the past 400 years (and perhaps longer) are likely the consequences of this CMB thermal heterogeneity. There are several studies on the impact of the thermal heterogeneity with numerical geodynamo simulations. However, direct correlation between the numerical results and the observations is found very difficult, except qualitative comparisons of certain features in the radial component of the magnetic field at the CMB. This makes it difficult to assess accurately the impact of thermal heterogeneity on the geodynamo and the geomagnetic secular variation. We revisit this problem with our MoSST_DAS system in which geomagnetic data are assimilated with our geodynamo model to predict geomagnetic secular variations. In this study, we implement a heterogeneous heat flux across the CMB that is chosen based on the seismic tomography of the lowermost mantle. The amplitude of the heat flux (relative to the mean heat flux across the CMB) varies in the simulation. With these assimilation studies, we will examine the influences of the heterogeneity on the forecast accuracies, e.g. the accuracies as functions of the heterogeneity amplitude. With these, we could be able to assess the model errors to the true core state, and thus the thermal heterogeneity in geodynamo modeling.

Kuang, Weijia↗
NASA NTRS2017

Macrosegregation Resulting from Directional Solidification Through an Abrupt Change in Cross-Sections

Simulations of the directional solidification of two hypoeutectic alloys (Al-7Si alloy and Al-19Cu) and resulting macrosegregation patterns are presented. The casting geometries include abrupt changes in cross-section from a larger width of 9.5 mm to a narrower 3.2 mm width then through an expansion back to a width of 9.5 mm. The alloys were chosen as model alloys because they have similar solidification shrinkages, but the effect of Cu on changing the density of the liquid alloy is about an order of magnitude greater than that of Si. The simulations compare well with experimental castings that were directionally solidified in a graphite mold in a Bridgman furnace. In addition to the simulations of the directional solidification in graphite molds, some simulations were effected for solidification in an alumina mold. This study showed that the mold must be included in numerical simulations of directional solidification because of its effect on the temperature field and solidification. For the model alloys used for the study, the simulations clearly show the interaction of the convection field with the solidifying alloys to produce a macrosegregation pattern known as "steepling" in sections with a uniform width. Details of the complex convection- and segregation-patterns at both the contraction and expansion of the cross-sectional area are revealed by the computer simulations. The convection and solidification through the expansions suggest a possible mechanism for the formation of stray grains. The computer simulations and the experimental castings have been part of on-going ground-based research with the goal of providing necessary background for eventual experiments aboard the ISS. For casting practitioners, the results of the simulations demonstrate that computer simulations should be applied to reveal interactions between alloy solidification properties, solidification conditions, and mold geometries on macrosegregation. The simulations also presents the possibility of engineering the mold-material to avoid, or mitigate, the effects of thermosolutal convection and macrosegregation by selecting a mold material with suitable thermal properties, especially its thermal conductivity.

Lauer, M.↗
NASA NTRS1994

Direct simulation Monte Carlo and Navier-Stokes simulations of blunt body wake flows

Numerical results obtained with direct simulation Monte Carlo and Navier-Stokes methods are presented for a Mach-20 nitrogen flow about a 70-deg blunted cone. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how the wake structure changes as a function of rare faction, what the afterbody levels of heating are, and to what limits the continuum models are realistic as rarefunction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the afterbody sting are emphasized in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

WAKE STRUCTURE↗
NASA NTRS1994

Direct simulation Monte Carlo and Navier-Stokes simulations of blunt body wake flows

Numerical results obtained with direct simulation Monte Carlo and Navier-Stokes methods are presented for a Mach-20 nitrogen flow about a 70-deg blunted cone. The flow conditions simuulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how the wake structure changes as a function of rarefaction, what the afterbody levels of heating are, and to what limits the continuum models are realistic as rarefaction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the after body sting are emphasizes in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

Moss, James N.↗