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 559 records · Page 31

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems Using MOOSE-based Multiphysics Simulations

As an increasing number of companies look toward commercial Inertial Fusion Energy (IFE) designs, there is a pressing need to understand the physics of thick-wall chamber gas dynamics. The thick liquid wall approach implements a renewable wall to mitigate the fusion target emissions, thereby reducing the radiation damage rate and significantly extending the lifetime of chamber structures, leading to increased plant availability and reduced waste streams in comparison to dry wall chamber designs. It is necessary, however, to assess the critical performance and safety aspects of these systems. For example, it is crucial to predict (1) where the mass ablated from the liquid walls will vent, which determines the placement of condensing surfaces; (2) debris propagation up the beam lines, which provides essential information for design and protection requirements; (3) peak pressures and impulse on chamber walls, which affect chamber structural design; and (4) momentum transfer to the liquid jets, which constrains the shape and positioning of the jets. In turn, the chamber design and its liquid walls affect shielding requirements, material activation, and tritium fuel cycle. Currently available simulation tools, however, are unable to accurately capture key thick-wall chamber dynamics. Significant assumptions are often made to simplify the system and reduce computational cost and modeling capability needs, but the impact of these assumptions on simulation predictions has not been evaluated. For example, no three-dimensional simulations can be found in the open literature to evaluate gas venting and momentum transfer to the jets with simulations using two-dimensional domains to represent complex three-dimensional geometries. Moreover, limited studies have been dedicated to jet breakup due to both turbulence and neutron heating, and no studies have been found that evaluate how jet breakup can impact shock-jet interaction. Furthermore, effects of radiative heat transfer have rarely been included for the hydrodynamic phase of shock propagation, and integration of proper equations of state in shock dynamics codes has been mostly exploratory. In this study, we use the flexible, high-fidelity Multiphysics Object-Oriented Simulation Environment (MOOSE) to model these complex phenomena and inform design and safety studies. Capabilities to model thick-wall chamber gas dynamics are being developed, and the impact of the assumptions listed above (i.e., two-dimensional vs three-dimensional, absence of jet breakout, no radiative heat transfer, and ideal gas behavior) are being quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems using MOOSE

As an increasing number of companies look toward commercial Inertial Fusion Energy (IFE) designs, there is a pressing need to understand the physics of thick-wall chamber gas dynamics. The thick liquid wall approach implements a renewable wall to mitigate the fusion target emissions, thereby reducing the radiation damage rate and significantly extending the lifetime of chamber structures, leading to increased plant availability and reduced waste streams in comparison to dry wall chamber designs. It is necessary, however, to assess the critical performance and safety aspects of these systems. For example, it is crucial to predict (1) where the mass ablated from the liquid walls will vent, which determines the placement of condensing surfaces; (2) debris propagation up the beam lines, which provides essential information for design and protection requirements; (3) peak pressures and impulse on chamber walls, which affect chamber structural design; and (4) momentum transfer to the liquid jets, which constrains the shape and positioning of the jets. In turn, the chamber design and its liquid walls affect shielding requirements, material activation, and tritium fuel cycle. Currently available simulation tools, however, are unable to accurately capture key thick-wall chamber dynamics. Significant assumptions are often made to simplify the system and reduce computational cost and modeling capability needs, but the impact of these assumptions on simulation predictions has not been evaluated. For example, no three-dimensional simulations can be found in the open literature to evaluate gas venting and momentum transfer to the jets with simulations using two-dimensional domains to represent complex three-dimensional geometries. Moreover, limited studies have been dedicated to jet breakup due to both turbulence and neutron heating, and no studies have been found that evaluate how jet breakup can impact shock-jet interaction. Furthermore, effects of radiative heat transfer have rarely been included for the hydrodynamic phase of shock propagation, and integration of proper equations of state in shock dynamics codes has been mostly exploratory. In this study, we use the flexible, high-fidelity Multiphysics Object-Oriented Simulation Environment (MOOSE) to model these complex phenomena and inform design and safety studies. Capabilities to model thick-wall chamber gas dynamics are being developed, and the impact of the assumptions listed above (i.e., two-dimensional vs three-dimensional, absence of jet breakout, no radiative heat transfer, and ideal gas behavior) are being quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗

Prevention of resistive wall tearing mode major disruptions with feedback

