Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density distribution”

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 1,351 records · Page 75

Bounds on intergalactic compact objects from observations of compact radio sources

The constraints that can be imposed on the cosmic density of uniformly distributed intergalactic compact objects in the mass range of 10 million to 1 billion solar masses, using VLBI observations of compact radio sources. A compact object at a sufficiently small angular separation from a compact radio source produces a deamplified and inverted secondary image. No such cases are known. Using published maps of compact radio sources and assuming an Einstein-de Sitter cosmology, a density parameter for mass in this range greater than 0.4 is excluded at the 99.7 percent confidence level. Observations using the VLBA should allow one to exclude or discover a density parameter greater than about 0.001 for compact objects in the mass range of 0.1-100 solar masses.

Kassiola, Aggeliki↗

High-Throughput Uniformity and Defect Monitoring in Low-Temperature Electrolysis Porous Transport Layers Using X-Ray Radiography

Effective quality control (QC) for manufacturing proton exchange membrane water electrolysis (PEMWE) components is critical to enabling widespread adoption of the technology for hydrogen generation. This study investigates X-ray radiography as a novel, high-throughput, potentially in-line QC technique for detecting defects and assessing material property distributions in titanium-based porous transport layers (PTLs) which constitute a crucial component of low temperature PEMWE stacks. We obtain radiographs of a set of fifteen PTLs and model their absorbance of the broadband radiation as a second-order polynomial to account for the non-monoenergetic radiation source used in this study. The resulting model serves as a basis for predicting the areal density and porosity distributions of the PTLs. We find radiography successful in detecting multiple instances of defects, including holes/depressions, cracks, and excess material on the surface or in the pores of the material, demonstrating its potential as a robust in-line QC tool for PTL manufacturing.

08 HYDROGEN↗

Study of auroral dynamics with combined spacecraft and incoherent-scatter radar data

We have examined Sondrestrom incoherent-scatter radar observations of ionospheric plasma density and temperature distributions, as well as measurements of F-region ion drifts that were made during a prenoon pass by the DMSP-F7 satellite through the radar field of view. The spacecraft traversed a region of intense electron precipitation with a characteristic energy below approximately 200 eV. Particles with such low characteristic energies are believed to originate, either directly or indirectly, in the magnetosheath. The precipitation region had a width of about 2 deg invariant latitude. The corotating radar observed a patch of enhanced electron density and elevated electron temperature in the F2 region between about 10.5 and 12 magnetic local time in the same invariant latitude range where DMSP-F7 detected the soft-electron flux. The ion drift pattern, also obtained by radar, shows that it is unlikely that the plasma patch was produced by solar radiation and advected into the radar field of view. We suggest that the radar observed modifications of the ionospheric plasma distribution, which resulted from direct entry of magnetosheath electrons into the magnetosphere and down to ionospheric altitudes. Model calculations of the ionospheric response to the observed electron flux support our interpretation.

Watermann, Juergen↗

Critical differences between electron beam melted and selective laser melted Ti-6Al-4 V

Effective optimization of the production of Ti-6Al-4 V using AM requires a fundamental understanding of the relative importance of different microstructural features to the deformation and failure mechanisms, particularly features that vary between production methods. In this study, the tensile response and deformation mechanisms of electron beam melted (EBM) AM Ti-6Al-4 V material loaded in different orientations and produced using various powder sizes were compared to those of selective laser melted (SLM) AM Ti-6Al-4 V material. The density and morphology of pores, phase fractions, prior-β grains, and defect microstructures were evaluated using scanning electron microscopy, X-ray computed tomography, electron backscatter diffraction, and transmission electron microscopy before and after deformation. The results were used to evaluate the relative importance of each feature on strengthening, deformation, and failure initiation mechanisms. Results focused primarily on coarse-powder EBM materials indicated that phase distribution and defect density were most influential for determining material yield strength as well as maximum possible strain to failure. Porosity was lower overall in EBM Ti-6Al-4 V than in SLM, allowing for occasional increases in part strain to failure, but remained a limiting factor determining overall part ductility.

36 MATERIALS SCIENCE↗

Direct local parametrization of nuclear state densities using the back-shifted Bethe formula

