SEARCH · Engineering Papers
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.
Electron acoustic waves present in a two temperature electron plasma can be driven unstable when ion beams propagate along the magnetic field. Both linear theory and numerical simulations indicate that this instability contributes to the generation of broadband electrostatic noise (BEN) in the geomagnetic tail.
Plasma processes dealing with ac and dc electric fields, the formation of ion beams and conics, and electron acceleration are considered, and similarities between simulation results and satellite-based observations are discussed. Electrostatic shock-type electric fields are found to occur near the current sheet edges, and double layers having upward electric fields form inside the sheet and are distinguishable from the large perpendicular electric fields only in wide sheets with thicknesses much greater than the ion Larmor radius. It is found that the most energetic ions have pitch angles near 90 deg, indicating a large perpendicular acceleration of the ions, and that the downward accelerating electrons inside the sheet are neither monoenergetic nor perfectly field aligned.
A new technique for vibration suppression in large space structures is investigated in laboratory experiments on a thin cantilever beam. This technique, called Positive Position Feedback, makes use of generalized displacement measurements to accomplish vibration suppression. Several features of Positive Position Feedback make it attractive for the large space structure control environment: The realization of the controller is simple and straightforward. Global stability conditions can be derived which are independent of the dynamical characteristics of the structure being controlled, i.e., all spillover is stabilizing. The method cannot be destabilized by finite actuator dynamics, and the technique is amenable to a strain-based sensing approach. The experiments control the first six bending modes of a cantilever beam, and make use of piezoelectric materials for actuators and sensors, simulating a piezoelectric active-member. The modal damping ratios are increased by factors ranging from 2 to 130.
The data for polarization phase difference Delta Phi between the HH- and VV-polarized backscattered waves from tree-covered fields were obtained with an airborne synthetic aperture radar at 1.225 GHz. The mean values over tree-covered fields were derived from the images of the phase difference and were examined as a function of incident beam angle from 15 to 55 deg. A theoretical model for simulating these data, based on the electromagnetic wave scatterings from the tree trunk and its branches, both of which are assumed as very long dielectric cylinders was developed. The radius and direction of a tree branch are taken as random variables and are chosen by a Monte Carlo method to encounter the incident waves in producing the scattering events. The Monte Carlo simulated results are in good agreement with the observations within experimental uncertainty.
Observational data from ISEE 1 are analyzed and a one-dimensional electromagnetic hybrid computer simulation is conducted for the heating of solar wind alpha particles in hot diamagnetic cavities (HDCs). In the simulation, which envisions alpha heating by ion-ion instabilities, low beam densities excite the proton/proton right-hand resonant instability that then pitch-angle scatters the beam without significantly heating the alphas. At greater beam densities, the proton/proton nonresonant instability undergoes saturation through a trapping of all three ion components. These results support the Thomsen et al. (1988) hypothesis that the nonresonant instability is the primary source of ion heating in hot diamagnetic cavities.
The primary research effort of this project is the development of analytical methods for the prediction of nonlinear random response of composite aerospace structures subjected to combined acoustic and thermal loads. The progress, accomplishments, and future plans of three random response research topics are discussed, namely acoustics-structure interactions using boundary/finite element methods, nonlinear vibrations of beams and composite plates under harmonic and random excitations, and numerical simulation of the nonlinear response of composite plates under combined thermal and acoustic loading.
One-dimensional electromagnetic particle simulations are used to investigate the characteristics of radiation from electron beams with a large temperature anisotropy, in order to identify the origin of the smooth and bursty radio emissions from Uranus and Neptune observed by the Voyager spacecraft. It is shown that these electron beams, which would typically originate from the sporadic or impulsive injections of energetic electrons, can generate electromagnetic radiation which should be able to escape into the solar wind despite the growth of the electrostatic instability. The amount of radiation with frequencies above the local x mode cutoff increases with the beam speed. It is also proposed that some of the radiation generated below the local x mode cutoff may also be able to escape the plasma and be detected remotely via mode conversion between regions where field-aligned currents produce local perturbations in the magnetic field.
A statistically generated weighting function for a second-order polynomial curve fit of residual functions has been developed. The residual flexibility test method, from which a residual function is generated, is a procedure for modal testing large structures in an external constraint-free environment to measure the effects of higher order modes and interface stiffness. This test method is applicable to structures with distinct degree-of-freedom interfaces to other system components. A theoretical residual function in the displacement/force domain has the characteristics of a relatively flat line in the lower frequencies and a slight upward curvature in the higher frequency range. In the test residual function, the above-mentioned characteristics can be seen in the data, but due to the present limitations in the modal parameter evaluation (natural frequencies and mode shapes) of test data, the residual function has regions of ragged data. A second order polynomial curve fit is required to obtain the residual flexibility term. A weighting function of the data is generated by examining the variances between neighboring data points. From a weighted second-order polynomial curve fit, an accurate residual flexibility value can be obtained. The residual flexibility value and free-free modes from testing are used to improve a mathematical model of the structure. The residual flexibility modal test method is applied to a straight beam with a trunnion appendage and a space shuttle payload pallet simulator.
The Imaging Calorimeter for ACCESS (ICA) utilizes a thin sampling calorimeter concept for direct measurements of high-energy cosmic rays. The ICA design uses arrays of small scintillating fibers to measure the energy and trajectory of the produced cascades. A test instrument has been developed to study the performance of this concept at accelerator energies and for comparison with simulations. Two test exposures have been completed using a CERN test beam. Some results from the accelerator tests are presented.
This paper describes a single-event upset experiment performed at the Los Alamos National Laboratory. A closed-loop control system consisting of a Quad-Redundant Flight Control Computer (FCC) and a B737 simulator was operated while the FCC was exposed to a neutron beam. The purpose of this test was to analyze the effects of neutron bombardment on avionics control systems operating at altitudes where neutron strikes are probable. The neutron energy spectrum produced at the Los Alamos National Laboratory is similar in shape to the spectrum of atmospheric neutrons but much more intense. The higher intensity results in accelerated life tests that are representative of the actual neutron radiation that a FCC may receive over a period of years.
A study is undertaken to develop a methodology for determining the suitability of various high-cycle fatigue models for metallic structures subjected to combined thermal-acoustic loadings. Two features of this problem differentiate it from the fatigue of structures subject to acoustic loading alone. Potentially large mean stresses associated with the thermally pre- and post-buckled states require models capable of handling those conditions. Snap-through motion between multiple post-buckled equilibrium positions introduces very high alternating stress. The thermal-acoustic time history response of a clamped aluminum beam structure with geometric and material nonlinearities is determined via numerical simulation. A cumulative damage model is employed using a rainflow cycle counting scheme and fatigue estimates are made for 2024-T3 aluminum using various non-zero mean fatigue models, including Walker, Morrow, Morrow with true fracture strength, and MMPDS. A baseline zero-mean model is additionally considered. It is shown that for this material, the Walker model produces the most conservative fatigue estimates when the stress response has a tensile mean introduced by geometric nonlinearity, but remains in the linear elastic range. However, when the loading level is sufficiently high to produce plasticity, the response becomes more fully reversed and the baseline, Morrow, and Morrow with true fracture strength models produce the most conservative fatigue estimates.
Muon g –2 experiment at Fermilab aims to measure the magnetic anomaly of the muon in a precision of 140 parts-per-billion (ppb). It is one of the most strongly motivated precision particle physics experiments to unveil the mystery of the long-standing 3σ-4σ tension between the previous measurement of the muon g –2 experiment at BNL and the Standard Model predictions. One of the main challenges of the experiment is to keep the total systematic uncertainty within 100 ppb budget. In this dissertation, a novel way of matching the phase space in the muon g –2 storage ring to reduce the beam dynamics systematic effects with RF electric fields is presented. Systematic simulation studies and hardware development of the RF matching system are covered, followed by the promising analysis results of the commissioning runs.
Satellite observations in recent years have confirmed that the plasma sheet boundary layer is a permanent feature of the earth's magnetotail located between the lobe and central plasma sheet during both quiet and active magnetic periods. Distinct features of the boundary layer include field aligned ion beams and intense electrostatic emissions known as broadband electrostatic noise. Since the plasma sheet boundary layer is a spatial feature of the magnetotail, within it will occur thermal mixing of the resident warm boundary layer plasma with inflowing (convecting) cold ionospheric plasma. A theoretical study involving linear theory and nonlinear numerical particle simulations is presented which examines ion beam instabilities in the presence of a thermally mixed hot and cold background plasma. It is found that the free energy in the ion beams can heat the cool ionospheric plasma to ambient plasma sheet boundary layer temperatures via broadband electrostatic noise. These results, along with recent observational reports that ionospheric outflow can account for measured plasma sheet densities, suggest that the ionospheric role in plasma sheet dynamics and content may be as large as the solar wind.
The inertial confinement fusion program relies upon detailed simulations with inertial confinement fusion (ICF) codes to design targets and to interpret the experimental results. These simulations treat as much physics from essential principles as is practical, including laser deposition, cross beam energy transfer, x-ray production and transport, nonlocal thermal equilibrium kinetics, thermal transport, hydrodynamic instabilities, thermonuclear burn, and transport of reaction products. Improvements in radiation hydrodynamic code capabilities and vast increases in computing power have enabled more realistic, accurate 3D simulations that treat all known asymmetry sources. We describe how numerical simulations helped to guide the program, assess the impediments to breakeven, and optimize every aspect of target design. A preshot simulation of the first National Ignition Facility experiment that surpassed breakeven predicted an increased yield that matches the experimental result, within the preshot predicted uncertainty, with a target gain of 1.5. We will cover the key developments in Lawrence Livermore National Laboratory ICF codes that enabled these simulations and give specific examples of how they helped to guide the program.
Abstract Energetic particle (EP) destabilized Alfvén eigenmode (AE) instabilities are simulated for a DIII-D experimental case with a pulsed neutral beam using a gyro-Landau moments model which introduces EP phase-mixing effects through closure relations. This provides a computationally efficient reduced model which is applied here in the nonlinear regime over timescales that would be difficult to address with more complete models. The long timescale nonlinear evolution and related collective transport losses are examined including the effects of zonal flow/current generation, nonlinear energy cascades, and EP profile flattening. The model predicts frequencies and mode structures that are consistent with experimental observations. These calculations address issues that have not been considered in previous modelling: the EP critical gradient profile evolution in the presence of zonal flows/currents, and the dynamical nature of the saturated state. A strong level of intermittency is present in the predicted instability-driven transport; this is connected to the zonal flow growth and decay cycles and nonlinear energy transfers. Simulation of intermittent AE-enhanced EP transport will be an important issue for the protection of plasma facing components in the next generation of fusion devices.
The interaction of Alfvén eigenmodes (AEs) and energetic particles is one of many important factors determining the success of future tokamaks. In JET, eight in-vessel antennas were installed to actively probe stable AEs with frequencies ranging 25–250 kHz and toroidal mode numbers |n| < 20. During the 2019–2020 deuterium campaign, almost 7500 resonances and their frequencies f 0 , net damping rates γ < 0, and toroidal mode numbers were measured in almost 800 plasma discharges. From a statistical analysis of this database, continuum and radiative damping are inferred to increase with edge safety factor, edge magnetic shear, and when including non-ideal effects. Both stable AE observations and their associated damping rates are found to decrease with |n|. Active antenna excitation is also found to be ineffective in H-mode as opposed to L-mode; this is likely due to the increased edge density gradient’s effect on accessibility and ELM-related noise’s impact on mode identification. A novel measurement is reported of a marginally stable, edge-localized ellipticity-induced AE probed by the antennas during high-power auxiliary heating (ion cyclotron resonance heating and neutral beam injection) up to 25 MW. NOVA-K kinetic-MHD simulations show good agreement with experimental measurements of f 0 , γ, and n, indicating the dominance of continuum and electron Landau damping in this case. Similar experimental and computational studies are planned for the recent hydrogen and ongoing tritium campaigns, in preparation for the upcoming DT campaign.
We use two-dimensional time-dependent hydrodynamical simulations to follow the growth of the Kelvin-Helmholtz (K-H) instability in cooling jets into the nonlinear regime. We focus primarily on asymmetric modes that give rise to transverse displacements of the jet beam. A variety of Mach numbers and two different cooling curves are studied. The growth rates of waves in the linear regime measured from the numerical simulations are in excellent agreement with the predictions of the linear stability analysis presented in the first paper in this series. In the nonlinear regime, the simulations show that asymmetric modes of the K-H instability can affect the structure and evolution of cooling jets in a number of ways. We find that jets in which the growth rate of the sinusoidal surface wave has a maximum at a so-called resonant frequency can be dominated by large-amplitude sinusoidal oscillations near this frequency. Eventually, growth of this wave can disrupt the jet. On the other hand, nonlinear body waves tend to produce low-amplitude wiggles in the shape of the jet but can result in strong shocks in the jet beam. In cooling jets, these shocks can produce dense knots and filaments of cooling gas within the jet. Ripples in the surface of the jet beam caused by both surface and body waves generate oblique shock "spurs" driven into the ambient gas. Our simulations show these shock "spurs" can accelerate ambient gas at large distances from the jet beam to low velocities, which represents a new mechanism by which low-velocity bipolar outflows may be driven by high-velocity jets. Rapid entrainment and acceleration of ambient gas may also occur if the jet is disrupted. For parameters typical of protostellar jets, the frequency at which K-H growth is a maximum (or highest frequency to which the entire jet can respond dynamically) will be associated with perturbations with a period of - 200 yr. Higher frequency (shorter period) perturbations excite waves associated with body modes that produce internal shocks and only small-amplitude wiggles within the jet. The fact that most observed systems show no evidence for large-amplitude sinusoidal oscillation leading to disruption is indicative that the perturbation frequencies are generally large, consistent with the suggestion that pro- tostellar jets arise from the inner regions (r less than 1 AU) of accretion disks.
To support the design of the Second Target Station, that aims to provide the world’s highest peak brightness of cold neutrons, studies that optimize the dimensions of the target and moderators are invaluable. In this work, we investigate the influence of the target dimensions and beam profile on the performance and optimal size of the moderators. We perform optimization runs with detailed MCNP6.2 simulations using high-fidelity unstructured mesh geometries generated from parametrized CAD models. We demonstrate that small changes in target height do not influence the moderator performance if the beam dimensions are chosen adequately. We quantify the effect of beam footprint and target width on the moderator performance and show that the optimal moderator dimensions are insensitive to limited changes in target and beam profile.