SEARCH · Engineering Papers

Results for “numerical modelling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Optimal gas cooling with minimal aero-opticaldistortion for next generation high-average-power lasers (Progress report, Year 2)

Next-generation lasers that combine high-peak-power ultrashort pulses with high repetition rates have the potential to efficiently deliver the high-average-power outputs required for inertial fusion energy generation, laser-based particle acceleration, and to generate compact secondary x- and 𝛾-ray radiation sources, among others. The combination of long gain length with high surface-to-volume ratio to remove residual heat makes gas-cooled multi-slab amplifier heads an inherently aperture- and average power scalable technology that enables this new class of solid-state lasers, and their thermal management becomes the most pressing technical challenge. As the index-of refraction depends on the gas density, heat removal from the slab surfaces unavoidably entails scattering of propagated light from regions of high temperature fluctuations. The requirement to maximize convective heat extraction from the gain medium under the constraint of minimum optical aberration of the beam renders amplifier slab cooling a technological challenge at the intersection of aero-optics, fluid-mechanics and heat transfer. The goal of this research is to leverage computational fluid dynamics and aero optics simulations of different levels of fidelity to devise a numerical modeling and optimization framework for gas-cooled multi-slab amplifiers. The project will be conducted in close collaboration with the Advanced Photon Technologies Group at Lawrence Livermore National Laboratory, where a gas-cooled test stand for laser materials and thermal management structures will be developed and installed. The final goal of the collaborative effort is to manufacture and benchmark a physical prototype of an numerically optimized amplifier head design. Measurement data on solid-state and thermo-optical distortions will be used for the validation of the numerical models. The key technical challenge in this multi-objective design problem is that of balancing optimal heat removal with minimal optical distortion. To facilitate computationally efficient design optimization and, at the same time, deepen our basic understanding of the interaction of light with turbulence, a complementary study that combines high and low-fidelity flow simulations and aero-optical models is proposed. Along the lines of these two objectives and solution strategies, the project is organized into two thrusts.

42 ENGINEERING↗

DOE OSTI2023

Effects of small marine energy deployments on oceanographic systems

The placement and operation of marine energy deployments in the ocean have the potential to change flow patterns, decrease wave heights, and/or remove energy from the oceanographic system. Changes in oceanographic systems resulting from harvesting marine energy, particularly tidal and wave energy, may be of concern. These changes include alterations in nearfield and farfield physical processes, as well as potential secondary environmental effects such as changes in sediment transport patterns, biological processes, or coastal erosion. Knowledge of changes in oceanographic systems associated with marine energy is primarily available from numerical modeling studies, informed by some laboratory tests and very few field measurements. A literature review was conducted using the Tethys knowledge base and other online sources, building on conclusions from the Ocean Energy Systems-Environmental State of the Science report. Potential changes in oceanographic systems that may be caused by marine energy differ between tidal and wave devices because of different extraction mechanisms and siting locations. Numerical models show that tidal extraction on the order of hundreds of megawatts or with significant channel blockage is required to create changes in oceanographic processes that exceed natural variability. Effects from wave energy extraction in arrays are localized and dependent on array spacing and proximity to the shore. Here the available evidence supports the conclusion that the risk of significant environmental effects from such changes could be retired (i.e., less investigation required for every project) for small deployments—those representative of the state of the industry in 2021. Determining changes in oceanographic systems to be low risk for small deployments can thereby streamline environmental consenting by reducing monitoring needs at this early stage in the industry.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Impurity leakage and radiative cooling in the first nitrogen and neon seeding study in the closed DIII-D SAS configuration

