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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

SUPIN: A Tool for the Aerodynamic Design and Analysis of Supersonic Inlets

The Supersonic Inlet Design and Analysis Tool (SUPIN) is a computational program that performs geometric modeling and aerodynamic design and analysis of a set of supersonic inlets whose types include the axisymmetric pitot, three-dimensional pitot, axisymmetric, two-dimensional, two-dimensional bifurcated-duct, and streamline-traced inlets. The inlets are modeled by dividing the inlet flowpath into components that start at the freestream and ends at the engine face. The primary components of the inlet include the external supersonic diffuser, throat section, and subsonic diffuser. Each component is characterized by a set of inputs that include geometric and aerodynamic factors. Examples of the geometric factors are angles, lengths, area ratios, and coordinates that describe the geometry of the components. Examples of the aerodynamic factors are Mach numbers, flow ratios, and pressure ratios that set constraints or desired aerodynamic properties of the component. These input factors are specified within a text-based input data file that is read into SUPIN. SUPIN is a Fortran 95 program that performs the inlet design and analysis as a batch or non-interactive process. SUPIN writes output data files describing the inlet geometry and aerodynamic performance. SUPIN uses analytic, empirical, and numerical methods to perform the inlet design and analysis. The inlet geometry is constructed using planar and three-dimensional geometry constructs and is represented by coordinates, angles, areas, profiles, and surfaces. Surfaces are output in the form of a Plot3D surface grid file. SUPIN does have some capability to automatically generate planar and three-dimensional, multi-block, structured grids for computational fluid dynamics (CFD) analysis. The aerodynamic analysis performed by SUPIN computes quasi-one-dimensional properties, such as Mach numbers, pressures, and temperatures, through the inlet flowpath. Planar aerodynamic properties, such as those for shock and Mach waves, are computed for the supersonic compression of the external supersonic diffuser and throat section. The inlet performance is characterized by the inlet flow rates, total pressure recovery, and inlet drag. The primary intent of this document is to describe the usage of SUPIN. The methods used within SUPIN are documented and may serve as a useful reference for understanding the aerodynamic design and analysis of supersonic inlets.

supersonic inlets↗
NASA NTRS1989

Linear instability of supersonic plane wakes

In this paper we present a theoretical and numerical study of the growth of linear disturbances in the high-Reynolds-number and laminar compressible wake behind a flat plate which is aligned with a uniform stream. No ad hoc assumptions are made as to the nature of the undisturbed flow (in contrast to previous investigations) but instead the theory is developed rationally by use of proper wake-profiles which satisfy the steady equations of motion. The initial growth of near wake perturbation is governed by the compressible Rayleigh equation which is studied analytically for long- and short-waves. These solutions emphasize the asymptotic structures involved and provide a rational basis for a nonlinear development. The evolution of arbitrary wavelength perturbations is addressed numerically and spatial stability solutions are presented that account for the relative importance of the different physical mechanisms present, such as three-dimensionality, increasing Mach numbers enough (subsonic) Mach numbers, there exists a region of absolute instability very close to the trailing-edge with the majority of the wake being convectively unstable. At higher Mach numbers (but still not large-hypersonic) the absolute instability region seems to disappear and the maximum available growth-rates decrease considerably. Three-dimensional perturbations provide the highest spatial growth-rates.

Papageorgiou, D. T.↗
DOE OSTI2024

Measuring thermal profiles in high explosives using neural networks

We present a new method for calculating the temperature profile of high explosive (HE) material using a Convolutional Neural Network (CNN). To train/test the CNN, we have developed a hybrid experiment/simulation method for collecting acoustic and temperature data. We experimentally heat cylindrical containers of HE material until detonation/deflagration, where we continuously measure the acoustic bursts through the HE using multiple acoustic transducers lined around the exterior container circumference. However, measuring the temperature profile in the HE in an experiment would require inserting a large number of thermal probes, which would disrupt the heating process. Thus, we use two thermal probes, one at the HE center and one at the wall. We then use numerical simulation of the heating process to calculate the temperature distribution and correct the simulated temperatures based on the experimental center and wall temperatures. We calculate temperature errors on the order of 15 °C, which is ∼12% of the range of temperatures in the experiment. We also investigate how the algorithm’s accuracy is affected by the number of acoustic receivers used to collect each measurement and the resolution of the temperature prediction. This work provides a means of assessing the safety status of HE material, which cannot be achieved using existing temperature measurement methods. In addition, it has implications for a range of other applications where internal temperature profile measurements would provide critical information. These applications include detecting chemical reactions, observing thermodynamic processes such as combustion, monitoring metal or plastic casting, determining the energy density in thermal storage capsules, and identifying abnormal battery operations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Measuring thermal profiles in high explosives using neural networks

