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 397 records · Page 22

Nonlinear power flow control for networked AC/DC microgrids

A method for designing feedforward and feedback controllers for integration of stochastic sources and loads into a nonlinear networked AC/DC microgrid system is provided. A reduced order model for general networked AC/DC microgrid systems is suitable for HSSPFC control design. A simple feedforward steady state solution is utilized for the feedforward controls block. Feedback control laws are provided for the energy storage systems. A HSSPFC controller design is implemented that incorporates energy storage systems that provides static and dynamic stability conditions for both the DC random stochastic input side and the AC random stochastic load side. Transient performance was investigated for the feedforward/feedback control case. Numerical simulations were performed and provided power and energy storage profile requirements for the networked AC/DC microgrid system overall performance. The HSSPFC design can be implemented in the Matlab/Simulink environment that is compatible with real time simulation/controllers.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Numerical Study and Validation of Melting and Solidification in PCM Embedded Heat Exchangers with Straight Tube

Latent heat thermal energy storage (LHTES) systems have shown great potential to enable reliable use of renewable energy and load shifting. LHTES offer high storage density and release energy at near constant temperature because of its use of phase change materials (PCMs). The cylindrical PCM heat exchangers (PCMHX) are one of the most used technologies due to their simplicity. Numerical models for such PCMHX enable engineers to estimate their performance for different design parameters and operating conditions without having to test them all. However, modeling the phase change phenomena can be challenging. To better understand the difficulties involving accurate modeling of PCMHX, a cylindrical latent storage unit filled with PCM and water as in-tube heat transfer fluid (HTF) is numerically investigated. This paper presents a study based on a 2D-axisymmetric model of a straight tube embedded in PCM in a cylindrical container. CFD is used to study the charging (melting) and discharging (solidification) phenomena. The models are validated against experimental and numerical data from the literature. The predicted local PCM temperature profile over time agrees within 2K compared to the experimental values. The paper also presents a simple method to estimate the melting and solidification phase change temperature range from limited data provided by PCM manufacturers.

PCM, melting, solidification, CFD, heat exchanger↗

A numerical study of three-dimensional diurnal variations within the thermosphere.

A thermosphere model with a realistic temperature profile is assumed. Heat conduction waves are introduced in addition to gravity waves. The temporal and spatial distribution of ion-neutral collisions is taken into account. However, the influence of viscosity waves is neglected. Viscosity-wave effects are simulated by an effective height-dependent collision number. Numerical calculations are conducted of the generation and propagation of two of the most important symmetric tidal waves at thermospheric heights. The influence of the solar EUV-heat upon the generation of the two tidal modes is investigated.

Volland, H.↗

Effect of radiometric errors on accuracy of temperature-profile measurement by spectral scanning using absorption-emission pyrometry

The spectral-scanning method may be used to determine the temperature profile of a jet- or rocket-engine exhaust stream by measurements of gas radiation and transmittance, at two or more wavelengths. A single, fixed line of sight is used, using immobile radiators outside of the gas stream, and there is no interference with the flow. At least two sets of measurements are made, each set consisting of the conventional three radiometric measurements of absorption-emission pyrometry, but each set is taken over a different spectral interval that gives different weight to the radiation from a different portion of the optical path. Thereby, discrimination is obtained with respect to location along the path. A given radiometric error causes an error in computed temperatures. The ratio between temperature error and radiometric error depends on profile shape, path length, temperature level, and strength of line absorption, and the absorption coefficient and its temperature dependency. These influence the choice of wavelengths, for any given gas. Conditions for minimum temperature error are derived. Numerical results are presented for a two-wavelength measurement on a family of profiles that may be expected in a practical case of hydrogen-oxygen combustion. Under favorable conditions, the fractional error in temperature approximates the fractional error in radiant-flux measurement.

Buchele, D. R.↗

Bayesian optimization of laser wakefield acceleration via spectral pulse shaping

In this paper, we investigate the effect of spectral pulse shaping of the laser driver on the performance of channel-guided, laser–plasma accelerators. The study was carried out with the assistance of Bayesian optimization using particle-in-cell simulations. We used a realistic plasma profile based on a novel optical-field-ionized channel technique with ionization injection and low, on-axis plasma densities to maximize the energy gain of the electron bunch trailing the laser. Spectral shaping allows us to modify the temporal profile of the laser driver while keeping the laser energy constant, affecting the acceleration and injection processes. In addition, we consider how modifying the plasma channel parameters may affect the target outputs. Given the complexity and breadth of the parameter space in question, we used numerical optimization to identify high-performing configurations. In particular, we found laser profiles with additional spectral content that, when used with optimal plasma channel parameters, result in charge content an order of magnitude higher than the baseline Gaussian case while also increasing the mean energy of the electron bunch.

