Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “V CODES”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Mixed Layer Sub-Mesoscale Parameterization - Part 1: Derivation and Assessment

Several studies have shown that sub-mesoscales (SM 1km horizontal scale) play an important role in mixed layer dynamics. In particular, high resolution simulations have shown that in the case of strong down-front wind, the re-stratification induced by the SM is of the same order of the de-stratification induced by small scale turbulence, as well as of that induced by the Ekman velocity. These studies have further concluded that it has become necessary to include SM in ocean global circulation models (OGCMs), especially those used in climate studies. The goal of our work is to derive and assess an analytic parameterization of the vertical tracer flux under baroclinic instabilities and wind of arbitrary directions and strength. To achieve this goal, we have divided the problem into two parts: first, in this work we derive and assess a parameterization of the SM vertical flux of an arbitrary tracer for ocean codes that resolve mesoscales, M, but not sub-mesoscales, SM. In Part 2, presented elsewhere, we have used the results of this work to derive a parameterization of SM fluxes for ocean codes that do not resolve either M or SM. To carry out the first part of our work, we solve the SM dynamic equations including the non-linear terms for which we employ a closure developed and assessed in previous work. We present a detailed analysis for down-front and up-front winds with the following results: (a) down-front wind (blowing in the direction of the surface geostrophic velocity) is the most favorable condition for generating vigorous SM eddies; the de-stratifying effect of the mean flow and re-stratifying effect of SM almost cancel each other out,

Canuto, V. M.↗
NASA NTRS2023

Role of Small-Scale Impulsive Events in Heating the X-Ray Bright Points of the Quiet Sun

Small-scale impulsive events, known as nanoflares, are thought to be one of the prime candidates that can keep the solar corona hot at its multimillion-Kelvin temperature. Individual nanoflares are difficult to detect with the current generation of instruments; however, their presence can be inferred through indirect techniques such as Differential Emission Measure (DEM) analysis. Here, we employ this technique to investigate the possibility of nanoflare heating of the quiet corona during the minimum of solar cycle 24. We estimate the DEM of disk-integrated quiet Sun and X-ray bright points (XBP) using the observations from XSM on board the Chandrayaan-2 orbiter and AIA on board the Solar Dynamic Observatory. XBPs are found to be the dominant contributor to disk-integrated X-rays, with a radiative flux of ∼2 × 10 5 erg cm −2 s −1 . XBPs consist of small-scale loops associated with bipolar magnetic fields. We simulate such XBP loops using the EBTEL hydrodynamic code. The lengths and magnetic field strengths of these loops are obtained through a potential field extrapolation of the photospheric magnetogram. Each loop is assumed to be heated by random nanoflares having an energy that depends on the loop properties. The composite nanoflare energy distribution for all the loops has a power-law slope close to −2.5. The simulation output is then used to obtain the integrated DEM. It agrees remarkably well with the observed DEM at temperatures above 1 MK, suggesting that the nanoflare distribution, as predicted by our model, can explain the XBP heating.

Solar coronal heating↗
DOE OSTI2024

Transcriptome-wide association analysis identifies candidate susceptibility genes for prostate-specific antigen levels in men without prostate cancer

Deciphering the genetic basis of prostate-specific antigen (PSA) levels may improve their utility for prostate cancer (PCa) screening. Using genome-wide association study (GWAS) summary statistics from 95,768 PCa-free men, we conducted a transcriptome-wide association study (TWAS) to examine impacts of genetically predicted gene expression on PSA. Analyses identified 41 statistically significant (p < 0.05/12,192 = 4.10 × 10 –6 ) associations in whole blood and 39 statistically significant (p < 0.05/13,844 = 3.61 × 10 –6 ) associations in prostate tissue, with 18 genes associated in both tissues. Cross-tissue analyses identified 155 statistically significantly (p < 0.05/22,249 = 2.25 × 10 –6 ) genes. Out of 173 unique PSA-associated genes across analyses, we replicated 151 (87.3%) in a TWAS of 209,318 PCa-free individuals from the Million Veteran Program. Based on conditional analyses, we found 20 genes (11 single tissue, nine cross-tissue) that were associated with PSA levels in the discovery TWAS that were not attributable to a lead variant from a GWAS. Ten of these 20 genes replicated, and two of the replicated genes had colocalization probability of >0.5: CCNA2 and HIST1H2BN. Six of the 20 identified genes are not known to impact PCa risk. Fine-mapping based on whole blood and prostate tissue revealed five protein-coding genes with evidence of causal relationships with PSA levels. Of these five genes, four exhibited evidence of colocalization and one was conditionally independent of previous GWAS findings. These results yield hypotheses that should be further explored to improve understanding of genetic factors underlying PSA levels.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1993

