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 559 records · Page 31

Numerical simulation of advection fog formation on multi-disperse aerosols due to combustion-related pollutants

The effects of multi-disperse distribution of the aerosol population are presented. Single component and multi-component aerosol species on the condensation/nucleation processes which affect the reduction in visibility are described. The aerosol population with a high particle concentration provided more favorable conditions for the formation of a denser fog than the aerosol population with a greater particle size distribution when the value of the mass concentration of the aerosols was kept constant. The results were used as numerical predictions of fog formation. Two dimensional observations in horizontal and vertical coordinates, together with time-dependent measurements were needed as initial values for the following physical parameters: (1)wind profiles; (2) temperature profiles; (3) humidity profiles; (4) mass concentration of aerosol particles; (5) particle size distribution of aerosols; and (6) chemical composition of aerosols. Formation and dissipation of advection fog, thus, can be forecasted numerically by introducing initial values obtained from the observations.

Hung, R. J.↗

Turbine rotor/stator flowfield analysis

A numerical study of the unsteady aerodynamic and thermal environment associated with axial turbine stages is presented. Computations are performed using a modification of the ROTOR1 rotor/stator interaction code. Two different turbine states are analyzed: the first state of the United Technologies Research Center large scale rotating rig and the first state of the Space Shuttle main engine (SSME) high pressure fuel turbopump. Time-averaged blade midspan pressure and heat transfer profiles are calculated using the following different surface boundary conditions: adiabatic wall, prescribed wall temperature, and prescribed heat flux. Numerical solutions for the large scale rotating rig are compared with experimental data. Unsteady pressure envelopes are also presented for each geometry. In addition, instantaneous contours are plotted for the SSME configuration.

Griffin, Lisa↗

The Vertical Error Characteristics of GOES-derived Winds: Description and Impact on Numerical Weather Prediction

Errors in the height assignment of some satellite-derived winds exist because the satellites sense radiation emitted from a finite layer of the atmosphere rather than a specific level. Potential problems in data assimilation may arise because the motion of a measured layer is often represented by a single-level value. In this research, cloud and water vapor motion winds that are derived from the Geostationary Operational Environmental Satellites (GOES winds) are compared to collocated rawinsonde observations (RAOBs). An important aspect of this work is that in addition to comparisons at each assigned height, the GOES winds are compared to the entire profile of the collocated RAOB data to determine the vertical error characteristics of the GOES winds. The impact of these results on numerical weather prediction is then investigated. The comparisons at individual vector height assignments indicate that the error of the GOES winds range from approx. 3 to 10 m/s and generally increase with height. However, if taken as a percentage of the total wind speed, accuracy is better at upper levels. As expected, comparisons with the entire profile of the collocated RAOBs indicate that clear-air water vapor winds represent deeper layers than do either infrared or water vapor cloud-tracked winds. This is because in cloud-free regions the signal from water vapor features may result from emittance over a thicker layer. To further investigate characteristics of the clear-air water vapor winds, they are stratified into two categories that are dependent on the depth of the layer represented by the vector. It is found that if the vertical gradient of moisture is smooth and uniform from near the height assignment upwards, the clear-air water vapor wind tends to represent a relatively deep layer. The information from the comparisons is then used in numerical model simulations of two separate events to determine the forecast impacts. Four simulations are performed for each case: 1) A control simulation that assimilates no satellite wind data, 2) assimilation of all GOES winds according to their assigned single level height, 3) assimilation of all GOES winds spread over multiple levels, and 4) assimilation of all GOES winds spread over multiple levels, but with variations in the vertical influence of clear-air water vapor winds based on the moisture profile in the model. In the first case, a strong mid-latitude cyclone is present and the use of the satellite data results in improved storm tracks during the initial approx. 36 h forecast period. This is because the satellite data improves the analysis of the environment into which the storm progresses. Statistics for mean wind vector and height differences show that, with the exception of the height field at later times in the first case, the use of GOES winds improves the simulation with time. The simulation results suggest that it is beneficial to spread the GOES wind information over multiple levels, particularly when the moisture profile is used to define the vertical influence.

