Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “higher order beam dynamics”

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 19 records

Review of Slow-Wave Structures

The majority of recent theoretical and experimental reports published in the literature dealing with helical slow-wave structures focus on the dispersion characteristics and their effects due to the finite helix wire thickness and attenuation, dielectric loading, metal loading, and the introduction of plasma. In many papers, an effective dielectric constant is used to take into account helix wire dimensions and conductivity losses, while the propagation constant of the signal and the interaction impedance of the structure are found to depend on the surface resistivity of the helix. Also, various dielectric supporting rods are simulated by one or several uniform cylinders having an effective dielectric constant, while metal vane loading and plasma effects are incorporated in the effective dielectric constant. The papers dealing with coupled cavities and folded or loaded wave guides describe equivalent circuit models, efficiency enhancement, and the prediction of instabilities for these structures. Equivalent circuit models of various structures are found using computer software programs SUPERFISH and TOUCHSTONE. Efficiency enhancement in tubes is achieved through dynamic velocity and phase adjusted tapers using computer techniques. The stability threshold of unwanted antisymmetric and higher order modes is predicted using SOS and MAGIC codes and the dependence of higher order modes on beam conductance, section length, and effective Q of a cavity is shown.

Wallett, Thomas M.↗

Nonlinear flexure and torsion of rotating beams, with application to helicopter rotor blades. I - Formulation. II - Response and stability results

The dynamic response and aeroelastic stability of rotating beams such as helicopter blades is investigated analytically. The Hamilton principle is used to formulate the equations of motion for extensional and inextensional beams with precone angles and variable pitch angles, taking higher-order nonlinearities into account. The derivation of the equations and their approximate solution by a Galerkin procedure are explained in detail, and numerical results of equilibrium solutions and stability analyses are presented graphically.

Crespo Da Silva, M. R. M.↗

Development of a composite tailoring procedure for airplane wing

The development of a composite wing box section using a higher order-theory is proposed for accurate and efficient estimation of both static and dynamic responses. The theory includes the effect of through-the-thickness transverse shear deformations which is important in laminated composites and is ignored in the classical approach. The box beam analysis is integrated with an aeroelastic analysis to investigate the effect of composite tailoring using a formal design optimization technique. A hybrid optimization procedure is proposed for addressing both continuous and discrete design variables.

Chattopadhyay, Aditi↗

Evaluation of several modal methods for transient structural analysis

The convergence properties of three modal methods for solving structural transient response (STR) problems are examined by comparing the accuracy of the total number of modes used to represent the responses. The goal is to identify a method which involves as few degrees of freedom as possible. Two mode superposition methods, the mode displacement method (MDM) and the mode acceleration method (MAM), are fitted with equivalent expressions for damped response. An error norm is employed to compare the accuracies of the methods in calculating displacements, moments, and shear forces for a viscously-damped cantilevered beam experiencing several dynamic loading conditions and levels of damping. The MAM converges quicker and uses fewer dof, mainly because it has a pseudo-static response term which accounts for some of the flexibility of higher modes. A third method, which integrates the convolution integral iteratively over time, generates a higher order approximation to the STR with a reduced number of modes. It converges fastest among the three, but requires definition of an appropriate forcing function. Finally, the MAM and MDM models are applied to study the effects of the forcing frequency and the initial conditions.

Camarda, Charles J.↗

A new pre-loaded beam geometric stiffness matrix with full rigid body capabilities

Space structures, such as the Space Station solar arrays, must be extremely light-weight, flexible structures. Accurate prediction of the natural frequencies and mode shapes is essential for determining the structural adequacy of components, and designing a controls system. The tension pre-load in the 'blanket' of photovoltaic solar collectors, and the free/free boundary conditions of a structure in space, causes serious reservations on the use of standard finite element techniques of solution. In particular, a phenomenon known as 'grounding', or false stiffening, of the stiffness matrix occurs during rigid body rotation. The authors have previously shown that the grounding phenomenon is caused by a lack of rigid body rotational capability, and is typical in beam geometric stiffness matrices formulated by others, including those which contain higher order effects. The cause of the problem was identified as the force imbalance inherent in the formulations. In this paper, the authors develop a beam geometric stiffness matrix for a directed force problem, and show that the resultant global stiffness matrix contains complete rigid body mode capabilities, and performs very well in the diagonalization methodology customarily used in dynamic analysis.

Bosela, P. A.↗

Structural dynamic model obtained from flight use with piloted simulation and handling qualities analysis

