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.

SourceResearch connectionFrom yearThrough year

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

At least 1,315 records · Page 73

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Cartesian Mesh Simulations for the Third AIAA Sonic Boom Prediction Workshop

Simulation results are presented for all cases from the Third AIAA Sonic Boom Prediction Workshop. An inviscid, embedded-boundary Cartesian-mesh flow solver is used in conjunction with adjoint-based mesh adaptation to compute nearfield pressure signatures. Specialized techniques are applied to maximize accuracy and minimize cost on Cartesian meshes. The Richardson-based error estimate highlights regions of the signatures most sensitive to mesh refinement. Timing results and coarse, medium, and fine mesh sizes for nearfield cases demonstrate that the parallel decomposition approach is efficient in both computational time and wall-clock. Pressure signals are propagated to the ground using an augmented Burgers’ equation solver to predict boom carpets. Ground signatures and loudness metrics are presented for a standard atmosphere as well as more realistic atmospheric profiles, which affect overall noise levels and can significantly widen the boom carpet. Mesh convergence studies show that high sampling frequencies, around 500 kHz, are required for propagation, and the sampling frequency increases at large off-track angles with longer acoustic ray paths and propagation times. The numerical methods yield accurate results for predicting low sonic boom signatures while being among the least computationally expensive of the workshop.

ARMD↗

NASA NTRS2005

Models for Water Isotopes Constrained with Data from Crystal Face

During the year covered by this proposal we conducted work on several different topics, as reflected by our publications listed below. One major activity was to work with a group of about 10 scientists from around the country to prepare a science-planning document (Tropical Composition, Cloud and Climate Coupling Experiment (TC4)) that outlined the rationale, locations, strategy to accomplish the goals, and possible payloads for a set of three tropical missions. We also prepared background materials for various NRAs being prepared at NASA Headquarters for missions in Costa Rica, Darwin and Guam. Unfortunately budgetary constraints prevented these missions from moving forward. In conjunction with the group NASA Ames we built a new numerical model for deep convection and have applied that model to simulate the CRYSTAL isotope data. Our goal in particular has been to better understanding how convection distributes water vapor isotopes. CRYSTAL observations of water isotopes are very different from those suggested by previous workers who assumed the isotopes would obey Rayleigh fractionation. The water isotope study has several implications. First it is a check on the realism of the deep convection model. Second, the isotopes are a measure of the precipitation removal in the atmosphere. Hence they provide a constraint on a parameter that is difficult to otherwise measure. Finally it has been suggested that isotopes may be the key to unraveling the water transport into the stratosphere and upper troposphere. Such transport is critical both for the radiation balance and for stratospheric chemistry. Ours is the first model that is able to treat this transport. Our initial results are now in press in Geophys. Res. Lett. Essentially we are able to explain the vertical profiles of isotopes in the tropical tropopause transition layer. We are also able to account for stratospheric humidity ana isotope abundances with this model. The data suggest that isotopes do not provide a clear constraint on the mechanism by which water enters the stratosphere- whether by convection, or by slow ascent. Our work is relevant for the water isotope comparison experiments recently done by NASA. We are conducting numerical experiments related to this project to help understand that data.

Toon, Owen B.↗

NASA NTRS1976

Airfoil section drag reduction at transonic speeds by numerical optimization

A practical procedure for the design of low drag, transonic airfoils is demonstrated. The procedure uses an optimization program, based on a gradient algorithm coupled with an aerodynamic analysis program, that solves the full, non-linear potential equation for transonic flow. The procedure is useful for the design of retrofit modifications for drag reduction of existing aircraft as well as for the design of low drag profiles for new aircraft. Results are presented for the modification of four different airfoils to decrease the drag at a given transonic Mach number.

Hicks, R. M.↗

NASA NTRS1976

Airfoil section drag reduction at transonic speeds by numerical optimization

Hicks, R. M.↗

DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗

NASA NTRS1995

