Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gaussian simulation”

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 793 records · Page 44

Discovering the Unknowns: A First Step

This article aims at discovering the unknown variables in the system through data analysis. The main idea is to use the time of data collection as a surrogate variable and try to identify the unknown variables by modeling gradual and sudden changes in the data. We use Gaussian process modeling and a sparse representation of the sudden changes to efficiently estimate the large number of parameters in the proposed statistical model. The method is tested on a realistic dataset generated using a one-dimensional implementation of a Magnetized Liner Inertial Fusion (MagLIF) simulation model, and encouraging results are obtained.

42 ENGINEERING↗

Utilizing HYSPLIT for Emergency Response Modeling at SRS

At SRS, emergency responders use a variety of tools to detect, track, and mitigate hazardous material releases into the atmosphere. Two models currently used at SRS are Puff-Plume and the Lagrangian Particle Dispersion Model (LPDM), a Gaussian and Lagrangian model, respectively. A decision has been made to replace LPDM with the more widely-supported Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model for evaluating inhalation and ingestion doses following a release. HYSPLIT is designed to compute complex dispersion and deposition simulation To achieve the implementation of HYSPLIT, we have developed a preliminary UI framework that will allow HYSPLIT to be run on ATG computers without the need for active network connections, thus avoiding the loss of capabilities in the event of a network outage during an emergency.

Earley, Ian↗

A new approach to active control of rotorcraft vibration

A state-variable feedback approach is utilized for active control of rotorcraft vibration. Fuselage accelerations are passed through undamped second-order filters with resonant frequencies at N/rev. The resulting outputs contain predominantly the N/rev vibration components, phase shifted by 180 deg, and are used to drive the blade pitch to cancel this component of fuselage vibration. The linear-quadratic-gaussian (LQG) method is used to design a feedback control system utilizing these filtered accelerations. The design is based on a nine-degree-of-freedom linear model of the Rotor System Research Aircraft (RSRA) in hover and is evaluated on a nonlinear blade-element simulation of the RSRA for this flight condition. The system is shown to essentially eliminate vibrations at N/rev in all axes. The required blade-pitch amplitude is within the capability of conventional actuators at the N/rev frequency.

Gupta, N. K.↗

Constraints on cosmological models from the IRAS QDOT survey results

The results of the recent IRAS survey are used to explore constraints on the nature of the Galaxy density field. Realistic simulations of the survey made from N-body simulations of the n = -1 power-law power spectrum and the Omega = 1 CDM models show that the power-law model is basically consistent with all QDOT variance results, while the CDM model has a smaller variance than the observation at the 20/h Mpc Gaussian smoothing scale. It is concluded that the QDOT results are best explained if the power spectrum of the Galaxy density field has a slope of about -1 near a wavelength of about 100/h Mpc, while the CDM model has a slope of about 0 at that scale. The possibility of marginal non-Gaussian behavior at the 20/h Mpc scale reported by Saunders et al. (1990) is examined for gravitational instability models with Gaussian initial conditions. It is argued that the observed (delta exp 3) versus (delta exp 2) in redshift space is not what these models would predict.

Park, Changbom↗

Measurement of Flaw Size From Thermographic Data

Simple methods for reducing the pulsed thermographic responses of delaminations tend to overestimate the size of the delamination, since the heat diffuses in the plane parallel to the surface. The result is a temperature profile over the delamination which is larger than the delamination size. A variational approach is presented for reducing the thermographic data to produce an estimated size for a flaw that is much closer to the true size of the delamination. The method is based on an estimate for the thermal response that is a convolution of a Gaussian kernel with the shape of the flaw. The size is determined from both the temporal and spatial thermal response of the exterior surface above the delamination and constraints on the length of the contour surrounding the delamination. Examples of the application of the technique to simulation and experimental data are presented to investigate the limitations of the technique.

Winfree, William P.↗

Performance of Variable Coded Modulations over a Nonlinear Channel for VCM Protocol Red Book

This presentation provides a summary of the results obtained from a study which evaluates the performance of variable coded modulations specified in the "Variable Coded Modulation Protocol" CCSDS Red Book, over a nonlinear channel that is currently on the CCSDS Blue Book track. The objectives are to perform simulation that identifies the operating signal-to-noise-power ratio required to achieve a codeword error rate of 1e-4 over the specified nonlinear channel for various modulations (i.e. BPSK, QPSK, 8-PSK) and LDPC coding schemes, and compare against performance over an ideal additive white Gaussian channel. We select a only subset of operating modes to perform the analysis which are based on a number of mitigation techniques to combat the nonlinear distortions. Specifically, transmitter centroidal pre-distortion, phase post-distortion and receiver mean phase estimation are analyzed in this study.

coding↗

Urban Flood Modeling: Uncertainty Quantification and Physics‐Informed Gaussian Processes Regression Forecasting

