SEARCH · Engineering Papers

Results for “viscous effects”

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 1,081 records · Page 60

All sources NASA NTRS DOE OSTI

NASA NTRS2005

Ice Forming Regions during Evolution of the Solar Nebula

The condensation/sublimation front using a two dimensional model of the evolving solar system is investigated based on combined viscous and radiative heating. The snow line is shown to be a two-branched curve reflecting the competing effects of solar heating in the photosphere and internal heating at the center plane. The evolution of the icy region is described from a limited region early in the disk evolution to final positions near 1 AU. The snow line evolution predicted using two surface density models, a Hayashi minimum mass power law and an analytical solution of the nebula evolution equation. Possible effects of this dynamic motion on disk chemistry and organic molecule formation is also described.

Davis, Sanford S.↗

NASA NTRS2017

Cause and Cure - Deterioration in Accuracy of CFD Simulations with Use of High-Aspect-Ratio Triangular/Tetrahedral Grids

Traditionally high-aspect ratio triangular/tetrahedral meshes are avoided by CFD researchers in the vicinity of a solid wall, as it is known to reduce the accuracy of gradient computations in those regions. Although for certain complex geometries, the use of high-aspect ratio triangular/tetrahedral elements in the vicinity of a solid wall can be replaced by quadrilateral/prismatic elements, ability to use triangular/tetrahedral elements in such regions without any degradation in accuracy can be beneficial from a mesh generation point of view. The benefits also carry over to numerical frameworks such as the space-time conservation element and solution element (CESE), where simplex elements are the mandatory building blocks. With the requirement of the CESE method in mind, a rigorous mathematical framework that clearly identifies the reason behind the difficulties in use of such high-aspect ratio simplex elements is formulated using two different approaches and presented here. Drawing insights from the analysis, a potential solution to avoid that pitfall is also provided as part of this work. Furthermore, through the use of numerical simulations of practical viscous problems involving high-Reynolds number flows, how the gradient evaluation procedures of the CESE framework can be effectively used to produce accurate and stable results on such high-aspect ratio simplex meshes is also showcased.

CESE Method↗

NASA NTRS1998

Deterministic Stress Modeling of Hot Gas Segregation in a Turbine

Simulation of unsteady viscous turbomachinery flowfields is presently impractical as a design tool due to the long run times required. Designers rely predominantly on steady-state simulations, but these simulations do not account for some of the important unsteady flow physics. Unsteady flow effects can be modeled as source terms in the steady flow equations. These source terms, referred to as Lumped Deterministic Stresses (LDS), can be used to drive steady flow solution procedures to reproduce the time-average of an unsteady flow solution. The goal of this work is to investigate the feasibility of using inviscid lumped deterministic stresses to model unsteady combustion hot streak migration effects on the turbine blade tip and outer air seal heat loads using a steady computational approach. The LDS model is obtained from an unsteady inviscid calculation. The LDS model is then used with a steady viscous computation to simulate the time-averaged viscous solution. Both two-dimensional and three-dimensional applications are examined. The inviscid LDS model produces good results for the two-dimensional case and requires less than 10% of the CPU time of the unsteady viscous run. For the three-dimensional case, the LDS model does a good job of reproducing the time-averaged viscous temperature migration and separation as well as heat load on the outer air seal at a CPU cost that is 25% of that of an unsteady viscous computation.

Busby, Judy↗

NASA NTRS1987

Computation of asymmetric supersonic flows around cones at large incidence

