Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “non-linear models”

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.

674 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

On the Synergy Between Numerics and Subgrid Scale Modeling in LES of Stratified Flows: Grid Convergence of a Stratocumulus-Topped Boundary Layer

The effectiveness of a linear upwinding scalar advection scheme to suppress numerical dispersion errors near sharp inversions in large-eddy simulations of a nocturnal stratocumulus-topped boundary layer is assessed. Linear upwinding is a trade-off between non-dissipative and non-linear positive definite advection schemes. It is shown that linear upwinding does not negatively impact the model's grid convergence properties and a sharp inversion free of numerical artifacts is maintained. Even though mean profiles and turbulence fluxes show good grid convergence characteristics the liquid water amount varies significantly with grid resolution. The entrainment rate is identical for all resolutions and independent of the liquid water amount. For the present stratocumulus case, the impact of cloud-top radiative cooling is negligible and turbulence is largely driven by convection emanating from the surface.

Matheou, Georgios↗
NASA NTRS2019

A Communication Channel Density Estimating Generative Adversarial Network

Autoencoder-based communication systems use neural network channel models to backwardly propagate message reconstruction error gradients across an approximation of the physical communication channel. In this work, we develop and test a new generative adversarial network (GAN) architecture for the purpose of training a stochastic channel approximating neural network. In previous research, investigators have focused on additive white Gaussian noise (AWGN) channels and/or simplified Rayleigh fading channels, both of which are linear and have well defined analytic solutions. Given that training a neural network is computationally expensive, channel approximation networks— and more generally the autoencoder systems—should be evaluated in communication environments that are traditionally difficult. To that end, our investigation focuses on channels that contain a combination of non-linear amplifier distortion, pulse shape filtering, intersymbol interference, frequency-dependent group delay, multipath, and non-Gaussian statistics. Each of our models are trained without any prior knowledge of the channel. We show that the trained models have learned to generalize over an arbitrary amplifier drive level and constellation alphabet. We demonstrate the versatility of our GAN architecture by comparing the marginal probability density function of several channel simulations with that of their corresponding neural network approximations

Smith, Aaron↗
NASA NTRS2019

Using solar eclipse events to validate VIIRS reflective solar band calibration at multiple radiance levels

The VIIRS instruments aboard the SNPP and NOAA-20 (N20) satellites have 14 reflective solar bands (RSB) covering a spectral range from 412 nm to 2250 nm. The radiance of each VIIRS RSB is calculated from a quadratic function of the background-subtracted digital count, 𝑑𝑛, with the quadratic coefficients determined during pre-launch testing. On orbit, calibration is maintained using observations of a sunlit solar diffuser (SD), which views the Sun every orbit through a fixed attenuation screen. The SD observations, which are at nearly the same radiance level every orbit, provide a time varying overall calibration adjustment factor, the F-factor. But there is no designed on-orbit mechanism for calibration at multiple radiance levels, so the relative strengths of the quadratic coefficients continue to be fixed at the pre-launch values. On a few rare occasions, the VIIRS instruments have passed through a partial solar eclipse during the part of the satellite orbit when the SD is illuminated by the Sun (near the South Pole). As of August 2019, the SNPP and N20 VIIRS SDs have observed five and three partial solar eclipses, respectively. While these events are rare, they offer a unique opportunity to test the RSB calibration using the SD at different radiance levels. In this paper, we compare the reduction in the measured SD signals during an eclipse to the predicted radiance reduction based on the Sun-Moon geometry and a solar radiance model. We find good agreement between the data and model for all events, indicating that the VIIRS RSB gain linearity has remained fairly stable on-orbit. The most significant deviation is for the N20 short-wave infrared bands, which had non-linearity concerns during pre-launch testing. We also investigate the SNPP results using different versions of the prelaunch quadratic gain coefficients.

solar eclipse↗
NASA NTRS2020

Evaluating impacts of snow, surface water, soil and vegetation on empirical vegetation and snow indices for the Utqiaġvik tundra ecosystem in Alaska with the LVS3 model