Helioseismological test of a new model for stellar convection

We compare the predictions of two solar models with the observed p-mode eigenfrequencies. The two models use the same input microphysics (nuclear reaction rates, opacity, and equation of a state) and the same numerical evolutionary code, but differ in the treatment of turbulent convection. The first model employs the standard mixing-length theory of convection (MLT), while the second model employs a new model of turbulent convection (CM) whose primary goal was that of accounting for the whole spectrum of turbulent eddies so as to avoid the MLT approximation that such a wide spectrum be represented by a single, large eddy. The main result of this paper is that the p-mode eigenfrequencies calculated with the CM model show an overall improvement with respect to those calculated with the standard MLT model.

Paterno, L.↗
NASA NTRS1997

Numerical modeling of HgCdTe solidification: Effects of phase diagram, double-diffusion convection and microgravity level

A numerical model of HgCdTe solidification was implemented using finite the element code FIDAP. Model verification was done using both experimental data and numerical test problems. The model was used to evaluate possible effects of double-diffusion convection in molten material, and microgravity level on concentration distribution in the solidified HgCdTe. Particular attention was paid to incorporation of HgCdTe phase diagram. It was found, that below a critical microgravity amplitude, the maximum convective velocity in the melt appears virtually independent on the microgravity vector orientation. Good agreement between predicted interface shape and an interface obtained experimentally by quenching was achieved. The results of numerical modeling are presented in the form of video film.

Bune, Andris V.↗
NASA NTRS2013

Aeromechanics Analysis of a Boundary Layer Ingesting Fan

Boundary layer ingesting propulsion systems have the potential to significantly reduce fuel burn but these systems must overcome the challe nges related to aeromechanics-fan flutter stability and forced response dynamic stresses. High-fidelity computational analysis of the fan a eromechanics is integral to the ongoing effort to design a boundary layer ingesting inlet and fan for fabrication and wind-tunnel test. A t hree-dimensional, time-accurate, Reynolds-averaged Navier Stokes computational fluid dynamics code is used to study aerothermodynamic and a eromechanical behavior of the fan in response to both clean and distorted inflows. The computational aeromechanics analyses performed in th is study show an intermediate design iteration of the fan to be flutter-free at the design conditions analyzed with both clean and distorte d in-flows. Dynamic stresses from forced response have been calculated for the design rotational speed. Additional work is ongoing to expan d the analyses to off-design conditions, and for on-resonance conditions.

Bakhle, Milind A.↗
NASA NTRS1959

Welded Titanium Case for Space-Probe Rocket Motor

Early in 1958, the Jet Propulsion Laboratory of the California Institute of Technology was requested to participate in a lunar-probe mission code-named Juno II which would place a 15-lb Instrumented payload (Pioneer IV) in the vicinity of the moon. The vehicle was to use the same high-speed upper-stage assembly as flown on the successful Jupiter-C configuration; however, the first-stage booster was to be a Jupiter rather than a Redstone. An analysis of the intended flight and payload configuration Indicated that the feasibility of accomplishing the mission was questionable and that additional performance would have to be obtained if the mission was to be feasible. Since the most efficient way of Increasing the performance of a staged vehicle is to increase the performance of the last stage, a study of possible ways of doing this was made.. Because of the time schedule placed on this effort It was decided to reduce the weight of the fourth-stage rocket-motor case by substituting the annealed 6Al--4V titanium alloy for the Type 410 stainless steel. Although this introduced an unfamiliar material, It reduced the changes in design and fabrication techniques. This particular titanium alloy was chosen on the basis of previous tests which proved the suitability of the alloy as a pressure-vessel material when used at an annealed yield strength of about 120, 000 psi. The titanium-case fourth stage of Juno U is shown with the payload and on the missile in Fig. 1; the stainless-steel motor cases used in the Jupiter-C vehicle are shown in Fig. 2. The fourth-stage motor case has a diameter of 6 in., a length of approximately 38 in. center dot and a nominal cylindrical wall thickness of 0.025 in. As shown in Fig. 1, the case serves as the structural support of the payload and is aligned to the upper stage assembly through an alignment ring. The nozzle is threaded into the end of the motor case, and is of the ceramic-coated steel design. Figure 3 shows a comparison of the components used to make the stainless steel and the 6A1--4V titanium alloy cases. The forward dome and aft fitting for the stainless steel assembly were fabricated from a combination of forged, spun and machined parts.. In order to facilitate the fabrication of the titanium alloy motor ) these components were machined from a large-diameter billet.

