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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Laser Beam Propagation Through Inhomogeneous Media with Shock-Like Profiles: Modeling and Computing

Wave propagation in inhomogeneous media has been studied for such diverse applications as propagation of radiowaves in atmosphere, light propagation through thin films and in inhomogeneous waveguides, flow visualization, and others. In recent years an increased interest has been developed in wave propagation through shocks in supersonic flows. Results of experiments conducted in the past few years has shown such interesting phenomena as a laser beam splitting and spreading. The paper describes a model constructed to propagate a laser beam through shock-like inhomogeneous media. Numerical techniques are presented to compute the beam through such media. The results of computation are presented, discussed, and compared with experimental data.

Adamovsky, Grigory↗
NASA NTRS1999

Real Time Optima Tracking Using Harvesting Models of the Genetic Algorithm

Tracking optima in real time propulsion control, particularly for non-stationary optimization problems is a challenging task. Several approaches have been put forward for such a study including the numerical method called the genetic algorithm. In brief, this approach is built upon Darwinian-style competition between numerical alternatives displayed in the form of binary strings, or by analogy to 'pseudogenes'. Breeding of improved solution is an often cited parallel to natural selection in.evolutionary or soft computing. In this report we present our results of applying a novel model of a genetic algorithm for tracking optima in propulsion engineering and in real time control. We specialize the algorithm to mission profiling and planning optimizations, both to select reduced propulsion needs through trajectory planning and to explore time or fuel conservation strategies.

Baskaran, Subbiah↗
NASA NTRS2011

Base Heating Sensitivity Study for a 4-Cluster Rocket Motor Configuration in Supersonic Freestream

In support of launch vehicle base heating and pressure prediction efforts using the Loci-CHEM Navier-Stokes computational fluid dynamics solver, 35 numerical simulations of the NASA TND-1093 wind tunnel test have been modeled and analyzed. This test article is composed of four JP-4/LOX 500 lbf rocket motors exhausting into a Mach 2 - 3.5 wind tunnel at various ambient pressure conditions. These water-cooled motors are attached to a base plate of a standard missile forebody. We explore the base heating profiles for fully coupled finite-rate chemistry simulations, one-way coupled RAMP (Reacting And Multiphase Program using Method of Characteristics)-BLIMPJ (Boundary Layer Integral Matrix Program - Jet Version) derived solutions and variable and constant specific heat ratio frozen flow simulations. Variations in turbulence models, temperature boundary conditions and thermodynamic properties of the plume have been investigated at two ambient pressure conditions: 255 lb/sq ft (simulated low altitude) and 35 lb/sq ft (simulated high altitude). It is observed that the convective base heat flux and base temperature are most sensitive to the nozzle inner wall thermal boundary layer profile which is dependent on the wall temperature, boundary layer s specific energy and chemical reactions. Recovery shock dynamics and afterburning significantly influences convective base heating. Turbulence models and external nozzle wall thermal boundary layer profiles show less sensitivity to base heating characteristics. Base heating rates are validated for the highest fidelity solutions which show an agreement within +/-10% with respect to test data.

Mehta, Manish↗
NASA NTRS2016

Observations and Simulations of the Na I D-1 Line profiles in an M-Class Solar Flare

We study the temporal evolution of the Na I D1 line profiles in the M3.9 flare SOL2014-06-11T21:03UT, using observations at high spectral resolution obtained with the Interferometric Bidimensional Spectrometer instrumentation the Dunn Solar Telescope combined with radiative hydrodynamic simulations. Our results show a significant increase in the intensities of the line core and wings during the flare. The analysis of the line profiles from the flare ribbons reveals that the Na I D1 line has a central reversal with excess emission in the blue wing (blue asymmetry).We combine RADYN and RH simulations to synthesize Na I D1 line profiles of the flaring atmosphere and find good agreement with the observations. Heating with a beam of electrons modifies the radiation field in the flaring atmosphere and excites electrons from the ground state 3s 2S to the first excited state 3p 2P, which in turn modifies the relative population of the two states. The change in temperature and the population density of the energy states make the sodium line profile revert from absorption into emission. Furthermore, the rapid changes in temperature break the pressure balance between the different layers of the lower atmosphere, generating upflow/downflow patterns. Analysis of the simulated spectra reveals that the asymmetries of the Na I D1 flare profile are produced by the velocity gradients in the lower solar atmosphere.