Bleed Roughness Mechanism

Suction applied through wall bleed is expected to reduce boundary-layer thickness and lead to a profile that can resist separation. However, experiments with rows of suction holes or slots have shown that the benefits of suction are not realized under certain circumstances. Lee, Sloan and Paynter attributed this to a 'bleed roughness', and modeled the phenomena using a modified inner-length scale with convective lag for the eddy viscosity. Test information was required to determine the roughness values at each Mach number and bleed rate. The present work used numerical calculation, subsonic and supersonic experimentation, and physical modeling were used to explore bleed roughness. It was found that the phenomena can exist in the absence of bleed, that it is a function of geometry and Mach number, and that it does not appear in low subsonic flow. The inner-length scale/convective lag approach was found to be inconsistent with recent physical observation. Alternate steady models for the roughness effect disagreed with mass flux trends and the choked flow limit, and suggested that non-steady mechanisms are present. While new supersonic experiments did exhibit low-frequency plenum oscillations, high-frequency unsteadiness associated with the orifice flow scale appears to be the cause of bleed roughness. In fact, boundary-layer calculations with an oscillatory orifice flow produced the same effect in the time-averaged downstream profile for incompressibile flow. A new empirical turbulence-model was proposed which increased the eddy viscosity throughout the entire boundary-layer by a constant which was calibrated using the zero-bleed profile. Preliminary calculations agreed with the experimental data. The model significantly reduced the test information required to determine roughness values for turbulence modeling, and provided prediction capability for supersonic bleed flows. An abstract of a thesis entitled 'Boundary-Layer Bleed Roughness' is included.

Russell, David A.↗

DOE OSTI2026

A buoyancy–shear–drag–scalar-based turbulence model for power-law acceleration-driven Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.

Buoyancy–drag↗

DOE OSTI2024

Development of a neural network model for peeling–ballooning stability analysis in the KSTAR tokamak pedestals

The neural network model, MISHKA-NN is developed to mitigate the computational burden associated with the linear ideal magnetohydrodynamic (MHD) stability analysis of the pedestal based on the peeling–ballooning (P–B) model. By utilizing both 1D plasma profiles (current density, pressure gradient, and safety factor) and 0D parameters (plasma geometry, total current, and toroidal mode number), the model predicts linear growth rate of edge-localized ideal MHD instability in a given equilibrium state. By enabling the prediction of each instability within a second, the model reduces the time required for plotting a pedestal P–B stability diagram (j - α diagram) from approximately 100 CPU hours to a few CPU minutes. Notably, even with the utilization of parametric pressure and current profiles and plasma boundary shapes for the training dataset, the model shows a satisfactory level of performance in benchmarking the j - α diagram for the reconstructed equilibrium from a KSTAR tokamak experiment. We anticipate the model to serve as a versatile alternative to 2D linear MHD stability codes, alleviating numerical costs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Control of plasmons in topological insulators via local perturbations

Here, we use a fully quantum mechanical approach to demonstrate control of plasmonic excitations in prototype models of topological insulators by molecule-scale perturbations. Strongly localized surface plasmons are present in the host systems, arising from the topologically nontrivial single-particle edge states. A numerical evaluation of the random phase approximation equations for the perturbed systems reveals how the positions and the internal electronic structure of the added molecules affect the degeneracy of the locally confined collective excitations, i.e., shifting the plasmonic energies of the host system and changing their spatial charge density profile. In particular, we identify conditions under which significant charge transfer from the host system to the added molecules occurs. Furthermore, the induced field energy density in the perturbed topological systems due to external electric fields is determined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Evaluation of a quasi-steady-state approximation of the cloud droplet growth equation (QDGE) scheme for aerosol activation in global models using multiple aircraft data over both continental and marine environments