Resistive wall tearing modes (RWTMs) can cause major disruptions. A signature of RWTMs is that the rational surface is sufficiently close to the wall to interact with it. For (m,n) = (2,1) modes, a RWTM requires normalized minor radius of the rational surface ρ q2 ≥ 0.75, which can also be expressed as q 75 ≤ 2. Major disruptions can occur when the criterion is satisfied. This is confirmed in simulations and theory and in a DIII-D locked mode disruption database. The q 75 < 2 criterion is valid at high β as well as at low β. A very important feature of RWTMs is that they can be feedback stabilized. If the ρ q2 criterion is not satisfied, or if the wall is ideally conducting, then the mode does not produce a major disruption, although it can produce a minor disruption. Feedback, or rotation of the mode at the wall by complex feedback, can emulate an ideal wall, preventing major disruptions. The ρ q2 criterion depends weakly on the wall radius. A simple geometric model of its dependence on wall radius is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of Channel Flow with Square Ribs for Blanket First-Wall Cooling: Geometry-Specific Tuning of k-ω Model Using Adjoint Method

Cooling of the plasma-facing first wall is challenging in the design of blanket components because of the high heat flux (on the order of 𝑀𝑊/𝑚2) from the plasma, especially when a low thermal mass medium like helium is chosen as the coolant. Therefore, heat transfer enhancement in which the convective heat transfer rate is augmented by the addition of turbulence-promoting structures becomes a key initiative for providing sufficient cooling capability with helium. Previously, computational fluid dynamics simulations had been performed on pipe flows with different transverse and longitudinal ribbed geometries at Oak Ridge National Laboratory to compare the enhancement performance among different ribbed geometries. Rib shape morphing had been conducted to obtain an optimized rib profile. In the work presented here, the adjoint method is adopted in the ANSYS Fluent solver for turbulence model augmentation, and the Generalized k-ω (GEKO) turbulence model is employed because of its ability of tuning the turbulence model. The Nusselt number and pressure drop obtained from the channel flow with bottom ribbed wall experiments are used as the targets. Sensitivity analysis provides information as guidance to improve the turbulence model accuracy. The augmented GEKO model is tuned for the studied ribbed channel geometry and flow conditions, providing improved predictive accuracy within this context. Extension to other configurations offers potential but may require additional tuning and validation.

Xu, Tracy [ORNL] (ORCID:0009000193700887)↗

Degradation of small mare surface features

Ring-moat structures are interpreted as remnants of flow topography associated with mare emplacement. Measurements of narrow moats characterizing these small mare features reveal a contradiction between their inferred preservation and the degree of surface degradation predicted by theoretical erosional models. Moat widths smaller than 50 m are recorded in areas where erosional models predict that walls of craters three to four times as large should be degraded to slopes of 1 deg. Crater statistics in such regions provide estimates of equilibrium conditions for impact crater formation that are consistent with the degree of preservation of the narrow moats. Although relative ages determined by observation of crater wall slopes near 20 deg, may be appropriate, it is concluded that current erosional models are not valid down to impact-degraded slopes of 1 deg. Relative age estimates based on the inferred degraded state of highly subdued craters are believed to be invalid, particularly for late Eratosthenian-age and Copernican-age mare units. In regions displaying ring-moat structures, highly subdued craters larger than the average moat widths are intepreted as either endogenic features associated with the emplacement of mare flow units or degraded secondary craters, which are typically subdued at the outset.

Schultz, P. H.↗

A model of filament-wound thin cylinders

A model was developed for simulating he manufacturing process of filament-wound cylinders made of a thermoset matrix composite. The model relates the process variables (winding speed, fiber tension, applied temperature) to the parameters characterizing the composite cylinder and the mandrel. The model is applicable to cylinders for which the diameter is large compared to the wall thickness. The model was implemented by a user-friendly computer code suitable for generating numerical results.

Calius, Emilio P.↗

Cascaded Perforates as One-Dimensional, Bulk Absorbers

Porous cell honeycomb liners for aircraft engine nacelles offer the possibility of exploiting extended reaction effects to improve liner attenuation bandwidth as generally attributed to the performance of bulk absorbers. This paper describes an analytical procedure, starting with an impedance prediction model for a single perforated plate, to estimate the bulk-absorber parameters for a cascade of such perforates - a first step to modeling a porous wall honeycomb structure. The objective is to build confidence in a lumped element impedance model, when applied to a uniformly-spaced set of porous plates to predict its .bulk. absorber properties. The model is based upon a modified version of the two-parameter flow resistance model of the form A + BV(sub inc), where A and B are physics-based, semi-empirical parameters that are adjusted to provide an optimum fit to a composite dataset from three plate porosities of 2.5, 5 and 10%. The composite dataset is achieved by reformulating the two-parameter flow resistance model into a .reduced pressure drop coefficient. dependency on perforate hole Reynolds number. The resulting impedance model is employed to calculate surface impedance spectra for N and 2N-layer perforate cascades. The well-known two-thickness method for experimental determination of bulk-absorber parameters is then applied to these .synthesized. data sets to predict the characteristic impedance and propagation constant for the perforate cascades. These results are then compared with experimental results reported in a companion paper.

