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 577 records · Page 32

The influence of line shape and band structure on temperatures in planetary atmospheres

Numerical experiments are performed to examine the effects of line shape and band structure on the radiative equilibrium temperature profile in planetary atmospheres. In order to accurately determine these effects, a method for calculating radiative terms is developed which avoids the usual approximations. It differs from the more commonly used methods in that it allows arbitrary dependence of the absorption coefficient on wave number, without requiring tedious line by line integration and without the constraints of band models. The present formulation is restricted to homogeneous atmospheres but the concept can be extended to the more general case. The numerical experiments reveal that the line shape and band structure of the absorbing gases have a large effect on temperatures in the higher layers of the atmosphere (corresponding to the stratosphere and mesosphere). The more nongrey the spectrum (that is, the higher the peaks and the deeper the troughs in the spectrum), the lower the temperature.

Arking, A.↗

High concentrations of ice: Investigations using polarimetric radar observations combined with in situ measurements and cloud modeling (Final Report)

This DOE-funded joint project had the over-arching goal of understanding the reasons for observed size distributions of ice particles in cold clouds. Its approach involves the use of cloud models and field observations by radar and aircraft of real storms. The project addresses three major research objectives: (1) Utilize a novel polarimetric radar technique to retrieve size distributions and amounts of ice from radar data collected during the previous DOE ARM field campaigns; (2) Evaluate the results of ice microphysical retrievals using available in situ aircraft measurements and other remote sensors (e.g., vertically pointing cloud radars and wind profilers); (3) Improve the treatment of microphysical processes leading to cloud glaciation and ice multiplication in numerical cloud models. The study was performed by three research organizations: University of Oklahoma, The Hebrew University of Jerusalem, Israel, and Lund University, Sweden. The Israel and Swedish partners were funded by the sub-awards from the University of Oklahoma. Cloud modeling studies have been performed by the research teams at The Hebrew University of Jerusalem (HUJ) and Lund University to identify the origins of high concentrations of cloud ice in areas of high ice water content (HIWC). The Hebrew University Cloud Model (HUCM) with full spectral bin microphysics and the Lund University aerosol-cloud (AC) model with a hybrid bin / bulk microphysics scheme complementing HUCM were utilized for simulations. Both research teams had particular focus on secondary ice production (SIP) as one of the possible sources of enhanced ice concentration. The HUJ group suggested a novel concept of ice multiplication during droplet freezing. It is assumed that splintering and droplet fragmentation during droplet freezing takes place because of dendritic growth within a supercooled drop. The resulting simulations of SIP generated small ice in concentrations exceeding hundreds per liter similar to what was observed in the HIWC regions of the tropical storms. The Lund team explored the SIP mechanisms such as breakup of ice particles due to ice-ice collisions and ice sublimation that are expected to dominate the continental storms. They also quantified the impact of homogeneous nucleation of cloud droplets on the total number concentration of ice at very low temperatures near the tops of the clouds. Additionally, the Lund AC model is able to simulate the effect of aerosols of various types (including biological) on the cloud life cycle and the corresponding ice production. The HUCM / AC model was used to simulate one of the “golden” cases of the DOE MC3E campaign on 20 May 2011. The model output was converted into the fields of polarimetric radar variables using the polarimetric radar forward operator developed at the University of Oklahoma and compared with radar observations and in situ microphysical measurements onboard research aircraft. It was demonstrated that specific differential phase KDP is the best radar parameter to identify the HIWC areas and quantify the corresponding ice parameters. For the first time, the shape of the vertical profile of KDP was realistically reproduced by the cloud model with the KDP maximum in the dendritic growth layer (DGL) centered at the -15°C isotherm. The University of Oklahoma team has developed a methodology for polarimetric radar retrievals of such microphysical parameters of ice as ice water content (IWC), mean volume diameter (Dm), and total number concentration (Nt) of ice particles. These retrievals have been validated using in situ aircraft measurements during 6 field campaigns and proved to be quite robust and reliable. This allowed to build the first climatology of the vertical profiles of polarimetric radar variables and retrieved microphysical parameters for the three types of weather systems: continental MCSs, maritime MCSs, and tropical cyclones / hurricanes (Hu and Ryzhkov 2022). The data were collected by a multitude of the WSR-88D radars in 13 continental and 10 maritime MCSs and 11 landfalling hurricanes. The HIWC areas were identified within the examined storms and the corresponding “HIWC statistics” was compared with the “background” one without HIWC. An overarching conclusion of the study is that maritime tropical storms (MCSs and hurricanes) are characterized by smaller size ice in higher concentration compared to the continental MCSs. High ice water content in the HIWC areas is primarily caused by a strong jump in a number concentration of ice particles rather than the increase of their size compared to the “background” environment. This may point to the homogeneous nucleation of excessive amounts of supercooled droplets and / or secondary ice production as the possible origins of HIWC. Such a climatology provides a good observational reference for the modelers to evaluate the performance of their models. As an example, the in-depth analysis of the 20 May 2011 MC3E case shows that the advanced cloud models developed in the course of this study still tend to underestimate the number concentration of ice in the HIWC areas although they succeed in reproducing realistically looking vertical profiles of IWC and Nt. The results of the project research are summarized in 13 journal papers.

