Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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,801 records · Page 100

Numerical Simulation of Ion Thruster Optics

A three-dimensional simulation code (ffx) designed to analyze ion thruster optics is described. It is an extension of an earlier code and includes special features like the ability to model a wide range of grid geometries, cusp details, and mis-aligned aperture pairs to name a few. However, the principle reason for advancing the code was in the study of ion optics erosion. Ground based testing of ion thruster optics, essential to the understanding of the processes of grid erosion, can be time consuming and costly. Simulation codes that can accurately predict grid lifetimes and the physical mechanisms of grid erosion can be of great utility in the development of future ion thruster optics designed for more ambitious applications. Results of simulations are presented that describe wear profiles for several standard and nonstandard aperture geometries, such as those grid sets with square- or slotted-hole layout patterns. The goal of this paper will be to introduce the methods employed in the ffx code and to briefly demonstrate their use.

Rawlin, Vincent K.↗

The scaling of kinetic and transport behaviors in the solution-phase chemistry of a plasma–liquid interface

The reactions at a plasma–liquid interface often involve species such as the solvated electron or the hydroxyl radical, which initiate the reduction or oxidation of solution-phase reactants (so-called scavengers) or are consumed by their own second-order recombination. Here, the mathematical scaling of the reaction–diffusion equations at the interface is used to obtain a characteristic time that can be used to determine the transition from highly efficient scavenger reduction or oxidation to lower efficiencies due to transport limitations. The characteristic time (tc) is validated using numerical solutions of the reaction–diffusion equations. When the scavenger kinetics are faster than second-order recombination, this characteristic transition time scales proportionally with the scavenger diffusivity (Ds) and the square of the scavenger bulk concentration (SB) and inversely proportional to the electron flux (J) squared; that is, t c = D s S B 2 F 2 /J 2 , where F is Faraday’s constant. However, when the scavenger kinetics are comparable or slower than second-order recombination, this scaling does not hold. Extending this analysis to three-dimensional systems shows that the profile of the electron flux on the surface affects the spatial location where reactions are most effective. Finally, the assessment of the implications of these behaviors for the reactor design highlights how effectively controlling the electron flux and solution transport may be necessary to improve the efficiency of scavenger reactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Large eddy simulation of turbulent channel flow: ILLIAC 4 calculation

The three-dimensional time dependent equations of motion were numerically integrated for fully-developed turbulent channel flow. A large scale flow field was obtained directly from the solution of these equations, and small scale field motions were simulated through an eddy viscosity model. The calculations were carried out on the ILLIAC 4 computer. The computed flow patterns show that the wall layer consists of coherent structures of low speed and high speed streaks alternating in the spanwise direction. These structures were absent in the regions away from the wall. Hot spots, small localized regions of very large turbulent shear stress, were frequently observed. The profiles of the pressure velocity-gradient correlations show a significant transfer of energy from the normal to the spanwise component of turbulent kinetic energy in the immediate neighborhood of the wall ('the splatting effect').

Kim, J.↗

Unsteady separated boundary layer in a transonic diffuser flow with self-excited oscillations

A numerical investigation of two-dimensional unsteady boundary layer in a transonic diffuser flow with self-excited oscillations and strong flow separation by solving the compressible, Reynolds-averaged, thin-layer Navier-Stokes equations with two-equations turbulence model is described. Three different meshes with constant streamwise mesh distribution and varying vertical mesh distribution were used. Results obtained indicate that a refinement of mesh studied here has minimal effect on the mean boundary layer flow but significantly increases the amplitude of oscillation of all flow variables. Comparisons of unsteady wall pressure, velocity profile, terminal shock, and separation pocket among computations and with experiment are presented.

Hsieh, T.↗

Numerical Solutions for the CAWAPI Configuration on Unstructured Grids at NASA LaRC, United States

In support of the Cranked Arrow Wing Aerodynamic Project International (CAWAPI) with its goal of improving the Technology Readiness Level of flow solvers by comparing results with measured F-16XL-1 flight data, NASA Langley employed the TetrUSS unstructured grid solver, USM3D, to obtain solutions for all seven flight conditions of interest. A newly available solver version that incorporates a number of turbulence models, including the two-equation linear and non-linear k- , was used in this study. As a first test, a choice was made to utilize only a single grid resolution with the solver for the simulation of the different flight conditions. Comparisons are presented with three turbulence models in USM3D, flight data for surface pressure, boundary-layer profiles, and skin-friction distribution, as well as limited predictions from other solvers. A result of these comparisons is that the USM3D solver can be used in an engineering environment to predict vortex-flow physics on a complex configuration at flight Reynolds numbers with a two-equation linear k- turbulence model.