Level densities are often parametrized using the back-shifted Bethe formula (BBF) for nuclei that possess experimental data for s-wave neutron resonance average spacings and a complete discrete level sequence at low excitation energies. However, these parametrizations require the additional modeling of the dependence of the spin-cutoff parameter on excitation energy. Here, in this work, we avoid the need to model the spin distribution of level densities by using the experimental data to parametrize directly the state densities, for which the BBF does not depend on the spin-cutoff parameter. This approach allows for a local parameterization of state densities that is independent of the spin-cutoff parameter. We provide these parameters in a tabulated form for applications in nuclear reaction calculations and for testing microscopic approaches to state densities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The free compressible viscous vortex

The present study investigates the effects of compressibility on free (unsteady) viscous heat-conducting vortices. Analytical solutions are found in the limit of large but finite Reynolds number and small but finite Mach number. It is shown that the spreading of the vortex causes a radial flow. This flow is given by the solution of an ordinary differential equation, which gives the dependence of the radial velocity on the tangential velocity, density, and temperature profiles of the vortex. Estimates of the radial velocity found by solving this equation are found to be in good agreement with numerical solutions of the full equations. The equations for the viscous evolution are expanded in powers of Mach number to obtain detailed analytical solutions. It is shown that swirling axisymmetric compressible flows generate negative radial velocities far from the vortex core owing to viscous effects, regardless of the initial distributions of vorticity, density, and entropy.

Colonius, Tim↗

Modeling of the high column density systems in the Lyman-Alpha forest

The Lyman-α forests observed in the spectra of high-redshift quasars can be used as a tracer of the cosmological matter density to study baryon acoustic oscillations (BAO) and the Alcock-Paczynski effect. Extraction of cosmological information from these studies requires modeling of the forest correlations. While the models depend most importantly on the bias parameters of the intergalactic medium (IGM), they also depend on the numbers and characteristics of high-column-density systems (HCDs) ranging from Lyman-limit systems with column densities log N HI /1cm -2 > 17 to damped Lyman-α systems (DLAs) with log N HI /1cm -2 > 20.2. These HCDs introduce broad damped absorption characteristic of a Voigt profile. Consequently they imprint a component on the power spectrum whose modes in the radial direction are suppressed, leading to a scale-dependent bias. Using mock data sets of known HCD content, we test a model that describes this effect in terms of the distribution of column densities of HCDs, the Fourier transforms of their Voigt profiles and the bias of the halos containing the HCDs. Our results show that this physically well-motivated model describes the effects of HCDs with an accuracy comparable to that of the ad-hoc models used in published forest analyses. We also discuss the problems of applying the model to real data, where the HCD content and their bias is uncertain.

Lyman alpha forest↗

Ionospheric footprint of magnetosheathlike particle precipitation observed by an incoherent scatter radar

We have examined Sondrestrom incoherent scatter radar observations of ionospheric plasma density and temperature distributions and measurements of F region ion drifts that were made during a prenoon pass of the Defense Meteorological Satellite Program (DMSP)-F7 satellite through the radar field of view. The spacecraft traversed a region of intense electron precipitation with a characteristic energy below approximately 200 eV. Particles with such low characteristic energies are believed to be directly or indirectly of magnetosheath origin. The precipitation region had a width about 2 deg invariant latitude and covered the low-latitude boundary layer (LLBL), the cusp, and the equatorward section of the plasma mantle (PM). The corotating radar observed a patch of enhanced electron density and elevated electron temperature in the F2 region between about 10.5 and 12 magnetic local time in the same invariant latitude range where DMSP-F7 detected the soft-electron flux. The ion drift pattern, also obtained by radar, shows that it is unlikely that the plasma patch was produced by solar radiation and advected into the radar field of view. We suggest that the radar observed modifications of the ionospheric plasma distribution, which resulted from direct entry of magnetosheath electrons into the magnetosphere and down to ionospheric altitudes. Model calculations of the ionospheric response to the observed electron precipitation support our interpretation. The spectral characteristics of the electron flux in the LLBL, cusp, and equatorward section of the PM were in this case too similar to allow to distinguish between them by using incoherent scatter radar measurements only.

Watermann, Jurgen↗

Origin and maintenance of the oxygen torus in Saturn's magnetosphere