A comparative study of nitrogen versus neon has been carried out to analyze the impact of the two radiative species on power dissipation, SOL impurity distribution, divertor and pedestal characteristics. The experimental results show that N remains compressed in the divertor, thereby providing high radiative losses without affecting the pedestal profiles and displacing carbon as dominant radiator. Neon, instead, radiates more upstream than N thus reducing the power flux through the separatrix leading to a reduced ELM frequency and compression in the divertor. A significant amount of neon is measured in the plasma core leading to a steeper density gradient. The different behavior between the two impurities is confirmed by SOLPS-ITER modeling which for the first time at DIII-D includes multiple impurity species and a treatment of full drifts, currents and neutral–neutral collisions. The impurity transport in the SOL is studied in terms of the parallel momentum balance showing that N is mostly retained in the divertor whereas Ne leaks out consistent with its higher ionization potential and longer mean free path. This is also in agreement with the enrichment factor calculations which indicate lower divertor enrichment for neon. The strong ionization source characterizing the SAS divertor causes a reversal of the main ions and impurity flows. The flow reversal together with plasma drifts and the effect of the thermal force contribute significantly in the shift of the impurity stagnation point affecting impurity leakage. This work provides a demonstration of the impurity leakage mechanism in a closed divertor structure and the consequent impact on pedestal. Since carbon is an intrinsic radiator at DIII-D, in this paper we have also demonstrated the different role of carbon in the N vs Ne seeded cases both in the experiments and in the numerical modeling. Here, carbon contributes more when neon seeding is injected compared to when nitrogen is used. Finally, the results highlight the importance of accompanying experimental studies with numerical modeling of plasma flows, drifts and ionization profile to determine the details of the SOL impurity transport as the latter may vary with changes in divertor regime and geometry. In the cases presented here, plasma drifts and flow reversal caused by high level of closure in the slot upper divertor at DIII-D play an important role in the underlined mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

KOMPASS: Compaction of crushed salt for the safe containment

In Germany, rock salt formations are a possible host rock taken into account for the safe disposal of heat-emitting radioactive waste. With respect to crushed salt will be used in the repository for backfilling of open cavitied (using dry material). With time, the crushed salt will be compacted by the convergence of the host rock and reaches porosities comparable with the rock salts. The compaction behaviour of crushed salt has been investigated within the last 40 years, however, its behaviour at low porosities and the resulting low permeabilities becomes relevant with the introduction of the approach of the containment providing rock zone. In the current state, the database and process understanding have some important gaps in knowledge referring the material behaviour, existing laboratory and numerical models, especially for the porosity range. The objective of this project was the development of methods and strategies for the reduction of deficits in the prediction of crushed salt compaction leading to an improvement of the prognosis quality. It includes the development of experimental methods for determining crushed salt properties in the range of low porosities, the enhancement of process understanding and the investigation and development of existing numerical models.

Crushed salt compaction↗

NASA NTRS2013

The Martian Dust Cycle: Observations and Modeling

The dust cycle is critically important for Mars' current climate system. Suspended atmospheric dust affects the radiative balance of the atmosphere, and thus greatly influences the thermal and dynamical state of the atmosphere. Evidence for the presence of dust in the Martian atmosphere can be traced back to yellow clouds telescopically observed as early as the early 19th century. The Mariner 9 orbiter arrived at Mars in November of 1971 to find a planet completely enshrouded in airborne dust. Since that time, the exchange of dust between the planet's surface and atmosphere and the role of airborne dust on Mars' weather and climate has been studied using observations and numerical models. The goal of this talk is to give an overview of the observations and to discuss the successes and challenges associated with modeling the dust cycle. Dust raising events on Mars range in size from meters to hundreds of kilometers. During some years, regional storms merge to produce hemispheric or planet encircling dust clouds that obscure the surface and raise atmospheric temperatures by tens of kelvin. The interannual variability of planet encircling dust storms is poorly understood. Although the occurrence and season of large regional and global dust storms are highly variable from one year to the next, there are many features of the dust cycle that occur year after year. A low-level dust haze is maintained during northern spring and summer, while elevated levels of atmospheric dust occur during northern autumn and winter. During years without global-scale dust storms, two peaks in total dust loading are generally observed: one peak occurs before northern winter solstice and one peak occurs after northern winter solstice. Numerical modeling studies attempting to interactively simulate the Martian dust cycle with general circulation models (GCMs) include the lifting, transport, and sedimentation of radiatively active dust. Two dust lifting processes are commonly represented in these models: wind-stress lifting (i.e., saltation) and dust devil lifting. Although the predicted patterns of dust lifting and atmospheric dust loading from these simulations capture some aspects of the observed dust cycle, there are many notable differences between the simulated and observed dust cycles. For example, it is common for models to predict one peak in global dust loading near northern winter solstice due to excessive dust lifting in the Hellas basin at this season. Additionally, it is difficult for models to realistically capture the observed interannual variability in global dust storms. New avenues of dust cycle modeling research include exploring the effects of finite surface dust reservoirs and the effects of coupling the dust and water cycles on the predicted dust cycle.