Parrott, T. L.↗

Prediction of compliant wall drag reduction, part 2

A numerical model of turbulent boundary layer flows over compliant walls was investigated. The model is based on Burton's observation that outer flow structures in turbulent boundary layers produce large scale pressure fluctuations near the wall. The results of calculations indicate that certain small wavelength wall motions can have a significant effect upon the stability of turbulent boundary layers.

Orszag, S. A.↗

CFD Modeling of an 81" Single Slope Barrier Wall Drainage Window for on Grade and Sag Conditions

The main goal of this study is to design a drainage window in an 81" single slope barrier wall for on-grade and sag conditions that could replace the standard ODOT I-3 barrier catch basin. The ODOT design guidelines are as follows: a) Examine configurations of the opening to the flume behind the wall under the wall section to maximize conveyance of flow to a roadside ditch section. b) Examine window lengths and configurations to maximize hydraulic efficiency on a grade up to an absolute maximum window length of 20 ft. c) Use a minimum opening height of 4 in. and consider tapering to an increased height at the outlet to the flume to reduce the potential for debris clogging beneath the barrier section. On-grade locations are analyzed for longitudinal slope up to 5% and spread up to 10 ft. Sag locations, represented by pavement at zero longitudinal slope, are analyzed for water depths up to 18 in. The roadway cross section consists of a 10-ft shoulder at a 4.0% cross slope and pavement at a 1.6% cross slope. All surfaces are represented with Manning’s n = 0.015. Computational fluid dynamics (CFD) is used to perform the analysis. The modeling approach follows previous studies by the authors. The new design is similar to a slotted inlet: a rectangular opening cut into the bottom of a barrier (e.g., a temporary barrier or single slope barrier wall), which typically has a uniform cross section. Findings from laboratory tests of flow under Florida DOT barrier walls were used for additional validation of the modeling approach used in this study.

42 ENGINEERING↗

Direct drag measurement on thin-element riblets with suction and blowing

Turbulent, direct drag has been measured for a series of riblet models in air with continuous, low-level suction or blowing applied through narrow streamwise slots located at the peak, base, or valley of the rectangularly shaped riblets. Riblet spacing was approximately that of low-speed wall streaks. The models were designed to test whether wall vortices which are known to occur along widely spaced riblets are integral to turbulence production by attempting to alter them with localized suction or blowing. The drag data show trends that are consistent with reduced or enhanced inflectional breakdown of wall vortices. The trends, however, may also be explained by more fundamental mean-flow effects.

Wilkinson, Stephen P.↗

A simple model to scope out parameter space for indirect drive designs on NIF

In this work, we present a simple model to scope out parameter space for indirect-drive, inertial confinement fusion designs for the National Ignition Facility laser. Because the parameter space is large, simple models can be used to identify regions of parameter space for further study with more sophisticated models and experiments. We include a model for Hohlraum radiation drive and symmetry—both based on empirical scalings from the data. The model for radiation drive is based on assuming that the high atomic number (Z) Hohlraum wall dominates the energy balance during the high power, peak of the pulse (≳ 300 TW). We find that the time-dependent radiation drive flux can be described by the running integral of the laser energy divided by the Hohlraum area multiplied by constant slopes in two distinct time periods. The first period is when the laser power rises rapidly, so the radiation temperature increases due to changes in laser power and wall albedo. The second period is during peak power—here, the laser power is typically held constant—so, the radiation temperature increases only due to changes in the wall albedo. This model is applied to several NIF designs with different Hohlraum sizes, laser pulse length durations, and peak powers and energies. Drive and symmetry models can be combined to find regions of parameter space that have high capsule absorbed energy while maintaining a symmetric implosion. We propose a new metric for evaluating designs based on minimizing the radius at which the maximum implosion kinetic energy is achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Finite element stress analysis of the human left ventricle whose irregular shape is developed from single plane cineangiocardiogram

