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 937 records · Page 52

Neural entropy-stable conservative flux form neural networks for learning hyperbolic conservation laws

We propose a neural entropy-stable conservative flux form neural network (NESCFN) for learning hyperbolic conservation laws and their associated entropy functions directly from solution trajectories, without requiring any predefined numerical discretization. While recent neural network architectures have successfully integrated classical numerical principles into learned models, most rely on prior knowledge of the governing equations or assume a fixed discretization. Our approach removes this dependency by embedding entropy-stable design principles into the learning process itself, enabling the discovery of physically consistent dynamics in a fully data-driven setting. By jointly learning both the flux function and a corresponding entropy, NESCFN promotes conservation and entropy dissipation, which is critical for long-term stability and fidelity in the system of hyperbolic conservation laws. Furthermore, numerical results demonstrate that the method achieves stability and conservation over extended time horizons and accurately captures shock propagation speeds, even without oracle access to future-time solution profiles in the training data.

Conservative flux form↗

Paradigm for global gyrokinetic turbulence

Using a novel wave-number-advection algorithm, we describe a method to add global profile curvature terms to the local gyrokinetic equations. This approach enables a high-precision global simulation capability without sacrificing the efficiency and spectral accuracy of local simulations. Preliminary numerical results suggest that there are two distinct types of global correction to the local result: (1) a trivial shift of the local result in radius, (2) a non-trivial decorrelation that is roughly equivalent to the well-known E×B shear effect. This interpretation provides a new, intuitive framework for understanding global gyrokinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗

Accretion process of the moon

Recent geochemical and geophysical data suggest that the initial temperature of the moon was strongly peaked toward the lunar surface. To explain such an initial temperature distribution, a simple model of accretion process of the moon is presented. The model assumes that the moon was formed from the accumulation of the solid particles or gases in the isolated, closed cloud. Two equations are derived to calculate the accretion rate and surface temperature of the accreting moon. Numerical calculations are made for a wide range of the parameters particle concentration and particle velocity in the cloud. A limited set of the parameters gives the initial temperature profiles as required by geochemical and geophysical data. These models of the proto-moon cloud indicate that the lunar outermost shell, about 400 km thick, was partially or completely molten just after the accretion of the moon and that the moon should have been formed in a period shorter than 1000 years. If the moon formed at a position nearer to the earth than its present one, the moon might have been formed in a period of less than one year.

Mizutani, H.↗

The interplanetary transport of solar cosmic rays

Numerical solutions are presented for the propagation of solar cosmic rays in interplanetary space, including the effects of pitch-angle scattering and adiabatic focusing. The intensity-time profiles can be well fitted by a simple radial spatial diffusion equation with scattering mean-free path lambda(fit). The radial mean-free path so obtained is significantly larger than the true scattering mean-free path for low-rigidity particles due to both adiabatic focusing and the inapplicability of the diffusive approximation early in the event. The well-known discrepancy between lambda(fit) and the theoretical predictions may be resolved by these calculations.

Gombosi, T. I.↗

Numerical study of solar flare particle propagation in the heliosphere

Numerical solutions are presented for the propagation of solar cosmic rays in interplanetary space, including the effects of pitch-angle scattering and adiabatic focusing. The intensity-time profiles can be well fitted by a simple radial spatial diffusion equation. For low-rigidity particles the radial mean free path so obtained is significantly larger than the mean free path calculated from the scattering coefficient due to the inapplicability of the diffusive approximation early in the event. The well-known discrepancy between the scattering mean free path and the theoretical predictions may be resolved by these calculations.

Gombosi, T. I.↗

Theoretical kinetic computations in complex reacting systems

Nasa Lewis' studies of complex reacting systems at high temperature are discussed. The changes which occur are the result of many different chemical reactions occurring at the same time. Both an experimental and a theoretical approach are needed to fully understand what happens in these systems. The latter approach is discussed. The differential equations which describe the chemical and thermodynamic changes are given. Their solution by numerical techniques using a detailed chemical mechanism is described. Several different comparisons of computed results with experimental measurements are also given. These include the computation of (1) species concentration profiles in batch and flow reactions, (2) rocket performance in nozzle expansions, and (3) pressure versus time profiles in hydrocarbon ignition processes. The examples illustrate the use of detailed kinetic computations to elucidate a chemical mechanism and to compute practical quantities such as rocket performance, ignition delay times, and ignition lengths in flow processes.

Bittker, David A.↗

Effect of gravity level fluctuations for rotating fluids in high and low rotating speeds

Time-dependent evolutions of the profile of the free surface (bubble shapes) for a cylindrical container partially filled with a Newtonian fluid of constant density, rotating about its axis of symmetry, have been studied. Numerical computations of the dynamics of bubble shapes have been carried out with sinusoidal-function vibration of the gravity environment at high and low cylinder speeds.

Hung, R. J.↗

Computer modeling of the dynamics of surface tension on rotating fluids in low and microgravity environments