Physics - Plasma physics↗

Instability of a supersonic shock free elliptic jet

This paper presents a comparison of the measured and the computed spatial stability properties of an aspect ratio 2 supersonic shock free elliptic jet. The shock free nature of the elliptic jet provides an ideal test of validity of modeling the large scale coherent structures in the initial mixing region of noncircular supersonic jets with linear hydrodynamic stability theory. Both aerodynamic and acoustic data were measured. The data are used to compute the mean velocity profiles and to provide a description of the spatial composition of pressure waves in the elliptic jet. A hybrid numerical scheme is applied to solve the Rayleigh problem governing the inviscid linear spatial stability of the jet. The measured mean velocity profiles are used to provide a qualitative model for the cross sectional geometry and the smooth velocity profiles used in the stability analysis. Computational results are presented for several modes of instability at two jet cross sections. The acoustic measurements show that a varicose instability is the jet's perferred mode of motion. The stability analysis predicts that the Strouhal number varies linearly as a function of axial distance in the jet's initial mixing region, which is in good qualitative agreement with previous measurements.

Baty, Roy S.↗

Special relativity effects for space-based coherent lidar experiments

There is a great need to develop a system that can measure accurately atmospheric wind profiles because an accurate data of wind profiles in the atmosphere constitutes single most input for reliable simulations of global climate numerical methods. Also such data helps us understand atmospheric circulation and climate dynamics better. Because of this need for accurate wind measurements, a space-based Laser Atmospheric Winds Sounder (LAWS) is being designed at MSFC to measure wind profiles in the lower atmosphere of the earth with an accuracy of 1 m/s at lower altitudes to 5m/s at higher altitudes. This system uses an orbiting spacecraft with a pulsed laser source and measures the Doppler shift between the transmitted and received frequencies to estimate the atmospheric wind velocities. If a significant return from the ground (sea) is possible, the spacecraft speed and height are estimated from it and these results and the Doppler shift are then used to estimate the wind velocities in the atmosphere. It is expected that at the proposed wavelengths, there will be enough backscatter from the aerosols but there may no be significant return from the ground. So a coherent (heterodyne) detection system is being proposed for signal processing because it can provide high signal to noise ratio and sensitivity and thus make the best use of low ground return. However, for a heterodyne detection scheme to provide the best results, it is important that the receiving aperture be aligned properly for the proposed wind sounder, this amounts to only a few microradians tolerance in alignment. It is suspected that the satellite motion relative to the ground may introduce errors in the order of a few microradians because of special relativity. Hence, the problem of laser scattering off a moving fixed target when the source and receiver are moving, which was not treated in the past in the literature, was analyzed in the following, using relativistic electrodynamics and applied to the case of the space-based coherent lidar, assuming flat ground. Here an interest in developing analytical expression for the location of the receiving point for the return with respect to the satellite, receiving angle and Doppler shift in frequency and amount of tip, all as measured in the satellite moving coordinate system and the diffuse scattering angle at the ground which does not require any compensation. All the three cases of retro-reflection, specular reflection and diffuse scattering by the ground should be treated though retro-reflection and diffuse scattering are more important.

Raogudimetla, V. S.↗

An analysis of a charring ablator with thermal nonequilibrium, chemical kinetics, and mass transfer

The differential equations governing the transient response of a one-dimensional ablative thermal protection system are presented for thermal nonequilibrium between the pyrolysis gases and the char layer and with finite rate chemical reactions occurring. The system consists of three layers (the char layer, the uncharred layer, and an optical insulation layer) with concentrated heat sinks at the back surface and between the second and third layers. The equations are solved numerically by using a modified implicit finite difference scheme to obtain solutions for the thickness of the charred and uncharred layers, surface recession and pyrolysis rates, solid temperatures, porosity profiles, and profiles of pyrolysis-gas temperature, pressure, composition, and flow rate. Good agreement is obtained between numerical results and exact solutions for a number of simplified cases. The complete numerical analysis is used to obtain solutions for an ablative system subjected to a constant heating environment. Effects of thermal, chemical, and mass transfer processes are shown.

Clark, R. K.↗

Design procedures for compressor blades