54 ENVIRONMENTAL SCIENCES↗

Cyclic oxidation behavior of beta+gamma overlay coatings on gamma and gamma+gamma-prime alloys

Detailed experimental studies of the cyclic oxidation behavior of low-pressure plasma sprayed beta+gamma coasting on gamma-phase Ni-Cr-Al alloys have shown the correlation of weight change, oxide type, and Cr and Al concentration-distance profiles as a function of oxidation time. Of special interest was the transition to breakway oxidation due to the loss of the Al flux to the oxide and the failure of the coated alloy to form an Al2O3-rich oxide scale. The experimental results on beta+gamma/gamma coating systems were used as the basis of a numerical model (ternary, semi-infinite, finite-difference analysis) which accurately predicted changes in Cr and Al concentration-distance profiles. The model was used to study parameters critical to enhancing the life of coatings which fail by a combination of Al loss in forming the oxide scale and Al loss via interdiffusion with the substrate alloy. Comparisons of beta+gamma/gamma coating behavior are made to the oxidation of coated gamma+gamma-prime substrates, both ternary Ni-Cr-Al alloys and Mar-M 247-type alloys.

Nesbitt, J. A.↗

A modeling study of the time-averaged electric currents in the vicinity of isolated thunderstorms

A thorough examination of the results of a time-dependent computer model of a dipole thunderstorm revealed that there are numerous similarities between the time-averaged electrical properties and the steady-state properties of an active thunderstorm. Thus, the electrical behavior of the atmosphere in the vicinity of a thunderstorm can be determined with a formulation similar to what was first described by Holzer and Saxon (1952). From the Maxwell continuity equation of electric current, a simple analytical equation was derived that expresses a thunderstorm's average current contribution to the global electric circuit in terms of the generator current within the thundercloud, the intracloud lightning current, the cloud-to-ground lightning current, the altitudes of the charge centers, and the conductivity profile of the atmosphere. This equation was found to be nearly as accurate as the more computationally expensive numerical model, even when it is applied to a thunderstorm with a reduced conductivity thundercloud, a time-varying generator current, a varying flash rate, and a changing lightning mix.

Driscoll, Kevin T.↗

A Methodology for Estimating the Energy and Moisture Budget of the Convective Boundary Layer Using Continuous Ground-Based Infrared Spectrometer Observations

