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 649 records · Page 36

Forced magnetic reconnection in a plasma sheet with localized resistivity profile excited by lower hybrid drift type instability

A forced magnetic reconnection process with a temporal evolution of resistivity is studied for a plasma sheet with a nonuniform resistivity profile based on the nonlocal mode structure of the lower hybrid drift type instability. The growth rate of the mode found is almost independent of the resistivity at the neutral sheet, but depends on the resistivity of the region of maximum density gradient away from the neutral sheet. This is studied by using both a nonlinear numerical MHD simulation and a linear theory. The mode may be relevant to the prevalent theoretical concept of MHD reconnection and the localized anomalous resistivity profile based on the lower hybrid drift instability.

Hoshino, M.↗

Fine Root and Soil Organic Carbon Depth Distributions are Inversely Related Across Fertility and Rainfall Gradients in Lowland Tropical Forests

Humid tropical forests contain some of the largest soil organic carbon (SOC) stocks on Earth. Much of this SOC is in subsoil, yet variation in the distribution of SOC through the soil profile remains poorly characterized across tropical forests. We used a correlative approach to quantify relationships among depth distributions of SOC, fine root biomass, nutrients and texture to 1 m depths across 43 lowland tropical forests in Panama. The sites span rainfall and soil fertility gradients, and these are largely uncorrelated for these sites. We used fitted b parameters to characterize depth distributions, where b is a numerical index based on an asymptotic relationship, such that larger b values indicate greater concentrations of root biomass or SOC at depth in the profile. Root b values ranged from 0.82 to 0.95 and were best predicted by soil pH and extractable potassium (K) stocks. For example, the three most acidic (pH < 4) and Kpoor (< 20 g K m -2 ) soils contained 76 ± 5% of fine root biomass from 0 to 10 cm depth, while the three least acidic (pH > 6.0) and most K-rich (> 50 g K m -2 ) soils contained only 41 ± 9% of fine root biomass at this depth. Root b and SOC b values were inversely related, such that a large fine root biomass in surface soils corresponded to large SOC stocks in subsoils (50–100 cm). SOC b values were best predicted by soil pH and base cation stocks, with the three most base-poor soils containing 34 ± 8% of SOC from 50 to 100 cm depth, and the three most base-rich soils containing just 9 ± 2% of SOC at this depth. Nutrient depth distributions were not related to Root b or SOC b values. Here, these data show that large surface root biomass stocks are associated with large subsoil C stocks in strongly weathered tropical soils. Further studies are required to evaluate why this occurs, and whether changes in surface root biomass, as may occur with global change, could in turn influence SOC storage in tropical forest subsoils.

54 ENVIRONMENTAL SCIENCES↗

Analysis and interpretation of lidar observations of the stratospheric aerosol

Data obtained with a 48 in. telescope lidar system are compared with results obtained using a one-dimensional stratospheric aerosol model to analyze various microphysical processes influencing the formation of this aerosol. Special attention is given to the following problems: (1) how lidar data can help determine the composition of the aerosol particles and (2) how the layer corresponds to temperature profile variations. The lidar record during the period 1974 to 1979 shows a considerable decrease of the peak value of the backscatter ratio. Seasonal variations in the aerosol layer and a gradual decrease in stratospheric loading are observed. The aerosol model simulates a background stratospheric aerosol layer, and it predicts stratospheric aerosol concentrations and compositions. Numerical experiments are carried out by using the model and by comparing the theoretical results with the experimentally obtained lidar record. Comparisons show that the backscatter profile is consistent with the composition when the particles are sulfuric acid and water; it is not consistent with an ammonium sulfate composition. It is shown that the backscatter ratio is not sensitive to the composition or stratospheric loading of condensation nuclei such as meteoritic debris.

Hamill, P.↗

Analytical and Numerical Solution for a Solidifying Liquid Alloy Slab