The conventional methods for the design of the blades in the case of axial turbomachines are considered, taking into account difficulties concerning the determination of optimal blade profiles. These difficulties have been partly overcome as a consequence of the introduction of new numerical methods during the last few years. It is pointed out that, in the case of the subsonic range, a new procedure is now available for the determination of the form of blade profile on the basis of a given velocity distribution on the profile surface. The search for a profile form with favorable characteristics is consequently transformed into a search for a favorable velocity or pressure distribution on the blade. The distribution of velocities depends to a large degree on the characteristics of the profile boundary layers. The considered concept is not new. However, its practical implementation has only recently become possible. The employment of the new design procedure is illustrated with the aid of an example involving a concrete design problem.

Starken, H.↗

Numerical study of boundary layer interaction with shocks: Method and code validation

A major problem in modeling of turbulent supersonic flows is the correct assessment of viscous-inviscid interaction problems. Of particular interest is the interaction of boundary layers with shocks. Present turbulence models give in most cases unsatisfactory results in the region of rapid distortion and in the separation region (if one is present) in particular with regard to mean flow profiles and turbulence quantities. The objective of the present work is the direct numerical simulation of shock boundary layer interaction. This report summarizes the first phase during which a numerical method suitable for this problem has been developed and a computer code has been written and tested.

Adams, Nikolaus A.↗

A prediction method for velocity and temperature profiles in a two-dimensional nominally steady turbulent boundary layer

This paper describes a recently developed boundary-layer prediction method for a variable property compressible flow, in which heat transfer takes place primarily by forced convection and for which the mainstream Mach number is small. The leading order terms, in asymptotic expansions for large Reynolds numbers, are obtained for the mean velocity and temperature distribution in both the inner and outer layer of the turbulent boundary layer. Closure in the inner layer is achieved using an analytical model for the mean profiles which is based on the observed coherent structure of the time-dependent inner layer flow. For the outer layer, simple eddy viscosity and conductivity models are developed without recourse to the Reynolds analogy. In the prediction method a numerical solution of the outer layer equations is matched to the analytical inner layer profiles as the computation procedes downstream. Calculations are presented for a range of adverse and favorable pressure gradient flows and the predicted results compare well with existing data.

Weigand, G. G.↗

Verifying raytracing/Fokker–Planck lower-hybrid current drive predictions with self-consistent full-wave/Fokker–Planck simulations

Raytracing/Fokker–Planck (FP) simulations used to model lower-hybrid current drive (LHCD) often fail to reproduce experimental results, particularly when LHCD is weakly damped. A proposed reason for this discrepancy is the lack of ‘full-wave’ effects, such as diffraction and interference, in raytracing simulations and the breakdown of the raytracing approximation. Previous studies of LHCD using non-Maxwellian full-wave/FP simulations have been performed, but these simulations were not self-consistent and enforced power conservation between the FP and full-wave code using a numerical rescaling factor. Here, we have created a fully self-consistent full-wave/FP model for LHCD that is automatically power conserving. This was accomplished by coupling an overhauled version of the non-Maxwellian TORLH full-wave solver and the CQL3D FP code using the Integrated Plasma Simulator. We performed converged full-wave/FP simulations of Alcator C-Mod discharges and compared them with raytracing. We found that excellent agreement in the power deposition profiles from raytracing and TORLH could be obtained, however, TORLH had somewhat lower current drive efficiency and broader power deposition profiles in some cases. This discrepancy appears to be a result of numerical limitations present in the TORLH model and a small amount of diffractional broadening of the TORLH wave spectrum. Our results suggest full-wave simulation of LHCD is likely not necessary as diffraction and interference represented only a small correction that could not account for the differences between simulations and experiment.

Physics↗

A numerical and experimental analysis of reactor performance and deposition rates for CVD on monofilaments

The computational fluid dynamics (CFD) code FLUENT is adopted to simulate a cylindrical upflow reactor designed for chemical vapor deposition (CVD) on monofilaments. Equilibrium temperature profiles along the fiber and quartz reactor wall are experimentally measured and used as boundary conditions in numerical simulations. Two-dimensional axisymmetric flow and temperature fields are calculated for hydrogen and argon; the effect of free convection is assessed. The gas and surface chemistry is included for predicting silicon deposition from silane. The model predictions are compared with experimentally measured silicon CVD rates. Inferences are made for optimum conditions to obtain uniformity.

Gokoglu, S. A.↗

Observations of neutral wind profiles between 115- and 175-km altitude in the dayside auroral oval