This research introduces a numerically efficient aerosol activation scheme and evaluates it by using stratus and stratocumulus cloud data sampled during multiple aircraft campaigns in Canada, Chile, Brazil, and China. The scheme employs a quasi-steady-state approximation of the cloud droplet growth equation (QDGE) to efficiently simulate aerosol activation, the vertical profile of supersaturation, and the activated cloud droplet number concentration (CDNC) near the cloud base. The calculated maximum supersaturation values using the QDGE scheme were compared with multiple parcel model simulations under various aerosol and environmental conditions. The differences are all below 0.18 %, indicating good performance and accuracy of the QDGE scheme. We evaluated the QDGE scheme by specifying observed environmental thermodynamic variables and aerosol information from 31 cloud cases as input and comparing the simulated CDNC with cloud observations. The average of mean relative error ($\overline{MRE}$) of the simulated CDNC for cloud cases in each campaign ranges from 17.30 % in Brazil to 25.90 % in China, indicating that the QDGE scheme successfully reproduces observed variations in CDNC over a wide range of different meteorological conditions and aerosol regimes. Additionally, we carried out an error analysis by calculating the maximum information coefficient (MIC) between the MRE and input variables for the individual campaigns and all cloud cases. MIC values were then sorted by aerosol properties, pollution level, environmental humidity, and dynamic condition according to their relative importance to MRE. Based on the error analysis, we found that the magnitude of MRE is more relevant to the specification of input aerosol pollution level in marine regions and aerosol hygroscopicity in continental regions than to other variables in the simulation.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1972

Assessment of a transitional boundary layer theory at low hypersonic Mach numbers

An investigation was carried out to assess the accuracy of a transitional boundary layer theory in the low hypersonic Mach number regime. The theory is based upon the simultaneous numerical solution of the boundary layer partial differential equations for the mean motion and an integral form of the turbulence kinetic energy equation which controls the magnitude and development of the Reynolds stress. Comparisions with experimental data show the theory is capable of accurately predicting heat transfer and velocity profiles through the transitional regime and correctly predicts the effects of Mach number and wall cooling on transition Reynolds number. The procedure shows promise of predicting the initiation of transition for given free stream disturbance levels. The effects on transition predictions of the pressure dilitation term and of direct absorption of acoustic energy by the boundary layer were evaluated.

Shamroth, S. J.↗

NASA NTRS1974

Thermodynamics and kinetics of reactions in protective coating systems

Investigations of fluoride activated packs with Al:Ni ratios greater than 50 a/o prove that the specimen surface is not in equilibrium with the pack at high Al:Ni ratios but that an activity gradient exists between pack and specimen. Therefore, gaseous diffusion and possibly surface reactions play a role in determining the overall rate of Al deposition in such packs. Noticeable differences in coating behavior have been obtained in packs activated with chloride and iodide, and it appears that poorest results are obtained with iodides, better with chlorides, and best with fluorides. A numerical method has been perfected for calculating rates of solid-state diffusion controlled coating formation, allowing for the variation of diffusivity with composition in the NiAl phase. Layer growth rates can now be accurately predicted from a knowledge of the surface and substrate compositions. Furthermore, the correct diffusion profiles are obtained by this method. These differ substantially from the profile obtained when the diffusivity is assumed constant.

Gupta, B.↗

NASA NTRS1977

Analytical study of laser supported combustion waves in hydrogen

A one-dimensional energy equation, with constant pressure and area, was used to model the LSC wave. This equation balances convection, conduction, laser energy absorption, radiation energy loss and radiation energy transport. Solutions of this energy equation were obtained to give profiles of temperature and other properties, as well as the relation between laser intensity and mass flux through the wave. The flow through the LSC wave was then conducted through a variable pressure, variable area streamtube to accelerate it to high speed, with the propulsion application in mind. A numerical method for coupling the LSC wave model to the streamtube flow was developed, and a sample calculation was performed. The result shows that 42% of the laser power has been radiated away by the time the gas reaches the throat. It was concluded that in the radially confined flows of interest for propulsion applications, transverse velocities would be less important than in the unconfined flows where air experiments have been conducted.

