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 883 records · Page 49

Coupled Target-Beam-Moderator Optimization for the Second Target Station

This report describes the results for a coupled target-beam-moderator optimization analysis for the Second Target Station (STS) at ORNL's Spallation Neutron Source. This study is a continuation of the optimization analysis for the moderators in the preliminary design of STS performed in 2022. In the 2022 analysis the dimensions of the moderators are parameterized, while the target and the proton beam profile are kept constant. In this analysis the target height and the proton beam profile are added as parameters. This allows to study the coupled effects of changing target, moderator and beam dimensions. Similar to the 2022 analysis, this work is performed with an automated optimization workflow that uses the optimization toolbox DAKOTA, parameterized geometries in CREO and SpaceClaim, the unstructured mesh generation in Attila4MC, and the particle transport code MCNP6.2©. This workflow enables an efficient optimization using high-fidelity geometries. The main conclusions of this analysis are the following: • Coupled beam-target-moderator optimization provides a few additional percent performance gain over stand-alone moderator optimization. • The moderator performance is not very sensitive to the target height (between ≈60 and ≈80 mm) as long as the beam profile is chosen adequately. • The moderator performance is sensitive to the choice of beam spatial standard deviations, even when the footprint is kept constant. • The optimal moderator radius is the same for a beam footprint of 30 cm 2 , 62.5 cm 2 , and 90 cm 2 . Also the slope of the super-gaussian beam profile does not significantly impact the optimal moderator radius. • The optimal parameters and sensitivities are very similar to the 2022 optimization analysis. These results only indicate a a difference in the optimal radius of the cylindrical moderator, however, this has been corrected in the final design moderator optimization. The main purpose of this report is to document the simulations, results and lessons learned. The most impactful results are summarized in. We also note that the target geometry used in this work is not the final design.

43 PARTICLE ACCELERATORS↗

Multi-Sensor Optimal Motion Planning for Radiological Contamination Surveys by Using Prediction-Difference Maps

Distributed and networked mobile sensor platforms using unmanned aerial and/or ground vehicles to survey areas of interest offer a safer and more efficient method for radiological contamination mapping; however, most applications rely on uniformly sweeping of the area in a raster-type motion without utilizing the information available in a dynamic sense. We have developed a fully autonomous optimal motion planning procedure for networks with two or more mobile sensors. The procedure utilizes well-established concepts of Gaussian processes in combination with control laws based on centroidal Voronoi tessellations to achieve optimal next-iteration sensor movements. A new method of informing optimal motion planning is proposed, whereby the absolute difference between the prior and current full-map prediction, referred to as the prediction-difference map, is used as the spatial density function within each Voronoi cell, providing immediate and iterative feedback for dynamic use of available information. The Gaussian process regression model used to estimate the contamination in unvisited locations also provides prediction uncertainties, and can be used as a quantitative metric to assess the confidence in the calculated contamination map; these estimates and prediction uncertainties are unavailable for standard uniform survey routines as they can only produce maps in the vicinity of observed locations. We present through simulation the achievable performance gains from using this new method by directly comparing to a uniform survey method. Results show that using the prediction-difference maps to inform motion planning procedures offers a faster rate of producing an accurate and convergent map relative to a uniform survey route.

47 OTHER INSTRUMENTATION↗

Decision-directed detector for overlapping PCM/NRZ signals.

A decision-directed (DD) technique for the detection of overlapping PCM/NRZ signals in the presence of white Gaussian noise is investigated. The performance of the DD detector is represented by probability of error Pe versus input signal-to-noise ratio (SNR). To examine how much improvement in performance can be achieved with this technique, Pe's with and without DD feedback are evaluated in parallel. Further, analytical results are compared with those found by Monte Carlo simulations. The results are in good agreement.

Wang, C. D.↗

Users manual for flight control design programs

Computer programs for the design of analog and digital flight control systems are documented. The program DIGADAPT uses linear-quadratic-gaussian synthesis algorithms in the design of command response controllers and state estimators, and it applies covariance propagation analysis to the selection of sampling intervals for digital systems. Program SCHED executes correlation and regression analyses for the development of gain and trim schedules to be used in open-loop explicit-adaptive control laws. A linear-time-varying simulation of aircraft motions is provided by the program TVHIS, which includes guidance and control logic, as well as models for control actuator dynamics. The programs are coded in FORTRAN and are compiled and executed on both IBM and CDC computers.