The ability to use flight data to determine an aircraft model with structural dynamic effects suitable for piloted simulation. and handling qualities analysis has been developed. This technique was demonstrated using SR-71 flight test data. For the SR-71 aircraft, the most significant structural response is the longitudinal first-bending mode. This mode was modeled as a second-order system, and the other higher order modes were modeled as a time delay. The distribution of the modal response at various fuselage locations was developed using a uniform beam solution, which can be calibrated using flight data. This approach was compared to the mode shape obtained from the ground vibration test, and the general form of the uniform beam solution was found to be a good representation of the mode shape in the areas of interest. To calibrate the solution, pitch-rate and normal-acceleration instrumentation is required for at least two locations. With the resulting structural model incorporated into the simulation, a good representation of the flight characteristics was provided for handling qualities analysis and piloted simulation.

Powers, Bruce G.↗

Derivation of a Tappered p-Version Beam Finite Element

A tapered p-version beam finite element suitable for dynamic applications is derived. The taper in the element is represented by allowing the area moments of inertia to vary as quartic polynomials along the length of the beam, and the cross-sectional area to vary as a quadratic polynomial. The p-version finite-element characteristics are implemented through a set of polynomial shape functions. The lower-order shape functions are identical to the classical cubic and linear shape functions normally associated with a beam element. The higher-order shape functions are a hierarchical set of polynomials that are integrals of orthogonal polynomials. Explicit expressions for the mass and stiffness matrices are presented for an arbitrary value of p. The element has been verified to be numerically stable using shape functions through 22nd order.

Hinnant, Howard E.↗

Rarefied gas dynamics; Proceedings of the Eighth International Symposium, Stanford University, Stanford, Calif., July 10-14, 1972

Recent studies in the theory of rarefied gases (kinetic theory, shock structure, molecular forces and interactions, evaporation and condensation, expansions and jets, flow past bodies, gas-solid interactions) are presented along with descriptions of numerous experiments in these areas. Some of the topics covered include: singular limits of solutions of Boltzmann's equations, the Boltzmann and Klimontovich formalisms with reference to turbulent resolution, higher-order kinetic theory for the Lorentz gas, intermolecular potentials from shock structure experiments, cooling of particles in aerosol beams, hypersonic-subsonic transition in spherically symmetric expanding flows, absolute velocity measurements in a rarefied gas flow by an ion time-of-flight technique, the product distribution of gas-solid reactions from absolute rate theory, and relaxation of a gas with internal degrees of freedom. Individual items are announced in this issue.

Karamcheti, K.↗

Current Flow and Pair Creation at Low Altitude in Rotation-Powered Pulsars' Force-Free Magnetospheres: Space Charge Limited Flow

We report the results of an investigation of particle acceleration and electron-positron plasma generation at low altitude in the polar magnetic flux tubes of rotation-powered pulsars, when the stellar surface is free to emit whatever charges and currents are demanded by the force-free magnetosphere. We apply a new 1D hybrid plasma simulation code to the dynamical problem, using Particle-in-Cell methods for the dynamics of the charged particles, including a determination of the collective electrostatic fluctuations in the plasma, combined with a Monte Carlo treatment of the high-energy gamma-rays that mediate the formation of the electron-positron pairs.We assume the electric current flowing through the pair creation zone is fixed by the much higher inductance magnetosphere, and adopt the results of force-free magnetosphere models to provide the currents which must be carried by the accelerator. The models are spatially one dimensional, and designed to explore the physics, although of practical relevance to young, high-voltage pulsars. We observe novel behaviour (a) When the current density j is less than the Goldreich-Julian value (0 < j/j(sub GJ) < 1), space charge limited acceleration of the current carrying beam is mild, with the full Goldreich-Julian charge density comprising the charge densities of the beam and a cloud of electrically trapped particles with the same sign of charge as the beam. The voltage drops are of the order of mc(sup 2)/e, and pair creation is absent. (b) When the current density exceeds the Goldreich-Julian value (j/j(sub GJ) > 1), the system develops high voltage drops (TV or greater), causing emission of curvature gamma-rays and intense bursts of pair creation. The bursts exhibit limit cycle behaviour, with characteristic time-scales somewhat longer than the relativistic fly-by time over distances comparable to the polar cap diameter (microseconds). (c) In return current regions, where j/j(sub GJ) < 0, the system develops similar bursts of pair creation. These discharges are similar to those encountered in previous calculations by Timokhin of pair creation when the surface has a high work function and cannot freely emit charge. In cases (b) and (c), the intermittently generated pairs allow the system to simultaneously carry the magnetospherically prescribed currents and adjust the charge density and average electric field to force-free conditions. We also elucidate the conditions for pair creating beam flow to be steady (stationary with small fluctuations in the rotating frame), finding that such steady flows can occupy only a small fraction of the current density parameter space exhibited by the force-free magnetospheric model. The generic polar flow dynamics and pair creation are strongly time dependent. The model has an essential difference from almost all previous quantitative studies, in that we sought the accelerating voltage (with pair creation, when the voltage drops are sufficiently large; without, when they are small) as a function of the applied current.