Abstract Estimating uncertainty in flood model predictions is important for many applications, including risk assessment and flood forecasting. We focus on uncertainty in physics‐based urban flooding models. We consider the effects of the model's complexity and uncertainty in key input parameters. The effect of rainfall intensity on the uncertainty in water depth predictions is also studied. As a test study, we choose the Interconnected Channel and Pond Routing (ICPR) model of a part of the city of Minneapolis. The uncertainty in the ICPR model's predictions of the floodwater depth is quantified in terms of the ensemble variance using the multilevel Monte Carlo (MC) simulation method. Our results show that uncertainties in the studied domain are highly localized. Model simplifications, such as disregarding the groundwater flow, lead to overly confident predictions, that is, predictions that are both less accurate and uncertain than those of the more complex model. We find that for the same number of uncertain parameters, increasing the model resolution reduces uncertainty in the model predictions (and increases the MC method's computational cost). We employ the multilevel MC method to reduce the cost of estimating uncertainty in a high‐resolution ICPR model. Finally, we use the ensemble estimates of the mean and covariance of the flood depth for real‐time flood depth forecasting using the physics‐informed Gaussian process regression method. We show that even with few measurements, the proposed framework results in a more accurate forecast than that provided by the mean prediction of the ICPR model.

Kohanpur, Amir H.↗

Connecting galaxy evolution in clusters with their radial profiles and phase space distribution: results from the IllustrisTNG hydrodynamical simulations

ABSTRACT We study the population of galaxies around galaxy clusters in the hydrodynamic simulation suite IllustrisTNG 300-1 to study the signatures of their evolutionary history on observable properties. We measure the radial number density profile, phase space distribution, and splashback radius for galaxies of different masses and colours over the redshift range z = 0−1. The three primary physical effects which shape the galaxy distribution within clusters are the galaxy quenching, angular momentum distribution, and dynamical friction. We find three distinct populations of galaxies by applying a Gaussian mixture model to their distribution in colour and mass. They have distinct evolutionary histories and leave distinct signatures on their distribution around cluster haloes. We find that low-mass red galaxies show the most concentrated distribution in clusters and the largest splashback radius, while high-mass red galaxies show a less concentrated distribution and a smaller splashback radius. Blue galaxies, which mostly quench into the low-mass red population, have the shallowest distribution within the clusters, with those on radial orbits quenched rapidly before reaching pericentre. Comparison with the distribution of galaxies from the Dark Energy Survey survey around Sunyaev–Zeldovich clusters from the Atacama Cosmology Telescope and South Pole Telescope surveys shows evidence for differences in galaxy evolution between simulations and data.

79 ASTRONOMY AND ASTROPHYSICS↗

Performance evaluation of land mobile satellite system under fading and interference using multiple TCM by Monte-Carlo simulation

The performance of the land mobile satellite system (LMSS) was evaluated by using two trellis coded modulation (TCM) schemes under multipath fading and interference. The results were also compared with uncoded QPSK. The trellis coding formats included a trellis code designed for optimum performance on the fading channel that typifies satellite to mobile communications, as well as one designed for optimum performance in the additive white Gaussian noise (AWGN) channel. The results show that, in Rayleigh fading, the TCM code designed for optimum performance in such an environment performs 4 dB better than the TCM code designed for optimum performance in AWGN, and 10 dB better than the uncoded format, at a bit error rate (BER) of 10 exp -4. Additional results on the performance degradation due to the nonlinearities in the satellite transponder and the adjacent and cochannel interference show that the TCM system is more sensitive than the uncoded system to the phase distortion caused by these impairments.

Kwa, S. C.↗

Nonlinear quantum behavior of ultrashort-pulse optical parametric oscillators

The quantum features of ultrashort-pulse optical parametric oscillators (OPOs) are investigated theoretically in the nonlinear regime near and above threshold. Viewing the pulsed OPO as a multimode open quantum system, we rigorously derive a general input-output model that features nonlinear coupling among many cavity (i.e., system) signal modes and a broadband single-pass (i.e., reservoir) pump field. Under appropriate assumptions, our model produces a Lindblad master equation with multimode nonlinear Lindblad operators describing two-photon dissipation and a multimode four-wave-mixing Hamiltonian describing a broadband, dispersive optical cascade, which we show is required to preserve causality. To simplify the multimode complexity of the model, we employ a supermode decomposition to perform numerical simulations in the regime where the pulsed supermodes experience strong single-photon nonlinearity. We find that the quantum nonlinear dynamics induces pump depletion as well as corrections to the below-threshold squeezing spectrum predicted by linearized models. We also observe the formation of non-Gaussian states with Wigner-function negativity and show that the multimode interactions with the pump, both dissipative and dispersive, can act as effective decoherence channels. Finally, we briefly discuss some experimental considerations for potentially observing such quantum nonlinear phenomena with ultrashort-pulse OPOs on nonlinear nanophotonic platforms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tree encoding of Gaussian sources

Tree codes are known to be capable of performing arbitrarily close to the rate-distortion function for any memoryless source and single-letter fidelity criterion. Tree coding and tree search strategies are investigated for the discrete-time memoryless Gaussian source encoded for a signal-power-to-mean-squared-error ratio of about 30 dB (about 5 binary digits per source output). Also, a theoretical lower bound on average search effort is derived. Two code search strategies (the Viterbi algorithm and the stack algorithm) were simulated in assembly language on a large digital computer. After suitable modifications, both strategies yielded encoding with a signal-to-distortion ratio about 1 dB below the limit set by the rate-distortion function. Although this performance is better than that of any previously known instrumentable scheme, it unfortunately requires search computation of the order of 100,000 machine cycles per source output encoded.

Dick, R. J.↗

An evaluation of thematic mapper simulator data for mapping forest cover

Computer-aided analysis techniques applied to Thematic Mapper Simulator (TMS) data were evaluated for the purpose of mapping forest cover types. Classification results obtained using a supervised set of training statistics and various combinations of three and four channel subsets of the seven available TMS channels are compared for the L2 (Minimum Euclidean Distance), GML (Gaussian Maximum Likelihood), and SECHO (Supervised Extraction and Classification of Homogeneous Objects) classification algorithms. SECHO performed significantly better than either of the two per-point classifiers for the untransformed data. Overall classification results of the Karhunen-Loeve transformation increased for the L2 algorithm, but decreased for both the GML and SECHO algorithms.

Dean, M. E.↗

Ridge Regression Signal Processing

The introduction of the Global Positioning System (GPS) into the National Airspace System (NAS) necessitates the development of Receiver Autonomous Integrity Monitoring (RAIM) techniques. In order to guarantee a certain level of integrity, a thorough understanding of modern estimation techniques applied to navigational problems is required. The extended Kalman filter (EKF) is derived and analyzed under poor geometry conditions. It was found that the performance of the EKF is difficult to predict, since the EKF is designed for a Gaussian environment. A novel approach is implemented which incorporates ridge regression to explain the behavior of an EKF in the presence of dynamics under poor geometry conditions. The basic principles of ridge regression theory are presented, followed by the derivation of a linearized recursive ridge estimator. Computer simulations are performed to confirm the underlying theory and to provide a comparative analysis of the EKF and the recursive ridge estimator.

Kuhl, Mark R.↗

Conditional H I Mass Functions and the H I -to-halo Mass Relation in the Local Universe

We present a new H I mass estimator that relates log 10 (M HI /M * ) to a linear combination of four galaxy properties: stellar surface mass density, color index u - r, stellar mass, and concentration index, with the scatter of individual galaxies around the mean H I mass modeled with a Gaussian distribution function. We calibrate the estimator using the xGASS sample, including both H I detection and nondetection, and constrain the model parameters through Bayesian inferences. Tests with mock catalogs demonstrate that our estimator provides unbiased H I masses for optical samples like SDSS. We apply our estimator to the SDSS spectroscopic sample to estimate the H I mass function (HIMF) of local galaxies, as well as the conditional H I mass function in galaxy groups and the H I –halo mass relation. Our HIMF agrees with the ALFALFA measurements at M HI ≳ 5 × 10 9 M ⊙ , but with higher amplitude and a steeper slope at lower masses. We show that this discrepancy is caused primarily by the cosmic variance, which is corrected for the SDSS sample but not for ALFALFA. The total CHIMFs for all halo masses can be described by a single Schechter function, while those of central galaxies show a double-Gaussian profile. The total H I mass in a group increases monotonically with halo mass, but for central galaxies, the H I mass shows weak dependence on halo mass when M h ≳ 10 12 M ⊙ . The observed H I –halo mass relation is not reproduced by current hydrodynamic simulations and semianalytic models of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗

Numerical simulations of collisions between rotating particles

Numerical simulations of Keplerian systems consisting of 200 mutually colliding rotating particles show that friction and surface irregularity reduce the equilibrium velocity dispersion and transfer some of the energy of random velocities to rotational velocities. The results confirm the theoretical predictions of Salo (1987). Simulations are performed for identical and different particle sizes, and for a power-law distribution of sizes. Taking rotation into account is not found to alter the previously observed Rayleigh distribution of eccentricities and inclinations, while the components of the rotational velocities are found to follow a Gaussian distribution. Application to the rarefied regions of Saturn's rings suggests that friction is able to reduce the equilibrium geometric thickness by about one-half.

Salo, H.↗

Vortex Beam Sidelobe Suppression via Bayliss Synthesis from a Rectangular Phased Array

We demonstrate the application of the Bayliss synthesis method for suppression of sidelobes in an OAM vortex beam, when generated by a commercial-of-the-shelf, rectangular phased array. The method is optimized for the current application and simulated for a 26.5 GHz signal in CST Microwave Studio. Sidelobes are notably reduced by 16.9 dB, with minimal reduction in gain and increase in beamwidth, as compared to a uniformed stimulated phased array. Further, the OAM modes are found to have greater purity that that of the uniformly stimulated array, as well as for other amplitude tapering methods. The Bayliss synthesis is then applied to a phased array in a test range, verifying the simulation results that demonstrate sidelobe suppression. These results indicate that OAM may be applied to existing phased arrays with minimal effort.

Orbital Angular Momentum (OAM)↗

A continuous mixing model for pdf simulations and its applications to combusting shear flows

The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in this work. A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models.

Hsu, A. T.↗