Search NASA⌕ Search

SEARCH · Search NASA

Results for “wall modeling”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Numerical modeling of runback water on ice protected aircraft surfaces

A numerical simulation for 'running wet' aircraft anti-icing systems is developed. The model includes breakup of the water film, which exists in regions of direct impingement, into individual rivulets. The wetness factor distribution resulting from the film breakup and the rivulet configuration on the surface are predicted in the numerical solution procedure. The solid wall is modeled as a multilayer structure and the anti-icing system used is of the thermal type utilizing hot air and/or electrical heating elements embedded with the layers. Details of the calculation procedure and the methods used are presented.

Al-Khalil, Kamel M.↗
NASA NTRS2002

The response of tissue-equivalent proportional counters to heavy ions

The paper presents a theoretical model for the response of a tissue-equivalent proportional counter (TEPC) irradiated with charged particles. Heavy ions and iron ions in particular constitute a significant part of radiation in space. TEPCs are used for all space shuttle and International Space Station (ISS) missions to estimate the dose and radiation quality (in terms of lineal energy) inside spacecraft. The response of the tissue-equivalent proportional counters shows distortions at the wall/cavity interface. In this paper, we present microdosimetric investigation using Monte Carlo track structure calculations to simulate the response of a TEPC to charged particles of various LET (1 MeV protons, 2.4 MeV alpha particles, 46 MeV/nucleon 20Ne, 55 MeV/nucleon 20Ne, 45 MeV/nucleon 40Ar, and 1.05 GeV/nucleon 56Fe). Data are presented for energy lost and energy absorbed in the counter cavity and wall. The model calculations are in good agreement with the results of Rademacher et al. (Radiat. Res. 149, 387-389, 1998), including the study of the interface between the wall and the sensitive region of the counter. It is shown that the anomalous response observed at large event sizes in the experiment is due to an enhanced entry of secondary electrons from the wall into the gas cavity.

NASA Discipline Radiation Health↗
NASA NTRS1982

A two-dimensional model of the hydrogen plasma for a laser powered rocket

A two-dimensional, closed-form model originally developed by Batteh and Keefer (1974) is modified and applied to the absorption of laser radiation by a hydrogen plasma. The model is used to predict the power absorbed by plasmas at one- and ten-atmosphere pressure as a function of laser beam radius. Predicted isotherms are given for one- and ten-atmosphere plasmas, together with thermal loading of the absorption chamber wall. The model is also used in predicting the laser power required to sustain a hydrogen plasma as a function of the absorption coefficient.

Keefer, D.↗
NASA NTRS1979

Modeling of thin, back-wall silicon solar cells

The performance of silicon solar cells with p-n junctions on the nonilluminated surface (i.e., upside-down or back-wall cells) was calculated. These structures consisted of a uniformly shaped p-type substrate layer, a p(+)-type field layer on the front (illuminated) surface, and a shallow, n-type junction on the back (nonilluminated) surface. A four-layer solar cell model was used to calculate efficiency, open-circuit voltage, and short-circuit current. The effect on performance of p-layer thickness and resistivity was determined. The diffusion length was varied to simulate the effect of radiation damage. The results show that peak initial efficiencies greater than 15 percent are possible for cell thicknesses or 100 micrometers or less. After 10 years of radiation damage in geosynchronous orbit, thin (25 to 50 micrometers thick) cells made from 10 to 100 ohm cm material show the smallest decrease (approximately 10 percent) in performance.

Baraona, C. R.↗
NASA NTRS1987

A finite element model for wave propagation in an inhomogeneous material including experimental validation

A finite element model was developed to solve for the acoustic pressure field in a nonhomogeneous region. The derivations from the governing equations assumed that the material properties could vary with position resulting in a nonhomogeneous variable property two-dimensional wave equation. This eliminated the necessity of finding the boundary conditions between the different materials. For a two media region consisting of part air (in the duct) and part bulk absorber (in the wall), a model was used to describe the bulk absorber properties in two directions. An experiment to verify the numerical theory was conducted in a rectangular duct with no flow and absorbing material mounted on one wall. Changes in the sound field, consisting of planar waves was measured on the wall opposite the absorbing material. As a function of distance along the duct, fairly good agreement was found in the standing wave pattern upstream of the absorber and in the decay of pressure level opposite the absorber.

Baumeister, Kenneth J.↗
NASA NTRS1987

A finite element model for wave propagation in an inhomogeneous material including experimental validation