Timokhin, A. N.↗

Higher Order Chemistry Models in the CFD Simulation of Laser-Ablated Carbon Plumes

Production of single-walled carbon nanotubes (SWNT) has taken place for a number of years and by a variety of methods such-as laser ablation, chemical vapor deposition, and arc-jet ablation. Yet, little is actually understood about the exact chemical kinetics and processes that occur in SWNT formation. In recent time, NASA Johnson Space Center has devoted a considerable effort to the experimental evaluation of the laser ablation production process for SWNT originally developed at Rice University. To fully understand the nature of the laser ablation process it is necessary to understand the development of the carbon plume dynamics within the laser ablation oven. The present work is a continuation of previous studies into the efforts to model plume dynamics using computational fluid dynamics (CFD). The ultimate goal of the work is to improve understanding of the laser ablation process, and through that improved understanding, refine the laser ablation production of SWNT. Fig. 1 shows a basic schematic of the laser-ablation oven at NASA-JSC. Construction of the facility is simple in concept. Two concentric quartz tubes of 1.5 mm thickness form the inner and outer tubes with inside diameters of 2.2 and 5.08 cm respectively. At one end of the inner tube are located two 60 Hz pulsed lasers operating at 1064 nm and 532 nm wavelength with beam diameters of 5 mm aligned coaxially with the longitudinal axis of the inner quartz tube. For standard nanotube production runs, a 10 ns 532 nm pulse is followed 50 ns later by a 10 ns 1064 nm pulse. Each pulse is of 300 mJ energy. A target of carbon graphite with approximately 1% nickel and cobalt catalysts is located at the other end of the inner quartz tube. In the ordinary processing of SWNT, a base flow of 100 sccm of argon is maintained from the laser location and exits past the carbon target at a pressure of 66.7 kPa. These conditions yield a baseline mass flow through the chamber of 2.723x10(exp -6)kg/s of argon. The whole oven facility is heated to a temperature of 1473 K prior to nanotube production runs. Upon laser irradiation, part of the carbon target ablates immediately and forms a carbon vapor plume that penetrates into the argon base flow towards the laser initially at supersonic velocities. In the previous studies either a single carbon species, C3, was used to model the plume development, or a simplified 11 species carbon model reduced down from the Krestinin and Moravsky model for full fullerene chemistry was used. While both of these options yielded significant results, it was felt that the actual chemistry occ&g in the carbon plume might have a greater affect on the plume than assumed. Indeed, in the earlier attempts at modeling the carbon plume, several thermophysical characteristics could never be matched to experimental observations of plume development - mainly the propagation distance of the plume itself. In the present study, two additional chemistry models will be used to duplicate the previous studies simulations of the carbon plume. The first chemistry model used in this study is again a reduced form of the Krestinin and Moravsky rates. However the highest order carbon species allowed has been increased from C6 to C30 - therefore allowing the simulation of up to a half of the standard C60 fullerene. The second chemistry model investigated is a reduced form of a full carbon nanotube model developed at NASA-JSC. The C30 studies have already been accomplished at the present time, and the reduced SWNT model studies are currently underway. To pursue the current study, one sacrifice had to be made in that the simulation grid spacing had to be increased from 0.5 mm spacing to 1 mm spacing for the sake of computational efficiency since computational effort is proportional to the square of the number of grid points multiplied by the number of species considered. propagation that is far more in line with the experimental results observed by Puretzky et al as shown in Fig. 2. e C6 studies had yielded a far greater propagation in previous studies. In addition, chemical species development with the C30 model indicates that many higher order carbon species are produced outside of the plume proper (indicated by plotting contours of the background argon concentrations in Fig. 3) - this result was not observed in previous studies. In fact, some species primarily occurred outside of the plume itself - as shown for C27 in Fig. 4 when it is compared to Fig. 3. It could be asserted that this has occurred because all of the C27 in the plume had already been consumed in the formation of C30, but this does not seem to be indicated over time. Several other factors that arose in the previous studies have also been made more clear by the use of the higher order chemical models - one being that the use of c6 as an indicator species was mistaken. C6 is the only carbon species in the previous studies that was not injected into the flowfield as a boundary condition; it was therefore hoped that this species would provide insight into the formation of higher order carbon species for comparison to full SWNT production. But, when the plot of total mass in the plume is examined on a species by species basis in Fig. 5, it is seen that Cg was a fairly insignificant contributor to the total carbon mass in the plume and would not provide information on higher order carbon formation. the thermophysical characteristics of the carbon plume as well as simulate the carbon plume using the reduced SWNT model to provide an even better simulation of full chemistry effects upon plume propagation.