We present a new method for calculating the temperature profile of high explosive (HE) material using a Convolutional Neural Network (CNN). To train/test the CNN, we have developed a hybrid experiment/simulation method for collecting acoustic and temperature data. We experimentally heat cylindrical containers of HE material until detonation/deflagration, where we continuously measure the acoustic bursts through the HE using multiple acoustic transducers lined around the exterior container circumference. However, measuring the temperature profile in the HE in an experiment would require inserting a large number of thermal probes, which would disrupt the heating process. Thus, we use two thermal probes, one at the HE center and one at the wall. We then use numerical simulation of the heating process to calculate the temperature distribution and correct the simulated temperatures based on the experimental center and wall temperatures. We calculate temperature errors on the order of 15 °C, which is ~12% of the range of temperatures in the experiment. We also investigate how the algorithm’s accuracy is affected by the number of acoustic receivers used to collect each measurement and the resolution of the temperature prediction. This work provides a means of assessing the safety status of HE material, which cannot be achieved using existing temperature measurement methods. In addition, it has implications for a range of other applications where internal temperature profile measurements would provide critical information. These applications include detecting chemical reactions, observing thermodynamic processes such as combustion, monitoring metal or plastic casting, determining the energy density in thermal storage capsules, and identifying abnormal battery operations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗
DOE OSTI2021

Comparison of Observations and Predictions of Daytime Planetary-Boundary-Layer Heights and Surface Meteorological Variables in the Columbia River Gorge and Basin During the Second Wind Forecast Improvement Project

The second Wind Forecast Improvement Project (WFIP2) is an 18-month field campaign in the Pacific Northwest U.S.A., whose goal is to improve the accuracy of numerical-weather-prediction forecasts in complex terrain. The WFIP2 campaign involved the deployment of a large suite of in situ and remote sensing instrumentation, including eight 915-MHz wind-profiling radars, and surface meteorological stations. The evolution and annual variability of the daytime convective planetary-boundary-layer (PBL) height is investigated using the wind-profiling radars. Three models with different horizontal grid spacing are evaluated: the Rapid Refresh, the High-Resolution Rapid Refresh, and its nested version. The results are used to assess errors in the prediction of PBL height within the experimental and control versions of the models, with the experimental versions including changes and additions to the model parametrizations developed during the field campaign, and the control version using the parametrizations present in the National Oceanic and Atmospheric Administration/National Centers for Environmental Prediction operational version of the models at the start of the project. Results show that the high-resolution models outperform the low-resolution versions, the experimental versions perform better compared with the control versions, model PBL height estimations are more accurate on cloud-free days, and model estimates of the PBL height growth rate are more accurate than model estimates of the rate of decay. Finally, using surface sensors, we assess surface meteorological variables, finding improved surface irradiance and, to a lesser extent, improved 2-m temperature in the experimental version of the model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

How machine learning can extend electroanalytical measurements beyond analytical interpretation

Electroanalytical measurements are routinely used to estimate material properties exhibiting current and voltage signatures. Analysis of such measurements relies on analytical expressions of material properties to describe the experiments. The need for analytical expressions limits the experiments that can be used to measure properties as well as the properties that can be estimated from a given experiment. Such analytical relations are essentially solutions of the physics-based differential equations (with properties as coefficients) describing the material behavior under certain specific conditions. In recent years, a new machine learning-based approach has been gaining popularity wherein the differential equations are numerically solved to interpret the electroanalytical experiments in terms of corresponding material properties. Since the physics-based differential equations are solved, one can additionally estimate underlying fields, e.g., concentration profile, using such an approach. To exemplify the characteristics of such a machine learning assisted interpretation of electroanalytical measurements, we use data from the Hebb–Wagner test on a magnesium spinel intercalation host. In conclusion, as compared to the traditional analytical expression-based interpretation, the emerging approach decreases experimental efforts to characterize relevant material properties as well as provides field information that was previously inaccessible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Design, Flow Field, and Heat Transfer Characterization of the Conjugate Aero-Thermal Test Facility at NETL