Kemp, N. H.↗

NASA NTRS1982

Collapse of accreting, rotating, isothermal, interstellar clouds

The evolutions of the envelopes of collapsing, accreting, isothermal clouds have been numerically calculated for both spherically symmetric and rotating (axisymmetric) clouds. The results provide a cohesive picture of isothermal collapse, and their relationship to previous numerical calculations and similarity solutions is discussed. Even with a large initial rotation rate, the majority of the cloud envelope is accreted, in one case leaving behind a large-scale circulation current. The calculations are performed for both initially uniform density and centrally condensed clouds. Density and velocity profiles for a wide variety of observed systems are compared with those obtained in this study, providing a preliminary assessment of the stage of evolution and initial structure for the observed systems.

Boss, A. P.↗

NASA NTRS1985

Numerical study of large-eddy breakup and its effect on the drag characteristics of boundary layers

The break-up of a field of eddies by a flat-plate obstacle embedded in a boundary layer is studied using numerical solutions to the two-dimensional Navier-Stokes equations. The flow is taken to be incompressible and unsteady. The flow field is initiated from rest. A train of eddies of predetermined size and strength are swept into the computational domain upstream of the plate. The undisturbed velocity profile is given by the Blasius solution. The disturbance vorticity generated at the plate and wall, plus that introduced with the eddies, mix with the background vorticity and is transported throughout the entire flow. All quantities are scaled by the plate length, the unidsturbed free-stream velocity, and the fluid kinematic viscosity. The Reynolds number is 1000, the Blasius boundary layer thickness is 2.0, and the plate is positioned a distance of 1.0 above the wall. The computational domain is four units high and sixteen units long.

Kinney, R. B.↗

NASA NTRS1984

Applications of VAS and TOVS to tropical cyclones

The results of the attempts to use VAS (visible IR spin scan radiometer atmospheric sounder) data for predicting tropical cyclones are reported. VAS, mounted on a GEO platform, has 14 detectors for upwelling radiance. The McIDAS software was, by the 1983 hurricane season, upgraded to providing imagery and mass, motion and moisture patterns using the satellite data, which also permitted determinations of accurate temperature and moisture profiles. Comparisons made between VAS and radiosonde data revealed fairly close agreement of satellite and in-situ conditions during the 1983 season. When combined with data from the Tiros operational vertical sounder (TOVS) to define wind streamlines and predict the path of a hurricane from data fed into numerical models, sufficiently good agreement was obtained to initiate plans to forecast a hurricane trajectory during the 1984 season.

Velden, C. S.↗

NASA NTRS1984

VAS data assimilation into a mesoscale model

In numerical weather prediction, one of the arising problems has always been related to insufficient data, especially in cases involving scales smaller than the synoptic rawinsonde network. The Visible Infrared Spin-Scan Radiometer (VISSR) Atmospheric Sounder (VAS) instrument on the Geostationary Operational Environmental Satellite (GOES) can provide mesoscale vertical profiles of temperature and moisture several times over a twelve hour period. However, as this data is only available over a relatively small area at any one time, the VAS data is more suited to data assimilation. A variational method has been developed to assimilate VAS data into a mesoscale model. In the present paper, the assimilation method and VAS data from one case are discussed. It is found that the greatest improvements to the forecast resulted when the VAS data resolved the mesoscale structure of the temperature and relative humidity fields.

Cram, J. M.↗

NASA NTRS1989

Plasma contacting - An enabling technology

An experimental study of plasma contacting with an emphasis on the electron collection mode of this process is described. Results illustrating variations in plasma property profiles and potential differences that develop at hollow cathode plasma contactors are presented. A model of the electron collection plasma contacting process that is consistent with experimentally measured results is reviewed. The shortcomings of laboratory results as direct predictors of contactor performance in space and their usefulness in validating numerical models of the contacting process, that can be used to predict such performance, are discussed.

Williams, John D.↗