Scott, C. D.↗

Relationship of Topside Ionospheric Ion Outflows to Auroral Forms and Precipitation, Plasma Waves, and Convection Observed by Polar

The POLAR satellite often observes upflowing ionospheric ions (UFIs) in and near the aurora] oval on southern perigee (approx. 5000 km altitude) passes. We present the UFI features observed by the thermal ion dynamics experiment (TIDE) and the toroidal imaging mass angle spectrograph (TIMAS) in the dusk-dawn sector under two different geomagnetic activity conditions in order to elicit their relationships with auroral forms, wave emissions, and convection pattern from additional POLAR instruments. During the active interval, the ultraviolet imager (UVI) observed a bright discrete aurora on the duskside after the substorm onset and then observed a small isolated aurora form and diffuse auroras on the dawnside during the recovery phase. The UFIs showed clear conic distributions when the plasma wave instrument (PWI) detected strong broadband wave emissions below approx. 10 kHz, while no significant auroral activities were observed by UVI. At higher latitudes, the low-energy UFI conics gradually changed to the polar wind component with decreasing intensity of the broadband emissions. V-shaped auroral kilometric radiation (AKR) signatures observed above -200 kHz by PWI coincided with the region where the discrete aurora and the UFI beams were detected. The latitude of these features was lower than that of the UFI conics. During the observations of the UFI beams and conics, the lower-frequency fluctuations observed by the electric field instrument were also enhanced, and the convection directions exhibited large fluctuations. It is evident that large electrostatic potential drops produced the precipitating electrons and discrete auroras, the UFI beams, and the AKR, which is also supported by the energetic plasma data from HYDRA. Since the intense broadband emissions were also observed with the UFIs, the ionospheric ions could be energized transversely before or during the parallel acceleration due to the potential drops.

Hirahara, M.↗

Relationship of Topside Ionospheric Ion Outflows to Auroral Forms and Precipitations, Plasma Waves, and Convection Observed by POLAR

The POLAR satellite often observes upflowing ionospheric ions (UFls) in and near the auroral oval on southern perigee (approximately 5000 km altitude) passes. We present the UFI features observed by the thermal ion dynamics experiment (TIDE) and the toroidal imaging mass-angle spectrograph (TIMAS) in the dusk-dawn sector under two different geomagnetic activity conditions in order to elicit their relationships with auroral forms, wave emissions, and convection pattern from additional POLAR instruments. During the active interval, the ultraviolet imager (UVI) observed a bright discrete aurora on the dusk side after the substorm onset and then observed a small isolated aurora form and diffuse auroras on the dawn side during the recovery phase. The UFls showed clear conic distributions when the plasma wave instrument (PWI) detected strong broadband wave emissions below approximately 10 kHz, while no significant auroral activities were observed by UVI. At higher latitudes, the low-energy UFI conics gradually changed to the polar wind component with decreasing intensity of the broadband emissions. V-shaped auroral kilometric radiation (AKR) signatures observed above approximately 200 kHz by PWI coincided with the region where the discrete aurora and the UFI beams were detected. The latitude of these features was lower than that of the UFI conics. During the observations of the UFI beams and conics, the lower-frequency fluctuations observed by the electric field instrument (EFI) were also enhanced, and the convection directions exhibited large fluctuations. It is evident that large electrostatic potential drops produced the precipitating electrons and discrete auroras, the UFI beams, and the AKR, which is also supported by the energetic plasma data from HYDRA. Since the intense broadband emissions were also observed with the UFIs. the ionospheric ions could be energized transversely before or during the parallel acceleration due to the potential drops.