Brothers, A. J.↗
NASA NTRS1973

System statistical reliability model and analysis

A digital computer code was developed to simulate the time-dependent behavior of the 5-kwe reactor thermoelectric system. The code was used to determine lifetime sensitivity coefficients for a number of system design parameters, such as thermoelectric module efficiency and degradation rate, radiator absorptivity and emissivity, fuel element barrier defect constant, beginning-of-life reactivity, etc. A probability distribution (mean and standard deviation) was estimated for each of these design parameters. Then, error analysis was used to obtain a probability distribution for the system lifetime (mean = 7.7 years, standard deviation = 1.1 years). From this, the probability that the system will achieve the design goal of 5 years lifetime is 0.993. This value represents an estimate of the degradation reliability of the system.

Lekach, V. S.↗
NASA NTRS1987

Performance and aerodynamic braking of a horizontal-axis wind turbine from small-scale wind tunnel tests

Wind tunnel tests of three 20" diameter, zero twist, zero pitch wind turbine rotor models were conducted in a 7' x 10' wind tunnel to determine the performance of such rotors with NACA 23024 and NACA 64 sub 3-621 airfoil sections. Aerodynamic braking characteristics of a 38% span, 30% chord, vented aileron configuration were measured on the NACA 23024 rotor. Surface flow patterns were observed using fluorescent mini-tufts attached to the suction side of the rotor blades. Experimental results with and without ailerons are compared to predictions using airfoil section data and a momentum performance code. Results of the performance studies show that the 64 sub 3-621 rotor produces higher peak power than the 23024 rotor for a given rotor speed. Analytical studies, however, indicate that the 23024 should produce higher power. Transition strip experiments show that the 23024 rotor is much more sensitive to roughness than the 64 sub 3-621 rotor. These trends agree with analytical predictions. Results of the aileron test show that this aileron, when deflected, produces a braking torque at all tip speed ratios. In free wheeling coastdowns the rotor blade stopped, then rotated backward at a tip speed ratio of -0.6.

Cao, H. V.↗
NASA NTRS1988

Turbulent magnetohydrodynamic density fluctuations

A spectral-method numerical code is used to compute mass-density fluctuation spectra in turbulent magnetofluids. The computations are used to test and extend the analytical theory of density variations in slightly compressible magnetofluids given by Montgomery, et al. (1987) and used to infer inertial-range density-fluctuation spectra for the nearby interstellar medium and solar wind. A local equation of state is assumed, relating density to pressure. Constant, scalar resistivities and viscosities are used. In the limit of low Mach numbers and high mechanical-to-magnetic pressure ratios, the fit of the computations to the analytical theory is seen to be close.

Shebalin, John V.↗
NASA NTRS1988

Overview of CFD methods and comparison with flight aerothermal data