A new aerothermal test facility was constructed for the purpose of studying film cooling performance in an environment that accurately simulates conjugate heat transfer characteristics that exist in engine operation. Here, this paper details the design of the facility and the plan for conducting steady state film cooling experiments to improve the understanding of conjugate heat transfer scaling from laboratory to engine conditions. The test facility consists of two separate low channels (hot gas/coolant) and each gas path has a flow conditioning section, a convergent nozzle and a test section/channel with viewports. Numerical simulations were conducted to predict flow field characteristics supporting the design of the flow loop facility. Preliminary experiments were conducted to characterize the flow field using velocity and temperature profile measurements. In addition, infrared (IR) thermography methods were developed to measure surface temperatures on the hot side of the test plate. The IR measurement methods including calibration of the IR camera is explained in detail. It was concluded that appropriate hot gas path flow conditioning could be achieved using a strainer-like tube, a perforated plate, and a honeycomb-mesh screen system upstream of the test section. Flow field measurements from preliminary experiments showed that the boundary layer profile follows the law of the wall.

47 OTHER INSTRUMENTATION↗
NASA NTRS1974

An experimental and numerical study of wave motion and upstream influence in a stratified fluid

A system consisting of two superimposed layers of liquid of different densities, with a thin transition layer at the interface, provides a good laboratory model of an ocean thermocline or of an atmospheric inversion layer. This research was to gain knowledge about the propagation of disturbances within these two geophysical systems. The technique used was to observe the propagation of internal waves and of upstream influence within the density-gradient region between the two layers of liquid. The disturbances created by the motion of a vertical flat plate, which was moved longitudinally through this region, were examined both experimentally and numerically. An upstream influence, which resulted from a balance of inertial and gravitational forces, was observed, and it was possible to predict the behavior of this influence with the numerical model. The prediction included a description of the propagation of the upstream influence to steadily increasing distances from the flat plate and the shapes and magnitudes of the velocity profiles.

Hurdis, D. A.↗
NASA NTRS1985

Ground based concept for time control of aircraft entering the terminal area

This paper summarizes current research for increasing the efficiency of traffic flow at hub airports by application of four dimensional (4D) guidance techniques. A method for generating 4D guidance commands to control the landing time of aircraft that are not equipped with on-board 4D guidance systems is described. In one possible implementation the commands are generated in an auxiliary processor linked to the ATC host computer and displayed on a controller's monitor in the form of profile descent advisories. Accurate time control is achieved by generating air traffic controller advisories in a ground-based algorithm that combines aerodynamic, thrust and atmospheric models with an efficient numerical integration method. The time accuracy and fuel efficiency achieved when the pilot responds to these advisories were evaluated in a piloted simulation of a transport aircraft.

Erzberger, H.↗
NASA NTRS1987

Numerical models of star formation in X-ray cluster cooling flows

A grid of numerical models is used here to study how accretion flows are affected by star formation. The observables associated with these models are investigated, presenting their X-ray surface brightness profiles. An X-ray diagnostic is searched for that, in the absence of spatially resolved X-ray spectral data, can be an unambiguous indicator for ongoing star formation in cooling flows. The spatial distributions of newly formed stars are presented and compared to profiles of the accreting galaxy in each model. Whether the accretion population can form the luminous parts of the accreting galaxy or whether it forms a dark component is considered.

White, Raymond E., III↗
NASA NTRS1991

Adjoint operator approach to shape design for internal incompressible flows

The problem of determining the profile of a channel or duct that provides the maximum static pressure rise is solved. Incompressible, laminar flow governed by the steady state Navier-Stokes equations is assumed. Recent advances in computational resources and algorithms have made it possible to solve the direct problem of determining such a flow through a body of known geometry. It is possible to obtain a set of adjoint equations, the solution to which permits the calculation of the direction and relative magnitude of change in the diffuser profile that leads to a higher pressure rise. The solution to the adjoint problem can be shown to represent an artificially constructed flow. This interpretation provides a means to construct numerical solutions to the adjoint equations that do not compromise the fully viscous nature of the problem. The algorithmic and computational aspects of solving the adjoint equations are addressed. The form of these set of equations is similar but not identical to the Navier-Stokes equations. In particular some issues related to boundary conditions and stability are discussed.