Observations of thermal ions in Saturn's inner magnetosphere suggest distributed local sources rather than diffusive mass loading from a source located further out. We suggest that the plasma is produced and maintained mainly by 'self-sputtering' of E ring dust. Sputtered particles are 'picked up' by the planetary magnetospheric field and accelerated to corotation energies (of the order of 8 eV/amu). The sputter yield for oxygen on ice at, for example, 120 eV is about 5, which implies that an avalanche of self-sputtering occurs. The plasma density is built up until it is balanced by local losses, presumably pitch angle scattering into the loss cone and absorption in the planet's ionosphere. The plasma density determines the distribution of dust in the E ring through plasma drag. Thus a feedback mechanism between the plasma and the E ring dust is established. The model accounts for the principal plasma observations and simultaneously the radial optical depth profile of the E ring.

Morfill, G. E.↗

Details of the neutral energy distribution and ionization source using spectrally resolved Balmer-alpha measurements on DIII-D

Spectrally resolved passive Balmer- α (D- α, H- α) measurements from the DIII-D 16 channel edge main-ion charge exchange recombination system confirm the presence of higher energy neutrals (“thermal” neutrals) in addition to the cold neutrals that recycle off the walls in the edge region of DIII-D plasmas. Charge exchange between thermal ions and edge neutrals transfers energy and momentum between the populations giving rise to thermal neutrals with energies approximating the ions in the pedestal region. Multiple charge exchange events in succession allow an electron to effectively take a random walk, transferring from ion to ion, providing a pathway of increasing energy and velocity, permitting a neutral to get deeper into the plasma before a final ionization event that contributes to the ion and electron particle fueling. Spectrally resolved measurements provide information about the density and velocity distribution of these neutrals, which has been historically valuable for validating Monte Carlo neutral models, which include the multi stage charge exchange dynamics. Here, in this study, a multi-channel set of such measurements is used to specifically isolate the details of the thermal neutrals that are responsible for fueling inside the pedestal top. Being able to separate the thermal from the cold emission overcomes several challenges associated with optical filter-based neutral density measurements. The neutral dynamics, deeper fueling by the thermal neutrals, and spectral measurement are modeled with the FIDASIM Monte Carlo collisional radiative code, which also produces synthetic spectra with a shape that is in close agreement with the measurements. By scaling the number of neutrals in the simulation to match the intensity of the thermal emission, we show it is possible to obtain local neutral densities and ionization source rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Consequences of Asteroid Characterization on the State of Knowledge about Inferred Physical Properties and Impact Risk

Physical characteristics of Near-Earth Objects (NEOs) are essential inputs to planetary defense assessments. The size, density, and strength of an NEO are critical inputs to modeling behavior during atmospheric entry as well as assessing the risk of impact. Similarly, knowledge of the physical characteristics of an object are necessary to evaluate the probable result of a mitigation mission. Usually, these attributes cannot be directly measured, but increasingly sophisticated methods have been developed to infer physical properties from related measurements of asteroids, meteors, and/or meteorites. Fortuitously, some of these measurements have been obtained for enough NEOs to elucidate the distribution of values across the sampled population. However, the situation becomes more challenging when considering a specific asteroid, since it is unlikely that all the relevant measurements have been made for any given object. We have developed a Bayesian network that can combine available information about a particular NEO with knowledge of the larger population to infer probabilistic values and uncertainties for physical characteristics of interest. Distributions of asteroid population albedos, taxonomic classes, and macroporosities, along with meteorite density distributions and associations between taxonomic classes and meteorite classes, provide the default distributions for the network’s parameter nodes. The inference network links parameters for each virtual asteroid either deterministically or probabilistically as appropriate, and eliminates any unphysical combinations of parameters. Within the context of planetary defense, our Bayesian network can be used to constrain the ranges of likely impactor properties, which can subsequently reduce the uncertainty in modelling of atmospheric entry, mitigation efficacy, and impact risk assessment. When additional measurements become available for a specific object, the network incorporates those measurements to generate virtual asteroids with property distributions that are consistent with the measurements. We will use the 2023 PDC scenario to demonstrate how the inference network can be combined with plausible characterization measurements to refine the state of knowledge about likely combinations of physical parameters and the resulting impact risk.

risk assessment↗

Consequences of Asteroid Characterization on the State of Knowledge about Inferred Physical Properties and Impact Risk