Abstract Land–atmosphere interactions play a critical role in both the atmospheric water and energy cycles. Changes in soil moisture and vegetation alter the partitioning of surface water and energy fluxes, influencing diurnal evolution of the planetary boundary layer (PBL). The mixing-diagram framework has proven useful in understanding the evolution of the heat and moisture budget within the convective boundary layer (CBL). We demonstrate that observations from the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains site provide all of the needed inputs needed for the mixing-diagram framework, allowing us to quantify the impact from the surface fluxes, advection, radiative heating, encroachment, and entrainment on the evolution of the CBL. Profiles of temperature and humidity retrieved from the ground-based infrared spectrometer [Atmospheric Emitted Radiance Interferometer (AERI)] are a critical component in this analysis. Large-eddy simulation results demonstrate that mean mixed-layer values derived are shown to be critical to close the energy and moisture budgets. A novel approach demonstrated here is the use of network of AERIs and Doppler lidars to quantify the advective fluxes of heat and moisture. The framework enables the estimation of the entrainment fluxes as a residual, providing a way to observe the entrainment fluxes without using multiple lidar systems. The high temporal resolution of the AERI observations enables the morning, midday, and afternoon evolution of the CBL to be quantified. This work provides a new way to use observations in this framework to evaluate weather and climate models. Significance Statement The energy and moisture budget of the planetary boundary layer (PBL) is influenced by multiple sources, and accurately representing this evolution in numerical models is critical for weather forecasts and climate predictions. The mixing-diagram approach, driven by profiling observations as illustrated here, provides a powerful way to quantify the contributions from each of these sources. In particular, the energy and moisture mixed into the PBL from above the PBL can be determined accurately from ground-based remote sensors using this approach.

54 ENVIRONMENTAL SCIENCES↗

Convolutional Neural Network for Transition Modeling Based on Linear Stability Theory

Transition prediction is an important aspect of aerodynamic design because of its impact on skin friction and potential coupling with flow separation characteristics. Traditionally, the modeling of transition has relied on correlation-based empirical formulas based on integral quantities such as the shape factor of the boundary layer. However, in many applications of computational fluid dynamics, the shape factor is not straightforwardly available or not well-defined. We propose using the complete velocity profile along with other quantities (e.g., frequency, Reynolds number) to predict the perturbation amplification factor. While this can be achieved with regression models based on a classical fully connected neural network, such a model can be computationally more demanding. We propose a novel convolutional neural network inspired by the underlying physics as described by the stability equations. Specifically, convolutional layers are first used to extract integral quantities from the velocity profiles, and then fully connected layers are used to map the extracted integral quantities, along with frequency and Reynolds number, to the output (amplification ratio). Numerical tests on classical boundary layers clearly demonstrate the merits of the proposed method. More importantly, we demonstrate that, for Tollmien-Schlichting instabilities in two-dimensional, low-speed boundary layers, the proposed network encodes information in the boundary layer profiles into an integral quantity that is strongly correlated to a well-known, physically defined parameter – the shape factor.

Laminar-turbulent transition↗

A computational study of turbulent separated flow over a wall-mounted cube at two different Reynolds numbers and incoming velocity profiles

Simulations of the canonical wall-mounted unit cube subjected to two distinct incident velocity profiles is performed at Reynolds numbers (Re) where the bulk flow characteristics are known to become relatively Re-insensitive using both direct numerical simulation (DNS) and wall-resolved large eddy simulation (WRLES). The aim of this work is to highlight the sensitivity of such bluff-body flows to mean shear and to provide a representative set of data for such flow scenarios where common turbulence modeling techniques often fail. Simple boundary conditions are selected to allow for easy comparisons for model development purposes and avoid confounding effects from freestream turbulence. In addition to mean velocity profiles, select turbulence statistics are presented in detail. We find that both mean velocity profiles and mean wake turbulence are Re-insensitive over the range examined here.

42 ENGINEERING↗

Human Exploration Missions Study: Space Surveillance Telescope Transfer to and Station at a Halo Orbit at the Earth-Sun Libration Point L2

This study was undertaken to determine mission profile and delta velocity requirements to place a telescope at the Earth-Sun libration point L2. The program, Integrated Mission Program (IMP), was selected to be used in the investigation. A description of IMP and its capabilities may be found in the Addenda. The Addenda also contains the libration halo equations, constants and other parameters. Comments regarding the chaotic nature of numerical integration near the libration points are also attached in the Addenda. A basic two stage S/C with a simple mission profile was selected. This profile is shown.

Dauro, Vincent A., Sr.↗

Recirculating and Sub-Atmospheric Rapid Cycle Amine Swing Bed Testing at Various Metabolic Profiles, Temperatures, and Half-Cycle Times