The Schiff-Steger parabolized Navier-Stokes (PNS) code has been modified to allow computation of conical flowfields around cones at high incidence. The improved algorithm of Degani and Schiff has been incorporated with the PNS code. This algorithm adds the cross derivative and circumferential viscous terms to the original PNS code and modifies the algebraic eddy viscosity turbulence model to take into account regions of so called cross-flow separation. Assuming the flowfield is conical (but not necessarily symmetric) a marching stepback procedure is used: the solution is marched one step downstream using improved PNS code and the flow variables are then scaled to place the solution back to the original station. The process is repeated until no change in the flow variables is observed with further marching. The flow variables are then constant along rays of the flowfield. The experiments obtained by Bannik and Nebbeling were chosen as a test case. In these experiments a cone of 7.5 deg. half angle at Mach number 2.94 and Reynolds number 1.372 x 10(7) was tested up 34 deg. angle of attack. At high angle of attack nonconical asymmetric leeward side vortex patterns were observed. In the first set of computations, using an earlier obtained solution of the above cone for angle of attack of 22.6 deg. and at station x=0.5 as a starting solution, the angle of attack was gradually increased up to 34 deg. During this procedure the grid was carfully adjusted to capture the bow shock. A stable, converged symmetric solution was obtained. Since the numerical code converged to a symmetric solution which is not the physical one, the stability was tested by a random perturbation at each point. The possible effect of surface roughness or non perfect body shape was also investigated. It was concluded that although the assumption of conical viscous flows can be very useful for certain cases, it can not be used for the present case. Thus the second part of the investigation attempted to obtain a marching (in space) solution with the PNS method using the conical solution as initial data. Finally, the solution of the full Navier-Stokes equations was carried out.

Degani, David↗

DOE OSTI2025

Bias-Variance Trade-Off in Physics-Informed Neural Networks with Randomized Smoothing for High-Dimensional PDEs

Physics-Informed Neural Networks (PINNs) have triggered a paradigm shift in scientific computing, leveraging mesh-free properties and robust approximation capabilities. While proving effective for low-dimensional partial differential equations (PDEs), the computational cost of PINNs remains a hurdle in high-dimensional scenarios. This is particularly pronounced when computing high-order and high-dimensional derivatives in the physics-informed loss. Randomized Smoothing PINN (RS-PINN) introduces Gaussian noise for stochastic smoothing of the original neural net model, enabling the use of Monte Carlo methods for derivative approximation, which eliminates the need for costly automatic differentiation. Despite its computational efficiency, especially in the approximation of high-dimensional derivatives, RS-PINN introduces biases in both loss and gradients, negatively impacting convergence, especially when coupled with stochastic gradient descent (SGD) algorithms. We present a comprehensive analysis of biases in RS-PINN, attributing them to the nonlinearity of the Mean Squared Error (MSE) loss as well as the intrinsic nonlinearity of the PDE itself. We propose tailored bias correction techniques, delineating their application based on the order of PDE nonlinearity. The derivation of an unbiased RS-PINN allows for a detailed examination of its advantages and disadvantages compared to the biased version. Specifically, the biased version has a lower variance and runs faster than the unbiased version, but it is less accurate due to the bias. To optimize the bias-variance trade-off, we combine the two approaches in a hybrid method that balances the rapid convergence of the biased version with the high accuracy of the unbiased version. In addition to methodological contributions, we present an enhanced implementation of RS-PINN. Extensive experiments on diverse high-dimensional PDEs, including Fokker-Planck, Hamilton-Jacobi-Bellman (HJB), viscous Burgers’, Allen-Cahn, and Sine-Gordon equations, illustrate the bias-variance trade-off and highlight the effectiveness of the hybrid RS-PINN. Empirical guidelines are provided for selecting biased, unbiased, or hybrid versions, depending on the dimensionality and nonlinearity of the specific PDE problem.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1975

Finite-rate chemistry effects upon convective and radiative heating of an atmospheric entry vehicle

A mathematical model of the aerothermochemical environment along the stagnation line of a planetary return spacecraft using an ablative thermal protection system was developed and solved for conditions typical of atmospheric entry from planetary missions. The model, implemented as a FORTRAN 4 computer program, was designed to predict viscous, reactive and radiative coupled shock layer structure and the resulting body heating rates. The analysis includes flow field coupling with the ablator surface, binary diffusion, coupled line and continuum radiative and equilibrium or finite rate chemistry effects. The gas model used includes thermodynamic, transport, kinetic and radiative properties of air and ablation product species, including 19 chemical species and 16 chemical reactions. Specifically, the impact of nonequilibrium chemistry effects upon stagnation line shock layer structure and body heating rates was investigated.