Hirahara, M.↗

The dynamics and control of large flexible space structures. Part B: Development of continuum model and computer simulation

The equations of motion of an arbitrary flexible body in orbit were derived. The model includes the effects of gravity with all its higher harmonics. As a specific example, the motion of a long, slender, uniform beam in circular orbit was modelled. The example considers both the inplane and three dimensional motion of the beam in orbit. In the case of planar motion with only flexible vibrations, the pitch motion is not influenced by the elastic motion of the beam. For large values of the square of the ratio of the structural modal frequency to the orbital angular rate the elastic motion was decoupled from the pitch motion. However, for small values of the ratio and small amplitude pitch motion, the elastic motion was governed by a Hill's 3 term equation. Numerical simulation of the equation indicates the possibilities of instability for very low values of the square of the ratio of the modal frequency to the orbit angular rate. Also numerical simulations of the first order nonlinear equations of motion for a long flexible beam in orbit were performed. The effect of varying the initial conditions and the number of modes was demonstrated.

Bainum, P. M.↗

Test Methodology Development for Experimental Structural Assessment of ASC Planar Spring Material for Long-Term Durability

A vibration-based testing methodology has been developed that will assess fatigue behavior of the metallic material of construction for the Advanced Stirling Convertor displacer (planar) spring component. To minimize the testing duration, the test setup is designed for base-excitation of a multiplespecimen arrangement, driven in a high-frequency resonant mode; this allows completion of fatigue testing in an accelerated period. A high performance electro-dynamic exciter (shaker) is used to generate harmonic oscillation of cantilever beam specimens, which are clasped on the shaker armature with specially-designed clamp fixtures. The shaker operates in closed-loop control with dynamic specimen response feedback provided by a scanning laser vibrometer. A test coordinator function synchronizes the shaker controller and the laser vibrometer to complete the closed-loop scheme. The test coordinator also monitors structural health of the test specimens throughout the test period, recognizing any change in specimen dynamic behavior. As this may be due to fatigue crack initiation, the test coordinator terminates test progression and then acquires test data in an orderly manner. Design of the specimen and fixture geometry was completed by finite element analysis such that peak stress does not occur at the clamping fixture attachment points. Experimental stress evaluation was conducted to verify the specimen stress predictions. A successful application of the experimental methodology was demonstrated by validation tests with carbon steel specimens subjected to fully-reversed bending stress; high-cycle fatigue failures were induced in such specimens using higher-than-prototypical stresses

fatigue testing machines↗

A viscoelastic higher-order beam finite element

A viscoelastic internal variable constitutive theory is applied to a higher-order elastic beam theory and finite element formulation. The behavior of the viscous material in the beam is approximately modeled as a Maxwell solid. The finite element formulation requires additional sets of nodal variables for each relaxation time constant needed by the Maxwell solid. Recent developments in modeling viscoelastic material behavior with strain variables that are conjugate to the elastic strain measures are combined with advances in modeling through-the-thickness stresses and strains in thick beams. The result is a viscous thick-beam finite element that possesses superior characteristics for transient analysis since its nodal viscous forces are not linearly dependent an the nodal velocities, which is the case when damping matrices are used. Instead, the nodal viscous forces are directly dependent on the material's relaxation spectrum and the history of the nodal variables through a differential form of the constitutive law for a Maxwell solid. The thick beam quasistatic analysis is explored herein as a first step towards developing more complex viscoelastic models for thick plates and shells, and for dynamic analyses. The internal variable constitutive theory is derived directly from the Boltzmann superposition theorem. The mechanical strains and the conjugate internal strains are shown to be related through a system of first-order, ordinary differential equations. The total time-dependent stress is the superposition of its elastic and viscous components. Equations of motion for the solid are derived from the virtual work principle using the total time-dependent stress. Numerical examples for the problems of relaxation, creep, and cyclic creep are carried out for a beam made from an orthotropic Maxwell solid.

Johnson, Arthur R.↗

Measurement and Prediction of the Thermomechanical Response of Shape Memory Alloy Hybrid Composite Beams