Sixty years since the first spacewalk, NASA’s Extravehicular Mobility Unit (xEMU) exploration space suit has undergone numerous design improvements and iterations. As commercial manufacturers develop next-generation systems to support increasingly complex missions, including Moon-to-Mars architecture prototyping and the design of the Martian Portable Life Support System (PLSS), material selection and individual subunit testing under realistic metabolic and pressure conditions are paramount. To determine whether a reusable carbon dioxide (CO 2 ) scrubber design and adsorbent can be utilized in new suit configurations, an all-encompassing test rig is required to compare existing and newly developed technologies. The test rig must be capable of sub-atmospheric testing, calibrated humid CO 2 dosing, and facilitation of the evaluation of sorbents across a wide range of operational requirements. Recently, XploSafe has developed an Extravehicular Mobility Unit (xEMU) testing rig specifically designed to explore CO 2 and humidity control materials for spacesuits. This closed-loop recirculating swing bed (Rapid Cycle Amine Test Rig) was utilized to compare adsorbent utility across a wide range of xEMU operating conditions. Following successful testing with Xplo-SA9T over multiple simulated eight-hour Extravehicular Activity (EVA) tests, improvements to the test stand were considered. Regeneration vacuum supply within the test rig was updated to more closely resemble space vacuum, theoretically increasing half-cycle timing. An improved water vapor delivery system was incorporated to dose a wider calibrated range of humidity over both short and long-duration EVAs. A range of metabolic profiles between 800 and 3000 BTU/h were tested, and the corresponding average CO 2 steady-states were compared. The effect of dosing temperature at the adsorbent bed inlet between 10 and 40 ̊C and the corresponding CO 2 removal efficiency was also evaluated.

John R Tidwell↗

Evaporation of LOX under supercritical and subcritical conditions

The evaporation of LOX under supercritical and subcritical conditions was studied experimentally and theoretically. In experiments, the evaporation rate and surface temperature were measured for LOX strand vaporizing in helium environments at pressures ranging from 5 to 68 atmospheres. Gas sampling and chromatography analysis were also employed to profile the gas composition above the LOX surface for the purpose of model validation. A comprehensive theoretical model was formulated and solved numerically to simulate the evaporation process of LOX at high pressures. The model was based on the conservation equations of mass, momentum, energy, and species concentrations for a multicomponent system, with consideration of gravitational body force, solubility of ambient gases in liquid, and variable thermophysical properties. Good agreement between predictions and measured oxygen mole fraction profiles was obtained. The effect of pressure on the distribution of the Lewis number, as well as the effect of variable diffusion coefficient, were further examined to elucidate the high-pressure transport behavior exhibited in the LOX vaporization process.

Yang, A. S.↗

Numerical Investigation of Wall-Cooling Effect on Aero-Optical Distortions for Hypersonic Boundary Layer

Compressible wall-modeled large-eddy simulations of Mach 8 turbulent boundary-layer flows over a flat plate were carried out for the conditions of the hypersonic wind tunnel at Sandia National Laboratories. The simulations provide new insight into the effect of wall cooling on the aero-optical path distortions for hypersonic turbulent boundary-layer flows. Four different wall-to-recovery temperature ratios, 0.3, 0.48, 0.71, and 0.89, are considered. Despite the much lower grid resolution, the mean velocity, temperature, and resolved Reynolds stress profiles from the simulation for a temperature ratio of 0.48 are in good agreement with those from a reference direct numerical simulation. The normalized root-mean-square optical path difference obtained from the present simulations is compared with that from reference direct numerical simulations, Sandia experiments, as well as predictions obtained with a semi-analytical model by Notre Dame University. Here the present analysis focuses on the effect of wall cooling on the wall-normal density correlations, on key underlying assumptions of the aforementioned model such as the strong Reynolds analogy, and on the elevation angle effect on the optical path difference. Wall cooling is found to increase the velocity fluctuations and decrease the density fluctuations, resulting in an overall reduction of the normalized optical path distortion. Compared to the simulations, the basic strong Reynolds analogy overpredicts the temperature fluctuations for cooled walls. Also different from the strong Reynolds analogy, the velocity and temperature fluctuations are not perfectly anticorrelated. Finally, as the wall temperature is raised, the density correlation length, away from the wall but inside the boundary layer, increases significantly for beam paths tilted in the downstream direction.