Lithium trail neutral wind measurements have been made in the morning dayside auroral oval as a function of altitude. During all four experiments the winds displayed a large altitude shear. In the most extreme case the wind at 135-km altitude was toward the south at 210 m/s, and at 185-km altitude toward the northwest at 170 m/s. A one-dimensional numerical solution to the Navier Stokes equation including an initial steady-state wind profile, viscosity, pressure gradient, Lorenz force and the Coriolis force was made in order to study the effect of auroral processes on the neutral atmosphere. Since electric fields and plasma density were measured, the Lorenz force could be modeled with reasonable accuracy. The pressure gradient was determined in an iterative way as a best fit to the observations. Within the limitations of a one-dimensional model the numerical solution is in reasonable agreement with one of the observations, which is applicable to such a modeling, and indicates that both Lorenz forces and pressure gradients due to auroral processes are required to explain the observation.

Pereira, E.↗

Optimal design of solidification processes

An optimal design algorithm is presented for the analysis of general solidification processes, and is demonstrated for the growth of GaAs crystals in a Bridgman furnace. The system is optimal in the sense that the prespecified temperature distribution in the solidifying materials is obtained to maximize product quality. The optimization uses traditional numerical programming techniques which require the evaluation of cost and constraint functions and their sensitivities. The finite element method is incorporated to analyze the crystal solidification problem, evaluate the cost and constraint functions, and compute the sensitivities. These techniques are demonstrated in the crystal growth application by determining an optimal furnace wall temperature distribution to obtain the desired temperature profile in the crystal, and hence to maximize the crystal's quality. Several numerical optimization algorithms are studied to determine the proper convergence criteria, effective 1-D search strategies, appropriate forms of the cost and constraint functions, etc. In particular, we incorporate the conjugate gradient and quasi-Newton methods for unconstrained problems. The efficiency and effectiveness of each algorithm is presented in the example problem.

Dantzig, Jonathan A.↗

Moisture Advection and Fresh Water Flux over Oceans

Moisture transport in the atmosphere is one of the most significant components in the hydrological cycle. Under stationary condition, ocean surface fresh water flux, which is the difference between precipitation (P) and evaporation (E), is balanced by the divergence of column-integrated moisture transport (IMT) in the atmosphere. Characterizing accurately a global picture of IMT from observation is a difficult task. It requires measurements of vertical profiles for wind vector and humidity. More specifically, IMT can be defined as the integration in pressure coordinates the product of specific humidity q and wind vector u, where g is the gravitational acceleration, and p, is the atmospheric pressure at ocean surface.In this study, a statistical relationship is derived between u, and u(sub)s using data from numerical weather prediction model. The relationship is then validated using surface and vertical profile from radiosonde data, before applied to spacebased measurements.

scatterometers↗

Large‐Area, High‐Numerical‐Aperture, Freeform Metasurfaces

Abstract Nanophotonic devices are optical platforms capable of unprecedented wavefront control. To push the limits of experimental device performance, scalable design methodologies that combine the simplicity and fabricability of conventional design paradigms with the extended capabilities of freeform optimization are required. A novel gradient‐based design framework for large‐area freeform metasurfaces is introduced in which nonlocal interactions between simply shaped nanostructures, placed on an irregular lattice, are tailored to produce high‐order hybridized modes that support customizable large‐angle scattering profiles. Utilizing this approach, multifunctional super‐dispersive metalenses are designed and experimentally demonstrated. The approach to high‐numerical‐aperture radial metalenses capable of diffraction limited focusing and the generation of donut‐shaped point spread functions is also extended. It is anticipated that these concepts will have utility in super‐resolution microscopy, particle trapping, additive manufacturing, and metrology applications that require ultra‐high numerical apertures.

Optics↗

Experiment and numerical simulation of RF heating in the Tandem Mirror plasma propulsion device

Recent results of experimental and theoretical studies of ICRF heating of the plasma in the Tendem Mirror rocket are presented. The radial and axial profile of the magnetic field of the wave have been measured, and the data agree with the results from numerical simulation. One very important new finding is that the wave damped when it approached the resonance plane. This is a strong indication of beach heating effect and that RF power is absorbed by the ions in the plasma as expected. This power absorption phenomenon was also further confirmed by the refined analytical study of the wave propagation in a slab model. The electron density in the central cell has been measured by a microwave interferometer. The existence of two types of discharge similar to those of H-alpha emission have also been observed from the density measurements. A preliminary design of a flight system for a 100-day Mars transit has been undertaken. The specific weight for a 10 MW rocket is only 0.04 kg/kW.

Yang, T. F.↗