Kahre, Melinda A.↗

NASA NTRS1998

NASA/NOAA/AMS Earth Science Electronic Theater

New method are demonstrated for visualizing, interpreting, comparing, organizing and analyzing immense HyperImage remote sensing datasets and three dimensional numerical model result. The traditional numerical spreadsheet paradigm has been extended to develop a scientific visualization approach for processing HyperImage datasets and three dimensional model results interactively.

Hasler, F.↗

NASA NTRS1991

Laboratory and theoretical models of planetary-scale instabilities and waves

Meteorologists and planetary astronomers interested in large-scale planetary and solar circulations recognize the importance of rotation and stratification in determining the character of these flows. The two outstanding problems of interest are: (1) the origins and nature of chaos in baroclinically unstable flows; and (2) the physical mechanisms responsible for high speed zonal winds and banding on the giant planets. The methods used to study these problems, and the insights gained, are useful in more general atmospheric and climate dynamic settings. Because the planetary curvature or beta-effect is crucial in the large scale nonlinear dynamics, the motions of rotating convecting liquids in spherical shells were studied using electrohydrodynamic polarization forces to generate radial gravity and centrally directed buoyancy forces in the laboratory. The Geophysical Fluid Flow Cell (GFFC) experiments performed on Spacelab 3 in 1985 were analyzed. The interpretation and extension of these results have led to the construction of efficient numerical models of rotating convection with an aim to understand the possible generation of zonal banding on Jupiter and the fate of banana cells in rapidly rotating convection as the heating is made strongly supercritical. Efforts to pose baroclinic wave experiments for future space missions using a modified version of the 1985 instrument have led us to develop theoretical and numerical models of baroclinic instability. Some surprising properties of both these models were discovered.

Hart, John E.↗

NASA NTRS1992

Prediction of high frequency combustion instability in liquid propellant rocket engines

The present use of a numerical model developed for the prediction of high-frequency combustion stabilities in liquid propellant rocket engines focuses on (1) the overall behavior of nonlinear combustion instabilities (2) the effects of acoustic oscillations on the fuel-droplet vaporization and combustion process in stable and unstable engine operating conditions, oscillating flowfields, and liquid-fuel trajectories during combustion instability, and (3) the effects of such design parameters as inlet boundary conditions, initial spray conditions, and baffle length. The numerical model has yielded predictions of the tangential-mode combustion instability; baffle length and droplet size variations are noted to have significant effects on engine stability.

Kim, Y. M.↗

NASA NTRS1984

Predictability of the Earth's polar motion