methods: numerical↗
DOE OSTI2025

EVSE Characterization: V2G EVSE Comparison

As part of the U.S. Department of Energy EVs@Scale consortium Next-Generation Profiles project, results and analysis from the characterization of high-power conductive and wireless charging infrastructure are presented. This characterization is conducted over a wide range of direct current (DC) current and DC voltage operation for nominal test conditions and off-nominal test conditions. Test plans and procedures were developed to define the test configurations and requirements, measurement parameters, and test procedures used throughout testing. Results from a 2024 study conducted on electric vehicle supply equipment (EVSE) characterization by the Idaho National Laboratory (INL) include two bi-directional vehicle-to-grid (V2G) capable EVSEs. These EVSE are referred to as V2G-EVSE9 and V2G-EVSE10. Laboratory testing is conducted at nominal test conditions to characterize the power transfer capabilities, efficiency, power factor, and other power quality metrics of the two DC EVSEs capable of V2G bi-directional power transfer. Results from testing show the performance is consistent for V2G-EVSE9 and V2G-EVSE10 when comparing charging to discharging performance, except for V2G-EVSE10 for power transfer when operating above 70% of the rated DC current. V2G-EVSE10 efficiency is >98% while charging and <91% while discharging at the same operating conditions, near maximum-rated current, at 300VDC. In contrast, V2G-EVSE9 results are consistent for charging and discharging. This EVSE is nearly 96% efficient while charging or discharging when operating over 50% of rated AC power. V2G EVSE performance is also characterized during off-nominal AC grid conditions involving AC voltage deviation (426 VAC to 518 VAC), AC frequency deviation of +2% (58.8 Hz to 61.2 Hz), and AC voltage harmonics injection. Many test conditions have little-to-no impact on performance characteristics of the two EVSEs; however, there are a few notable findings with significant power transfer capability impacts. AC voltage harmonics injection resulted in negative impacts on power quality attributes for both EVSEs, but with no impact on power transfer capability. Off-nominal AC voltage and frequency conditions resulted in unstable or lack of power transfer capability for both EVSEs. V2G-EVSE9 is unable to transfer power when AC voltage is >300V L-N. V2G-EVSE10 is unable to transfer power when AC frequency deviation exceeds +0.8%. V2G energy management system transient response and latency are quantified during laboratory testing. V2G-EVSE9 and V2G-EVSE10 utilize cloud-based V2G energy management systems that command the power transfer level between the EVSE and EV. The latency and response characteristics of the entire systems (web-based user interface, V2G energy management system, cellular communications, and EVSE response) are quantified through laboratory testing for V2G-EVSE9 and V2G-EVSE10. V2G-EVSE9 latency ranges from 0.8 to 1.8 seconds, whereas V2G-EVSE10 latency ranges from 3.4 to 8.8 seconds. The ramp rate to a change in power transfer request also differs between the two EVSEs. V2G-EVSE10 ramp rate ranges from 50% to -250% of rated AC power per second, whereas V2G-EVSE9 rate ranges from 95% to -95% of rated AC power per second. At the highest rate of change in power transfer, V2G-EVSE10 can change from full charge power to full discharge power in less than one second. The V2G EVSE characterization presented in this report provides valuable insights and results for use by numerous entities. This includes modeling and simulation organizations, decision makers, fleet planning, industry stakeholders, and many others involved with the development and deployment of electrified transportation technologies. Additional high-power DC chargers, bidirectional chargers, and inductive power transfer EVSE characterization results are anticipated from additional EVSE brands and models, which will be detailed in future publications in support of the U.S. Department of Energy EVs@Scale consortium Next-Gen Profiles project.

25 ENERGY STORAGE↗
NASA NTRS2020

Clouds and Convective Self-Aggregation in a Multi-Model Ensemble of Radiative-Convective Equilibrium Simulations

The Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative-convective56equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climatecsensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique amongst intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud-resolving models (CRMs), large eddy simulations (LES), and global cloud-resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self-aggregation in large domains and agree that self-aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self-aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than un-aggregated simulations. Plain Language Summary This study investigates tropical clouds and climate using results from more than 30 different numerical models set up in a simplified framework. The dataset of model simulations is unique in that it includes a wide range of model types configured in a consistent manner. We address some of the biggest open questions in climate science, including how cloud properties change with warming and the role that the tendency of clouds to form clusters plays in determining the average climate and how climate changes. While there are large differences in how the different models simulate average temperature, humidity, and cloudiness, in a majority of models, the amount of high clouds decreases as climate warms. Nearly all models simulate a tendency for clouds to cluster together. There is agreement that when the clouds are clustered, the atmosphere is drier with fewer clouds overall. We don’t find a conclusive result for how cloud clustering changes as the climate warms.

Allison A. Wing↗
DOE OSTI2023

Machine learning methods for particle stress development in suspension Poiseuille flows

Numerical simulations are used to study the dynamics of a developing suspension Poiseuille flow with monodispersed and bidispersed neutrally buoyant particles in a planar channel, and machine learning is applied to learn the evolving stresses of the developing suspension. The particle stresses and pressure develop on a slower time scale than the volume fraction, indicating that once the particles reach a steady volume fraction profile, they rearrange to minimize the contact pressure on each particle. Here we consider how the stress development leads to particle migration, time scales for stress development, and present a new physics-informed Galerkin neural network that allows for learning the particle stresses when direct measurements are not possible. The particle fluxes are compared with the Suspension Balance Model with good agreement. We show that when stress measurements are possible, the MOR-physics operator learning method can also capture the particle stresses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Collisionless adiabatic afterglow

In this work, we study, by numerical and analytical means, the evolution of a collisionless plasma initiated between absorbing walls. The ensuing flow is described by rarefaction waves that travel inward from the boundaries, interact, and eventually vanish after crossing through, leading up to the asymptotic stage of the decay. Particle simulations indicate that the kinetic evolution strongly resembles one found in isentropic gas dynamics. Namely, a very gradual density profile forms in the expanding central region where the rarefaction waves interact, with an accompanying linear velocity profile. Asymptotically, the density falls off as 1/t. The density and the flux at the boundary show little variation over the period when rarefaction waves still exist. Plasma potential, on the other hand, drops quite rapidly (on the underlying ion-acoustic timescale) to less than initial electron temperature Te when over 70% of the particles still remain in the system. This is due to electron kinetics being governed by conservation of adiabatic invariant in a slowly varying potential well. Analytical model of the velocity distribution is presented to explain the simulations. The results have implications for afterglow plasmas used in material processing and also for ion-extraction devices. One property of potential interest is good uniformity of the decaying plasma that occurs after approximately one ion-acoustic time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Turbulence model reduction by deep learning

A defining problem of turbulence theory is to produce a predictive model for turbulent fluxes. These have profound implications for virtually all aspects of the turbulence dynamics. In magnetic confinement devices, drift-wave turbulence produces anomalous fluxes via cross-correlations between fluctuations. In this work, we introduce an alternative, data-driven method for parametrizing these fluxes. The method uses deep supervised learning to infer a reduced mean-field model from a set of numerical simulations. We apply the method to a simple drift-wave turbulence system and find a significant new effect which couples the particle flux to the local gradient of vorticity. Notably, here, this effect is much stronger than the oft-invoked shear suppression effect. We also recover the result via a simple calculation. The vorticity gradient effect tends to modulate the density profile. In addition, our method recovers a model for spontaneous zonal flow generation by negative viscosity, stabilized by nonlinear and hyperviscous terms. We highlight the important role of symmetry to implementation of the new method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effect of Internal Structure on Warpage in a Large-Scale Additive Manufacturing Process With Bio-Derived Composites

Extrusion-based large-format additive manufacturing (LFAM) often results in unintended deformation and failures due to thermal residual stress between layers. This LFAM involves a hot molten polymer deposited on a previously deposited and cooled layer, generating a temperature mismatch between the layers. The temperature discrepancy causes a thermal contraction mismatch, generating thermal residual stress. Due to the residual stress, printed structures experience warpage and delamination Therefore, understanding the heat distribution behavior is essential to preventing undesired deformation and failures. The goal of this study is to investigate the effect of infill patterns on heat distribution and deformation in an LFAM system. Wood fiber-reinforced polylactic acid was used as a sustainable and low-cost reinforcement material. A numerical model was developed to predict the temperature and deformation field with thermo-mechanical properties, which were measured by a dynamic mechanical analysis and a thermal expansion test. Simulations were performed on a box geometry with three different infill patterns. The simulation model was verified with temperature data gathered from an infrared camera. The results show good agreement between predicted and measured temperature profiles. The developed simulation model was applied to a roof tray with different infill patterns for a case study. We found that the different thermal mass distributions resulting from the various infill patterns affected the heat distribution and deformation. These findings contribute to a better understanding of the thermal and mechanical behavior of LFAM.