Numerical and analytical solutions are presented for the temperature and concentration distributions during the solidification of a binary liquid alloy slab. The slab is taken to be of a finite depth but infinite in the horizontal direction. The solidification process is started by withdrawing a fixed amount of heat from the lower surface of the slab. The upper surface of the slab is subjected to both radiation and convective conditions. The solution gives the concentration and temperature profiles and the interface position as a function of time. Due to the smallness of the mass diffusion coefficient in the solid, the numerical solution method breaks down whenever the ratio of the diffusivities in the solid and the liquid falls below a certain value. An analytical method is developed which gives accurate solution for any value of the diffusivity ratio.

Antar, B. N.↗

Radiative Heating Rate Profiles over the Southeast Atlantic Ocean during the 2016 and 2017 Biomass Burning Seasons

Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. Further uncertainties in the cloud structure arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean where biomass burning aerosol is transported from the African continent. As the aerosol plume progresses across the southeast Atlantic Ocean, radiative heating within the aerosol layer has the potential to alter the thermodynamic environment and therefore the cloud structure; however, limited work has been done to quantify this along the trajectory of the aerosol plume in the region. The deployment of the First Atmospheric Radiation Measurement Mobile Facility in support of the Layered Atlantic Smoke Interactions with Clouds field campaign provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96 °S, 14.35 °W). Using observed profiles of temperature, humidity, and clouds from the field campaign, alongside aerosol optical properties from MERRA-2 as input for the Rapid Radiation Transfer Model, profiles of the radiative heating rate due to aerosols and clouds were computed. Radiative heating is also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory Model. Idealized experiments using the Rapid Radiation Transfer Model with and without aerosols and a range of values for the single scattering albedo demonstrate that shortwave heating within the aerosol layer above Ascension Island can locally range between 2 and 8 K per day depending on the aerosol optical properties, though impacts of the aerosol can be felt elsewhere in the atmospheric column. When considered under clear conditions, the aerosol has a cooling effect at the TOA and based on the observed cloud properties at Ascension Island, the cloud albedo is not large enough to overcome this. Shortwave radiative heating due to biomass burning aerosol is not balanced by additional longwave cooling, and the net radiative impact results in a stabilization of the lower troposphere. However, these results are extremely sensitive to the single scattering albedo assumptions in models.

Allison B. Marquardt Collow↗

An Observational Evaluation of RKW Theory over the U.S. Southern Great Plains

The theory of Rotunno et al. (“RKW” theory) addresses the behavior of squall-line cold pools in vertically sheared flows. It predicts that, within a given thermodynamic environment, a balance between baroclinic vorticity generation by the cold pool and low-level environmental vertical wind shear induces an upright updraft along the gust front that maximizes the initiation of new convective cells. Although this theory has been evaluated numerically, its applicability to observed systems remains unclear and is limited by a lack of critical measurements, including high-frequency thermodynamic and wind profiles across the gust front. Herein, observations from the Atmospheric Radiation Measurement Southern Great Plains (ARM-SGP) observatory near Lamont, Oklahoma, are used to evaluate RKW theory for 10 well-observed squall lines over a 11-yr period. For this evaluation, RKW parameters including cold-pool intensity (c), low-level ambient, line-normal vertical shear (ΔV n ), subcloud and cloud-layer updraft tilts, and multiple measures of system intensity are estimated. Furthermore, the c estimates rely on thermodynamic retrievals from the Atmosphere Emitted Radiance Interferometer (AERI), which are uncertain but verify reasonably well against independent observations. As predicted by the theory, for c/ΔV n ≥ 1, c/ΔV n correlates positively with updraft tilt and negatively with system intensity, but these results are not always statistically significant and are also sensitive to the method by which ΔV n is evaluated. Specifically, ΔV n evaluations that extend above the cold-pool top yield greater consistency with RKW predictions. Also, some measures of intensity correlate more strongly with standard moist instability metrics than with RKW parameters.

Cold pools↗