Nalbandian, J. Y.↗

Viscoelastic relaxation and topological fluctuations in glass-forming liquids

In this study, a method for characterizing the topological fluctuations in liquids is proposed. This approach exploits the concept of the weighted gyration tensor of a collection of particles and permits the definition of a local configurational unit (LCU). The first principal axis of the gyration tensor serves as the director of the LCU, which can be tracked and analyzed by molecular dynamics simulations. Analysis of moderately supercooled Kob–Andersen mixtures suggests that orientational relaxation of the LCU closely follows viscoelastic relaxation and exhibits a two-stage behavior. The slow relaxing component of the LCU corresponds to the structural, Maxwellian mechanical relaxation. Additionally, it is found that the mean curvature of the LCUs is approximately zero at the Maxwell relaxation time with the Gaussian curvature being negative. This observation implies that structural relaxation occurs when the configurationally stable and destabilized regions interpenetrate each other in a bicontinuous manner. Finally, the mean and Gaussian curvatures of the LCUs can serve as reduced variables for the shear stress correlation, providing a compelling proof of the close connection between viscoelastic relaxation and topological fluctuations in glass-forming liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A novel scheme to aid coherent detection of GMSK signals in fast Rayleigh fading channels

A novel scheme to insert carrier pilot to Gaussian Minimum Shift Keying (GMSK) signal using Binary Block Code (BBC) and a highpass filter in baseband is proposed. This allows the signal to be coherently demodulated even in a fast Rayleigh fading environment. As an illustrative example, the scheme is applied to a 16 kb/s GMSK signal, and its performance over a fast Rayleigh fading channel is investigated using computer simulation. This modem's 'irreducible error rate' is found to be Pe = 5.5 x 10(exp -5) which is more than that of differential detection. The modem's performance in Rician fading channel is currently under investigation.

Leung, Patrick S. K.↗

Prognostics of Power MOSFET

This paper demonstrates how to apply prognostics to power MOSFETs (metal oxide field effect transistor). The methodology uses thermal cycling to age devices and Gaussian process regression to perform prognostics. The approach is validated with experiments on 100V power MOSFETs. The failure mechanism for the stress conditions is determined to be die-attachment degradation. Change in ON-state resistance is used as a precursor of failure due to its dependence on junction temperature. The experimental data is augmented with a finite element analysis simulation that is based on a two-transistor model. The simulation assists in the interpretation of the degradation phenomena and SOA (safe operation area) change.

Celaya, Jose Ramon↗

First Constraints on the Epoch of Reionization Using the Non-Gaussianity of the Kinematic Sunyaev-Zel’dovich Effect from the South Pole Telescope and Herschel-SPIRE Observations

We report results from an analysis aimed at detecting the trispectrum of the kinematic Sunyaev-Zel’dovich (kSZ) effect by combining data from the South Pole Telescope (SPT) and Herschel-SPIRE experiments over a 100 deg 2 field. The SPT observations combine data from the previous and current surveys, namely SPTpol and SPT-3G, to achieve depths of 4.5, 3, and 16 μ K - arcmin in bands centered at 95, 150, and 220 GHz. For SPIRE, we include data from the 600 and 857 GHz bands. We reconstruct the velocity-induced large-scale correlation of the small-scale kSZ signal with a quadratic estimator that uses two cosmic microwave background (CMB) temperature maps, constructed by optimally combining data from all the frequency bands. We reject the null hypothesis of a zero trispectrum at 10.3 σ level. However, the measured trispectrum contains contributions from both the kSZ and other undesired components, such as CMB lensing and astrophysical foregrounds, with kSZ being sub-dominant. We use the agora simulations to estimate the expected signal from CMB lensing and astrophysical foregrounds. After accounting for the contributions from CMB lensing and foreground signals, we do not detect an excess kSZ-only trispectrum and use this nondetection to set constraints on reionization. By applying a prior based on observations of the Gunn-Peterson trough, we obtain an upper limit on the duration of reionization of Δ z re , 50 < 4.5 (95% confidence level). We find these constraints are fairly robust to foregrounds assumptions. This trispectrum measurement is independent of, but consistent with, Planck’s optical depth measurement. This result is the first constraint on the epoch of reionization using the non-Gaussian nature of the kSZ signal.

79 ASTRONOMY AND ASTROPHYSICS↗