Jo, Eonyeon↗
DOE OSTI2021

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗
DOE OSTI2025

Image-Driven Hybrid Structural Analysis Based on Continuum Point Cloud Method with Boundary Capturing Technique

Conventional approaches for the structural health monitoring of infrastructures often rely on physical sensors or targets attached to structural members, which require considerable preparation, maintenance, and operational effort, including continuous on-site adjustments. This paper presents an image-driven hybrid structural analysis technique that combines digital image processing (DIP) and regression analysis with a continuum point cloud method (CPCM) built on a particle-based strong formulation. Polynomial regressions capture the boundary shape change due to the structural loading and precisely identify the edge and corner coordinates of the deformed structure. The captured edge profiles are transformed into essential boundary conditions. This allows the construction of a strongly formulated boundary value problem (BVP), classified as the Dirichlet problem. Capturing boundary conditions from the digital image is novel, although a similar approach was applied to the point cloud data. It was shown that the CPCM is more efficient in this hybrid simulation framework than the weak-form-based numerical schemes. Unlike the finite element method (FEM), it can avoid aligning boundary nodes with regression points. A three-point bending test of a rubber beam was simulated to validate the developed technique. The simulation results were benchmarked against numerical results by ANSYS and various relevant numerical schemes. The technique can effectively solve the Dirichlet-type BVP, yielding accurate deformation, stress, and strain values across the entire problem domain when employing a linear strain model and increasing the number of CPCM nodes. In addition, comparative analysis with conventional displacement tracking techniques verifies the developed technique’s robustness. The proposed technique effectively circumvents the inherent limitations of traditional monitoring methods resulting from the reliance on physical gauges or target markers so that a robust and non-contact solution for remote structural health monitoring in real-scale infrastructures can be provided, even in unfavorable experimental environments.

Chemistry↗
DOE OSTI2025

Irradiation Testing of Ultrasonic Transducers

Ultrasonic technologies offer the potential for high accuracy and resolution in-pile measurement of numerous parameters, including geometry changes, temperature, crack initiation and growth, gas pressure and composition, and microstructural changes. Many Department of Energy-Office of Nuclear Energy (DOE-NE) programs are exploring the use of ultrasonic technologies to provide enhanced sensors for in-pile instrumentation during irradiation testing. For example, the ability of single, small diameter ultrasonic thermometers (UTs) to provide a temperature profile in candidate metallic and oxide fuel would provide much needed data for validating new fuel performance models. Other efforts include an ultrasonic technique to detect morphology changes (such as crack initiation and growth) and acoustic techniques to evaluate fission gas composition and pressure. These efforts are limited by the lack of existing knowledge of ultrasonic transducer material survivability under irradiation conditions. To address this need, the Pennsylvania State University (PSU) was awarded an Advanced Test Reactor National Scientific User Facility (ATR NSUF) project to evaluate promising magnetostrictive and piezoelectric transducer performance in the Massachusetts Institute of Technology Research Reactor (MITR) up to a fast fluence of at least 1021 n/cm2 (E> 0.1 MeV). This test will be an instrumented lead test; and real-time transducer performance data will be collected along with temperature and neutron and gamma flux data. By characterizing magnetostrictive and piezoelectric transducer survivability during irradiation, test results will enable the development of novel radiation tolerant ultrasonic sensors for use in Material and Test Reactors (MTRs). The current work bridges the gap between proven out-of-pile ultrasonic techniques and in-pile deployment of ultrasonic sensors by acquiring the data necessary to demonstrate the performance of ultrasonic transducers

Daw, J.↗
NASA NTRS1985

A comparison of observed and forecast energetics over North America

