SEARCH · Search NASA
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.
RIPPLE computer program models transient, two-dimensional flows of incompressible fluids with surface tension on free surfaces of general shape. Surface tension modeled as volume force derived from continuum-surface-force model, giving RIPPLE both robustness and accuracy in modeling surface-tension effects at free surface. Also models wall adhesion effects. Written in FORTRAN 77.
This paper highlights the accomplishments on a joint effort between NASA - Marshall Space Flight Center and Texas A and M University to develop accurate seal analysis software for use in rocket turbopump design, design audits and trouble shooting. Results for arbitrary clearance profile, transient simulation, thermal effects solution and flexible seal wall model are presented. A new solution for eccentric seals based on cubic spline interpolation and ordinary differential equation integration is also presented.
Observations have been made and reported that the experimental normal force coefficients at a constant angle of attack were constant with a variation of more than 2 orders of magnitude of Reynolds number at a free-stream Mach number M(sub infinity) of 8.00 and more than 1 order of magnitude variation at M(sub infinity) = 6.00 on the same body-wing hypersonic cruise configuration. These data were recorded under laminar, transitional, and turbulent boundary layer conditions with both hot-wall and cold-wall models. This report presents experimental data on 25 configurations of 17 models of both simple and complex geometry taken at M(sub infinity) = 6.00, 6.86, and 8.00 in 4 different hypersonic facilities. Aerodynamic calculations were made by computational fluid dynamics (CID) and engineering methods to analyze these data. The conclusions were that the normal force coefficients at a given altitude are constant with Reynolds numbers at hypersonic speeds and that the axial force coefficients recorded under laminar boundary-layer conditions at several Reynolds numbers may be plotted against the laminar parameter (the reciprocal of the Reynolds number to the one-half power) and extrapolated to the ordinate axis to determine the inviscid-wave-drag coefficient at the intercept.
Structural skeleton of high Reynolds number flow gives insight into smooth and rough walls. (modeling implications not yet clear) Dynamic roughness as a model of "designer" roughness: a) Impose length scale(s) and dominant frequency. b) Irregular roughness well-represented by first "few" POD modes (Christensen, 2009). Experiments and simple model demonstrate: a) Harmonics associated with forcing (and w) important. b) Change to the mean profile (skin friction).
The 2013-2022 Decaedal survey for planetary exploration has identified probe missions to Uranus and Saturn as high priorities. This work endeavors to examine the uncertainty for determining aeroheating in such entry environments. Representative entry trajectories are constructed using the TRAJ software. Flowfields at selected points on the trajectories are then computed using the Data Parallel Line Relaxation (DPLR) Computational Fluid Dynamics Code. A Monte Carlo study is performed on the DPLR input parameters to determine the uncertainty in the predicted aeroheating, and correlation coefficients are examined to identify which input parameters show the most influence on the uncertainty. A review of the present best practices for input parameters (e.g. transport coefficient and vibrational relaxation time) is also conducted. It is found that the 2(sigma) - uncertainty for heating on Uranus entry is no more than 2.1%, assuming an equilibrium catalytic wall, with the uncertainty being determined primarily by diffusion and H(sub 2) recombination rate within the boundary layer. However, if the wall is assumed to be partially or non-catalytic, this uncertainty may increase to as large as 18%. The catalytic wall model can contribute over 3x change in heat flux and a 20% variation in film coefficient. Therefore, coupled material response/fluid dynamic models are recommended for this problem. It was also found that much of this variability is artificially suppressed when a constant Schmidt number approach is implemented. Because the boundary layer is reacting, it is necessary to employ self-consistent effective binary diffusion to obtain a correct thermal transport solution. For Saturn entries, the 2(sigma) - uncertainty for convective heating was less than 3.7%. The major uncertainty driver was dependent on shock temperature/velocity, changing from boundary layer thermal conductivity to diffusivity and then to shock layer ionization rate as velocity increases. While radiative heating for Uranus entry was negligible, the nominal solution for Saturn computed up to 20% radiative heating at the highest velocity examined. The radiative heating followed a non-normal distribution, with up to a 3x variation in magnitude. This uncertainty is driven by the H(sub 2) dissociation rate, as H(sub 2) that persists in the hot non-equilibrium zone contributes significantly to radiation.
The design of a modular multi-physics high-order space-time finite-element framework is presented together with its extension to allow monolithic coupling of different physics. One of the main objectives of the framework is to perform efficient high- fidelity simulations of capsule/parachute systems. This problem requires simulating multiple physics including, but not limited to, the compressible Navier-Stokes equations, the dynamics of a moving body with mesh deformations and adaptation, the linear shell equations, non-re effective boundary conditions and wall modeling. The solver is based on high-order space-time - finite element methods. Continuous, discontinuous and C1-discontinuous Galerkin methods are implemented, allowing one to discretize various physical models. Tangent and adjoint sensitivity analysis are also targeted in order to conduct gradient-based optimization, error estimation, mesh adaptation, and flow control, adding another layer of complexity to the framework. The decisions made to tackle these challenges are presented. The discussion focuses first on the "single-physics" solver and later on its extension to the monolithic coupling of different physics. The implementation of different physics modules, relevant to the capsule/parachute system, are also presented. Finally, examples of coupled computations are presented, paving the way to the simulation of the full capsule/parachute system.
eddy is a collection of tools - nonlinear solvers, meshing, post-processing, visualization, optimization, etc. - for performing scale-resolving simulations of multi-physics applications. The framework is designed to enable advanced R&D on a variety of topics by leveraging a mature capability for scale resolving simulations, and simultaneously be an appropriate tool for application analysis and support. Currently, eddy is at a relatively low technical readiness level (TRL), and users and developers should maintain appropriate expectations. The technical details behind eddy are outlined in several publications which can be consulted for more information [1–10]. The solvers are built around an unstructured high-order capability, and heavily utilize the tensor product sum-factorization approach for efficiency. The unsteady formulation utilizes a fully implicit space-time approach with a matrix-free Newton- Krylov method. A primitive steady-state solver is available for testing purposes, but is not expected to converge for all but simple verification cases. The Navier-Stokes fluid solvers do not support either RANS or hybrid-RANS capability, only LES and wall-modeled LES approaches. All of the solvers within eddy support three modes of operation: a primal solve of the full nonlinear problem, and two linearization approaches of the primal solve - the ad joint and the tangent solution. Details on how to select and use these three modes are outlined in Sec. 3.
Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES),and wall-modeled large-eddy simulation (WMLES)are employed to simulate transonic flow over an NACA-0012 airfoil at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75,and the Reynolds number based on the chord length is 10million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985).1The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model.The unsteady RANS simulations,only with the compressibility correction,predict the pre- and post-buffet characteristics,which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment.The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment
Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.
Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.
This paper presents results from ongoing research on jet noise prediction using wall-modeled large eddy simulations (WMLES) conducted within the Launch Ascent and Vehicle Aerodynamics (LAVA) computational framework. The study primarily focuses on the aeroacoustic implications of multi-stream nozzle configurations with internal mixing and both internal and external plugs, at a Reynolds number of 1×10 6 based on the nozzle exit diameter. While internal mixing nozzles have long been considered for jet noise reduction, the complexities of their impact on overall noise levels remain insufficiently understood. This research applies established best practices for WMLES to these complex nozzle configurations, aiming to assess their efficacy and identify limitations in accurately predicting noise behaviors. Through detailed comparisons with experimental data obtained from NASA’s Glenn Research Center, initial findings underscore the challenges inherent in current simulation practices when confronted with intricate geometric and operational conditions. Responding to these challenges, the study explores innovative computational approaches to rectify discrepancies noted between experimental outcomes and CFD predictions. Once example is the introduction of a resonating sound source in the simulation environment to mimic potential unobserved acoustic phenomena.A second example is the strategic modifications to the geometry of external plugs to account for real-world deformations caused by heating and gravity. These novel strategies aim to enhance the accuracy and reliability of noise predictions from supersonic jets, advancing our understanding and capability to effectively reduce jet noise in commercial supersonic aircraft
The behavior of various terms in the Reynolds-stress transport equations is evaluated in view results of near-wall turbulence characterizations which show that all components of the velocity-pressure gradient correlation vanish at the wall. A methodology for deriving an asymptotically correct model of the velocity-pressure gradient correlation is proposed, and a model capable of approaching the high-Reynolds number model for pressure distribution far from the wall is derived. A parametric study of the model constants introduced by the near-wall closure reveals that one constant in the dissipation-rate equation is Reynolds number-dependent; with this modification, excellent agreement is obtained with measured and simulated near-wall turbulence statistics.
Numerical simulations and wind tunnel experiments were performed to investigate flow separation and its control over a wall-mounted hump model. Three-dimensional, unsteady fluid dynamic simulations were used to capture the relevant flow physics of unsteady flow separation and its control via active and passive flow control methods. Surface pressure measurements and surface oilflow visualization were obtained in the experiments. Current measurements documented a considerably thinner incoming boundary layer, but similar pressure distribution, compared to the reference data reported in the CFD (Computational Fluid Dynamics) Validation (CFDVAL 2004) Experiments. Consistent with the data available in the literature, numerical simulations predicted a longer separation bubble, which manifested itself as a shift in the pressure distributions. Streamwise vortices generated by passive vortex generators increased the suction pressure and provided substantial pressure recovery. Numerical simulations of vortex generators qualitatively predicted the effects of control, and revealed patterns of attached/separated flow regions. On the other hand, active flow control using steady discrete jets reduced, but was not able to eliminate, flow separation for the tested configuration. Although the major disagreement between the numerical and experimental results was in the pressure recovery region due to the inaccurate prediction of separation bubbles, the relative pressure recoveries from the baselines were comparable. Surface oilflow visualization and simulated surface streamlines were qualitatively in good agreement revealing the key flow features.
A k-epsilon model is proposed for wall bonded turbulent flows. In this model, the eddy viscosity is characterized by a turbulent velocity scale and a turbulent time scale. The time scale is bounded from below by the Kolmogorov time scale. The dissipation equation is reformulated using this time scale and no singularity exists at the wall. The damping function used in the eddy viscosity is chosen to be a function of R(sub y) = (k(sup 1/2)y)/v instead of y(+). Hence, the model could be used for flows with separation. The model constants used are the same as in the high Reynolds number standard k-epsilon model. Thus, the proposed model will be also suitable for flows far from the wall. Turbulent channel flows at different Reynolds numbers and turbulent boundary layer flows with and without pressure gradient are calculated. Results show that the model predictions are in good agreement with direct numerical simulation and experimental data.
In this paper we present a new thin-wall eddy current modeling code, ThinCurr, for studying inductively-coupled currents in 3D conducting structures -- with primary application focused on the interaction between currents flowing in coils, plasma, and conducting structures of magnetically-confined plasma devices. The code utilizes a boundary finite element method on an unstructured, triangular grid to accurately capture device structures. The new code, part of the broader Open FUSION Toolkit, is open-source and designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. Scalability to large models is enabled through use of hierarchical off-diagonal low-rank compression of the inductance matrix, which is otherwise dense. Ease of handling large models of complicated geometry is further supported by automatic determination of supplemental elements through a greedy homology approach. Here, a detailed description of the numerical methods of the code and verification of the implementation of those methods using cross-code comparisons against the VALEN code and Ansys commercial analysis software is shown.
Plant secondary cell walls constitute the dominant reservoir of renewable biomass, comprising tightly packed cellulose, hemicellulose, and lignin at the nanoscale. Recent advances in solid-state NMR spectroscopy and the availability of small-angle X-ray scattering for biomass characterization have led to an accumulation of experimental data on cell wall organization, yet no explicit structure model has simultaneously satisfied both Xray and NMR observations. Using wheat straw as a model system, we propose a structural framework consistent with current knowledge of cellulose biosynthesis, X-ray scattering data, and one- and two-dimensional 13 C solid-state NMR spectra. In this model, 18-chain elementary fibrils align in parallel and populate the cross-section at random. Arabinose-substituted xylan shows no conformational dependence for cellulose-binding in wheat, and only a minor fraction of 2-fold xylan appears in close proximity to cellulose, unlike in Arabidopsis, where xylan is more tightly attached to the cellulose surface. While NMR data cannot unambiguously resolve the internal arrangement of the 18 glucan chains, X-ray scattering profiles uniquely constrain the fibril size and exclude the possibility of tight bundling in the intact walls. The specific interaction between the matrix polymers and the cellulose elementary fibrils must be reconsidered in light of the small interfibril spaces, which bring the matrix components into spatial proximity with cellulose even in the absence of attractive interactions. These findings provide fundamental molecular-level insight into cellulose fibril architecture and matrix−polymer interactions, resolving longstanding discrepancies between spectroscopic and scattering data and advancing our understanding of biopolymer assembly into structurally and functionally versatile lignocellulosic biomaterials.
A wall-injected, hydrogen-fueled scramjet is modeled numerically using a parabolic Navier-Stokes computer code with a k-epsilon turbulence model and finite-rate chemistry. The numerical results are compared to experimental scramjet data taken in a shock tunnel and are found to be in reasonable agreement. Under the conditions studied, the numerical results show that combustion in the scramjet is mainly mixing limited.
This is the final report for the DOE award DE-SC0021657, titled “Validation of Simulation Models for Wall-Connected Scrape off Layer Currents during MHD Activity in the HBT-EP Tokamak” for the period September 1, 2020 – April 30, 2025. This award is part of a broader collaborative project, which supported training of a PhD student at Columbia University under the joint supervision of Dr. Levesque and Dr. Hansen (on award DE-SC0021325). The award supported research into computational models for plasma and 3D conducting structures in the vicinity relevant to tokamak disruptions and other mode activity (eg. RWMs) using the NIMROD and PSI-Tet codes. Additional detail on progress in support of these tasks is provided below.