A reactive electrochemomechanical theory for growth and remodeling of polyelectrolyte hydrogels and application to dynamic polymerization of $\text{DNA}$ hydrogels

Here this study develops a framework for growth and remodeling of active polyelectrolyte hydrogels that accounts for effects of compositional changes on the mechanical response. By developing a reactive electrochemomechanical theory, thermodynamical constraints upon reactive and remodeling processes are elucidated within a general framework that allows any number of chemical reactions to evolve the response of the gel and transfer mass and charge between constituents. Fully coupled, nonlinear constitutive relations are adopted for molar fluxes, allowing exploration of effects including cross-diffusion, electrophoresis, and electro-osmosis. A robust finite element implementation is developed in the open source FEBio software (febio.org) by exploiting an equivalence between electrochemomechanics and mixture theory. The implementation is verified against analytical solutions for free swelling, and a proper reduction to a prior chemomechanical theory is demonstrated for neutral gels swollen only by a solvent with no solutes. The theory and implementation are then applied to model the tunable large swelling achieved through dynamic polymerization of DNA crosslinkers seen in our recently developed experimental hydrogel system (Cangialosi et al., 2017). A novel constitutive model for reaction-driven evolution of the locking stretch λ L in a non-Gaussian mechanical free energy was developed, where the increasing concentration of DNA crosslinkers makes further swelling energetically favorable. With a single free parameter, excellent agreement was found between measured and predicted equilibrium swelling ratios. This study demonstrated the ability to extend the electrochemomechanical framework to include chemical reactions and composition-aware constitutive models, and showed that development of reactive models allows simulation of complex dynamic polymerization phenomena not treated before. The theoretical frame here can be further expanded in scope to incorporate additional non-ideal and nonlinear phenomena.

42 ENGINEERING↗

2D-FFTLog: efficient computation of real-space covariance matrices for galaxy clustering and weak lensing

ABSTRACT Accurate covariance matrices for two-point functions are critical for inferring cosmological parameters in likelihood analyses of large-scale structure surveys. Among various approaches to obtaining the covariance, analytic computation is much faster and less noisy than estimation from data or simulations. However, the transform of covariances from Fourier space to real space involves integrals with two Bessel integrals, which are numerically slow and easily affected by numerical uncertainties. Inaccurate covariances may lead to significant errors in the inference of the cosmological parameters. In this paper, we introduce a 2D-FFTLog algorithm for efficient, accurate, and numerically stable computation of non-Gaussian real-space covariances for both 3D and projected statistics. The 2D-FFTLog algorithm is easily extended to perform real-space bin-averaging. We apply the algorithm to the covariances for galaxy clustering and weak lensing for a Dark Energy Survey Year 3-like and a Rubin Observatory’s Legacy Survey of Space and Time Year 1-like survey, and demonstrate that for both surveys, our algorithm can produce numerically stable angular bin-averaged covariances with the flat sky approximation, which are sufficiently accurate for inferring cosmological parameters. The code CosmoCov for computing the real-space covariances with or without the flat-sky approximation is released along with this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Final Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing: Development of Microwave Cavity Flow Meter for Advanced Reactor High Temperature Fluids

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with one of the flat walls of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the microwave resonant frequency, which can be detected with a spectrum analyzer. We have developed a continuum electromechanics model of the microwave resonant cavity transducer performance. The deflection of the membrane is calculated using equations from the theory of plates and shells for deflection of radially constrained thin circular disk under uniform load. The microwave frequency shift is calculated using the equations from the microwave cavity perturbation theory. The transducer model was validated with proof-of-concept experimental data of water flow sensing with a Brass 360 right circular cylinder with 22.2mm diameter and 203µm thick wall. The cavity was excited through a subwavelength hole in TE 011 mode with resonant frequency f ≈ 17.8GHz. Subsequently, we have performed a preliminary proof-of-principle test of flow sensing in high temperature liquid sodium in environment. For this test, we have developed a cylindrical resonator with the same dimensions as for the Brass 360 cavity. The resonator for high temperature liquid sodium test was machined from stainless steel 316 and electroplated with silver on the interior surfaces. We have also developed and insertion probe consisting of a 50cm WR-42 brass waveguide enclosed in a protective SS316 tube. The cavity was excited through a subwavelength hole on the side of the wall of the cylinder in the TE 011 mode with resonant frequency f ≈ 17.8GHz. The liquid sodium setup consists of a cylindrical vessel with a center feed line, where a transducer inserted through the top cap of the vessel measures velocity of the impinging liquid jet. The sodium flow sensing study was performed in a liquid sodium vessel at 340°C temperature and ambient pressure. The flow rate was changed by varying sodium pump power, and a frequency shift of several hundred kHz was observed. The transducer was removed after 70 days of testing. No visual damage was observed, and the electromagnetic response remained the same. As a calibration experiment, we have assembled a water vessel with a center feed, and with dimensions similar to those of the liquid sodium setup. Water flow sensing was performed before and after liquid sodium test. The responses in both cases are similar and agree with COMSOL computer simulations. Because long-term (multi-year) experimental tests of transducer resilience to harsh environment are not practical, we have developed a probabilistic model of creep to estimate transducer resilience to the harsh environment. The probabilistic model considers diffusion creep under the condition of high temperature and low stress, where the stress and temperature are allowed to be random variables with Gaussian distributions. Using the probabilistic model, we estimate inelastic membrane deflections due to creep for several temperature ranges. We conclude that for temperatures less than 650°C, creep has negligible long-term effect on the transducer performance.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Creation and testing of an artificial neural network based carbonate detector for Mars rovers