Cabuk, H.↗
NASA NTRS1992

High-resolution processing of long-pulse-lidar data

The purpose of this work is to demonstrate the performance of the Fourier-deconvolution technique developed to improve the resolution of the long-pulse coherent lidar power profile, and to take into account the multiplicative fluctuations. The behavior of the error due to the power fluctuations and a method to reduce it are analyzed theoretically and simulated numerically. A processing of real data obtained from the National Oceanic and Atmospheric Administration (NOAA) ground-based Doppler lidar is also presented. Similar problems were investigated using a different approach based on introducing a correction function to the lidar equation.

Gurdev, L. L.↗
NASA NTRS1991

IMP: A performance code

IMP (Integrated Mission Program) is a simulation language and code used to model present and future Earth, Moon, or Mars missions. The profile is user controlled through selection from a large menu of events and maneuvers. A Fehlberg 7/13 Runge-Kutta integrator with error and step size control is used to numerically integrate the differential equations of motion (DEQ) of three spacecraft, a main, a target, and an observer. Through selection, the DEQ's include guided thrust, oblate gravity, atmosphere drag, solar pressure, and Moon gravity effects. Guide parameters for thrust events and performance parameters of velocity changes (Delta-V) and propellant usage (maximum of five systems) are developed as needed. Print, plot, summary, and debug files are output.

Dauro, Vincent A., Sr.↗
NASA NTRS1960

The Blunt Plate In Hypersonic Flow

The sonic-wedge characteristics method has been used to obtain the shock shapes and surface pressure distributions on several blunt two-dimensional shapes in a hypersonic stream for several values of the ratio of specific heats. These shapes include the blunt slab at angle of attack and power profiles of the form yb = a)P, where 0 les than m less than 1, Yb and x are coordinates of the body surface, and a is a constant. These numerical results have been compared with the results of blast-wave theory, and methods of predicting the pressure distributions and shock shapes are proposed in each case. The effects of a free-stream conical-flow gradient on the pressure distribution on a blunt slab in hypersonic flow were investigated by the sonic-wedge characteristics method and were found to be sizable in many cases. Procedures which are satisfactory for reducing pressure data obtained in conical flows with small gradients are presented.

Baradell, Donald L.↗
NASA NTRS1999

Directional Solidification of Pure Succinonitrile and a Succinonitrile-Acetone Alloy