Physical characteristics of Near-Earth Objects (NEOs) are essential inputs to planetary defense assessments. The size, density, and strength of an NEO are critical inputs to modeling behavior during atmospheric entry as well as assessing the risk of impact. Similarly, knowledge of the physical characteristics of an object are necessary to evaluate the probable result of a mitigation mission. Usually, these attributes cannot be directly measured, but increasingly sophisticated methods have been developed to infer physical properties from related measurements of asteroids, meteors, and/or meteorites. Fortuitously, some of these measurements have been obtained for enough NEOs to elucidate the distribution of values across the sampled population. However, the situation becomes more challenging when considering a specific asteroid, since it is unlikely that all the relevant measurements have been made for any given object. We have developed a Bayesian network that can combine available information about a particular NEO with knowledge of the larger population to infer probabilistic values and uncertainties for physical characteristics of interest. Distributions of asteroid population albedos, taxonomic classes, and macroporosities, along with meteorite density distributions and associations between taxonomic classes and meteorite classes, provide the default distributions for the network’s parameter nodes. The inference network links parameters for each virtual asteroid either deterministically or probabilistically as appropriate, and eliminates any unphysical combinations of parameters. Within the context of planetary defense, our Bayesian network can be used to constrain the ranges of likely impactor properties, which can subsequently reduce the uncertainty in modelling of atmospheric entry, mitigation efficacy, and impact risk assessment. When additional measurements become available for a specific object, the network incorporates those measurements to generate virtual asteroids with property distributions that are consistent with the measurements. We will use the 2023 PDC scenario to demonstrate how the inference network can be combined with plausible characterization measurements to refine the state of knowledge about likely combinations of physical parameters and the resulting impact risk.

risk assessment↗

Conditional diffusion machine-learning framework for mapping valence electron distribution from convergent beam electron diffraction

Quantitative convergent beam electron diffraction (CBED) enables determination of aspherical valence electron distributions through refinement of low-order structure factors, which are highly sensitive to chemical bonding and charge density variations. However, conventional quantitative CBED (QCBED) requires solving a highly nonlinear inverse problem with many coupled parameters, and computationally intensive dynamical diffraction calculations, making it time-consuming and difficult to apply to complex systems. More broadly, reconstructing charge density and orbital electron distribution from diffraction data has long been a central challenge in both x-ray and electron crystallography. Here, in this study, we introduce an artificial-intelligence (AI)-based framework that replaces traditional refinement with a data-driven inverse solver. Using a large synthetic CBED dataset generated by Bloch-wave simulations, we train a conditional diffusion model to directly infer crystal structural parameters and multipole density formalism parameters, and hence valence electron distributions, from CBED patterns alone. By learning from forward simulations across realistic parameter space, the model effectively solves the inverse problem. Compared with direct regression approaches, the diffusion-based framework provides posterior parameter distributions for rigorous uncertainty quantification while preserving quantitative fidelity and reducing analysis time by orders of magnitude. By eliminating the need for external single-crystal x-ray diffraction data and complex nonlinear refinement, this approach enables practical, high-throughput, and in situ quantitative CBED, enabling real-time mapping of valence electron distributions and their correlation with functional responses in quantum and energy materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Probing Electrolyte Influence on CO 2 Reduction in Aprotic Solvents

Selective CO 2 capture and electrochemical conversion are important tools in the fight against climate change. Industrially, CO 2 is captured using a variety of aprotic solvents due to their high CO 2 solubility. However, most research efforts on electrochemical CO 2 conversion use aqueous media and are plagued by competing hydrogen evolution reaction (HER) from water breakdown. Fortunately, aprotic solvents can circumvent HER, making it important to develop strategies that enable integrated CO 2 capture and conversion. However, the influence of ion solvation and solvent selection within nonaqueous electrolytes for efficient and selective CO 2 reduction is unclear. In this work, we show that the bulk solvation behavior within the nonaqueous electrolyte can control the CO2 reduction reaction and product distribution occurring at the catalyst-electrolyte interface. We study different tetrabutylammonium (TBA) salts in two electrolyte systems with glyme ethers (e.g., 1,2 dimethoxyethane or DME) and dimethyl sulfoxide (DMSO) as a low and high dielectric constant medium, respectively. Using spectroscopic tools, we quantify the fraction of ion pairs that forms within the electrolyte. Also, we show how ion pair formation is prevalent in DME and is dependent on the anion type. More importantly, we show that as ion pair formation decreases within the electrolyte, CO 2 current densities increase, and a higher CO Faradaic efficiency is observed at low overpotentials. Meanwhile, in an electrolyte medium where the ion pair fraction does not change with the anion type (such as in DMSO), a smaller influence of solvation is observed on CO 2 current densities and product distribution. By directly coupling bulk solvation to interfacial reactions and product distribution, we showcase the importance and utility of controlling the reaction microenvironment in tuning the electrocatalytic reaction pathways. Insights gained from this work will enable novel electrolyte designs for efficient and selective CO 2 conversion to desired fuels and chemicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Initial energy-momentum to final flow: A general framework for heavy-ion collisions