Time-dependent evolutions of the profile of the free surface (bubble shapes) for a cylindrical container partially filled with a Newtonian fluid of constant density, rotating about its axis of symmetry, have been studied. Numerical computations have been carried out with the following situations: (1) linear functions of spin-up and spin-down in low- and microgravity environments, (2) linear functions of increasing and decreasing gravity environments at high- and low-rotating cylinder speeds, and (3) step functions of spin-up and spin-down in a low-gravity environment.

Hung, R. J.↗

Nonlinear effects in a plain journal bearing. I - Analytical study. II - Results

In the first part of this work, a numerical model is presented which couples the variable-property Reynolds equation with a rotor-dynamics model for the calculation of a plain journal bearing's nonlinear characteristics when working with a cryogenic fluid, LOX. The effects of load on the linear/nonlinear plain journal bearing characteristics are analyzed and presented in a parametric form. The second part of this work presents numerical results obtained for specific parametric-study input variables (lubricant inlet temperature, external load, angular rotational speed, and axial misalignment). Attention is given to the interrelations between pressure profiles and bearing linear and nonlinear characteristics.

Choy, F. K.↗

An investigation of DTNS2D for use as an incompressible turbulence modelling test-bed

This paper documents an investigation of a two dimensional, incompressible Navier-Stokes solver for use as a test-bed for turbulence modelling. DTNS2D is the code under consideration for use at the Center for Modelling of Turbulence and Transition (CMOTT). This code was created by Gorski at the David Taylor Research Center and incorporates the pseudo compressibility method. Two laminar benchmark flows are used to measure the performance and implementation of the method. The classical solution of the Blasius boundary layer is used for validating the flat plate flow, while experimental data is incorporated in the validation of backward facing step flow. Velocity profiles, convergence histories, and reattachment lengths are used to quantify these calculations. The organization and adaptability of the code are also examined in light of the role as a numerical test-bed.

Steffen, Christopher J., Jr.↗

Numerical models for the collapse and fragmentation of centrally condensed molecular cloud cores

The gravitational collapse and fragmentation of centrally condensed molecular cloud cores are investigated using a new hydrodynamical code. The numerical scheme is second-order accurate and uses explicit finite difference methods to advance the fluid variables on a 3D Cartesian grid. Two initial power-law density profiles, rho varies as r exp -1 and rho varies as r exp -2, are considered, as well as two initial density perturbations in the azimuthal coordinate theta, rho-i - rho(1 + a cos 2theta) where a = 0.1 and 0.5. Fragmentation is found to be possible in these centrally condensed cores if the initial conditions also include differential rotation. Models which collapse with initial uniform rotation do not produce fragments. If molecular cloud cores are indeed centrally condensed, as suggested by observations of star-forming regions and by studies of ambipolar diffusion, then differential rotation may be a mechanism for producing binary protostars during gravitational collapse.

Myhill, Elizabeth A.↗

Computational/experimental investigation of staged injection into a Mach 2 flow

The stage normal injection of two jets of air into a Mach 2 freestream behind a rearward-facing step has been investigated using laser-induced iodine fluorescence and laser doppler anemometry techniques. A detailed data set has been compiled that includes profiles of pressure, temperature, two components of velocity, and mole fraction of the injectant at 44 locations, plus planar surveys of pressure, temperature, velocity, and mole fraction. A companion numerical study was performed using the SPARK Navier-Stokes code. Good overall agreement was observed for this complex, highly three-dimensional flowfield. Discrepancies were observed in the computed and measured total temperatures, turbulent mixing, and shock strengths. The effect of grid resolution was investigated by calculating solutions on grids of 60,000 points, 200,000 points, and 450,000 points. Differences in the solutions on the two finer grids were small. The effect of turbulence modeling was investigated by calculating solutions with three different algebraic models for the jet turbulence. Overall, the turbulence models were found to have the greatest effect on the numerical solutions, followed by the grid resolution and the injectant Mach number.

Eklund, D. R.↗

Fundamental Phenomena on Fuel Decomposition and Boundary-Layer Combustion Precesses with Applications to Hybrid Rocket Motors: Experimental Investigation - Part 1

This final report summarizes the major findings on the subject of 'Fundamental Phenomena on Fuel Decomposition and Boundary-Layer Combustion Processes with Applications to Hybrid Rocket Motors', performed from 1 April 1994 to 30 June 1996. Both experimental results from Task 1 and theoretical/numerical results from Task 2 are reported here in two parts. Part 1 covers the experimental work performed and describes the test facility setup, data reduction techniques employed, and results of the test firings, including effects of operating conditions and fuel additives on solid fuel regression rate and thermal profiles of the condensed phase. Part 2 concerns the theoretical/numerical work. It covers physical modeling of the combustion processes including gas/surface coupling, and radiation effect on regression rate. The numerical solution of the flowfield structure and condensed phase regression behavior are presented. Experimental data from the test firings were used for numerical model validation.

Kuo, Kenneth K.↗

Experimental and Computational Study fo CH, CH*, and OH* in an Axisymmetric Laminar Diffusion Flame