Radiative heating rate profiles over the southeast Atlantic Ocean during the 2016 and 2017 biomass burning seasons

Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. Further uncertainties in the cloud structure arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean, where biomass burning aerosol is transported from the African continent. As the aerosol plume progresses across the southeast Atlantic Ocean, radiative heating within the aerosol layer has the potential to alter the thermodynamic environment and therefore the cloud structure; however, limited work has been done to quantify this along the trajectory of the aerosol plume in the region. The deployment of the first Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF1) in support of the Layered Atlantic Smoke Interactions with Clouds field campaign provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96° S, 14.35° W). Using observed profiles of temperature, humidity, and clouds from the field campaign alongside aerosol optical properties from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM), profiles of the radiative heating rate due to aerosols and clouds were computed. Radiative heating is also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. Idealized experiments using the RRTM with and without aerosols and a range of values for the single-scattering albedo (SSA) demonstrate that shortwave (SW) heating within the aerosol layer above Ascension Island can locally range between 2 and 8 K d -1 depending on the aerosol optical properties, though impacts of the aerosol can be felt elsewhere in the atmospheric column. When considered under clear conditions, the aerosol has a cooling effect at the TOA, and based on the observed cloud properties at Ascension Island, the cloud albedo is not large enough to overcome this. Shortwave radiative heating due to biomass burning aerosol is not balanced by additional longwave (LW) cooling, and the net radiative impact results in a stabilization of the lower troposphere. However, these results are extremely sensitive to the single-scattering albedo assumptions in models.

54 ENVIRONMENTAL SCIENCES↗

Beam Propagation Through Atmospheric Turbulence Using an Altitude-Dependent Structure Profile with Non-Uniformly Distributed Phase Screens

Modeling the effects of atmospheric turbulence on optical beam propagation is a key element in the design and analysis of free-space optical communication systems. Numerical wave optics simulations provide a particularly useful technique for understanding the degradation of the optical field in the receiver plane when the analytical theory is insufficient for characterizing the atmospheric channel. Motivated by such an application, we use a split-step method modeling the turbulence along the propagation path as a series of thin random phase screens with modified von Karman refractive index statistics using the Hufnagel-Valley turbulence profile to determine the effective structure constant for each screen. In this work, we employ a space-to-ground case study to examine the irradiance and phase statistics for both uniformly and non-uniformly spaced screens along the propagation path and compare to analytical results. We find that better agreement with the analytical theory is obtained using a non-uniform spacing with the effective structure constant for each screen chosen to minimize its contribution to the scintillation in the receiver plane. We evaluate this method as a flexible alternative to other standard layered models used in astronomical imaging applications.

Chahine, Yousef K.↗

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of negative triangularity on tearing mode stability in tokamak plasmas

Abstract The influence of negative triangularity (NT) of the plasma shape on the n = 1 (n is the toroidal mode number) tearing mode (TM) stability has been numerically investigated, with results compared to that of the positive triangularity (PT) counterpart. By matching the safety factor profile for a series of toroidal equilibria, several important plasma parameters, including the triangularity, the plasma equilibrium pressure, the plasma resistivity as well as the toroidal rotation, have been varied. The TM localized near the plasma edge is found to be more unstable in the NT plasmas as compared to the PT counterpart. The fundamental reason for this difference is the lack of favorable average curvature stabilization in NT configurations. Direct comparison of the Mercier index corroborates this conclusion. For the core-localized mode, where the difference in the local triangularity between NT and PT becomes small and the curvature stabilization is significantly reduced, larger Shafranov shift in the plasma core associated with the NT configuration results in more stable TM. The plasma toroidal flow generally stabilizes the TM in plasmas with both NTs and PTs. The flow stabilization is however weaker in the case of negative triangularity with finite plasma pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Beam-Based Diagnostics of Electric Guide Fields and Lattice Parameters for Run-1 of the Muon g-2 Storage Ring at Fermilab

A portion of the Muon g-2 Storage Ring electric system, which provides vertical beam focusing, exhibited an unexpected time dependence that produced a characteristic evolution of the stored beam during Run-1 of the Muon g-2 Experiment at Fermilab (E989). A method to reconstruct the Run-1 electric guide fields has been developed, which is based on a numerical model of the muon storage ring and optimization algorithms supported by COSY INFINITY. This method takes beam profile measurements from the Muon g-2 straw tracking detectors as input, and it produces a full reconstruction of the time-dependent fields. The fields can then be used for the reproduction of detailed beam tracking simulations and the calculation of ring lattice parameters for acceptance studies and systematic error evaluations.