Previous work at NASA Langley Research Center (LaRC) involved fabrication and testing of composite beams with embedded, pre-strained shape memory alloy (SMA) ribbons within the beam structures. That study also provided comparison of experimental results with numerical predictions from a research code making use of a new thermoelastic model for shape memory alloy hybrid composite (SMAHC) structures. The previous work showed qualitative validation of the numerical model. However, deficiencies in the experimental-numerical correlation were noted and hypotheses for the discrepancies were given for further investigation. The goal of this work is to refine the experimental measurement and numerical modeling approaches in order to better understand the discrepancies, improve the correlation between prediction and measurement, and provide rigorous quantitative validation of the numerical analysis/design tool. The experimental investigation is refined by a more thorough test procedure and incorporation of higher fidelity measurements such as infrared thermography and projection moire interferometry. The numerical results are produced by a recently commercialized version of the constitutive model as implemented in ABAQUS and are refined by incorporation of additional measured parameters such as geometric imperfection. Thermal buckling, post-buckling, and random responses to thermal and inertial (base acceleration) loads are studied. The results demonstrate the effectiveness of SMAHC structures in controlling static and dynamic responses by adaptive stiffening. Excellent agreement is achieved between the predicted and measured results of the static and dynamic thermomechanical response, thereby providing quantitative validation of the numerical tool.

Davis, Brian↗

Precision in Assembled Discrete Lattice Space Structures for Next-Generation ISAM Applications

Robust autonomous robotic assembly of large-scale space structures has been a long-term and challenging goal to enable higher quality in-space communication and science instrumentation. For optical observatory support structures and antenna structures, the challenge is the strict dimensional precision requirements (generally RMS surface error) driven by the operational radiation wavelength. Early theoretical work in the area linked RMS surface error of a reflector support truss plate to the natural passive dynamic mode frequencies of the structure, as well as to the error distribution in dimensions of constituent structural elements. Prior robotic truss assembly demonstrations focused on designing ultra-high precision structural elements, joints, and robotic actuators. Recently, an alternative space structure assembly strategy based on a programmable matter approach (NASA ARMADAS) was demonstrated. This approach uses lattice building blocks (voxels) that are reversibly, mechanically joined into a bulk lattice structure by robots that locomote in and on the structure itself. Such a system can achieve high-level autonomy with low computation, robustly assemble utilizing inexpensive and imprecise robots, and efficiently build structures several orders of magnitude larger than the assembly robots. However, for instrumentation support structure applications, the resulting precision of these building block-based lattice structures is not well studied. Since they are demonstrated with many more assembly units than prior art trusses, it is unclear whether the same precision design approaches apply. In this work, we study in simulation the effects of voxel geometry, error in distribution, assembly resolution (module size), and assembly geometry on the error of both beam and plate lattice structures. While average RMS error of a plate increases with plate size, increasing plate thickness quickly collapses RMS error towards a limit that is on the order of the error of the constituent parts. We validate our models against previously published precision measurements of built systems. Results from this study will guide manufacturing precision requirements, as well as designs, for future robotically assembled structural applications and establish feasibility for different applications.

Christine E Gregg↗

Traveling-Wave Tube Amplifier Model to Predict High-Order Modulation Intersymbol Interference

Demands for increased data rates in satellite communications necessitate higher order modulation schemes, larger system bandwidth, and minimum distortion of the modulated signal as it is passed through the traveling wave tube amplifier (TWTA). One type of distortion that the TWTA contributes to is intersymbol interference (ISI), and this becomes particularly disruptive with wide-band, complex modulation schemes. It is suspected that in addition to the dispersion of the TWT, frequency dependent reflections due to mismatches within the TWT are a significant contributor to ISI. To experimentally investigate the effect of these mismatches within the physical TWT on ISI would be prohibitively expensive, as it would require manufacturing numerous amplifiers in addition to the acquisition of the required digital hardware. In an attempt to develop a more accurate model to correlate IS1 with the TWTA and the operational signal, a fully three-dimensional (3D), time-dependent, TWT interaction model has been developed using the electromagnetic particle-in-cell (PIC) code MAFIA (solution of Maxwell's equations by the Finite-Integration-Algorithm). The model includes a user defined slow-wave circuit with a spatially tapered region of loss to implement a sever, and spatially varied geometry (such as helical pitch) to implement a phase velocity taper. The model also includes user defined input/output coupling and an electron beam contained by solenoidal, electrostatic, or periodic permanent magnet (PPM) focusing allowing standard or novel TWTs to be investigated. This model comprehensively takes into account the effects of frequency dependent nonlinear distortions (MAM and AMPM); gain ripple due to frequency dependent reflections at the input/output coupling, severs, and mismatches from dynamic pitch variations; drive induced oscillations; harmonic generation; intermodulation products; and backward waves.

Kory, Carol L.↗