The three-dimensional left ventricular chamber geometrical model is developed from single plane cineangiocardiogram. This left ventricular model is loaded by an internal pressure monitored by cardiac catheterization. The resulting stresses in the left ventricular model chamber's wall are determined by computerized finite element procedure. For the discretization of this left ventricular model structure, a 20-node, isoparametric finite element is employed. The analysis and formulation of the computerised procedure is presented in the paper, along with the detailed algorithms and computer programs. The procedure is applied to determine the stresses in a left ventricle at an instant, during systole. Next, a portion (represented by a finite element) of this left ventricular chamber is simulated as being infarcted by making its active-state modulus value equal to its passive-state value; the neighbouring elements are shown to relieve the 'infarcted' element of stress by themselves taking on more stress.

Ghista, D. N.↗

Validation of a wall-interference assessment/correction procedure for airfoil tests in the Langley 0.3-meter Transonic Cryogenic Tunnel

Validation of a wall-interference assessment/correction (WIAC) procedure for a given facility requires its application to several sets of test data. It is necessary therefore to create a somewhat automated method for processing data through the various steps in the procedure. An automated procedure is also welcomed by the test engineer or eventual user in order to reduce the required effort and opportunity for error. Such a procedure has been developed for the Langley 0.3-m Transonic Cryogenic Tunnel using the TWINTN4 WIAC code. This code provides a four-wall, 2-D, transonic correction; that is, it accounts for sidewall boundary-layer effects, as well as for top and bottom wall effects on the airfoil tests. The TWINTN4 code utilizes measured pressure data at the tunnel walls and on the airfoil model; thus, classical homogeneous-wall boundary conditions are not used in the correction procedure.

Gumbert, C. R.↗

Monte Carlo Technique Used to Model the Degradation of Internal Spacecraft Surfaces by Atomic Oxygen

Atomic oxygen is one of the predominant constituents of Earth's upper atmosphere. It is created by the photodissociation of molecular oxygen (O2) into single O atoms by ultraviolet radiation. It is chemically very reactive because a single O atom readily combines with another O atom or with other atoms or molecules that can form a stable oxide. The effects of atomic oxygen on the external surfaces of spacecraft in low Earth orbit can have dire consequences for spacecraft life, and this is a well-known and much studied problem. Much less information is known about the effects of atomic oxygen on the internal surfaces of spacecraft. This degradation can occur when openings in components of the spacecraft exterior exist that allow the entry of atomic oxygen into regions that may not have direct atomic oxygen attack but rather scattered attack. Openings can exist because of spacecraft venting, microwave cavities, and apertures for Earth viewing, Sun sensors, or star trackers. The effects of atomic oxygen erosion of polymers interior to an aperture on a spacecraft were simulated at the NASA Glenn Research Center by using Monte Carlo computational techniques. A two-dimensional model was used to provide quantitative indications of the attenuation of atomic oxygen flux as a function of the distance into a parallel-walled cavity. The model allows the atomic oxygen arrival direction, the Maxwell Boltzman temperature, and the ram energy to be varied along with the interaction parameters of the degree of recombination upon impact with polymer or nonreactive surfaces, the initial reaction probability, the reaction probability dependence upon energy and angle of attack, degree of specularity of scattering of reactive and nonreactive surfaces, and the degree of thermal accommodation upon impact with reactive and non-reactive surfaces to be varied to allow the model to produce atomic oxygen erosion geometries that replicate actual experimental results from space. The degree of erosion of various interior locations was compared with the erosion that would occur external to the spacecraft. Results of one cavity model indicate that, at depths into a two-dimensional cavity that are equal to 10 cavity widths, the erosion on the walls of the cavity is less than that on the top surface by over 2 orders of magnitude. Wall erosion near the surface of a cavity depends on which wall is receiving direct atomic oxygen attack. However, deep in the cavity little difference is present. Testing of various cavity models such as these gives spacecraft designers an indication of the level of threat to sensitive interior surfaces for different geometries. Even though the Monte Carlo model is two-dimensional, it can be used to provide qualitative information about spacecraft openings that are three-dimensional by offering reasonable insight as to the nature of the attenuation of damage that occurs within a spacecraft in low Earth orbit. As shown, there is more erosion on the side seeing direct atomic oxygen attack until a depth of approximately 5 times the width of the opening, where the erosion is the same on both sides.

Banks, Bruce A.↗

Aeroacoustic Study of a Subscale Large Civil Transport (STAR) Model – Part 1: Simulations