A finite element model was developed to solve for the acoustic pressure field in a nonhomogeneous region. The derivations from the governing equations assumed that the material properties could vary with position resulting in a nonhomogeneous variable property two-dimensional wave equation. This eliminated the necessity of finding the boundary conditions between the different materials. For a two media region consisting of part air (in the duct) and part bulk absorber (in the wall), a model was used to describe the bulk absorber properties in two directions. An experiment to verify the numerical theory was conducted in a rectangular duct with no flow and absorbing material mounted on one wall. Changes in the sound field, consisting of planar waves, was measured on the wall opposite the absorbing material. As a function of distance along the duct, fairly good agreement was found in the standing wave pattern upstream of the absorber and in the decay of pressure level opposite the absorber.

Baumeister, Kenneth J.↗
DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

The effect of heavy ions on the formation and structure of cometary bow shocks

A hybrid simulation model is used to investigate the effects of heavy cometary ions on the formation and structure of a cometary bow shock. The calculations are carried out over various Mach numbers and heavy ion velocity distribution functions. The model is based on previous formulations for phenomena in the solar wind and at the earth's bow shock. The generation of the shock is described in terms of particles injected from one side of the simulation field and reflected from the other end of the field, i.e., a solid wall boundary model. This technique permits the steep buildup of the ion density near the cometary nucleus, followed by coupling of the incident and reflected ion streams to produce a shock. It is shown that at low Mach numbers (up to Mach 2) the shock is transitory and periodically formed by protons, then destroyed by heavy ions (O+). Slightly higher Mach numbers lead to a true stationary shock. An examination of coupling effects between the solar wind and the heavy ions at low Mach numbers by using the Rankine-Hugoniot relations reveals that the ions and the solar wind protons cannot be treated as a single fluid calculating the shock characteristics.

Omidi, N.↗
NASA NTRS2015

Uncertainty Quantification of Turbulence Model Closure Coefficients for Transonic Wall-Bounded Flows

The goal of this work was to quantify the uncertainty and sensitivity of commonly used turbulence models in Reynolds-Averaged Navier-Stokes codes due to uncertainty in the values of closure coefficients for transonic, wall-bounded flows and to rank the contribution of each coefficient to uncertainty in various output flow quantities of interest. Specifically, uncertainty quantification of turbulence model closure coefficients was performed for transonic flow over an axisymmetric bump at zero degrees angle of attack and the RAE 2822 transonic airfoil at a lift coefficient of 0.744. Three turbulence models were considered: the Spalart-Allmaras Model, Wilcox (2006) k-w Model, and the Menter Shear-Stress Trans- port Model. The FUN3D code developed by NASA Langley Research Center was used as the flow solver. The uncertainty quantification analysis employed stochastic expansions based on non-intrusive polynomial chaos as an efficient means of uncertainty propagation. Several integrated and point-quantities are considered as uncertain outputs for both CFD problems. All closure coefficients were treated as epistemic uncertain variables represented with intervals. Sobol indices were used to rank the relative contributions of each closure coefficient to the total uncertainty in the output quantities of interest. This study identified a number of closure coefficients for each turbulence model for which more information will reduce the amount of uncertainty in the output significantly for transonic, wall-bounded flows.

Schaefer, John↗
DOE OSTI2022

Improvement of two-phase closure models in CTF using Bayesian inference

Under the Consortium for Advanced Simulation of Light Water Reactors (CASL) program, extensive capabilities have been developed in CTF to analyze light-water reactors (LWRs) for normal operating conditions, departure from nucleate boiling (DNB), and system transients. However, further improvements are required in the modeling and simulation of boiling water reactors (BWRs), which is a focus of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this work, CTF validation results were used to optimize selected modeling coefficients by calibrating to experimental data using a Bayesian inference approach. Here, calibration studies were conducted to improve (vapor) void fraction prediction without worsening the two-phase pressure drop prediction, as well as to improve the two-phase pressure drop prediction. Calibration was performed for interfacial drag and wall shear models. Surrogates were developed to alleviate the computational expense required for sampling the parameter space using Markov chain Monte Carlo (MCMC). An assessment performed with calibrated models demonstrated an improvement of CTF in its prediction of key parameters such as void fraction and two-phase pressure drop.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1979

Three component velocity measurements in an upper surface blowing configuration

In four recent papers, results of an experimental investigation of a freely expanding coflowing jet as well as a three dimensional wall jet have been presented. A flat plate and a curved wall surface modelling a wing-flap combination in a high lift V/STOL configuration have been documented as to their effects on the flow field. The ratio of the jet exit plane velocity to the free stream velocity was varied from 5.1 to 10.9 to 1. This report continues the documentation of the three dimensional wall jet with special emphasis on the vertical and lateral mean and fluctuating velocities. Additionally, boundary layer profiles over both the flat plate and the flap are presented.