Guillermo, P.↗

DOE OSTI2025

Polarization of the 𝜙 meson in the hadronic phase with nucleon scatterings and a viscous hydrodynamic background

We extend our previous work on the spin alignment of the ϕ vector meson in the hadronic phase of the quark gluon plasma, by including effects of nucleon scatterings. The emission rates are calculated in a realistic hydrodynamic background simulated with the code FLUIDuM , for different beam energies. We find that all the effects taken into account cannot explain the out-of-plane spin alignment of the ϕ meson observed experimentally.

Polarization phenomena↗

NASA NTRS1993

The effect of unsteadiness on the time-mean thermal loads in a turbine stage

Two steady numerical analysis methods and one unsteady method are used to study the viscous three-dimensional flow in the middle stage of the Pratt & Whitney alternate design Space Shuttle Main Engine fuel turbine. The principal characteristic of this flow is that the secondary flows generated in the rotor blade reconfigure a radial inlet total temperature distortion into one with a pitchwise exit hot streak distortion. Secondary flows in the following vane redistribute the radial variation while unsteadiness causes a segregation of hot and cold flow from the hot streak within the vane. Such redistribution and segregation can lead to unexpected thermal loads and reduced durability. The physical phenomena and the ability of a steady analysis to capture them are investigated by performing a numerical experiment whereby the results of the two steady analysis methods are compared to the time-mean of the unsteady simulation. The flow physics related to the segregation and mixing of total temperature are discussed.

Kirtley, K. R.↗

NASA NTRS1991

The effect of viscosity on steady transonic flow with a nodal solution topology

The effect of viscosity on a steady, transonic flow for which the inviscid limit has a nodal solution topology near the critical point is investigated. For the accelerating case, viscous solutions tend to repel each other, so that a very delicate choice of initial conditions is required to prevent them from diverging. Only the two critical solutions extend to arbitrarily large distances into both the subsonic and supersonic flows. For the decelerating case, the solutions tend to attract, and so an entire two-parameter family of solutions now extends over large distances. The general effect of viscosity on the solution degeneracy of a nodal topology is thus to reduce or limit it for the accelerating case and to enhance it for the decelerating case. The astrophysical implications of these findings are addressed.

Owocki, Stanley P.↗

NASA NTRS1981

Calculation of Unsteady Transonic Flow over an Airfoil

An implicit finite-difference solver for either the Euler equations or the "thin-layer" Navier-Stokes equations was used to calculate a transonic flow over the NACA 64A010 airfoil pitching about its one-quarter chord. An unsteady automatic grid-generation procedure that will improve significantly the computational efficiency of various unsteady flow problems is described. The calculated results for both inviscid and viscous flows at Much number 0.8 over the airfoil oscillating with reduced frequency referenced to one-half chord, 0.2, are compared with experimental data measured in the Ames 11 x 11-ft Transonic Wind Tunnel. Nonlinear, unsteady effects of the flow on the surface pressure variations, shock-wave excursions, and overall airloads are examined. Good agreements between the results of computations and experiments were obtained. In the shock-wave region, however, the results of the viscous-flow computations showed closer agreement with the experimental data.

Chyu, W. J.↗

NASA NTRS1981

Calculation of Unsteady Transonic Flow Over an Airfoil

An implicit finite-difference solver for either the Euler equations or the "thin-layer" Navier-Stokes equations was used to calculate a transonic flow over the NACA 64A010 airfoil pitching about its one-quarter chord. An unsteady automatic grid-generation procedure that will improve significantly the computational efficiency of various unsteady flow problems is described. The calculated results for both inviscid and viscous flows at Mach number 0.8 over the airfoil oscillating with reduced frequency referenced to one-half chord, 0.2, are compared with experimental data measured in the Ames 11 x 11 ft Transonic Wind Tunnel. Nonlinear, unsteady effects of the flow on the surface pressure variations, shock-wave excursions, and overall airloads are examined. Good agreements between the results of computations and experiments were obtained. In the shock-wave region, however, the results of the viscous-flow computations showed closer agreement with the experimental data.