Satellite observations for the Arctic and boreal region may contain information of vegetation, soil, snow, snowmelt, and/or other surface water bodies. We investigated the impacts of vegetation, soil, snow and surface water on empirical vegetation/snow indices on a tundra ecosystem area located around Utqiaġvik (formerly Barrow) of Alaska with the Moderate Resolution Imaging Spectrometer (MODIS) images in 2001–2014. Empirical vegetation indices, such as normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), the index of near infrared of vegetation (NIRv), and modified EVI (EVI2), have been used to monitor vegetation. Normalized difference snow index (NDSI) has been widely applied to monitor snow. The vegetation cover fraction (VGCF), the soil cover fraction (SOILCF), the snow cover fraction (SNOWCF), the surface water body cover fraction (WaterBodyCF), the fractional absorption of photosynthetically active radiation (PAR) by vegetation chlorophyll (fAPARchl), the fractional absorption of PAR by non-chlorophyll components of the vegetation (fAPARnon-chl), and the fractional absorption of PAR by the entire canopy (fAPARcanopy) are retrieved with the MODIS images and a coupled Leaf-Vegetation-Soil-Snow-Surface water body radiative transfer model, LVS3. The vegetation indices (NDVI, EVI, EVI2 and NIRv) differ from VGCF, fAPARchl, fAPARnon-chl, and fAPARcanopy. In addition to vegetation, we find that soil, snow and surface water also have impacts on vegetation indices NDVI, EVI (EVI2), and NIRv. Presence of snow makes lower the observed values of NDVI, EVI2 and NIRv. After snowmelt is gone, the vegetation indices (NDVI, EVI, EVI2 and NIRv) linearly decrease with SOILCF and WaterBodyCF, and WaterBodyCF has stronger impacts on these vegetation indices than SOILCF. The relationship between EVI and snow is complicated. NDSI non-linearly increases with SNOWCF, but linearly increases with sum of SNOWCF and WaterBodyCF (sum = 0.5893 × NDSI +0.4342, R2 = 0.976). NDSI linearly decreases with VGCF, and the relationship between NDSI and SOILCF is complex. Retrievals of VGCF, fAPARchl, fAPARnon-chl and fAPARcanopy with the LVS3 model provide alternatives for vegetation monitoring and ecological modeling.

Qingyuan Zhang↗
NASA NTRSDate not supplied

Progress on Inverse Estimation Technique of Non-Linear Pitch Damping Coefficient Curves Using Free-Flight CFD Generated Trajectories

Characterization of entry vehicle pitch damping coefficient curves is crucial to ensure appropriate re-entry and overall mission success. The pitch damping coefficient (C_(m_q )+C_(m_α ̇ )) is used to encapsulate the oscillatory growth or decay of a body during a trajectory. The inverse estimation technique utilizes an existing Free-Flight CFD (FF-CFD) dataset and wraps a reconstruction algorithm in an optimizer. The reconstruction integrates the planar equations of motion derived by Schoenenberger, Queen [1] using Python’s scipy.integrate.solve_ivp. The optimizer’s objective function is the normalized 𝐿2 residual of the angle of attack peaks between the reconstructed trajectory and the original data produced with FF-CFD. Inclusion of the peak times in this residual calculation allows for simultaneous optimization of the pitch moment coefficient, C_(m_α ). This residual equation is shown below in Eq. 1. The optimizer scipy.optimize.minimize was used with the gradient-based Powell method for the analysis presented, however the differential evolution method was investigated as means of comparison, and was found to produce marginally lower residual values with prohibitively longer run times. Further, the pitch damping curve is found by fitting a cubic interpolation function to a set of (α, (C_(m_q )+C_(m_α ̇ ))) control points, where the α points are held constant and the (C_(m_q )+C_(m_α ̇ )) values are the optimized parameters. The pitch moment curve uses a linear interpolation between the minimum and maximum α in the dataset. FF-CFD generated trajectories using the Dragonfly capsule geometry with the Genesis ballistic range model parameters were simulated and used for this analysis. These FF-CFD trajectories simulate planar motion, as restricted by the reconstructing the equations of motion, of three different cases: 1-DoF (free-to-pitch), 2-DoF (free-to-pitch and heave), and 3-DoF (free-to-pitch, heave, and decelerate). Pitch damping coefficient curves generated using this inverse estimation curve technique with FF-CFD 1-DoF Dragonfly data are found in Fig. 1. Preliminary results reconstructing ballistic range shots using these FF-CFD derived predictions of the pitch damping curve (Fig. 1) are shown in Fig. 2. It should be noted that the ballistic range shot used a Genesis model whereas the FF-CFD data used a Dragonfly geometry, however these geometries are similar.