A comprehensive, experimental study of the predictability of the polar motion using a homogeneous BIH (Bureau International de l'Heure) data set is presented. Based on knowledge of the physics of the annual and the Chandler wobbles, the numerical model for the polar motion is constructed by allowing the wobble periods to vary. Using an optimum base length of 6 years for prediction, this floating-period model, equipped with a non-linear least-squares estimator, is found to yield polar motion predictions accurate from 0.012 to 0.024 inches depending on the prediction length up to one year, corresponding to a predictability of 91-83%. This represents a considerable improvement over the conventional fixed-period predictor, which does not respond to variations in the apparent wobble periods. The superiority of the floating-period predictor to other predictors based on critically different numerical models is also demonstrated.

Chao, B. F.↗

NASA NTRS1988

Cage stability analysis for SSME HPOTP bearings

A numerical model of cage motion (CAGEDYN) was used to analyze the stability of bearing cages in the Space Shuttle main engine (SSME) high pressure oxygen turbopump (HPOTP). The stability of existing bearing geometries, as well as perturbations of these geometries, was analyzed for various operating conditions. Results of the analyses show that some combinations of operating parameters, exacerbated by the sparse lubrication that exist in the HPOTP bearings, can cause unstable cage oscillations. Frequencies of cage oscillations were predicted by the CAGEDYN numerical model by Fourier analysis of predicted cage motions. Under conditions that cause unstable cage motion, high frequency oscillations were predicted that could cause premature cage failures.

Merriman, T. L.↗

DOE OSTI2020

Pre-Ionization Controlled Laser Plasma Formation for Ignition Applications (Final Scientific/Technical Report)

This report consists of three grand funded papers. Abstracts for each follow. 1) In this work, we present experimental and modeling studies of air pre-ionization using ultraviolet (UV) laser pulses and its effect on laser breakdown of an overlapped near-infrared (NIR) pulse. Experimental studies are conducted with a 266 nm beam (fourth harmonic of Nd:YAG) for UV pre-ionization and an overlapped 1064 nm NIR beam (fundamental of Nd:YAG), both having pulse duration of 10 ns. Results show that the UV beam produces a pre-ionized volume which assists in breakdown of the NIR beam, leading to reduction in NIR breakdown threshold by factor of >2. Numerical modeling is performed to examine the ionization and breakdown of both beams. The modeled breakdown threshold of the NIR, including assist by pre-ionization, is in reasonable agreement with the experimental results. 2) The present contribution compares the energy absorption, optical emission, temperature, and fluid dynamics of ultraviolet (UV) λ = 266 nm and near infrared (NIR) λ = 1064 nm nanosecond laser induced plasmas in ambient air. For UV pulses at the conditions studied, energy absorption by the plasmas increases relatively gradually with laser pulse energy starting at delivered energy of E~8 mJ. Corresponding measurements of plasma luminosity show that the absorption of UV radiation does not necessarily result in visible plasma emission. For the NIR induced plasmas, the energy absorption profile is far more abrupt and begins at ~55 mJ. In contrast with UV, the absorption of NIR radiation is always accompanied by intense optical emission. The temperatures of both types of plasma have been measured with Rayleigh scattering thermometry (at times after the Thomson signal sufficiently diminishes). The UV plasmas can attain a wider range of temperatures, including lower temperatures, depending on the pulse energy (e.g., T ~400–2000K for E ~7–35 mJ at Δt = 10 μs after the pulse) while the NIR plasmas show only hotter temperatures (e.g., T ~ 12 000 K for E = 75 mJ at Δt = 10 μs after the pulse) as is consistent with the literature. Differences in the fluid dynamics for UV versus NIR pulses are shown with Schlieren imaging. In this work, the contrast in the UV and NIR plasma threshold behavior is attributed to differing roles of avalanche ionization and multiphoton ionization as is also illustrated by a simple numerical model. 3) The present contribution examines the impact of plasma dynamics and plasma-driven fluid dynamics on the flame growth of laser ignited mixtures and shows that a new dual-pulse scheme can be used to control the kernel formation process in ways that extend the lean ignition limit. We perform a comparative study between (conventional) single-pulse laser ignition (λ = 1064 nm) and a novel dual-pulse method based on combining an ultraviolet (UV) pre-ionization pulse (λ = 266 nm) with an overlapped near-infrared (NIR) energy addition pulse (λ = 1064 nm). We employ OH* chemiluminescence to visualize the evolution of the early flame kernel. For single-pulse laser ignition at lean conditions, the flame kernel separates through third lobe detachment, corresponding to high strain rates that extinguish the flame. In this work, we investigate the capabilities of the dual-pulse to control the plasma-driven fluid dynamics by adjusting the axial offset of the two focal points. In particular, we find there exists a beam waist offset whereby the resulting vorticity suppresses formation of the third lobe, consequently reducing flame stretch. With this approach, we demonstrate that the dual-pulse method enables reduced flame speeds (at early times), an extended lean limit, increased combustion efficiency, and decreased laser energy requirements.