Catalano, G. D.↗
DOE OSTI2021

Modeling active control of resistive wall mode with power saturation and sensor noise on HL-2M

The resistive wall mode (RWM) control on the HL-2M tokamak is simulated with the MARS-F code (Liu et al 2000 Phys. Plasmas 7 3681), aiming at quantifying control current and voltage requirements when more realistic issues are taken into account, i.e. the control power saturation and the sensor signal noise. The fluid model predicts a narrow stability region for the n = 1 RWM without magnetic feedback, in the 2D parameter space of the plasma pressure versus the toroidal flow speed. Magnetic feedback can fully stabilize the RWM on HL-2M. Without considering the voltage limitation and the sensor signal noise, it is found that plasma flow helps active control of the mode, by reducing the required critical feedback gain for both flux-to-current and flux-to-voltage control schemes. In the absence of the sensor signal noise, the lowest control voltage saturation level, below which the RWM control is lost, is found to roughly satisfy a linear relation to the plasma flow frequency, indicating that subsonic plasma flow is effective in relaxing the control power requirement for the RWM feedback stabilization. The presence of the sensor signal noise substantially modifies the feedback results. A statistical study finds that the sensor signal noise, with the standard deviation of 0.1 G on HL-2M, roughly doubles the required control voltage for successful mode control. The synergistic stabilization effect due to plasma flow is somewhat weakened by the presence of the sensor signal noise. At a given rotation, the tolerable voltage limit generally increases with increasing feedback gain due to the sensor signal noise.

Physics↗
NASA NTRS2007

Modeling Creep-Induced Stress Relaxation at the Leading Edge of SiC/SiC Airfoils

Anticipating the implementation of advanced SiC/SiC composites into internally cooled airfoil components within the turbine section of future aero-propulsion engines, the primary objective of this study was to develop physics-based analytical and finite-element modeling tools to predict the effects of composite creep and stress relaxation at the airfoil leading edges, which will generally experience large thermal gradients at high temperatures. A second objective was to examine how some advanced NASA-developed SiC/SiC systems coated with typical EBC materials would behave as leading edge materials in terms of long-term steady-state operating temperatures. Because of the complexities introduced by mechanical stresses inherent in internally cooled airfoils, a simple cylindrical thin-walled tube model subjected to thermal stresses only is employed for the leading edge, thereby obtaining a best-case scenario for the material behavior. In addition, the SiC/SiC composite materials are assumed to behave as isotropic materials with temperature-dependent viscoelastic creep behavior as measured in-plane on thin-walled panels. Key findings include: (1) without mechanical stresses and for typical airfoil geometries, as heat flux is increased through the leading edge, life-limiting tensile crack formation will occur first in the hoop direction on the inside wall of the leading edge; (2) thermal gradients through all current SiC/SiC systems should be kept below approx.300 F at high temperatures to avoid this cracking; (3) at temperatures near the maximum operating temperatures of advanced SiC/SiC systems, thermal stresses induced by the thermal gradients will beneficially relax with time due to creep; (4) although stress relaxation occurs, the maximum gradient should still not exceed 300oF because of residual tensile stress buildup on the airfoil outer wall during cool-down; and (5) without film cooling and mechanical stresses, the NASA-developed N26 SiC/SiC system with thru-thickness Sylramic-iBN fiber reinforcement and a typical EBC coating has the potential of offering a maximum long-term steady-state operating temperature of approx.3100 F at the surface of the EBC.

Lang, Jerry↗
DOE OSTI2021

A Scoping Study to Develop a Computational Fluid Dynamics Based Model to Predict Radiological Materials Packaging Temperatures within a Generic Staging Building

The objective of this work is to perform computational fluid dynamics simulations of a ventilated radiological-material-package staging building to determine the effect of including or excluding a variety of physical effects and computational methods on both the predicted package temperatures and required computational resources. The generic building contains 640 drum-packages containing heat-generating, radiological material supported on four racks that are eight levels tall. Additionally, there is a forced ventilation system, lighting, and insulated walls. Computational models were constructed that included or excluded (a) shelving, (b) effects of unsteadiness, and (c) radiation heat transfer. Simulations with each drum modeled separately were compared to simpler simulations with sets of four drums represented by an equivalent box-package. These models employed between 106 to 107 elements. Steady state simulations predicted package temperatures that were within 0.1°C of simulations that included transient effects and required only one-eighth the computational resources. Calculations that excluded shelving predicted temperatures within 0.6°C of simulations that included shelving and required one-fourth the computational resources. Excluding radiation heat transfer systematically increased temperatures by around 1.5°C but reduced computational resources by a factor of four. The box-package model reduced the computational resources by a factor of 3, but systematically predicted higher temperatures by around 1.1°C. These results will be used to develop an efficient computational fluid dynamics model to assess the ability of different staging building designs to prevent the temperature of package components from exceeding specified limits.