We have developed an artificial neural network (ANN) based carbonate detector capable of running on current and future rover hardware. The detector can identify calcite in visible/NIR (350-2500 nm) spectra of both laboratory specimens covered by ferric dust and rocks in Mars analogue field environments. The ANN was trained using the Backpropagation algorithm with sigmoid activation neurons. For the training dataset, we chose nine carbonate and eight non-carbonate representative mineral spectra from the USGS spectral library. Using these spectra as seeds, we generated 10,000 variants with up to 2% Gaussian noise in each reflectance measurement. We cross-validated several ANN architectures, training on 9,900 spectra and testing on the remaining 100. The best performing ANN correctly detected, with perfect accuracy, the presence (or absence) of carbonate in spectral data taken on field samples from the Mojave desert and clean, pure marbles from CT. Sensitivity experiments with JSC Mars-1 simulant dust suggest the carbonate detector would perform well in aeolian Martian environments.

neural networks↗

Uncertainty Propagation in Pre-Flight Prediction of Unmanned Aerial Vehicle Separation Violation

Current forecasts on the future of aeronautics suggest an in- creasing number of unmanned aerial vehicles entering the low- altitude airspace in the next decades (FAA, 2018; Kopardekar et al., 2016). Small vehicles for package delivery as well as larger vehicles for urban air mobility will change the airspace drastically, increasing density of operations both in time, i.e. high number of take-off and landings per unit time, and in space, operating in dense urban environment. This scenario poses challenges to the current approach to air traffic control, and large efforts from academia, industry and regulatory bodies are dedicated to the development of new traffic management strategies that leverage higher computing and simulating capabilities available today. In this paper, we propose a simple look-ahead approach to predict potential minimum separation violations at the strategic level, that is before vehicles start flying, depending on the predefined 4D trajectories and uncertainty affecting the wind acting along those routes. The wind field is extracted from the NOAA North America Mesoscale Forecast System and interpolated using Gaussian process regression, while uncertainty affecting the expected cruise airspeed is propagated through error intervals. The approach allows the prediction of aircraft separation as a function of time, highlighting potential safety violations that would go undetected if uncertainty affecting the expected 4D trajectories is not considered. The paper will also discuss issues related to accuracy and scalability of the approach to multiple vehicle operations.

Trajectory Prediction↗

A Self-learning Neural Network Approach for Radio Frequency Interference Detection and Removal in Radio Astronomy

In this work, we present a novel neural network (NN) method for the detection and removal of radio frequency interference (RFI) from the raw digitized signal in the signal processing chain of a typical radio astronomy experiment. The main advantage of our method is that it does not require a training set. Instead, our method relies on the fact that the true signal of interest coming from astronomical sources is thermal and therefore described as a Gaussian random process, which cannot be compressed. We employ a variational encoder/decoder network to find the compressible information in the data stream that can explain the most variance with the fewest degrees of freedom. We demonstrate it on a set of toy problems and stored ring buffers from the Baryon Mapping eXperiment prototype. We find that the RFI subtraction is effective at cleaning simulated time streams: While we find that the power spectra of the RFI-cleaned time streams output by the NN suffer from extra signal consistent with additive noise, we find that it is generally around percent level across the band and sub 10% in contaminated spectral channels even when RFI power is an order of magnitude larger than the signal. We discuss the advantages and limitations of this method and possible implementation in the front end of future radio experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermal Overloading Risk Mitigation With a Semi-Analytical Probabilistic Model on Branch Current