Results from a computational aeroacoustic study of the Subsonic Transport Aeroacoustic Research (STAR) model are presented. The model, a 26%-scale semispan Boeing 777-200 aircraft, was an inceptive configuration in a comprehensive simulation campaign to accurately predict the airframe noise of a full-scale, large civil transport during landing. The lengthy process of obtaining a high-fidelity, CFD-ready, vetted, digital representation of the model is described in detail. The simulations were performed with the STAR model installed in the NASA Ames 40- by 80-foot full-scale wind tunnel to enable direct comparisons with aerodynamic and acoustic measurements of the model, presented separately in Part 2 of the study. Two versions of the wind tunnel test section were considered. The entire test section was used to simulate the aerodynamic behavior of the model; an open-wall representation of the section was adopted during the acoustic simulations to eliminate wall reflections. To preserve the aerodynamic equivalency of the model for the two simulated tunnel configurations, the open-wall computations were executed at an angle-of-attack one degree higher than that of the corresponding closed-wall tunnel. The effect on aerodynamic behavior of several minor geometry alterations necessary to con-form the simulated to the tested model was assessed – the modifications considerably altered the strength and extent of a flow separation zone on the inboard segment of the flap. Also evaluated was the impact of spatial resolution on aerodynamic and acoustic characteristics of the model – in general, medium resolution was sufficient to establish convergence trends for averaged surface pressure behavior and proper localization of noise sources. After the aeroacoustic characteristics of the baseline configuration were predicted, a toboggan-shaped fairing designed to reduce main landing gear noise was evaluated computationally. Far-field noise spectra for the baseline and toboggan-equipped STAR model were calculated via a Ffowcs-Williams and Hawkings integral approach, with flow quantities on a permeable data surface enclosing the source regions used as input. Results from the simulations indicated that, while the fairing reduced gear noise in the near-field, it had a negligible effect on the far-field noise signature of the model. This trend agreed with previously published experimental results obtained during sub-scale, isolated gear tests and full-scale flight tests

airframe noise↗

Earthquake-induced impact of base-isolated buildings: theory, numerical modeling, and design solutions

Earthquake shaking more intense than that used to size the horizontal clearance between a base–isolated building and near–rigid perimeter moat wall will result in hard impact, producing high–frequency, high–amplitude acceleration response in the structure and supported equipment. Herein this paper provides a design solution for the damaging effects of hard impact by installing a compliant engineered element in the load path between the base–isolated building and the moat wall, resulting in soft impact and a much smaller acceleration response. The engineered element assumed herein is a commercial–off–the–shelf marine fender with mechanical properties determined by physical testing. The attachment of a flexible engineered element, with well–defined stiffness and damping, to a near–rigid moat wall, simplifies the numerical modeling of the building–moat wall system and eliminates the need to bound the lateral stiffness of the wall for impact calculations. The simple model of the engineered element can be implemented in commercial finite element codes. Theory is developed for two–sided impact of a single–degree–of–freedom oscillator. Analytical solutions are derived for the shifted first–mode frequency of the impacted oscillator and for its free–vibration response. The shifted first–mode frequency is a function of the composite lateral stiffness of the isolator–engineered element assembly and its earthquake–induced displacement. Local peaks in the spectral response of the impacted oscillator form at odd integer multiples of the shifted first–mode frequency. The analytical solutions can be used to verify, in part, the numerical model used for impact analysis.

42 ENGINEERING↗

QCD-collapsed domain walls: QCD phase transition and gravitational wave spectroscopy

Abstract For a discrete symmetry that is anomalous under QCD, the domain walls produced in the early universe from its spontaneous breaking can naturally annihilate due to QCD instanton effects. The gravitational waves generated from wall annihilation have their amplitude and frequency determined by both the discrete symmetry breaking scale and the QCD scale. The evidence of stochastic gravitational waves at nanohertz observed by pulsar timing array experiments suggests that the discrete-symmetry-breaking scale is around 100 TeV, assuming the domain-wall explanation. The annihilation temperature is about 100 MeV, which could naturally be below the QCD phase transition temperature. We point out that the QCD phase transition within some domains with an effective large QCDθangle could be a first-order one. To derive the phase diagram inθand temperature, we adopt a phenomenological linear sigma model with three quark flavors. The domain-wall explanation for the NANOGrav, EPTA, PPTA and CPTA results hints at a first-order QCD phase transition, which predicts additional gravitational waves at higher frequencies. If the initial formation of domain walls is also a first-order process, this class of domain-wall models predicts an interesting gravitational wave spectroscopy with frequencies spanning more than ten orders of magnitude, from nanohertz to 100 Hz.

Physics↗