42 ENGINEERING↗
NASA NTRS2003

Numerical Simulations of STOVL Hot Gas Ingestion in Ground Proximity Using a Multigrid Solution Procedure

A multi grid solution procedure for the numerical simulation of turbulent flows in complex geometries has been developed. A Full Multigrid-Full Approximation Scheme (FMG-FAS) is incorporated into the continuity and momentum equations, while the scalars are decoupled from the multi grid V-cycle. A standard kappa-Epsilon turbulence model with wall functions has been used to close the governing equations. The numerical solution is accomplished by solving for the Cartesian velocity components either with a traditional grid staggering arrangement or with a multiple velocity grid staggering arrangement. The two solution methodologies are evaluated for relative computational efficiency. The solution procedure with traditional staggering arrangement is subsequently applied to calculate the flow and temperature fields around a model Short Take-off and Vertical Landing (STOVL) aircraft hovering in ground proximity.

Wang, Gang↗
DOE OSTI2021

Sheath collapse at critical shallow angle due to kinetic effects

The Debye sheath is known to vanish completely in magnetised plasmas for a sufficiently small electron gyroradius and small angle between the magnetic field and the wall. This angle depends on the current onto the wall. When the Debye sheath vanishes, there is still a potential drop between the wall and the plasma across the magnetic presheath. The magnetic field angle corresponding to the predicted sheath collapse is shown to be much smaller than previous estimates, scaling with the electron-ion mass ratio and not with the square root of the mass ratio. This is shown to be a consequence of the kinetic electron and finite ion orbit width effects, which are not captured by fluid models. The wall potential with respect to the bulk plasma at which the Debye sheath vanishes is calculated. Above this wall potential, it is possible that the Debye sheath will invert.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Development of a near-wall Reynolds-stress closure based on the SSG model for the pressure strain

In this research, a near-wall second-order closure based on the Speziable et al.(1991) or SSG model for the pressure-strain term is proposed. Unlike the LRR model, the SSG model is quasi-nonlinear and yields better results when applied to calculate rotating homogeneous turbulent flows. An asymptotic analysis near the wall is applied to both the exact and modeled, equations so that appropriate near-wall corrections to the SSG model and the modeled dissipation-rate equation can be derived to satisfy the physical wall boundary conditions as well as the asymptotic near-wall behavior of the exact equations. Two additional model constants are introduced and they are determined by calibrating against one set of near-wall channel flow data. Once determined, their values are found to remain constant irrespective of the type of flow examined. The resultant model is used to calculate simple turbulent flows, near separating turbulent flows, complex turbulent flows and compressible turbulent flows with a freestream Mach number as high as 10. In all the flow cases investigated, the calculated results are in good agreement with data. This new near-wall model is less ad hoc, physically and mathematically more sound and eliminates the empiricism introduced by Zhang. Therefore, it is quite general, as demonstrated by the good agreement achieved with measurements covering a wide range of Reynolds numbers and Mach numbers.

So, R. M. C.↗
DOE OSTI2025

Optimization of spray breakup model parameters for predicting fuel spray and film characteristics in gasoline direct injection engines

This study investigated the behavior of gasoline direct injection (GDI) sprays using computational fluid dynamics (CFD). The authors developed an approach to identify optimal spray breakup model parameters by evaluating an error function across numerous simulations, with the goal of minimizing discrepancies from experimental data. Using the optimal setup, the simulated spray matched well with projected liquid volume distributions, liquid penetration, and spray width measured in a constant-pressure continuous-flow chamber. To further validate the approach, the same setup was tested across various fuels, injectors, and operating conditions. Subsequently, the optimal setup, along with a recently developed spray-wall interaction model, were applied to a direct-injection spark-ignited engine under late-injection conditions to predict and evaluate fuel film formation and evolution at varying engine coolant temperatures. Here, with the centrally mounted injector directing the spray toward the piston, simulations indicated that the spray tends to impinge on the piston surface. The proposed simulation framework also accurately captured the aggregate film area on the piston surface, aligning with previously published experimental results. Moreover, simulations showed that increasing the coolant temperature from cold start conditions (333 K) to warm conditions (363 K) reduced the fuel mass deposited on the piston by roughly 50%. Furthermore, for the spray-guided engine configuration studied in this work, the CFD model predicted minimal film deposition on the spark plug electrodes regardless of the coolant temperatures due to a relatively weak in-cylinder flow during the compression phase.

Computational fluid dynamics (CFD)↗