A semi-analytical formulation is presented in this paper for the probability computation of branch current in multiphase systems. The developed formula is derived based on the linear power flow model in rectangular coordinates. The system uncertainty injections can be renewable energy resources or loads and are modeled using a Gaussian mixture model (GMM). The developed formula can be used to compute the line current violation probability as well as integrate into optimal power flow problem as chance-constraint relaxation. Here, the proposed formula is first compared with the Matlab embedded numerical integration function to show its performance. Besides, the semi-analytical formula is validated and compared with the Monte Carlo simulation method using the IEEE 123-bus system, EPRI Ckt5, and Ckt7 systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Logarithmic Bayesian Approach to Quantum Error Detection

We consider the problem of continuous quantum error correction from a Bayesian perspective, proposing a pair of digital filters using logarithmic probabilities that are able to achieve near-optimal performance on a three-qubit bit-flip code, while still being reasonable to implement on low-latency hardware. These practical filters are approximations of an optimal filter that we derive explicitly for finite time steps, in contrast with previous work that has relied on stochastic differential equations such as the Wonham filter. By utilizing logarithmic probabilities, we are able to eliminate the need for explicit normalization and can reduce the Gaussian noise distribution to a simple quadratic expression. The state transitions induced by the bit-flip errors are modeled using a Markov chain, which for log-probabilities must be evaluated using a LogSumExp function. We develop the two versions of our filter by constraining this LogSumExp to have either one or two inputs, which favors either simplicity or accuracy, respectively. Using simulated data, we demonstrate that the single-term and two-term filters are able to significantly outperform both a double threshold scheme and a linearized version of the Wonham filter in tests of error detection under a wide variety of error rates and time steps.

97 MATHEMATICS AND COMPUTING↗

VTOL controls for shipboard landing

The problem of landing a VTOL aircraft on a small ship in rough seas using an automatic controller is examined. The controller design uses the linear quadratic Gaussian results of modern control theory. Linear time invariant dynamic models are developed for the aircraft, ship, and wave motions. A hover controller commands the aircraft to track position and orientation of the ship deck using only low levels of control power. Commands for this task are generated by the solution of the steady state linear quadratic gaussian regulator problem. Analytical performance and control requirement tradeoffs are obtained. A landing controller commands the aircraft from stationary hover along a smooth, low control effort trajectory, to a touchdown on a predicted crest of ship motion. The design problem is formulated and solved as an approximate finite-time linear quadratic stochastic regulator. Performance and control results are found by Monte Carlo simulations.

Mcmuldroch, C. G.↗

Time-Dependent Satellite Explosion Probabilities for Long-Term Orbital Debris Environment Modeling

On-orbit accidental explosions are a significant contributor to the growth of the orbital debris population. Certain types of satellites have been seen to exhibit different behaviors and timelines of explosions. Recent assessments of on-orbit explosions for three categories of satellites – Russian Sistema Obespecheniya Zapuska (SOZ, Proton 4th stage attitude and ullage motors) units, spacecraft, and rocket bodies – suggest that on-orbit explosions can be modeled using a continuous, time-dependent probability of explosion. This paper discusses the methodology for developing a time-dependent probability of explosion as a function of the satellite’s orbital lifetime based on historical explosions. The SOZ units are found to be best modeled by a modified Gaussian probability distribution, with the peak probability of explosion occurring around 10.4 years on-orbit. Spacecraft and rocket bodies both exhibit exponentially decaying explosion probabilities, with rocket bodies showing fast-, medium-, and slow-time modes of decay. Total cumulative probabilities of explosion are approximately 57% for SOZ units, 4% for spacecraft, and 2% for rocket bodies. Examples of implementation in a simulation of the orbital debris environment using NASA's long-term evolutionary model are given, showing the overall explosion behavior over two centuries is significantly different from that using a previously implemented constant-probability type-dependent explosion rate model and provides a better representation of the historical explosion record.

Alyssa P Manis↗