The evolution of a relativistic heavy-ion collision is typically understood as a process that transmutes the initial geometry of the system into the final momentum distribution of observed hadrons, which can be described via a cumulant expansion of the initial distribution of energy density and is represented at leading order as the well-known eccentricity scaling of anisotropic flow. We extend this framework to include the contribution from initial momentum-space properties, as encoded in other components of the energy-momentum tensor. Furthermore, we confirm the validity of the framework in state-of-the-art hydrodynamic simulations of large and small systems. With this framework, it is possible to separate the effects of early time dynamics from those of final-state evolution, even in the case when the distribution of energy does not fully determine subsequent evolution, as for example, in small systems. Specifically, we answer the question of when and how azimuthal correlations from the initial state survive to the final state. In very small systems such as 𝑝−𝑝, for example, initial momentum degrees of freedom dominate over energy. Thus, even if the system forms a quark-gluon plasma that is well described by hydrodynamics, the usual hydrodynamic picture of the transmutation of initial geometry to final momentum anisotropy is broken. Nevertheless, we show that the hydrodynamic response to the full energy-momentum tensor can be well understood in a similar manner as larger systems. Additionally, this framework elucidates the generic features of the system's evolution that are responsible for the impressive success of hydrodynamic simulations, but which may still hold even in cases where hydrodynamics is not applicable.

Relativistic heavy-ion collisions↗

Photospheric Current Spikes and Their Possible Association with Flares - Results from an HMI Data Driven Model

A data driven, near photospheric magnetohydrodynamic model predicts spikes in the horizontal current density, and associated resistive heating rate per unit volume Q. The spikes appear as increases by orders of magnitude above background values in neutral line regions (NLRs) of active regions (ARs). The largest spikes typically occur a few hours to a few days prior to M or X flares. The spikes correspond to large vertical derivatives of the horizontal magnetic field. The model takes as input the photospheric magnetic field observed by the Helioseismic & Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO) satellite. This 2.5 D field is used to determine an analytic expression for a 3 D magnetic field, from which the current density, vector potential, and electric field are computed in every AR pixel for 14 ARs. The field is not assumed to be force-free. The spurious 6, 12, and 24 hour Doppler periods due to SDO orbital motion are filtered out of the time series of the HMI magnetic field for each pixel using a band pass filter. The subset of spikes analyzed at the pixel level are found to occur on HMI and granulation scales of 1 arcsec and 12 minutes. Spikes are found in ARs with and without M or X flares, and outside as well as inside NLRs, but the largest spikes are localized in the NLRs of ARs with M or X flares. The energy to drive the heating associated with the largest current spikes comes from bulk flow kinetic energy, not the electromagnetic field, and the current density is highly non-force free. The results suggest that, in combination with the model, HMI is revealing strong, convection driven, non-force free heating events on granulation scales, and that it is plausible these events are correlated with subsequent M or X flares. More and longer time series need to be analyzed to determine if such a correlation exists. Above an AR dependent threshold value of Q, the number of events N(Q) with heating rates greater than or equal to Q obeys a scale invariant power law distribution for each AR given by N(Q) varies Q(sup -s), where 0.40 less than or equal to S less than or equal to 0.53, with a mean and standard deviation across the 14 ARs of 0.47 and 0.045, showing there is little variation of S from one AR to another. These properties of N(Q) are in close agreement with those of the distribution N(E) for the total energy E of solar flares, determined from observations to be N(E) = constant x E(sup -alpha). From observations of nanoflares in the 0.7 to 4 MK range, and from observations of flares in hard X-rays, it is found that 0.51 less than or equal to alpha less than or equal to 0.57, and 0.4 less than or equal to alpha less than or equal to 0.6, respectively (Crosby et al. 1993, Sol. Phys., 143, 275; Aschwanden & Parnell 2002, ApJ, 572, 1048). Observations also show that, as is found here for the exponent S, there is little variation of alpha with AR (Wheatland 2000, ApJ, 532, 1209), indicating N(E) and N(Q) are largely independent of individual properties of ARs such as area, total magnetic flux, and distribution of current density (i.e. non-potentiality). Therefore the power law scaling of the photospheric heating rate Q computed here on granulation scales is essentially identical to that found for coronal observations of flare energies on scales 1-2 orders of magnitude larger. This suggests the physical mechanisms that cause Q and coronal flares are closely related. It seems likely that Q is the signature of a magnetic reconnection process in an energy range and volume orders of magnitude smaller than those of flares. In this context, at least the larger spikes in Q might be signatures of UV photospheric or lower chromospheric bombs in which plasma is heated to temperatures approximately 10(exp -5) K (Peter et al. 2014, Science 346, 1255726; Judge 2015, ApJ, 808, 116). In addition, lattice based avalanche simulations of flare energy release predict 0.4 less than or equal to alpha less than or equal to 0.5, while analytic, fractal-diffusive self-organized criticality models predict 0.4 less than or equal to alpha less than or equal to 0.67, in excellent agreement with observations, and the results presented here (Aschwanden & Parnell 2002, ApJ, 572, 1048; Aschwanden 2012, A&A, 539, A2; Aschwanden 2013, in "Self Organized Criticality Systems"; Aschwanden et al. 2016, SSR, 198, 47).