Rao, P. Anil↗

Jet thermalization in QCD kinetic theory

We perform numerical studies in the framework of QCD kinetic theory to investigate the energy and angular profiles of a high energy parton — as a proxy for a jet produced in heavy ion collisions — passing through a Quark-Gluon Plasma (QGP). We find that the fast parton loses energy to the plasma mainly via a radiative turbulent quark and gluon cascade that transports energy locally from the jet down to the temperature scale where dissipation takes place. In this first stage of the system time evolution, the angular structure of the turbulent cascade is found to be relatively collimated. However, when the lost energy reaches the plasma temperature it is rapidly transported to large angles w.r.t. the jet axis and thermalizes. We investigate the contribution of the soft jet constituents to the total jet energy. We show that for jet opening angles of about 0.3 rad or smaller, the effect is negligible. Conversely, larger opening angles become more and more sensitive to the thermal component of the jet and thus to medium response. Our result showcases the importance of the jet cone size in mitigating or enhancing the details of dissipation in jet quenching observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

FreeMHD: Validation and verification of the open-source, multi-domain, multi-phase solver for electrically conductive flows

The extreme heat fluxes in the divertor region of tokamaks may require an alternative to solid plasma-facing components, for the extraction of heat and the protection of the surrounding walls. Flowing liquid metals are proposed as an alternative, but raise additional challenges that require investigation and numerical simulations. Free surface designs are desirable for plasma-facing components, but steady flow profiles and surface stability must be ensured to limit undesirable interactions with the plasma. Previous studies have mainly used steady-state, 2D, or simplified models for internal flows and have not been able to adequately model free-surface liquid metal (LM) experiments. Therefore, FreeMHD has been recently developed as an open-source magnetohydrodynamics (MHD) solver for free-surface electrically conductive flows subject to a strong external magnetic field. The FreeMHD solver computes incompressible free-surface flows with multi-region coupling for the investigation of MHD phenomena involving fluid and solid domains. The model utilizes the finite-volume OpenFOAM framework under the low magnetic Reynolds number approximation. FreeMHD is validated using analytical solutions for the velocity profiles of closed channel flows with various Hartmann numbers and wall conductance ratios. Next, experimental measurements are then used to verify FreeMHD, through a series of cases involving dam breaking, 3D magnetic fields, and free-surface LM flows. These results demonstrate that FreeMHD is a reliable tool for the design of LM systems under free surface conditions at the reactor scale. Furthermore, it is flexible, computationally inexpensive, and can be used to solve fully 3D transient MHD flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Regime Characterization of Offshore Wind Resource Using Unsupervised Learning

Predictability of wind resource conditions is critical for offshore wind design and operations. While many studies of extreme wind conditions focus on specific events such as low-level jets or ramps, these rely on threshold definitions that limit generality. Here we present a data-driven framework that combines principal component analysis (PCA), self-organizing maps (SOM), and k-means clustering to classify wind resource conditions as typical and anomalous from climatological data. Anomalies are defined not by fixed thresholds but by flagging samples located far from SOM node centers inside the baseline SOM structure. This reframes extremes as rare ebents and hence, likely difficult to anticipate by numerical weather prediction models. We applied this approach to 23 years (2000–2022) of hourly profiles from the NOW-23 hindcast model at the Humboldt Wind Energy Area. Classification is conducted on a feature space consisting of 10 m wind speed and direction, bulk shear and veer across 30–270 m, and a low-level jet index. Dimensionality reduction is achieved through PC. A 2 × 3 OM lattice trained on the PCA vectors identified six baseline regimes spanning weak to strong flow states. High quantization-error profiles are identified and re-clustered into four anomalous regimes. The baseline regimes exhibited clear seasonal and diurnal cycles. Meanwhile, the anomalous regimes represented <10 % of all hours but showed distinct combinations of speed, shear, and veer, when compared to the baseline regimes. Anomalous regimes are typically short-lived (~few hours), yet their transitions can lead to hub-height wind changes of −18 to +9 m s -1 . For a representative 15 MW turbine, these shifts imply rapid swings in capacity factor from near-full output to negligible generation. Validation with lidar buoy data showed 51% agreement in SOM labels across ~6,000 overlapping hours, with most mismatches confined to adjacent speed classes. HRRR comparisons further revealed that anomalous regimes were disproportionately associated with forecast biases exceeding 5 m s -1 . Together, these results reframe extremes in offshore wind from absolute maxima or minima to weather states that are difficult to anticipate from models.