Lamar, John E.↗

Experimental measurements of gas-particle flows in large-scale strippers

Measurements of gas–solid particle flow in a large-scale stripper unit are reported. The experiments target industrial (pilot) scale measurements of gas-particle flow that can be used to validate numerical methods for modeling multiphase flows, such as computational fluid dynamics coupled to the discrete element method (CFD-DEM). Specifically, experiments were performed with Geldart Group B (533 µm) glass beads in a stripper unit. In general, previous CFD-DEM validation datasets were performed in benchtop experimental systems and limited to O(10 5 ) particles. Here, in this work, the number of particles in the isolated section of the stripper is O(10 7 ) for all tests. Measurements are reported for a ~2 m tall section of the ~1 m diameter stripper. Radial mass flux profiles measured at five axial positions and three axial pressure drop measurements are reported for six operating conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Controlling pollutant emissions in a high-pressure combustor with fuel-diluent blending

This paper investigates the formation of nitric oxide (NO) and carbon monoxide (CO) emissions as a function of nitrogen (N 2 ) or carbon dioxide (CO 2 ) diluent content, added to a premixed reacting jet in crossflow. Here, reaction characteristics of a rich methane-air jet injected into a lean vitiated crossflow were analyzed at an elevated pressure of 5 atm. The jet was pre-heated and enriched with 0%, 15%, 30%, and 50% mass fraction to quantify the effect on pollutant emissions. Simulated results of the full chemistry Star-CCM+CFD model were verified with data taken in the experimental high-pressure combustion facility, which provides pressure, temperature and velocity profiles, as well as line-of-sight chemiluminescence and exit emission measurements. The analysis revealed the significant influence of the diluent to delay axial combustion; this mixing delay increased the axial flame lift-off and controlled any thermal ignition limitations. Increasing the diluent content increased the timescale for the flame to stabilize, which allowed for greater entrainment of crossflow oxidizer into the axial jet stream, and led to decreased pollutant emissions. Hence, crossflow entrainment is a critical driving force at high diluent content. Local nitric oxide (NO) emission formation in the axial stage was predicted numerically, showing the correlation between diluent addition, axial heat release, and the formation of nitric oxide pollution. The high diluent levels leaned out the jet mixture and delayed axial combustion, while reducing prompt and thermal NOx by mitigating flame hotspots and minimizing the timescale that the products remain at high temperatures.

03 NATURAL GAS↗

Thermodynamic stability of a spin microemulsion in Rashba spin-orbit-coupled bosons

Recent finite-temperature numerical simulations have unveiled a quantum “spin” microemulsion analog, found by raising the temperature of a stripe supersolid phase in a Rashba spin-orbit coupled Bose gas. This microemulsion state is a highly correlated, isotropic normal fluid where atoms self-arrange based on their internal pseudospin into patterns that resemble bicontinuous microemulsions. This finding leaves several open questions regarding the broader accessibility of this phase in experiments. Here, we use equilibrium finite-temperature numerical simulations based on a coherent-state path integral representation to perform a computational investigation into the thermodynamic stability of the spin microemulsion state. Numerical simulations emphasize the requirement of a nearly, but not perfectly, isotropic spin-orbit coupling in order to achieve the microemulsion phase in cold-atom experiments. Moreover, the microemulsion state exists independent of miscibility of the pseudospin components and for a wide range of pseudospin population imbalance, suggesting a high degree of flexibility in choosing the atom and hyperfine states in an experimental realization. Lastly, we demonstrate this feasibility by mimicking a Rashba spin-orbit-coupled 87 Rb experiment in an isotropic harmonic trap, where we confirm the microemulsion's existence via its density profile and equilibrium quasimomentum distribution.

Complex Langevin dynamics↗

Physics Guided Simulation of Electrostatic Discharge: Technical Report