An overview of previously published aerothermal investigations which demonstrate the capabilities of detailed computational fluid dynamics and engineering codes to predict the aerothermal environment about an entry vehicle is presented. The overview consists of a brief discussion of the computational methods and experimental data and includes comparisons between the computed results and data. The overview focuses primarily on analyses of flight data since these data provide the unique capability to assess the real-gas chemistry options in the codes. The computed results are based on a series of codes which are employed by the Aerothermodynamics Branch of the Space System Division at the Langley Research Center. The flight data, which were measured on the Reentry F, the Space Shuttle, and the Fire II vehicles, represent a wide range of vehicle configurations and freestream conditions. Also, results of one recent set of ground tests are included since the tests provide data on a model of a pending flight project. The comparisons of the predicted results and data demonstrate the adequacy of the present computational fluid dynamics capabilities and indicate the potential to predict the aerothermal environment about future flight vehicles.

Sutton, Kenneth↗
NASA NTRS1993

Averaging techniques for steady and unsteady calculations of a transonic fan stage

It is often desirable to characterize a turbomachinery flow field with a few lumped parameters such as total pressure ratio or stage efficiency. Various averaging schemes may be used to compute these parameters. The momentum, energy, and area averaging schemes are described and compared. The schemes were compared for two computed solutions of the midspan section of a transonic fan stage: a steady averaging-plane solution in which average rotor outflow conditions were used as stator inflow conditions, and an unsteady rotor-stator interaction solution. The solutions were computed on identical grids using similar Navier-Stokes codes and an algebraic turbulence model. The unsteady solution is described, some unsteady flow phenomena are discussed, and the steady pressure distributions are compared. Despite large unsteady pressure fluctuations on the stator surface, the steady pressure distribution matched the average unsteady distribution almost exactly. Stator wake profiles, stator loss coefficient, and stage efficiency were computed for the two solutions with the three averaging schemes and are compared. In general, the energy averaging scheme gave good agreement between the averaging-plane solution and the time-averaged unsteady solution, even though certain phenomena due to unsteady wake migration were neglected.

Wyss, M. L.↗
NASA NTRS1993

Formulation for Simultaneous Aerodynamic Analysis and Design Optimization

An efficient approach for simultaneous aerodynamic analysis and design optimization is presented. This approach does not require the performance of many flow analyses at each design optimization step, which can be an expensive procedure. Thus, this approach brings us one step closer to meeting the challenge of incorporating computational fluid dynamic codes into gradient-based optimization techniques for aerodynamic design. An adjoint-variable method is introduced to nullify the effect of the increased number of design variables in the problem formulation. The method has been successfully tested on one-dimensional nozzle flow problems, including a sample problem with a normal shock. Implementations of the above algorithm are also presented that incorporate Newton iterations to secure a high-quality flow solution at the end of the design process. Implementations with iterative flow solvers are possible and will be required for large, multidimensional flow problems.

Hou, G. W.↗
NASA NTRS2011

Computational Analysis of the Large Scale Low-Boom Supersonic Inlet

This presentation describes two computational fluid dynamic (CFD) analyses done in support of a supersonic inlet test performed at NASA Glenn Research Center in the fall of 2010. The large-scale-low-boom supersonic inlet was designed for a small supersonic aircraft that would cruise at a Mach number of 1.6. It uses an axisymmetric, external compression spike to reduce the Mach number to 0.65 at the fan face. The inlet was tested in the 8x6 supersonic wind tunnel at NASA GRC using conventional pressure probes, pressure sensitive paint, and high-speed schlieren. Two CFD analyses of the inlet were performed before the test, and compared to the experimental data afterwards. Both analyses used the WIND-US code. First, an axisymmetric analysis of the inlet, diffuser, cold pipe, and mass flow plug was performed to predict the performance of the entire system in the wind tunnel. Then a 3-D analysis of the inlet with all its interior struts was performed to predict details of the flow field and effects of angle of attack. Test results showed that the inlet had excellent performance, with a peak total pressure recovery of 96 percent, and a buzz point far outside the engine operating range. The computations agreed very well with the data, with predicted recoveries within 0.3 - 0.5 points of the measurements.

Chima, ROdrick V.↗
NASA NTRSDate not supplied

Toward Verification of the γ-Reθt Transition Model in OVERFLOW and FUN3D