17 WIND ENERGY↗

Model-Free Dynamic Voltage Control of Distributed Energy Resource (DER)-Based Microgrids

In this paper, we present a new control technique for sustaining dynamic voltage stability by effective reactive power control and coordination of distributed energy resources (DERs) in microgrids. The proposed control technique is based on model-free control (MFC), which has shown successful operation and improved performance in different domains and applications. This paper presents its first use in the voltage stability of a microgrid setting employing multiple synchronous generator (SG)-based and power electronic (PE)-based DERs. MFC is a computationally efficient, data-driven control technique that does not require modelling of the different components and disturbances in the power system. It is utilized as an online controller to achieve the dynamic voltage stability of a microgrid system under different disturbances and fault conditions. A 21-bus microgrid system fed by multiple DERs is considered as a case study and the overall dynamic voltage stability is investigated using time-domain dynamic simulations. Numerical results show that the proposed MFC provides improvements on the dynamic load bus voltage profiles and requires less computational time as compared to the traditional enhanced microgrid voltage stabilizer (EMGVS) scheme. Due to its simplicity and low computational requirement, MFC can be easily implemented in resource-constrained computing devices such as smart inverters.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Current conduction in junction gate field effect transistors

The internal physical mechanism that governs the current conduction in junction-gate field effect transistors is studied. A numerical method of analyzing the devices with different length-to-width ratios and doping profiles is developed. This method takes into account the two dimensional character of the electric field and the field dependent mobility. Application of the method to various device models shows that the channel width and the carrier concentration in the conductive channel decrease with increasing drain-to-source voltage for conventional devices. It also shows larger differential drain conductances for shorter devices when the drift velocity is not saturated. The interaction of the source and the drain gives the carrier accumulation in the channel which leads to the space-charge-limited current flow. The important parameters for the space-charge-limited current flow are found to be the L/L sub DE ratio and the crossover voltage.

Kim, C.↗

Removal of spurious modes encountered in solving stability problems by spectral methods

A technique based on the Galerkin approximation is developed to remove spurious roots arising when Chebyshev spectral methods are used to solve eigenvalue problems in hydrodynamic stability. The derivation of Galerkin-Chebyshev approximations is explained, and numerical results for the Orr-Sommerfeld equations of plane Poiseuille flow and a Blasius profile are presented in tables and compared with those obtained by the method of Zebib (1984). It is pointed out that the present method does not increase the size of the algebraic system to be solved.

Zebib, Abdelfattah↗

Numerical methods for one-dimensional reaction-diffusion equations arising in combustion theory

A review of numerical methods for one-dimensional reaction-diffusion equations arising in combustion theory is presented. The methods reviewed include explicit, implicit, quasi-linearization, time linearization, operator-splitting, random walk and finite-element techniques and methods of lines. Adaptive and nonadaptive procedures are also reviewed. These techniques are applied first to solve two model problems which have exact traveling wave solutions with which the numerical results can be compared. This comparison is performed in terms of both the wave profile and computed wave speed. It is shown that the computed wave speed is not a good indicator of the accuracy of a particular method. A fourth-order time-linearized, Hermitian compact operator technique is found to be the most accurate method for a variety of time and space sizes.

Ramos, J. I.↗

Axial strain insensitivity of weakly guiding optical fibers