Chyu, W. J.↗

DOE OSTI2021

Jet wake from linearized hydrodynamics

We explore how to improve the hybrid model description of the particles originating from the wake that a jet produced in a heavy ion collision leaves in the droplet of quark-gluon plasma (QGP) through which it propagates, using linearized hydrodynamics on a background Bjorken flow. Jet energy and momentum loss described by the hybrid model become currents sourcing linearized hydrodynamics. By solving the linearized hydrodynamic equations numerically, we investigate the development of the wake in the dynamically evolving droplet of QGP, study the effect of viscosity, scrutinize energy-momentum conservation, and check the validity of the linear approximation. We find that linearized hydrodynamics works better in the viscous case because diffusive modes damp the energy-momentum perturbation produced by the jet. We calculate the distribution of particles produced from the jet wake by using the Cooper-Frye prescription and find that both the transverse momentum spectrum and the distribution of particles in azimuthal angle are similar in shape in linearized hydrodynamics and in the hybrid model. Their normalizations are different because the momentum-rapidity distribution in the linearized hydrodynamics analysis is more spread out, due to sound modes. Since the Bjorken flow has no transverse expansion, we explore the effect of transverse flow by using local boosts to add it into the Cooper-Frye formula. After including the effects of transverse flow in this way, the transverse momentum spectrum becomes harder: more particles with transverse momenta bigger than 2 GeV are produced than in the hybrid model. Although we defer implementing this analysis in a jet Monte Carlo, as would be needed to make quantitative comparisons to data, we gain a qualitative sense of how the jet wake may modify jet observables by computing proxies for two example observables: the lost energy recovered in a cone of varying open angle, and the fragmentation function. We find that linearized hydrodynamics with transverse flow effects added improves the description of the jet wake in the hybrid model in just the way that comparison to data indicates is needed. Our study illuminates a path to improving the description of the wake in the hybrid model, highlighting the need to take into account the effects of both transverse flow and the broadening of the energy-momentum perturbation in spacetime rapidity on particle production.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Bayesian calibration and uncertainty quantification of a rate-dependent cohesive zone model for polymer interfaces

In this work we present a rate-dependent cohesive zone model for the fracture of polymeric interfaces and performs a Bayesian calibration, an uncertainty quantification, and a sensitivity analysis for the model. The proposed cohesive zone model accounts for both reversible elastic and irreversible rate-dependent separation sliding deformation at the interface. The viscous dissipation due to the irreversible opening at the interface is modeled using elastic-viscoplastic kinematics that incorporates the effects of strain rate. Inverse calibration of parameters for such complex models through trial and error is challenging due to the large number of parameters of the model. Moreover, the calibrated parameter values are often non-unique and uncertain when the available experimental data is limited. To tackle this challenge, we employ a Bayesian calibration approach to identify parameters from experimental data, the resulting parameters significantly enhance the accuracy of the model. To quantify the uncertainty associated with the inverse parameter estimation, a modular Bayesian approach is employed to calibrate the unknown model parameters, accounting for the parameter uncertainty of the cohesive zone model. The advantages of the Bayesian calibration over a deterministic parameter fit are demonstrated. Further, to quantify the model uncertainties, such as incorrect assumptions or missing physics, a discrepancy function is introduced, which significantly improves the model’s prediction. Finally, the total uncertainty of the model is quantified in a predictive setting. A sensitivity analysis is performed to assess how changes in the input variables of the model affect the peak load, facilitating the identification of a concise set of highly influential parameters. The present approach can be used for calibration and uncertainty quantification for other complex computational mechanics models. It should also facilitate the designing of interface materials under uncertainty.