42 ENGINEERING↗

Pair plasma instability in homogeneous magnetic guide fields

Pair plasmas, collections of both matter and antimatter particles of equal mass, represent a paradigm for the study of basic plasma science, and many open questions exist regarding these unique systems. They are found in many astrophysical settings, such as gamma-ray bursts, and have recently also been produced in carefully designed laboratory experiments. A central research topic in plasma physics is instability; however, unlike their more common ion–electron siblings, pair plasmas are generally thought to be stable to cross field pressure gradients in homogeneous magnetic fields. It is shown here by means of kinetic full-f simulations that, when a pressure gradient is first established, the Gradient-driven Drift Coupling mode is destabilized and becomes turbulent. Force balance is eventually achieved by a combination of flattened pressure profiles due to turbulent transport and establishment of a magnetic field gradient, saturating the growth. During the unstable phase, key physics can be captured by a δf gyrokinetic description, where it is shown analytically and numerically that parallel particle motion results in a coupling of all electromagnetic field components. A fluid model derived therefrom accurately predicts linear eigenmodes and is used to resolve global profile effects. For laser-based electron–positron plasma experiments, prompt instability is predicted with growth times much shorter than plasma lifetimes. Similarly, growth rates are calculated for the planned APEX experiment as well as gamma-ray burst scenarios, suggesting that the instability may contribute to the early evolution of these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Large-eddy simulation of a backward facing step flow using a least-squares spectral element method

We report preliminary results obtained from the large eddy simulation of a backward facing step at a Reynolds number of 5100. The numerical platform is based on a high order Legendre spectral element spatial discretization and a least squares time integration scheme. A non-reflective outflow boundary condition is in place to minimize the effect of downstream influence. Smagorinsky model with Van Driest near wall damping is used for sub-grid scale modeling. Comparisons of mean velocity profiles and wall pressure show good agreement with benchmark data. More studies are needed to evaluate the sensitivity of this method on numerical parameters before it is applied to complex engineering problems.

Chan, Daniel C.↗

Shock compression of silver to 300 GPa: Wave profile measurements and melting transition

Recent in-situ x-ray diffraction (XRD) measurements on laser-shock compressed Ag foils demonstrated a face-centered-cubic to body-centered-cubic transformation at ~150 GPa and melting between 172 and 197 GPa [Phys. Rev. Lett. 124, 235701 (2020)]. As a complement to the XRD work, we conducted plate impact experiments to obtain shock velocity and wave profile measurements on thicker Ag samples shock-compressed to peak stresses between 30 and 300 GPa. The shock velocity–particle velocity results were fitted very well by a linear relation over the entire stress range, providing an accurate determination of the Ag Hugoniot (locus of the stress-volume states achieved under shock compression). For peak stresses below 187 GPa and above 210 GPa—corresponding to the solid and liquid phases, respectively—the wave profiles show clean single waves. No wave profile features related to the fcc-bcc transformation at ~150 GPa were observed, implying minimal volume change for the transformation. For stresses between 187 and 210 GPa, an initial jump was followed by a time-dependent increase in the particle velocity (20–80 ns risetime) to the peak state—corresponding to the solid-liquid mixed phase response. Unlike the solid and liquid response, the mixed-phase response cannot be readily analyzed analytically. Instead, numerical simulations incorporating an accurate multiphase equation of state for Ag—not currently available—are required to analyze the wave profiles measured at 187–210 GPa stresses. Furthermore, the present work shows the potential for using wave profile measurements to examine the melting transition under shock compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Directed energy deposition of radially functionally graded nickel based superalloys IN625 and IN738: Microstructure and mechanical properties