The results of an ongoing assessment of the transition modeling capability in NASA's OVERFLOW and FUN3D programs are presented, with a focus on the Langtry-Menter γ-Reθt transition model in combination with the shear-stress transport (SST) turbulence model. While the effect of numerics and boundary conditions on the accuracy and iterative convergence of RANS solutions for fully turbulent flows has been well documented, especially in the context of canonical flow configurations, the same cannot be said for transport-equation-based transition models coupled with RANS-based turbulence models. Given the criticality of transition modeling for new aircraft design and optimization, there has been a renewed focus on the accuracy of such transition models and their inconsistent implementation across different flow solvers as seen from the AIAA and NATO-AVT workshops. In this work, we aim to establish the verification of such models via these two well established CFD codes with different numerics. The goal is to produce high-quality data, such as grids, solutions, and other auxiliary data, that may be utilized for code verification by others in the CFD community. As a first step, the work reported here is focused on the Langtry-Menter transition model as applied to two simple 2D configurations, namely the flat plate, and the NLF-0416, respectively. The paper also highlights how the boundary conditions and baseline turbulence model can affect the solutions from the SST-based Langtry-Menter transition model. A preliminary evaluation of the automatic mesh adaption capabilities of these solvers and potential benefits for flow configurations involving a mix of laminar, transitional, and fully turbulent flows is also reported.

CFD Modeling↗
DOE OSTI2021

Integrated modeling methodology for ash agglomeration in poly-disperse fluidized beds using particle population framework

This article discusses a unique agglomeration modeling methodology developed based on binary collisions to combine the effects of heterogeneity in ash chemistry and granular physics. A simple population balance is defined to find changes in the particle size distribution (PSD) of a fluidized bed. Thermodynamic equilibrium calculations and a computational fluid dynamics (CFD) code, are used to obtain hydrodynamic parameters. A method to calculate and use a distribution of collision frequencies in fluidized beds, corresponding to a poly-disperse particle size distribution, was developed in order to incorporate the particle level heterogeneities. A distribution of collision frequencies obtained for poly-dispersed particles showed a three orders of magnitude higher collision frequency amongst the smaller particles than the coarser ones, at the initiation of agglomeration. Ash agglomeration occurred when the slag amount (binder) was less than 10 wt% at temperatures <850 °C for Pittsburgh No. 8 coal.

42 ENGINEERING↗
DOE OSTI2021

Dispersive shell pellet injection modeling and validation for DIII-D disruption mitigation

Dispersive shell pellet injection for disruption mitigation on the DIII-D tokamak is modeled with the 3D MHD code NIMROD. The simulations attempt to reproduce as closely as possible a pellet velocity scan that was performed on the DIII-D tokamak [Hollmann, et al., Phys. Rev. Lett. 122, 65001 (2019)]. This scan involved injecting a hollow diamond shell filled with boron dust that is released in the core of the plasma. Trends that are reproduced in the simulation include improved mitigation, reduced plasma current spike amplitude, and increased likelihood of hot-tail runaway electron production with increased pellet speed. The simulations provide additional physics insight into these findings. Limitations of the present model and requirements for a predictive model are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Resolution requirements for energy conservation in kinetic plasma simulations

The kinetic codes used to model the coupled dynamics of electromagnetic fields and charged particle transport have requirements for spatial, temporal, and charge resolution. These requirements may vary by the solution technique and scope of the problem. In this report, we investigate the resolution limits in the energy-conserving implicit particle-in-cell code CHICAGO. This report has the narrow aim of determining the maximum acceptable grid spacing for the dense plasmas generated in models of z-pinch target gases and power-flow electrode plasmas. In the 2D sample problem, the plasma drifts without external forces with velocity of 10 cm/µs. Simulations are scaled by plasma density to maintain uniform strides across the plasma and from the plasma to the boundaries. Additionally, the cloud-in-cell technique is used with 400 particles per cell and Δ t = 0.85× the Courant limit. For the linear cloud distribution, the criterion for conserving energy is Δ E/E tot < 0.01 for 50,000 time steps. The grid resolution is determined to crudely be Δx ≲ 3 l s , where l s is the electron collisionless skin depth. For the second-order cloud distribution the criterion is Δ E/E tot < 0.005 yielding Δx ≤ 15 l s . These scalings are functions of the chosen v d , Δ t , particles-per-cell, and number of steps.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