Triboelectrically-charged objects may create threshold sparks, electrostatic discharge (ESD) events, to equilibrate charge between themselves and other relatively charged objects. ESD events exhibit many complex physical phenomena. They are a nexus of several fields of physics with disparate characteristic scales: plasma physics, chemical kinetics, hydrodynamics, circuit models, etc. These scales can span many orders of magnitude from the varied collisions thermalizing information within a plasma on the $\mathcal{O}(fs/ps)$ to the physical size of the plasma channel on the $\mathcal{O}(100µm)$, to the speed of a nonlinear hydrodynamic wave propagating at $\mathcal{O}(µm, ns)$. These threshold ESD events may occur in situations of programmatic importance, delivering energy and power profiles to a “victim load” generating deleterious consequences. To predict and mitigate these consequences we must answer questions about the spark’s energy budget: how much energy goes into producing the spark channel; how much gets radiated away; how much energy is advected away into the hydrodynamics; and how much energy is delivered to a victim load. An ESD simulation toolset has been created and evolved in order to answer these questions. An appropriate, physics guided implementation for simulation can be done by gaining insight into its constituent physics and leveraging that intuition to choose a suitable numerical operator. We examine in detail the chemical kinetics, circuit discharge, and hydrodynamics to deter mine dominant regimes, values, timescales, and interactions to uncover the underlying physical dynamics. We also examine and propose model reduction schemes for high-dimensional chemical kinetics. We use past and current work with experimentally validated and theoretically-verified hydrodynamics to calculate applicability limits of the non-ionizing strong shock limit. We quantify the energy budget from a hydrodynamic perspective and demonstrate that a significant fraction of the stored energy is “earmarked” for hydrodynamic advection as an energy terminus. Lastly, we combine the constituent physics of an ESD event (chemical kinetics, circuit model, and hydrodynamics) into a cohesive, actionable toolset and obtain promising results from an isothermal test case. We then propose a viable, modular evolution of the ESD toolset based upon the performed examination of the physics uncovering dominant physical scales and the stiffness of the compositional differential system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Acute Stent-Induced Endothelial Denudation: Biomechanical Predictors of Vascular Injury

Recent concern for local drug delivery and withdrawal of the first Food and Drug Administration-approved bioresorbable scaffold emphasizes the need to optimize the relationships between stent design and drug release with imposed arterial injury and observed pharmacodynamics. In this study, we examine the hypothesis that vascular injury is predictable from stent design and that the expanding force of stent deployment results in increased circumferential stress in the arterial tissue, which may explain acute injury poststent deployment. Using both numerical simulations and ex vivo experiments on three different stent designs (slotted tube, corrugated ring, and delta wing), arterial injury due to device deployment was examined. Furthermore, using numerical simulations, the consequence of changing stent strut radial thickness on arterial wall shear stress and arterial circumferential stress distributions was examined. Regions with predicted arterial circumferential stress exceeding a threshold of 49.5 kPa compared favorably with observed ex vivo endothelial denudation for the three considered stent designs. In addition, increasing strut thickness was predicted to result in more areas of denudation and larger areas exposed to low wall shear stress. We conclude that the acute arterial injury, observed immediately following stent expansion, is caused by high circumferential hoop stresses in the interstrut region, and denuded area profiles are dependent on unit cell geometric features. Such findings when coupled with where drugs move might explain the drug–device interactions.

Conway, Claire↗

Control theory based airfoil design using the Euler equations

This paper describes the implementation of optimization techniques based on control theory for airfoil design. In our previous work it was shown that control theory could be employed to devise effective optimization procedures for two-dimensional profiles by using the potential flow equation with either a conformal mapping or a general coordinate system. The goal of our present work is to extend the development to treat the Euler equations in two-dimensions by procedures that can readily be generalized to treat complex shapes in three-dimensions. Therefore, we have developed methods which can address airfoil design through either an analytic mapping or an arbitrary grid perturbation method applied to a finite volume discretization of the Euler equations. Here the control law serves to provide computationally inexpensive gradient information to a standard numerical optimization method. Results are presented for both the inverse problem and drag minimization problem.

Jameson, Antony↗

Crystal Growth of ZnSe by Physical Vapor Transport: A Modeling Study