The observed kinetic energy balance is calculated over North America and compared with that computed from forecast fields for the 13-15 January 1979 cyclone. The FGGE upper-air rawinsonde network serves as the observational database while the forecast energetics are derived from a numerical integration with the GLAS fourth-order general circulation model initialized at 00 GMT 13 January. Maps of the observed and predicted kinetic energy and eddy conversion are in good qualitative agreement, although the model eddy conversion tends to be 2 to 3 times stronger than the observed values. Both the forecast and observations exhibit the lower and upper tropospheric maxima in vertical profiles of kinetic energy generation and dissipation typically found in cyclonic disturbances. An interesting time lag is noted in the observational analysis with the maximum observed kinetic energy occurring 12 h later than the maximum eddy conversion over the same region.

Baker, W. E.↗
NASA NTRS1992

BLSTA: A boundary layer code for stability analysis

A computer program is developed to solve the compressible, laminar boundary-layer equations for two-dimensional flow, axisymmetric flow, and quasi-three-dimensional flows including the flow along the plane of symmetry, flow along the leading-edge attachment line, and swept-wing flows with a conical flow approximation. The finite-difference numerical procedure used to solve the governing equations is second-order accurate. The flow over a wide range of speed, from subsonic to hypersonic speed with perfect gas assumption, can be calculated. Various wall boundary conditions, such as wall suction or blowing and hot or cold walls, can be applied. The results indicate that this boundary-layer code gives velocity and temperature profiles which are accurate, smooth, and continuous through the first and second normal derivatives. The code presented herein can be coupled with a stability analysis code and used to predict the onset of the boundary-layer transition which enables the assessment of the laminar flow control techniques. A user's manual is also included.

Wie, Yong-Sun↗
NASA NTRS2001

Quasi 1-D Analysis of a Circular, Compressible, Turbulent Jet Laden with Water Droplets

Recent experimental studies indicate that presence of small amount of liquid droplets reduces the Overall Sound Pressure Level (OASPL) of a jet. Present study is aimed at numerically investigating the effect of liquid particles on the overall flow quantities of a heated, compressible round jet. The jet is assumed perfectly expanded. A quasi-1D model was developed for this purpose which uses area-averaged quantities that satisfy integral conservation equations. Special attention is given to represent the early development region since it is acoustically important. Approximate velocity and temperature profiles were assumed in this region to evaluate entrainment rate. Experimental correlations were used to obtain spreading rate of shear layer. The base flow thus obtained is then laden with water droplets at the exit of the nozzle. Mass, momentum and energy coupling between the two phases is represented using empirical relations. Droplet size and mass loading are varied to observe their effect on flow variables.

Source record↗
NASA NTRS2000

Calcaneal loading during walking and running

PURPOSE: This study of the foot uses experimentally measured kinematic and kinetic data with a numerical model to evaluate in vivo calcaneal stresses during walking and running. METHODS: External ground reaction forces (GRF) and kinematic data were measured during walking and running using cineradiography and force plate measurements. A contact-coupled finite element model of the foot was developed to assess the forces acting on the calcaneus during gait. RESULTS: We found that the calculated force-time profiles of the joint contact, ligament, and Achilles tendon forces varied with the time-history curve of the moment about the ankle joint. The model predicted peak talocalcaneal and calcaneocuboid joint loads of 5.4 and 4.2 body weights (BW) during walking and 11.1 and 7.9 BW during running. The maximum predicted Achilles tendon forces were 3.9 and 7.7 BW for walking and running. CONCLUSIONS: Large magnitude forces and calcaneal stresses are generated late in the stance phase, with maximum loads occurring at approximately 70% of the stance phase during walking and at approximately 60% of the stance phase during running, for the gait velocities analyzed. The trajectories of the principal stresses, during both walking and running, corresponded to each other and qualitatively to the calcaneal trabecular architecture.

NASA Center ARC↗
NASA NTRS2005

Progress Towards Computational Method for Circulation Control Airfoils

The compressible Reynolds-averaged Navier-Stokes equations are solved for circulation control airfoil flows. Numerical solutions are computed with both structured and unstructured grid solvers. Several turbulence models are considered, including the Spalart-Allmaras model with and without curvature corrections, the shear stress transport model of Menter, and the k-enstrophy model. Circulation control flows with jet momentum coefficients of 0.03, 0.10, and 0.226 are considered. Comparisons are made between computed and experimental pressure distributions, velocity profiles, Reynolds stress profiles, and streamline patterns. Including curvature effects yields the closest agreement with the measured data.

Swanson, R. C.↗