In this study, we extend the results of previous combined numerical and experimental investigations of an axisymmetric laminar diffusion flame in which difference Raman spectroscopy, laser-induced fluorescence (LIF), and a multidimensional flame model were used to generate profiles of the temperature and major and minor species. A procedure is outlined by which the number densities of ground-state CH (X(sup 2)II) excited-state CH (A(sup 2)Delta, denoted CH*), and excited-state OH (A(sup 2)Sigma, denoted OH*) are measured and modeled. CH* and OH* number densities are deconvoluted from line-of-sight flame-emission measurements. Ground-state CH is measured using linear LIF. The computations are done with GRI Mech 2.11 as well as an alternate hydrocarbon mechanism. In both cases, additional reactions for the production and consumption of CH* and OH* are added from recent kinetic studies. Collisional quenching and spontaneous emission are responsible for the de-excitation of the excited-state radicals. As with our previous investigations, GRI Mech 2.11 continues to produce very good agreement with the overall flame length observed in the experiments, while significantly under predicting the flame lift-off height. The alternate kinetic scheme is much more accurate in predicting lift-off height but overpredicts the over-all flame length. Ground-state CH profiles predicted with GRI Mech 2.11 are in excellent agreement with the corresponding measurements, regarding both spatial distribution and absolute concentration (measured at 4 ppm) of the CH radical. Calculations of the excited-state species show reasonable agreement with the measurements as far as spatial distribution and overall characteristics are concerned. For OH*, the measured peak mole fraction, 1.3 x 10(exp -8), compared well with computed peaks, while the measured peak level for CH*, 2 x 10(exp -9), was severely underpredicted by both kinetic schemes, indicating that the formation and destruction kinetics associated with excited-state species in flames require further research.

Walsh, K. T.↗

Application of an Orbital GPR Model to Detecting Martian Polar Subsurface Features

There are numerous challenges in successfully implementing and interpreting planetary ground penetrating radar (GPR) measurements. Many are due to substantial uncertainties in the target ground parameters and the intervening medium (i.e., the ionosphere). These uncertainties generate a compelling need for meaningful quantitative simulation of the planetary GPR problem. An accurate numerical model would enable realistic numerical GPR simulations using parameter regimes much broader than are possible in laboratory or field experiments. Parameters such as source bandwidth and power, surface and subsurface features, and ionospheric profiles could be rapidly iterated to understand their impact on GPR performance and the reliable interpretation of GPR data.

Xu, Y.↗

The Effect of the Saharan Air Layer on the Formation of Hurricane Isabel (2003) Simulated with AIRS Data

The crucial physics of how the atmosphere really accomplishes the tropical cyclogenesis process is still poorly understood. The presence of the Saharan Air Layer (SAL), an elevated mixed layer of warm and dry air that extends from Africa to the tropical Atlantic and contains a substantial amount of mineral dust, adds more complexity to the tropical cyclogenesis process in the Atlantic basin. The impact of the SAL on tropical cyclogenesis is still uncertain. Karyampudi and Carlson (1988) conclude that a strong SAL can potentially aid tropical cyclone development while Dunion and Velden (2004) argue that the SAL generally inhibits tropical cyclogenesis and intensification. Advancing our understanding of the physical mechanisms of tropical cyclogenesis and the associated roles of the SAL strongly depends on the improvement in the observations over the data-sparse ocean areas. After the Atmospheric Infrared Sounder (AIRS), the Advanced Microwave Sounding Unit (AMSU), and the microwave Humidity Sounder of Brazil (HSB) were launched with the NASA Aqua satellite in 2002, new data products retrieved from the AIRS suite became available for studying the effect of the warm, dry air mass associated with the SAL (referred to as the thermodynamic effect). The vertical profiles of the AIRS retrieved temperature and humidity provide an unprecedented opportunity to examine the thermodynamic effect of the SAL. The observational data can be analyzed and assimilated into numerical models, in which the model thermodynamic state is allowed to relax to the observed state from AIRS data. The objective of this study is to numerically demonstrate that the thermodynamic effect of the SAL on the formation of Hurricane Isabel (2003) can be largely simulated through nudging of the AIRS data.

Wu, iguang↗

Submm-Wave Radiometry for Cloud/Humidity/Precipitation Sciences

Although active sensors can provide cloud profiles at good vertical resolution, clouds are often coupled with dynamics to form fast and organized structures. Lack of understanding of these organized systems leads to great challenge for numerical models. The deficiency is partly reflected, for example, in poorly modeled intraseasonal variations (e.g., MJD). Remote sensing clouds in the middle and upper troposphere has been challenging from space. Vis/IR sensors are sensitive to the topmost cloud layers whereas low-frequency MW techniques are sensitivity to liquid and precipitation at the bottom of cloud layers. The middle-level clouds, mostly in the ice phase, require a sensor that has moderate penetration and sensitivity to cloud scattering, in order to measure cloud water content. Sensors at submm wavelengths provide promising sensitivity and coverage with the spatial resolution needed to measure cloud water content floating in the upper air. In addition, submm-wave sensors are able to provide better measurements of upper-tropospheric humidity than traditional microwave instruments.

Wu, Dong L.↗