A numerical model has been developed to calculate the modal phase shift of circular step index profile weakly guiding fibers under axial strain. Whenever an optical fiber is under stress, the optical path length, the index of refraction, and the propagation constants of each mode change. In consequence, the phase of each mode is also modified. A relationship for the modal phase shift is presented. This relation is applied to both single mode and two-mode fibers in order to determine the sensitivity characteristics of strained fibers. It was found that the phase shift is strongly dependent on the core refractive index, n(co). It was also found that it is possible to design fibers which are insensitive to axial strain. Practical applications of strain insensitive fibers are discussed.

Egalon, Claudio O.↗

Radiation Shielding of the Stagnation Region by Transpiration of an Opaque Gas

The laminar compressible boundary layer in the two-dimensional and axisymmetric stagnation regions has been analyzed to show the effects of the injection of a radiation absorbing foreign gas on an incident radiation field, and on the enthalpy profiles across the boundary layer. Total heat transfer to the stagnation region is evaluated for numerous cases and the results are compared with the no shielding case. Required absorption properties of the foreign gas are determined and compared with properties of known gases.

Howe, John Thomas↗

Lunar and Planetary Science XXXV: Martian Aeolian and Mass Wasting Processes: Blowing and Flowing

The session Martian Aeolian and Mass Wasting Processes: BLowing and Flowing included the following topics: 1) Three Decades of Martian Surface Changes; 2) Thermophysical Properties of Isidis Basin, Mars; 3) Intracrater Material in Eastern Arabia Terra: THEMIS, MOC, and MOLA Analysis of Wind-blown Deposits and Possible High-Inertia Source Material; 4) Thermal Properties of Sand from TES and THEMIS: Do Martian Dunes Make a Good Control for Thermal Inertia Calculations? 5) A Comparative Analysis of Barchan Dunes in the Intra-Crater Dune Fields and the North Polar Sand Sea; 6) Diluvial Dunes in Athabasca Valles, Mars: Morphology, Modeling and Implications; 7) Surface Profiling of Natural Dust Devils; 8) Martian Dust Devil Tracks: Inferred Directions of Movement; 9) Numerical Simulations of Anastomosing Slope Streaks on Mars; 10) Young Fans in an Equatorial Crater in Xanthe Terra, Mars; 11) Large Well-exposed Alluvual Fans in Deep Late-Noachian Craters; 12) New Evidence for the Formation of Large Landslides on Mars; and 13) What Can We Learn from the Ages of Valles Marineris Landslides on Martian Impact History?

Source record↗

The DUV Stability of Superlattice-Doped CMOS Detector Arrays

In this paper, we present experimental results and band structure calculations that illuminate the unique properties of superlattice-doped detectors. Numerical band structure calculations are presented to analyze the dependencies of surface passivation on dopant profiles and interface trap densities (Figure 3). Experiments and calculations show that quantum-engineered surfaces, grown at JPL by low temperature molecular beam epitaxy, achieve a qualitative as well as quantitative uniqueness in their near-immunity to high densities of surface and interface traps.

doped silicon surfaces↗

Analysis of Summer-Time Ozone and Precursor Species in the Southeast United States