42 ENGINEERING↗

DOE OSTI2023

Particle size and shape effect of Crumbler® rotary shear-milled granular woody biomass on the performance of Acrison® screw feeder: A computational and experimental investigation

Physical experiments and discrete element model (DEM) simulations are conducted to evaluate particle characteristics and operation parameter effects on screw feeding performance for rotary shear-milled Douglas fir. Three performance metrics are used: mass flow rate, shaft driving torque, and specific energy consumption. The impact of particle size, particle size distribution (PSD), shaft rotational speed (rpm), and hopper dimensions on the performance are investigated. All employed performance metrics reveal the superior flowability of the 2-mm particles in contrast to the larger 6-mm counterpart. Remarkably, wider PSD results in poorer flowability than the two mono-sized particles, proving the flowability enhancement achieved by narrower PSD of the woody feedstock. More importantly, DEM simulations unveil PSD-induced degradation in flowability is attributed to mechanical interlocking and particle segregation effects. Furthermore, higher shaft rpm causes higher mass flow rate at the cost of higher specific energy consumption due to viscous dissipation and changes in flow pattern.

09 BIOMASS FUELS↗

NASA NTRS1972

Factors which influence the analysis and design of ejector nozzles

A theoretical analysis of the viscous interaction between the primary and secondary streams of ejector nozzles was developed. The analysis accounts for real sonic line effects and the streamwise variation in stream mixing and boundary layer within the ejector. The aspects of the analysis are explained and illustrated by applying the theory to a variety of ejector configurations including cylindrical shroud, contoured flap and plug nozzles. Extensive comparisons are made between theory and data to show the importance of various analytical assumptions and such design variables as diameter ratio, spacing ratio, total temperature ratio, and primary nozzle geometry.

Anderson, B. H.↗

NASA NTRS1973

The analysis of the transient pressure response of the shuttle EPS-ECS cryogenic tanks with external pressurization systems

An analysis of transient pressures in externally pressurized cryogenic hydrogen and oxygen tanks was conducted and the effects of design variables on pressure response determined. The analysis was conducted with a computer program which solves the compressible viscous flow equations in two-dimensional regions representing the tank and external loop. The external loop volume, thermal mass, and heat leak were the dominant design variables affecting the system pressure response. No significant temperature stratification occurred in the fluid contained in the tank.

Barton, J. E.↗

NASA NTRS1974

Supersonic laminar base pressure, heat-transfer, and upstream influence correlations for small steps

Some results of the application of a unified theory of small disturbances in nonuniform high-speed boundary layers to the small-step problem for the case of supersonic laminar nonadiabatic flow are described. Comparisons with experimental data are given which show good agreement in a variety of physical features, including useful new base pressure and heat-transfer correlation relations. Using a small perturbation approach, an analytical theory is developed that predicts the steady-state pressure, temperature, and heat-transfer disturbance fields including suction or blowing through the step. Particular emphasis is placed on the important effects of the highly nonuniform flow across the boundary layer, lateral pressure gradients, and the upstream influence associated with viscous-inviscid interaction. Approximate analytical expressions for the wall pressure, skin friction, and heat-transfer perturbations and upstream influence distance are derived by a Fourier transformation approach.

Inger, G. R.↗

NASA NTRS1978

A numerical solution of the Navier-Stokes equations for chemically nonequilibrium, merged stagnation shock layers on spheres and two-dimensional cylinders in air

Results of solving the Navier-Stokes equations for chemically nonequilibrium, merged stagnation shock layers on spheres and two-dimensional cylinders are presented. The effects of wall catalysis and slip are also examined. The thin shock layer assumption is not made, and the thick viscous shock is allowed to develop within the computational domain. The results show good comparison with existing data. Due to the more pronounced merging of shock layer and boundary layer for the sphere, the heating rates for spheres become higher than those for cylinders as the altitude is increased.

Johnston, K. D.↗