This study investigates the microstructure and mechanical properties of radially functionally graded (RFG) materials fabricated from nickel-based superalloys IN625 and IN738 using directed energy deposition (DED). RFG cylindrical deposits with radial gradients of IN625-IN738 and IN738-IN625 were successfully fabricated and confirmed to be crack-free through tailored composition grading, scan strategy and process parameters that minimized thermal gradients and solidification velocities. Microstructural analysis revealed distinct variations across the gradient, characterized by differences in dendrite arm spacing, microsegregation, and the formation of distinct γ' precipitates and carbide/boride phases (MC, M 6 C, M 2 C, and MB 2 ) in as-deposited (AD) conditions, which evolved into M 23 C 6 and M 2 B upon subsequent heat treatment (HT). EBSD showed predominant columnar grains with a [100] texture, while microhardness profiles demonstrated consistent variation across the gradients. Mechanical testing highlighted the synergy of the graded structure, combining IN738's strength with IN625's ductility, resulting in enhanced properties such as a yield strength of 1131 MPa, ultimate tensile strength of 1382 MPa, and elongation of 12 % for HT IN625-IN738 grading. Numerical simulations employing CALPHAD thermodynamic calculations and finite element thermal modeling are used to establish phase evolution and microstructural feature correlation with thermal profiles, aiding in the optimization of process-structure-property relationship.

Additive manufacturing↗

Utah FORGE: LBNL Reports on VEMP Electromagnetic Data Collection and Processing - 2024

This archive contains reports related to Vertical Electromagnetic Profiling (VEMP) tool data collection and processing at Utah FORGE in 2024. The first report describes LBNL's effort to collect electromagnetic geophysical data with the tool in well 78-32B and a downhole electrode in well 16A. The second report describes the final data acquisition and processing of the VEMP electromagnetic data collected at the Utah FORGE site in May of 2024. Also included are a noise analysis as well as a comparison of the data to numerical models. This was originally presented as a paper at the 2025 Stanford Geothermal Workshop.

15 GEOTHERMAL ENERGY↗

Three-dimensional full-field velocity measurements in shock compression experiments using stereo digital image correlation

Shock compression plate impact experiments conventionally rely on point-wise velocimetry measurements based on laser-based interferometric techniques. This study presents an experimental methodology to measure the free surface full-field particle velocity in shock compression experiments using high-speed imaging and three-dimensional (3D) digital image correlation (DIC). The experimental setup has a temporal resolution of 100 ns with a spatial resolution varying from 90 to 200 μm/pixel. Experiments were conducted under three different plate impact configurations to measure spatially resolved free surface velocity and validate the experimental technique. First, a normal impact experiment was conducted on polycarbonate to measure the macroscopic full-field normal free surface velocity. Second, an isentropic compression experiment on Y-cut quartz–tungsten carbide assembly is performed to measure the particle velocity for experiments involving ramp compression waves. To explore the capability of the technique in multiaxial loading conditions, a pressure shear plate impact experiment was conducted to measure both the normal and transverse free surface velocities under combined normal and shear loading. The velocities measured in the experiments using digital image correlation are validated against previous data obtained from laser interferometry. Numerical simulations were also performed using established material models to compare and validate the experimental velocity profiles for these different impact configurations. The novel ability of the employed experimental setup to measure full-field free surface velocities with high spatial resolutions in shock compression experiments is demonstrated for the first time in this work.

47 OTHER INSTRUMENTATION↗

Ideal internal kink stability in presence of plasma flow and neoclassical toroidal viscosity due to energetic particles

Abstract The influence of energetic particles (EPs) on the ideal internal kink mode, in rotating tokamak plasmas, is numerically investigated by simultaneously solving MHD-kinetic hybrid equations together with a toroidal momentum balance equation utilizing the MARS-Q code (Liu et al 2013 Phys. Plasmas 20 042503). The neoclassical toroidal viscous (NTV) torque, induced by precessional drift resonances of trapped energetic particles, acts as the momentum sink term to damp the plasma flow. Quasi-linear initial value simulations show local reduction of the flow amplitude and enhancement of the flow shear near the q = 1 rational surface ( q is the safety factor) due to EP induced NTV. Both effects in turn destabilize the internal kink mode. These numerical findings are robust against the initial linear stability of internal kink, the initial plasma flow profile, as well as the equilibrium distribution model for EPs.

Physics↗