entry↗
NASA NTRS2023

Global Climate Impacts of Greenland and Antarctic Meltwater: A Comparative Study

Both the Greenland and Antarctic ice sheets have been melting at an accelerating rate over recent decades. Meltwater from Greenland might be expected to initiate a climate response which is distinct, and perhaps different from, that associated with Antarctic meltwater. Which one might elicit a greater climate response, and what mechanisms are involved? To explore these questions, we apply “Climate Response Functions (CRFs)” to guide a series of meltwater perturbation experiments using a fully-coupled climate model. In both hemispheres, meltwater drives atmospheric cooling, sea-ice expansion, and strengthened Hadley and Ferrel cells. Greenland meltwater induces a slowdown of the Atlantic Meridional Overturning Circulation (AMOC) and a cooling of the subsurface ocean in the northern high-latitudes. Antarctic meltwater, instead, induces a slowdown of the Antarctic Bottom Water formation and a warming of the subsurface ocean around Antarctica. For melt-rates up to 2000 Gt yr −1 , the climate response is rather linear. However, as melt-rates increase to 5000 Gt yr −1 , the climate response becomes non-linear. Due to a collapsed AMOC, the climate response is super-linear at high Greenland melt-rates. Instead, the climate response is sub-linear at high Antarctic melt-rates, due to the halting of the northward expansion of Antarctic sea ice by warm surface waters. Finally, in the linear limit, we use CRFs and linear convolution theory to make projections of important climate parameters in response to meltwater scenarios, which suggest that Antarctic meltwater will become a major driver of climate change, dominating that of Greenland meltwater, as the current century proceeds.

Greenland ice sheet↗
NASA NTRS2024

MSFC-TVC | TB-03: Derivation of Thrust Vector Control (TVC) Actuator-Force / Gimbal-Torque Transformation Matrix

Thrust vector control (TVC) systems for rocket engine propulsion traditionally use a simple linear relationship to convert between actuator forces and torques about the engine gimbal’s center-of-rotation (COR). As shown in Equation (1), the torque about the gimbal COR is proportional to the applied actuator force and the TVC moment arm (MA)—the perpendicular distance between the TVC actuator’s line-of-action (LOA) and the engine gimbal’s COR. While this fundamental relationship remains valid and accurate in a two-dimensional (2D), one-degree-of-freedom (1-DOF) context—particularly in its non-linear formulation as described in the ER63 Technical Bulletin TB-02 (Derivation of Thrust Vector Control (TVC) Engine-Gimbal / Actuator Moment-Arm Geometry)—it becomes limited when extended to three-dimensional (3D), two-degree-of-freedom (2-DOF) analyses. In 3D, out-of-plane angular displacements can arise, causing the gimbaled engine plane-of-motion to become non-coplanar with the respective actuator plane-of-motion. Such misalignments occur due to the rod-end (RE) and/or tail-stock (TE) mounting geometry of individual TVC actuators. These geometric complexities lead to inaccuracies in calculating engine torques and corresponding actuator forces if only the traditional TVC MA relationship is employed. For higher gimbal angular displacements (e.g., greater than 5 degrees), these inaccuracies become more pronounced, necessitating a more robust mathematical framework. This document presents a comprehensive 3D (2-DOF) derivation of a transformation matrix that accurately converts between gimbaled engine torques and TVC actuator forces. By incorporating the necessary geometric and rotational considerations, this new approach corrects the limitations of the traditional TVC MA method. Subsequent sections compare the newly formulated approach to the traditional method, demonstrating its enhanced accuracy and reliability for a broad range of gimbaled engine conditions.

Thrust Vector Control↗
NASA NTRSDate not supplied

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics↗