43 PARTICLE ACCELERATORS↗

Thermospheric winds and exospheric temperatures from incoherent scatter radar measurements in four seasons

The time-dependent equations of neutral air motion are solved subject to three constraints: two of them are the usual upper and lower boundary conditions and the third is the value of the wind-induced ion drift at any given height. Using incoherent radar data, this procedure leads to a fast, direct numerical integration of the two coupled differential equations describing the horizontal wind components and yields time dependent wind profiles and meridional exospheric neutral temperature gradients. The diurnal behavior of the neutral wind system and of the exospheric temperature is presented for two solstice and two equinox days. The data used were obtained by the St. Santin and the Millstone Hill incoherent scatter radars. The derived geographic distributions of the exospheric temperatures are compared with those predicted by the OGO-6 empirical thermospheric model.

Antoniadis, D. A.↗

A note of sound radiation from distributed sources

The power output from a normally vibrating strip radiator is expressed in alternative general forms, one of these being chosen to refine and correct some particular estimates given by Heckl for different numerical ratios of strip width to wave length. An exact and explicit calculation is effected for sinusoidal velocity profiles when the strip width equals an integer number of half wave lengths.

Levine, H.↗

A note on sound radiation from distributed sources

The power output from a normally vibrating strip radiator is expressed in alternative general forms, one of these being chosen to refine and correct some particular estimates given by Heckl for different numerical ratios of strip width to wave length. An exact and explicit calculation is effected for sinusoidal velocity profiles when the strip width equals an integer number of half wave lengths.

Levine, H.↗

Measurements and predictions of turbulence in sound excited jets

The modifications to the structure of turbulence in a round jet excited by acoustic tones upstream of the jet exit are measured and compared with predictions. The prediction scheme assumes that the acoustic waves excite instability waves or large turbulent structures in the jet. The presence of these waves at finite amplitude causes an increase in the levels of the random turbulence. Models are presented for the interaction between the instability waves and the turbulence. The numerical scheme permits the enhanced turbulence and large turbulent structures to modify the shape of the mean velocity profile. The changes in the turbulence spectra with excitation are measured. The jet has a Reynolds number of 3.75 x 10 to the 5th and the initial boundary layer is turbulent.

Morris, P. J.↗

Flux-vector splitting for the Euler equations

Flux-vector splitting is shown to be the simplest way of implementing upwind differencing for the full or isenthalpic Euler equations combined with the ideal-gas law. The scheme presented produces steady shock profiles with two interior zones. A disadvantage in using flux-vector splitting is that it leads to numerical diffusion of a contact discontinuity at rest. This diffusion, however, can be removed, and present research is aimed at achieving this with minimal computational effort.

Van Leer, B.↗

Studies on the interference of wings and propeller slipstreams

The small disturbance potential flow theory is applied to determine the lift of an airfoil in a nonuniform parallel stream. The given stream is replaced by an equivalent stream with a certain number of velocity discontinuities, and the influence of these discontinuities is obtained by the method of images. Next, this method is extended to the problem of an airfoil in a nonuniform stream of smooth velocity profile. This model allows perturbation velocity potential in a rotational undisturbed stream. A comparison of these results with numerical solutions of Euler equations indicates that, although approximate, the present method provides useful information about the interaction problem while avoiding the need to solve the Euler equations.

Prabhu, R. K.↗

Budget study of a mesoscale convective system - Model simulation

A tropical squall-type cloud cluster is examined as an example of a mesoscale convective complex. The precipitation growth processes and air circulations that develop in the convective and stratiform regions are studied using a data set generated from a time-dependent numerical cloud model. The relationship of the stratiform region of the mesoscale convective complex is discussed. The vertical profiles of heating and drying are calculated. Comparisons are made between simulations with and without ice-phase microphysical processes and a simulation with forcing by weaker lifting at middle and upper levels.

Tao, Wei-Kuo↗