Mechanical engineering↗

NASA NTRS1991

A study of baroclinic instability in a cylindrical annulus with the temperature gradient imposed on the lower surface

Results are presented from experimental and numerical-modeling studies carried out on a rotating thermally driven fluid system in a cylindrical annulus with the horizontal gradient imposed upon the lower horizontal surface. The rotation rate and the temperature difference were varied to construct a regime diagram in thermal Rossby-number/Taylor-number space. It is shown that the curve separating the axisymmetric flow from the wave flow is 'knee-shaped', similar to the conventional side-heated and side-cooled baroclinic annulus. The experimentally observed transition curve agrees well with numerical calculations and the agreement between the predicted and measured wavenumbers is good both near the transition and within the wave regime. Away from the transition curve and well within the wave regime, there is an evolution from larger to smaller wavenumbers as the flow equilibrates, which is well captured by the numerical model for these cases.

Miller, Timothy L.↗

NASA NTRS1979

Time optimal control of a jet engine using a quasi-Hermite interpolation model

This work made preliminary efforts to generate nonlinear numerical models of a two-spooled turbofan jet engine, and subject these models to a known method of generating global, nonlinear, time optimal control laws. The models were derived numerically, directly from empirical data, as a first step in developing an automatic modelling procedure.

Comiskey, J. G.↗

NASA NTRS2003

Mechanisms of Rotating Instability in Axial Compressors Investigated

Rotating instability is a phenomenon that occurs in the tip flow region of axial compressor stages during stable operation. It can be observed in highly staggered rotors with significant tip clearance and is strongest at high-load operating points where the characteristic levels off. In this condition, the single-stage fan under investigation radiates an audible, whistling tone, and wall pressure spectra in the vicinity of the rotor disk exhibit nonrotational components. The graph shows the spectrum of static pressure at a point on the endwall near the leading edge. A hump appears at roughly half of the blade passing frequency (BPF) and is characteristic of rotating instability. A computational model was developed at the NASA Glenn Research Center to investigate the mechanism behind this phenomenon. A three-dimensional steady Navier-Stokes code that has been successfully tested for a wide range of turbomachinery flows was modified to execute a time-accurate simulation of the full annulus of the compressor. At the inlet of the computational domain, the total pressure, total temperature, and two velocity components are specified. Since no unsteady measurements of static pressure or other flow variables were available downstream of the rotor, circumferentially averaged static pressure was specified on the shroud at the outlet of the computational domain. A three-dimensional view of the vortex from the numerical model is shown. Particle traces released near the leading edge tip have rolled up to illustrate the tip clearance vortex. Flow near the trailing edge is pushed forward by the axially reversed flow. It then interacts with the tip clearance flow and the incoming flow and results in the rotating instability vortex, the core of which is illustrated by total pressure shading on planes located successively downstream. The rotating instability vortex is formed periodically midway between the blades and moves toward the pressure side of the passage. The unsteady behavior of this vortex structure is the main mechanism of the rotating instability is shown. The numerical model can be used to detect any possible occurrence of rotating instability when the tip clearance increases during engine service.

Hah, Chunill↗

NASA NTRS1988

Fluid dynamic modeling and numerical simulation of low-density hypersonic flow