Crystal growth from the vapor phase has various advantages over melt growth. The main advantage is from a lower processing temperature which makes the process more amenable in instances where the melting temperature of the crystal is high. Other benefits stem from the inherent purification mechanism in the process due to differences in the vapor pressures of the native elements and impurities, and the enhanced interfacial morphological stability during the growth process. Further, the implementation of Physical Vapor Transport (PVT) growth in closed ampoules affords experimental simplicity with minimal needs for complex process control which makes it an ideal candidate for space investigations in systems where gravity tends to have undesirable effects on the growth process. Bulk growth of wide band gap II-VI semiconductors by physical vapor transport has been developed and refined over the past several years at NASA MSFC. Results from a modeling study of PVT crystal growth of ZnSe arc reported in this paper. The PVI process is numerically investigated using both two-dimensional and fully three-dimensional formulation of the governing equations and associated boundary conditions. Both the incompressible Boussinesq approximation and the compressible model are tested to determine the influence of gravity on the process and to discern the differences between the two approaches. The influence of a residual gas is included in the models. The preliminary results show that both the incompressible and compressible approximations provide comparable results and the presence of a residual gas tends to measurably reduce the mass flux in the system. Detailed flow, thermal and concentration profiles will be provided in the final manuscript along with computed heat and mass transfer rates. Comparisons with the 1-D model will also be provided.

Ramachandran, Narayanan↗

Assessment Study of Small Space Debris Removal by Laser Satellites

Space debris in Earth orbit poses significant danger to satellites, humans in space, and future space exploration activities. In particular, the increasing number of unidentifiable objects, smaller than 10 cm, presents a serious hazard. Numerous technologies have been studied for removing unwanted objects in space. Our approach uses a short wavelength laser stationed in orbit to vaporize these small objects. This paper discusses the power requirements for space debris removal using lasers. A short wavelength laser pumped directly or indirectly by solar energy can scan, identify, position, and illuminate the target, which will then be vaporized or slow down the orbital speed of debris by laser detonation until it re-enters the atmosphere. The laser-induced plasma plume has a dispersive motion of approximately 105 m/sec with a Lambertian profile in the direction of the incoming beam [1-2]. The resulting fast ejecting jet plume of vaporized material should prevent matter recombination and condensation. If it allows any condensation of vaporized material, the size of condensed material will be no more than a nanoscale level [3]. Lasers for this purpose can be indirectly pumped by power from an array of solar cells or directly pumped by the solar spectrum [4]. The energy required for vaporization and ionization of a 10 cm cube (~ 2700 gm) of aluminum is 87,160 kJ. To remove this amount of aluminum in 3 minutes requires a continuous laser beam power of at least 5.38 MW under the consideration of 9% laser absorption by aluminum [5] and 5% laser pumping efficiency. The power needed for pumping 5.38 MW laser is approximately 108 MW, which can be obtained from a large solar array with 40% efficiency solar cells and a minimal area of 450 meters by 450 meters. This solar array would collect approximately 108 MW. The power required for system operation and maneuvering can be obtained by increasing solar panel size. This feasibility assessment covers roughly the power requirement, laser system, and a potential operational scenario.

Choi, Sang H.↗

Effect of thermal environments on fast charging Li-ion batteries

Battery thermal management systems (BTMSs) are expected to keep the battery temperature at a moderate level (~30 °C) to minimize the thermally exacerbated degradation. However, during fast charging, a strong cooling system is required to restrict the temperature rise of Li-ion batteries (LiBs), which significantly increases the cost and weight of battery packs, and induces a large temperature variation inside the battery. Here, in this work, we find that all these drawbacks could be relieved by allowing LiBs to charge at higher temperatures. Since the fast charging of a LiB only takes a tiny fraction of its lifetime, the aging rate is limited even at a charging temperature of 60 °C. Three types of thermal environments are proposed: kept constant at 30 °C, preheated to 60 °C, and adiabatic fast charging. With an experimentally validated electrochemical-thermal (ECT) coupled model, we explore the interplay between thermal management and the fast-charging performance. It is found that a gradually increasing temperature profile is the best option to balance the lithium plating and thermal management of the battery. Combining adiabatic fast charging with a preheating step, we can achieve minimal cooling need, perfect temperature uniformity within a battery, and fast-charging capability simultaneously.

