SEARCH · Search NASA
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.
Computer program for determining properties of high speed turbulent boundary layer with heat transfer and arbitrary pressure gradient
The outer surface of the high-Z inner shell in the double shell configuration of inertial confinement fusion experiments experiences Rayleigh–Taylor instability growth during the implosion process due to inverted density and pressure gradients between a highly compressed foam interstitial layer and the accelerating dense inner shell. Graded density layers have long been known to reduce instability growth rates. In this study, we employ high-fidelity radiation hydrodynamic simulations to demonstrate this improved stability when grading beryllium into tungsten. We first characterize the response to L-band preheat of these layers using a newly calibrated radiation drive. While graded layer capsules suffer reduced performance (here, measured as DD neutron yield from a CD foam fuel) in 1D simulations due to reduced kinetic energy coupling and reduced fuel compression, they suffer less of a performance drop when 2D instabilities are accounted for. With the improved stability of graded layers, we explore the performance of capsules with larger fuel radii and thinner shells as a preliminary study to find new designs in which graded layers produce the highest yields.
Compressibility transformation properties and high speed turbulent boundary layer flow with mass transfer and arbitrary pressure gradient
One of the most important outcomes of partial drop coalescence is the ratio of the secondary drop radius to the primary drop radius, known as the drop ratio, ri. This ratio is thought to be approximately constant and independent of physical parameters of the fluids involved. However, here we reveal that ambient fluid viscosity can alter the size of the secondary drop and the drop ratio consequently. Using scaling analysis, we derive a model that predicts the behavior of the drop ratio as a function of the Ohnesorge number, a dimensionless ratio of viscous to inertial forces. In addition, we present our experimental results of coalescing drops on a planar interface under the influence of surface tension gradients. A high-speed digital camera is used to observe the evolution of drops as they coalesce with a bulk liquid. We show that this process is influenced by the surface tension gradient between the drop and the bulk liquid. The ratio of the secondary drop to the primary drop in partial coalescence is smaller than the reported values for coalescence without a surface tension gradient. The analytical model derived through this study is based on a new modified Ohnesorge number that includes surface tension gradients. Our analytical model is compared against other models and the results illustrate good agreement with our experimental findings and experimental data in the literature.
One concept for high lift aerodynamics is the trapped-vortex. Recently it has been suggested that several spoilers located at different chordwise locations on an airfoil could be deployed to form several cavities in the chordwise direction. This may provide a means of increasing upper surface camber and thereby increase wing lift. It is envisioned that the cavities would be located on an airfoil in its pressure recovery region. Investigations of the effect of one or more cavities on the adjacent boundary layer in an adverse boundary layer were not found in a literature survey. Since flow separation can be caused by flow over cavities, boundary layer measurements were made in the vicinity of cavities which were located in an adverse pressure gradient to study the associated flow in an effort to identify any benefits. The experimental data from the present investigation indicate a minimal impact on the boundary layer profiles is caused by the presence of cavities in either a constant pressure freestream or an adverse pressure gradient. The data from the experimental investigation has been obtained and analyzed. The computational investigation using an incompressible Navier-Stokes code show that reasonable results are being obtained when they are compared with the experimental data.
Photovoltaic characteristics have been measured on a large number of crucible-grown lithium-containing solar cells irradiated by 1-MeV electrons to fluences ranging from 3 x 10 to the 13th power to 3 x 10 to the 15th power electrons per sq cm. These measurements have established empirical relationships between cell photovoltaic parameters and lithium donor density gradient. Short-circuit current and maximum power measured immediately after irradiation decrease logarithmically with lithium gradient. Open-circuit voltage increases logarithmically with lithium gradient both immediately after irradiation and after recovery, the degree of recovery being strongly gradient-dependent at high fluence. As a result, the maximum power and the power at 0.43 V after recovery from 3 x 10 to the 15th power electrons per sq cm increase with increasing lithium gradient.
We describe in detail the instrumentation and calibrations used in the ALE, GAGE and AGAGE experiments and present a history of the majority of the anthropogenic ozone- depleting and climate-forcing gases in air based on these experiments. Beginning in 1978, these three successive automated high frequency in-situ experiments have documented the long-term behavior of the measured concentrations of these gases over the past twenty years, and show both the evolution of latitudinal gradients and the high frequency variability due to sources and circulation. We provide estimates of the long-term trends in total chlorine contained in long- lived halocarbons involved in ozone depletion. We summarize interpretations of these measurements using inverse methods to determine trace gas lifetimes and emissions. Finally, we provide a combined observational and modeled reconstruction of the evolution of chlorocarbons by latitude in the atmosphere over the past sixty years which can be used as boundary conditions for interpreting trapped air in glaciers and oceanic measurements of chlorocarbon tracers of the deep oceanic circulation. Some specific conclusions are: (a) International compliance with the Montreal Protocol is so far resulting in chlorofluorocarbon and chlorocarbon mole fractions comparable to target levels, (b) Mole fractions of total chlorine contained in long-lived halocarbons (CCl2F2, CCl3F, CH3CCl3, CCl4, CHClF2, CCl2FCClF2, CH3Cl, CH2Cl2, CHCl3, CCl2=CCl2) in the lower troposphere reached maximum values of about 3.6 ppb in 1993 and are beginning to slowly decrease in the global lower atmosphere, (c) The chlorofluorocarbons have atmospheric lifetimes consistent with destruction in the stratosphere being their principal removal mechanism, (d) Multi-annual variations in chlorofluorocarbon and chlorocarbon emissions deduced from ALUGAGWAGAGE data are consistent approximately with variations estimated independently from industrial production and sales data where available (CCl2F2 (CFC-12) and CCl2FCClF2 (CFC-113) show the greatest discrepancies), (e) The mole fractions of the hydrochlorofluorocarbons and hydrofluorocarbons, which are replacing the regulated halocarbons, are rising very rapidly in the atmosphere but, with the exception of the much longer manufactured CHClF2 (HCFC-22), they are not yet at levels sufficient to contribute significantly to atmospheric chlorine loading. These replacement species could in the future provide independent estimates of the global weighted-average OH concentration provided their industrial emissions are accurately documented, (f) In the future, analysis of pollution events measured using high frequency in-situ measurements of chlorofluorocarbons and their replacements may enable emission estimates at the regional level which, together with industrial end-use data, are of sufficient accuracy to be capable of identifying regional non-compliance with the Montreal Protocol.
An oscillating, weakly ionized surface plasma has been investigated for use in turbulent boundary layer viscous drag reduction. The study was based on reports showing that mechanical spanwise oscillations of a wall can reduce viscous drag due to a turbulent boundary layer by up to 40%. It was hypothesized that the plasma induced body force in high electric field gradients of a surface plasma along strip electrodes could also be configured to oscillate the flow. Thin dielectric panels with millimeter-scale, flush- mounted, triad electrode arrays with one and two-phase high voltage excitation were tested. Results showed that while a small oscillation could be obtained, the effect was lost at a low frequency (less than 100Hz). Furthermore, a mean flow was generated during the oscillation that complicates the effect. Hot-wire and pitot probe diagnostics are presented along with phase-averaged images revealing plasma structure.
The focus of this research is on supercritical C7H16/N2 and O2/H2 mixing layers undergoing transitions to turbulence. The C7H16/N2 system serves as a simplified model of hydrocarbon/air systems in gas-turbine and diesel engines; the O2/H2 system is representative of liquid rocket engines. One goal of this research is to identify ways of controlling area production to increase disintegration of fluids and enhance combustion in such engines. As used in this research, "area production" signifies the fractional rate of change of surface area oriented perpendicular to the mass-fraction gradient of a mixing layer. In the study, a database of transitional states obtained from direct numerical simulations of the aforementioned mixing layers was analyzed to investigate global layer characteristics, phenomena in regions of high density-gradient magnitude (HDGM), irreversible entropy production and its relationship to the HDGM regions, and mechanisms leading to area production.
An analysis is given of some results of the LIMS experiment, launched on the Nimbus 7 satellite on October 24, 1978 for the purpose of sounding the middle atmosphere composition and structure. One of the LIMS channels was centered spectrally in the 6 micron region to measure vertical profiles and the global distribution of NO2 mixing ratio. The data reveal that NO2 is highly variable with altitude, latitude, longitude and time. Steep latitudinal gradients are observed in high latitude winter ('the NO2 cliff') when a wave number one pattern exists in the geopotential height field. The daytime gradient under these conditions is less than that for the night. The summer hemisphere column amount is considerably greater (by a factor of about 2) than in the winter hemisphere. The largest mixing ratios occur at 4 mb at night (about 20 ppbv) and 9 mb in the day (about 7 ppbv) and the latitude region of the peak mixing ratio is skewed toward the Southern Hemisphere. The mixing ratio variability is greater at night (by a factor of about 2 to 3) than in the day and there are significant long-term changes revealed by mixing ratio time series analyses on a pressure vs time grid.
Taking advantage of the improved radiometric resolution of Landsat-8 OLI which, unlike previous Landsat sensors, does not saturate over snow, the progress of fire recovery progress at the landscape scale (less than 100 m) is examined. High quality Landsat-8 albedo retrievals can now capture the true reflective and layered character of snow cover over a full range of land surface conditions and vegetation densities. This new capability particularly improves the assessment of post-fire vegetation dynamics across low- to high-burn severity gradients in Arctic and boreal regions in the early spring, when the albedos during recovery show the greatest variation. We use 30 m resolution Landsat-8 surface reflectances with concurrent coarser resolution (500 m) MODIS high quality full inversion surface Bidirectional Reflectance Distribution Functions (BRDF) products to produce higher resolution values of surface albedo. The high resolution full expression shortwave blue sky albedo product performs well with an overall RMSE of 0.0267 between tower and satellite measures under both snow-free and snow-covered conditions. While the importance of post-fire albedo recovery can be discerned from the MODIS albedo product at regional and global scales, our study addresses the particular importance of early spring post-fire albedo recovery at the landscape scale by considering the significant spatial heterogeneity of burn severity, and the impact of snow on the early spring albedo of various vegetation recovery types. We found that variations in early spring albedo within a single MODIS gridded pixel can be larger than 0.6. Since the frequency and severity of wildfires in Arctic and boreal systems is expected to increase in the coming decades, the dynamics of albedo in response to these rapid surface changes will increasingly impact the energy balance and contribute to other climate processes and physical feedback mechanisms. Surface radiation products derived from Landsat-8 data will thus play an important role in characterizing the carbon cycle and ecosystem processes of high latitude systems.
Mesoscale eddies play an important role in transport of heat and biogeochemical tracers in the global ocean circulation. Resolving these energetic eddies, however, is challenging in ocean general circulation models (OGCM) because it requires a horizontal grid spacing of ≲1/8° that is computationally expensive. As a result, we are required to parameterize mesoscale eddy effects on large-scale ocean flows. In this work, we introduce a new subgrid-scale (SGS) model that is developed based on a Taylor series expansion of resolved variables to parameterize subgrid mesoscale eddy transports and momentum fluxes in OGCM. We have performed an a priori study to evaluate the performance of our new gradient model using high-resolution ocean simulations. Our results show that the gradient model well predicts the actual SGS thickness fluxes in the zonal and meridional directions in coarse-resolution simulations with the grid spacing ≳1/4°. The unresolved kinetic energy at the ocean surface is also skillfully estimated. More importantly, unlike current mesoscale eddy parameterizations, which are mainly developed based on an assumption of flat bottom ocean, our new SGS model can capture the structure of unresolved standing meanders at the ocean surface. We have also developed a dynamic procedure for setting in non-dimensional parameters in our new parameterization through a non-ad hoc and tuning-free method. Overall, this work suggests that implementing the gradient model in OGCM can improve the model accuracy with an affordable computational cost in eddy-permitting and non-eddying simulations.
This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.
The objective of MULTI-LEADER is to cut design costs by 80% while generating more energy-efficient designs of multi-stage compressors by developing and implementing novel machine learning (ML) techniques, which enable faster and fewer design iterations, improved solver performance, and concurrent multi-disciplinary design. Current industrial practices for the design of multi-stage compressors involve simulation-based design optimization with successive levels of model fidelity, iteratively evaluated between distinct disciplines, one stage at a time to tackle the high dimensional design variations. This project addresses these key design challenges: (1) concurrent optimization of multiple stages under many non-linear constraints; (2) multitude of evaluation of high-fidelity and expensive solvers and their gradients during optimization convergence in high-dimensional design; (3) multi-disciplinary design to maximize aerodynamic performance while guaranteeing structural integrity and additive manufacturability; (4) utilization of multiple fidelity of solvers with disparate parameterization and modeling assumptions. MULTI-LEADER achieved more than 5x speed up in detailed design of more energy-efficient compressors via these machine learning (ML) innovations: (i) rapid design surrogates by multi-source learning from diverse fidelities across multiple disciplines, (ii) physics-constrained data-augmented modeling for improved empiricism, (iii) generative manifold embedding for high dimensional concurrent design without gradient information; (iv) budget-constrained fidelity-adaptive sampling towards fewer design iterations.
The high solidification rates during additive manufacturing cause highly localized thermal and strain gradients. The effect of these gradients on the evolution of local orientation misorientations within a grain is not well understood. In this study, stainless steel 316H parts were fabricated via laser powder bed fusion using three different energy densities: 43, 71, and 135 J/mm 3 . Electron backscatter diffraction showed that the maximum misorientations of the grains can be up to 25° along the build direction. Misorientation gradients (RM g ) within grains are process-dependent and can change from 0.036°/μm to 0.015°/μm with increased volumetric energy densities. In conclusion, the characterized misorientation gradients are an indication of the level of dislocations and, to an extent, the plastic deformation resulting from the rapid solidification during laser powder bed fusion.
Laser-driven (LD) ion acceleration has been explored in a newly constructed short focal length laser beamline at the BELLA petawatt facility (interaction point 2, iP2). For applications utilizing such LD ion beams, a beam transport system is required, which for reasons of compactness be ideally contained within 3 m. While they are generated from a micron-scale source, large divergence and energy spread of LD ion beams present a unique challenge to transporting them compared to beams from conventional accelerators. This study gives an overview of proposed compact transport designs using permanent magnets satisfying different requirements depending on the application for the iP2 laser beamline such as radiation biology, material science, and high-energy density science. These designs are optimized for different parameters such as energy spread and peak proton density according to the application’s need. The various designs consist solely of permanent magnet elements, which can provide high magnetic field gradients on a small footprint. While the field strengths are fixed, we have shown that the beam size is able to be tuned effectively by varying the placement of the magnets. The performance of each design was evaluated based on high-order particle tracking simulations of typical LD proton beams. We also examine the ability of certain configurations to tune and select beam energies, critical for specific applications. A more detailed investigation was carried out for a design to deliver 10 MeV LD accelerated ions for radiation biology applications. With these transport system designs, the iP2 laser beamline is ready to house various application experiments.
The influence of a constant acceleration on a turbulent premixed flame is studied by direct numerical simulation. This acceleration induces a mean pressure gradient across the flame brush, leading to a modification of the turbulent flame structure due to differential buoyancy mechanisms between heavy cold fresh and light hot burnt gases. Such a pressure gradient may be encountered in practical applications in ducted flames. A favorable pressure gradient, i.e. the pressure decreases from unburnt to burnt gases, is found to decrease the flame wrinkling, the flame brush thickness, and the turbulent flame speed. A favorable pressure gradient also promotes counter-gradient turbulent transport. On the other hand, adverse pressure gradients tend to increase the flame brush thickness and turbulent flame speed, and promote classical gradient turbulent transport. The balance equation for the turbulent flux of the Favre averaged progress variable is also analyzed. The first results show that the fluctuating pressure term, cannot be neglected as generally assumed in models. Simple models assuming that a high mean pressure gradient may only be balanced by the cross-dissipation term seem too approximate. This analysis has to be continued to compare simulation data and closure schemes proposed for the transport equation. The analysis developed by Veynante et al.(1995) has been extended to imposed acceleration and mean pressure gradients. A simple model for the turbulent flux is proposed and validated from simulation data. Then, a modified criterion is derived to delineate between counter-gradient and gradient turbulent diffusion. In fact, counter-gradient diffusion may occur in most practical applications, especially for ducted flames.
Retinal threshold along horizontal meridian for dark and light adapted eyes for stray light from small foveally fixated high luminance target