The concept of a viscous shock-layer and several related versions of continuum theories/methods are examined for their adequacy as a viable framework to study flow physics and aerothermodynamics of relevance to sustained hypersonic flights. Considering the flat plate at angle of attack, or the wedge, as a generic example for the major aerodynamic component of a hypersonic vehicle, the relative importance of the molecular-transport effects behind the shock (in the form of the 'shock slip') and the wall-slip effects are studied. In the flow regime where the shock-transition-zone thickness remains small compared to the shock radius of curvature, a quasi-one-dimensional shock structure under the Burnett/thirteen-moment approximation, as well as particulate/collisional models, can be consistently developed. The fully viscous version of the shock-layer model is shown to provide the crucial boundary condition downstream the shock in this case. The gas-kinetic basis of the continuum description for the flow behind the bow shock, and certain features affecting the non-equilibrium flow chemistry, are also discussed.

Cheng, H. K.↗

NASA NTRS1981

Heat storage in alloy transformations

The feasibility of using metal alloys as thermal energy storage media was determined. The following major elements were studied: (1) identification of congruently transforming alloys and thermochemical property measurements; (2) development of a precise and convenient method for measuring volume change during phase transformation and thermal expansion coefficients; (3) development of a numerical modeling routine for calculating heat flow in cylindrical heat exchangers containing phase change materials; and (4) identification of materials that could be used to contain the metal alloys. Several eutectic alloys and ternary intermetallic phases were determined. A method employing X-ray absorption techniques was developed to determine the coefficients of thermal expansion of both the solid and liquid phases and the volume change during phase transformation from data obtained during one continuous experimental test. The method and apparatus are discussed and the experimental results are presented. The development of the numerical modeling method is presented and results are discussed for both salt and metal alloy phase change media.

Birchenall, C. E.↗

DOE OSTI2020

Ionization wave propagation in a He plasma jet in a controlled gas environment

Characterizing ionization wave propagation in low temperature plasma jets is critical to predicting production of reactive species and plasma–surface interactions for biomedical applications and surface functionalization. In this work, results from optical emission and laser induced fluorescence measurements of the ionization wave in a He plasma jet operating in a controlled gas environment are discussed and used for comparison with numerical modeling. Additionally, the ionization wave was observed using ICCD (Intensified Charge Coupled Device) imaging and characterized by time and spatially resolved electron density measurements using laser-collision-induced fluorescence. The plasma jet was initially characterized using pure He (nominally at 200 Torr), while varying pressure and voltage. When operating in pure He, the ionization wave broadly expands exiting the plasma tube. Increasing the operating pressure reduces the speed and isotropic expansion of the ionization wave. Furthermore, the jet operated with a humid He shroud was also studied. The humid He shroud results in the electron density increasing and having an annular profile due to the lower ionization potential of H2O compared to He and localized photoionization in the mixing region. Numerical modeling highlighted the importance of resonance radiation emitted by excited states of He, photoelectron emission from the quartz tube, and the kinetic behavior of the electrons produced by photoionization ahead of the ionization front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1984

Parametric analysis of a passive cyclic control device for helicopters

A parametric study of a passive device which provides a cyclic longitudinal control moment for a helicopter rotor was performed. It utilizes a rotor blade tip which is structurally decoupled from the blade inboard section. This rotor configuration is generally called the Free-Tip Rotor. A two dimensional numerical model was used to review the Constant Lift Tip Rotor, a predecessor of the current configuration, and then the same model was applied to the Passive Cyclic Control Device. The Constant Lift Tip was proven to have the ability to suppress the vibratory lift loading on the tip around the azimuth and to eliminate a significant negative lift peak on the advancing tip. The Passive Cyclic Control Device showed a once-per-revolution lift oscillation with a large amplitude, while minimizing the higher harmonic terms of the lift oscillation. This once-per-revolution oscillation results in the cyclic moment to trim the rotor longitudinally. A rotor performance analysis was performed by a three dimensional numerical model. It indicated that the vortices shed from the junction between the tip and the inboard section has a strong influence on the tip, and it may severely limit the tip performance. It was also shown that the Free-Tip allows the inboard section to have a larger twist, which results in a better performance.

Kumagai, H.↗