25 ENERGY STORAGE↗

Analysis of thermal radiation of a gas-particulate-cloud in a pressurized oxy-coal combustor

Pressurized oxy-combustion (POC) aims to reduce carbon emissions from coal-fired power plants. In the POC process, pulverized coal is burned at elevated pressure in an O2-CO2 environment. The elevated pressure and high CO2 and H2O concentrations strongly impact thermal radiation, thereby significantly distinguishing POC system from conventional atmospheric-pressure air-fired and oxy-fuel combustion. Consequently, a thorough understanding of thermal radiation in POC is required to promote the development of novel combustors. The present computational work analyzes thermal radiation of the gaseous and particle phases in the pilot-scale POC combustor in Washington University in St. Louis. Particularly, the radiation characteristics of the gaseous and particulate cloud were estimated by employing the statistical narrow band model and the large-particle model. It is found that thermal radiation of a particulate cloud dominates in the combustor under the conditions of a furnace temperature of 1500 K and no substantial particle loss. Another important result is that the gas-and-particulate cloud can be approximately treated as a graybody under the same conditions. To be more specific, the results of spectral radiation intensity of a gas comprised of 40% (vol) H2O and 60% CO2 show that when the pressure is increased to 15 bar, and the radiation pathlength to 100 cm, the spectral radiation profile of the gas phase approaches that of a blackbody at the respective temperature. In addition, the emissivity of the particulate cloud has been evaluated as a function of the particle concentration and diameter by means of the large particle model. It is shown that the emissivity grows with the particle concentration but decreases with the particle size for the same mass of a particulate cloud. The emissivity of a particulate cloud in oxy- combustion exceeds that in air-fired combustion. This work could be used to validate the assumption of the gray gas model adopted in numerical simulations.

Li, Lei↗

The role of forced planetary waves in the annual cycle of the zonal mean circulation of the middle atmosphere

A severely truncated semispectral numerical model is used to simulate the annual cycle of the zonally averaged circulation in the middle atmosphere (16-96 km). The model includes only a single zonal harmonic wave component which interacts with the mean flow; the circulation is driven by diabatic heating and by a specified perturbation in the topography of the lower boundary, which is taken to be the 100 mb surface. A comparison of the annual cycle simulated by this model with the results of an analogous two-dimensional model indicates that planetary waves have relatively little influence on the zonal mean temperature profiles and on the solstice mean zonal winds at high latitudes. The primary effects of the forced waves are in decelerating the mean winds at low latitudes in the winter hemisphere to produce a region of weak westerlies, and in generating final warmings at the spring equinoxes.

Holton, J. R.↗

A comparison of radiative transfer models for predicting the microwave emission from soils

Noncoherent and coherent numerical models for predicting emission from soils are compared. Coherent models use the boundary conditions on the electric fields across the layer boundaries to calculate the radiation intensity, and noncoherent models consider radiation intensities directly. Interference may cause different results in the two approaches when coupling between soil layers in coherent models causes greater soil moisture sampling depths. Calculations performed at frequencies of 1.4 and 19.4 GHz show little difference between the models at 19.4 GHz, although differences are apparent at the lower frequency. A definition for an effective emissivity is also given for when a nonuniform temperature profile is present, and measurements made from a tower show good agreement with calculations from the coherent model.

Schmugge, T. J.↗

Evolution of rotationally and tidally distorted low-mass, close binary systems - Implications for the minimum orbital period of cataclysmic variables

A (1 + 0.4) solar mass close binary system consisting of a compact primary and a red dwarf secondary has been evolved numerically. Such a binary system should effectively model cataclysmic variables for which a minimum orbital period cutoff of about 81 minutes has been observed. The influence of gravitational radiation losses which drive Roche lobe overflow has been studied, and the effects of rotational and tidal distortion have also been incorporated in the calculation. The evolution of the He-3 abundance, which has been suggested as a possible explanation for the upper limit of the apparent orbital period gap exhibited by cataclysmic variables, is also considered. It is found that both the distortional effects and the He-3 chemical profile can play an important role in determining the subsequent evolution of these systems. Specifically, when distortion is included, the theoretical minimum orbital period is increased by about 10 percent, yielding better agreement with observations.

Nelson, L. A.↗

Refine your search to explore more results.