An experimental study of the horizontal Bridgman growth of pure succinonitrile (SCN) and of a succinonitrile-1.0 mol.% acetone alloy (SCN-1.0 mol.% ACE) has been performed. Experiments involving both a stationary thermal field (no-growth case) and a translating thermal field (growth case) were conducted. Growth rates of 2 and 40 micrometers/s were investigated. For the pure SCN experiments, the velocity field in the melt was estimated using video images of seed particles in the melt. Observations of the seed particles indicate that a primary longitudinal convective cell is formed. The maximum velocity of two different particles which traveled along similar paths was the same and equal to 1.49 +/- 0.01 mm/s. The general accuracy of velocity measurements is estimated to be +/- 0.08 mm/s, though the data shows consistency to within +/- 0.02 mm/s. The shape of the solid/liquid interface was also quantitatively determined. The solid/liquid interface was stable (non-dendritic and non-cellular) but not flat; rather it was significantly distorted by the influence of convection in (he melt and, for the growth case, by the moving temperature boundary conditions along the ampoule. It was found that the interface shape and position were highly dependent on the alignment of the ampoule in the apparatus. Consequently, the ampoule was carefully aligned for all experiments. The values for front location agree with those determined in previous experiments. For the alloy experiments, the solid/liquid interface was determined to be unstable at growth rates greater than 2.8 micrometers/s, but stable for the cases of no-growth and growth at 2 micrometers/s. When compared to the shape of the pure SCN interface, the alloy interface forms closer to the cold zone, indicating that the melting temperature decreased due to the alloying element. Extensive temperature measurements were performed on the outside of the ampoule containing pure SCN. The resulting thermal profiles are presented in detail in the results section. It is intended that the interface shape, thermal boundary condition and velocity data presented in this paper be used to test numerical simulations.

Simpson, James E.↗
NASA NTRS2000

Directional Solidification of Pure Succinonitrile and a Succinonitrile-Acetone Alloy

An experimental study of the horizontal Bridgman growth of pure succinonitrile (SCN) and of a succinonitrile-1.0 mol% acetone alloy (SCN-1.0 mol.% ACE) has been performed. Experiments involving both a stationary thermal field (no-growth case) and a translating thermal field (growth case) were conducted. Growth rates of 2 and 40 micrometers/sec were investigated. For the pure SCN experiments, the velocity field in the melt was estimated using video images of seed particles in the melt. Observations of the seed particles indicate that a primary longitudinal convective cell is formed. The maximum velocity of two different particles which traveled along similar paths was the same and equal to 1.49 +/- 0.01 mm/s. The general accuracy of velocity measurements is estimated to be +/-0.08 mm/s, though the data shows consistency to within +/- 0.02 mm/s. The shape of the solid/liquid interface was also quantitatively determined. The solid/liquid interface was stable (non-dendritic and non-cellular) but not flat: rather it was significantly distorted by the influence of connection in the melt and, for the growth case, by the moving temperature boundary conditions along the ampoule. It was found that the interface shape and position were highly dependent on the alignment of the ampoule in the apparatus. Consequently, the ampoule was carefully aligned for all experiments. The values for front location agree with those determined in previous experiments. For the alloy experiments, the solid/liquid interface was determined to be unstable at growth rates greater than 2.8 micrometers/sec, but stable for the cases of no-growth and growth at 2 micrometers/sec. When compared to the shape of the pure SCN interface, the alloy interface forms closer to the cold zone, indicating that the melting temperature decreased due to the alloying element. Extensive temperature measurements were performed on the outside of the ampoule containing pure SCN. The resulting thermal profiles are presented in detail in the results section. It is intended that the interface shape, thermal boundary condition, and velocity data presented in this paper be used to test numerical simulations.

Simpson, James E.↗
NASA NTRS2000

Towards LES Models of Jets and Plumes

As pointed out by Rodi standard integral solutions for jets and plumes developed for discharge into infinite, quiescent ambient are difficult to extend to complex situations, particularly in the presence of boundaries such as the sea floor or ocean surface. In such cases the assumption of similarity breaks down and it is impossible to find a suitable entrainment coefficient. The models are also incapable of describing any but the most slowly varying unsteady motions. There is therefore a need for full time-dependent modeling of the flow field for which there are three main approaches: (1) Reynolds averaged numerical simulation (RANS), (2) large eddy simulation (LES), and (3) direct numerical simulation (DNS). Rodi applied RANS modeling to both jets and plumes with considerable success, the test being a match with experimental data for time-averaged velocity and temperature profiles as well as turbulent kinetic energy and rms axial turbulent velocity fluctuations. This model still relies on empirical constants, some eleven in the case of the buoyant jet, and so would not be applicable to a partly laminar plume, may have limited use in the presence of boundaries, and would also be unsuitable if one is after details of the unsteady component of the flow (the turbulent eddies). At the other end of the scale DNS modeling includes all motions down to the viscous scales. Boersma et al. have built such a model for the non-buoyant case which also compares well with measured data for mean and turbulent velocity components. The model demonstrates its versatility by application to a laminar flow case. As its name implies, DNS directly models the Navier-Stokes equations without recourse to subgrid modeling so for flows with a broad spectrum of motions (high Re) the cost can be prohibitive - the number of required grid points scaling with Re(exp 9/4) and the number of time steps with Re(exp 3/4). The middle road is provided by LES whereby the Navier-Stokes equations are formally filtered with the filter chosen to only exclude the smallest turbulent motions. If successful, LES should provide much of the detail available to DNS but at more bearable cost. Fatica et al. in comparing LES with DNS for a low Reynolds number jet showed that the LES could simulate the temporally evolving behavior including growth of the jet thickness. It is the intention of this report to explore the application of an LES model to jets and plumes. As always, before tackling complex situations, the model must be tested for the simplest of cases and so we address only two, a non-buoyant axisymmetric jet issuing steadily from an orifice into a semi-infinite stationary environment and a buoyant jet in the same environment. The work is a continuation of Basu and Mansour.

Webb, A. T.↗