Ozone (O3) is a greenhouse gas and toxic pollutant which plays a major role in air quality and atmospheric chemistry. The understanding and ability to model the horizontal and vertical structure of O3 mixing ratios is difficult due to the complex formation/destruction processes and transport pathways that cause large variability of O3. The Environmental Protection Agency has National Ambient Air Quality Standards for O3 set at 75 ppb with future standards proposed to be as low as 65 ppb. These lower values emphasize the need to better understand/simulate the transport processes, emission sources, and chemical processes controlling precursor species (e.g., NOx, VOCs, and CO) which influence O3 mixing ratios. The uncertainty of these controlling variables is particularly large in the southeast United States (US) which is a region impacted by multiple different emission sources of precursor species (anthropogenic and biogenic) and transport processes resulting in complex spatio-temporal O3 patterns. During this work we will evaluate O3 and precursor species in the southeast US applying models, ground-based and airborne in situ data, and lidar observations. In the summer of 2013, the UAH O3 Differential Absorption Lidar (DIAL) (part of the Tropospheric Ozone Lidar Network (TOLNet)) measured vertical O3 profiles from the surface up to approximately 12 km. During this period, the lidar observed numerous periods of dynamic temporal and vertical O3 structures. In order to determine the sources/processes impacting these O3 mixing ratios we will apply the CTM GEOS-Chem (v9-02) at a 0.25 deg x 0.3125 deg resolution. Using in situ ground-based (e.g., SEARCH Network, CASTNET), airborne (e.g., NOAA WP-3D - SENEX 2013, DC-8 - SEAC4RS), and TOLNet lidar data we will first evaluate the model to determine the capability of GEOS-Chem to simulate the spatio-temporal variability of O3 in the southeast US. Secondly, we will perform model sensitivity studies in order to quantify which emission sources (e.g., anthropogenic, biogenic, lighting, wildfire) and transport processes (e.g., stratospheric, long-range, local scale) are contributing to these TOLNet-observed dynamic O3 patterns. Results from the evaluation of the model and the study of sources/processes impacting observed O3 mixing ratios will be presented.

Optical Radar↗

Coordinated Optimization and Control of Residential Solid State Power Substations in Electrical Distribution Network

This paper proposes a computationally efficient mixed integer linear programming (MILP) model for the coordinated optimization of solid-state power substations (SSPSs) in a feeder considering the full unbalanced three-phase structure of the distribution grid and its characteristics. The proposed model determines the optimal real and reactive power of each SSPS at the point of common coupling (PCC) that minimizes the operating cost and maximizes the system performance, e.g., voltage regulation and phase balancing. To improve the computational efficiency, an inscribed octagon is introduced to approximate the quadratic capacity constraints of components. Numerical simulation results show the effectiveness of the proposed model and significant improvements in voltage profiles and power imbalance between phases.

Liu, Guodong↗

Turbulent mixing in supercritical jets: effect of compressibility factor and inflow condition

Fuel injection and turbulent mixing at supercritical pressures determines ignition and combustion in numerous engineering applications. Flow evolution under such conditions is characterized by strong non-linear coupling between dynamics, transport coefficients, and thermodynamics. Experimental studies observe that the jets injected at supercritical pressures exhibit significantly different dynamics from the jets at subcritical conditions, owing to the lack of distinct liquid and gas phases in supercritical state. Thus, the averaged flow quantities such as the potential core length, jet spatial growth rate and velocity decay profiles differ in the two conditions, resulting in different mixed-fluid distributions. In this study, turbulent jet direct numerical simulations (DNS) are performed to examine the variations in flow statistics between injection of Nitrogen (N₂) in Nitrogen (N₂) at both subcritical (perfect-gas) and supercritical conditions. In all cases, isothermal round jets at Reynolds number (Re_{D}), based on jet diameter (D) and jet orifice velocity (U₀), of 5000 are considered. For mixing analyses, a passive scalar transported with the flow is examined.

Sharan, Nek↗

Numerical studies of motion of vortex filaments - Implementing the asymptotic analysis

A computational code is developed for the integro-differential equations governing the motion of the centerlines of vortex filaments submerged in a background potential flow. These equations, which are derived from the method of matched asymptotic analysis, include the effect of the decaying large-magnitude circumferential and axial velocity components in the vortical cores. Numerical examples are presented to assess the effect of a large axial velocity and that of nonsimilar initial profiles in the vortical cores. The initial configurations of the filaments are chosen so as to fulfill the basic assumption of the asymptotic analysis, which is that the effective vortical core size is much smaller than all the other length scales in the flowfield, e.g., the radius of curvature and the interfilament distance. The computations are continued until the basic assumption is no longer valid, that is when the merging or intersection of filaments has begun. A classification of the various types of local or global merging or intersection of filaments is made and demonstrated by numerical examples. It is then shown that the asymptotic solution not only provides the initial data but also can be used to formulate the appropriate boundary conditions for the numerical solution of a merged region.

Liu, C. H.↗