Goodman, Michael↗

Organic matter distribution in the icy environments of Taylor Valley, Antarctica

Glaciers can accumulate and release organic matter affecting the structure and function of associated terrestrial and aquatic ecosystems. Here we analyzed 18 ice cores collected from six locations in Taylor Valley (McMurdo Dry Valleys), Antarctica to determine the spatial abundance and quality of organic matter, and the spatial distribution of bacterial density and community structure from the terminus of the Taylor Glacier to the coast (McMurdo Sound). Our results showed that dissolved and particulate organic carbon (DOC and POC) concentrations in the ice core samples increased from the Taylor Glacier to McMurdo Sound, a pattern also shown by bacterial cell density. Fluorescence Excitation Emission Matrices Spectroscopy (EEMs) and multivariate parallel factor (PARAFAC) modeling identified one humic-like (C1) and one protein-like (C2) component in ice cores whose fluorescent intensities all increased from the Polar Plateau to the coast. The fluorescence index showed that the bioavailability of dissolved organic matter (DOM) also decreased from the Polar Plateau to the coast. Partial least squares path modeling analysis revealed that bacterial abundance was the main positive biotic factor influencing both the quantity and quality of organic matter. Marine aerosol influenced the spatial distribution of DOC more than katabatic winds in the ice cores. Certain bacterial taxa showed significant correlations with DOC and POC concentrations. Collectively, our results show the tight connectivity among organic matter spatial distribution, bacterial abundance and meteorology in the McMurdo Dry Valley ecosystem.

54 ENVIRONMENTAL SCIENCES↗

A non-gaussian model of continuous atmospheric turbulence for use in aircraft design

A non-Gaussian model of atmospheric turbulence is presented and analyzed. The model is restricted to the regions of the atmosphere where the turbulence is steady or continuous, and the assumptions of homogeneity and stationarity are justified. Also spatial distribution of turbulence is neglected, so the model consists of three independent, stationary stochastic processes which represent the vertical, lateral, and longitudinal gust components. The non-Gaussian and Gaussian models are compared with experimental data, and it is shown that the Gaussian model underestimates the number of high velocity gusts which occur in the atmosphere, while the non-Gaussian model can be adjusted to match the observed high velocity gusts more satisfactorily. Application of the proposed model to aircraft response is investigated, with particular attention to the response power spectral density, the probability distribution, and the level crossing frequency. A numerical example is presented which illustrates the application of the non-Gaussian model to the study of an aircraft autopilot system. Listings and sample results of a number of computer programs used in working with the model are included.

Reeves, P. M.↗