Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smooth approximation”

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 703 records · Page 39

Wrinkle Ridges on Mercury and the Moon Within and Outside of Mascons

Found on all terrestrial planets, wrinkle ridges are anticlines formed by thrust faulting and folding resulting from crustal shortening. The MErcury Surface, Space Environment, Geochemistry, and Ranging (MESSENGER) spacecraft's orbital phase returned high resolution images and topographic data of the previously unimaged northern high latitudes of Mercury where there are large expanses of smooth plains deformed by wrinkle ridges. Concurrently, the Lunar Reconnaissance Orbiter (LRO) is obtaining high resolution images and topographic data covering lunar mare wrinkle ridges. These data allow quantitative comparison of the scale of wrinkle ridges in smooth plains volcanic units on Mercury with mare wrinkle ridges. We evaluate the topographic relief of 300 wrinkle ridges within and outside of mascon basins on the Moon and Mercury. Measured wrinkle ridges range from ~112 to 776 m in relief with a mean of ~350 m (median = ~340 m, n = 150) on Mercury and from ~47 to 678 m in relief with a mean of ~198 m (median = ~168 m, n = 150) on the Moon. Wrinkle ridges on Mercury thus are approximately twice as large in mean relief compared to their counterparts on the Moon. The larger scale of Mercury's wrinkle ridges suggests that their formation can be attributed, in part, to global contraction. As global contraction on the Moon is estimated to be an order of magnitude smaller than on Mercury, the smaller scale of lunar wrinkle ridges suggests they most likely form primarily by load induced subsidence of the mare basalt. Wrinkle ridges located in lunar mascon basins and in the Caloris mascon on Mercury are not statistically significantly different in relief than ridges in non-mascon regions, suggesting comparable levels of contractional strain. The fact that mascon basins do not host wrinkle ridges with greater structural relief relative to non-mascon units may indicate the critical role lithospheric thickness plays in controlling subsidence and contraction of thick volcanic sequences on the Moon and Mercury.

Global contraction↗

Coal to Carbon Fiber (C2CF) Continuous Processing for High Value Composites (Final Report)

Coal tar is a condensed and recovered by-product of the coking of metallurgical coal for steel production. The heaviest fraction of distilled coal tar is an isotropic pitch largely used as a binder in the manufacturing of carbonaceous electrodes for primary aluminum smelting and in electric arc furnaces. Coal tar pitch offers high carbon yield upon carbonization. In this project, a process to convert the domestically sourced isotropic coal tar pitch, containing very low particulates (QI = 0.32 wt.%), to form flow-domain mesophase pitch amenable for melt spinning into precursor (green) fibers for carbon fiber, was developed and optimized. The final reproducible processing developed is reviewed in this report, along with several characterizations of the mesophase pitch. A final definition of the characteristics of a ‘spinnable’ mesophase pitch is presented. With this mesophase pitch, reproducible and stable multifilament melt spinning was developed, producing green fiber tows with filament diameters of approximately 20 m. The multifilament melt spinning was the most challenging aspect of the project and required the most effort. The best practices learned from this project for melt spinning are reviewed in this report. Once spun, the green fibers were oxidatively stabilized, carbonized and graphitized under inert gas atmosphere to form the final carbon fibers. Given the relatively high softening point of the mesophase pitch, no issues of interfilament fusion were observed during batch oxidation, and subsequent batch carbonization and graphitization went smoothly in all cases. An ~ 80 wt.% conversion of the green fiber mass to final carbon fiber was achieved. After graphitization, the carbon fibers showed high tensile moduli (most were ~ 600 GPa, or 87 Msi) consistent with commercially-available high-performance pitch-based carbon fiber. However, tensile strength and strain to failure were comparatively low. Further work to reduce defects in and on the fiber surfaces would increase these properties. SEM imaging of the graphitic fiber textures is presented herein. Rudimentary composites were fabricated from the carbon fibers and characterized showing similar modulus to baseline composites fabricated with commercial carbon fiber. Finally, a basic economic analysis was done showing the potential to increase the value of the isotropic coal tar pitch by up to 13.6 to 136 times based on a carbon fiber value of $\$$5/lb to $\$$50/lb, respectively. Moreover, the site case study suggests that the coal tar from the single integrated steel mill could supply production of up to 16 kt/yr of carbon fiber. Finally, a technological gap analysis was done which shed light on remaining technical challenges. These challenges included: recovery and utilization of condensates from the mesophase pitch processing, further advancing and increasing stability of the multifilament melt spinning processing, optimization of the oxidation processing, defect reduction for increased carbon fiber strength, and the development of a weaving process towards carbon fiber fabrics. To maximize the coal value chain, the primary objectives of this project were to (a) develop and scale efficient processing technology for producing melt-spinnable mesophase pitch from isotropic coal tar pitch, (b) clarify and simplify tedious continuous fiber processing technologies (particularly multifilament melt spinning of mesophase pitch) towards the efficient production of high performance carbon fiber, and (c) demonstrate and characterize representative composite parts derived from the final carbon fiber. Immense progress was made on all 3 objectives and is detailed in this report.

01 COAL, LIGNITE, AND PEAT↗

Analysis of the EM scattering from arbitrary open-ended waveguide cavities using axial Gaussian Beam tracking

The electromagnetic (EM) scattering from a planar termination located inside relatively arbitrarily shaped open-ended waveguide cavities with smoothly curved interior walls is analyzed using a Gaussian Beam (GB) expansion of the incident plane wave fields in the open end. The cavities under consideration may contain perfectly-conducting interior walls with or without a thin layer of material coating, or the walls may be characterized by an impedance boundary condition. In the present approach, the GB's are tracked only to the termination of the waveguide cavity via beam reflections from interior waveguide cavity walls. The Gaussian beams are tracked approximately only along their beam axes; this approximation which remains valid for relatively well focussed beams assumes that an incident GB gives rise to a reflected GB with parameters related to the incident beam and the radius of curvature of the wall. It is found that this approximation breaks down for GB's which come close to grazing a convex surface and when the width of the incident beam is comparable to the radius of curvature of the surface. The expansion of the fields at the open end depend on the incidence angle only through the expansion coefficients, so the GB's need to be tracked through the waveguide cavity only once for a wide range of incidence angles. At the termination, the sum of all the GB's are integrated using a result developed from a generalized reciprocity principle, to give the fields scattered from the interior of the cavity. The rim edge at the open end of the cavity is assumed to be sharp and the external scattering from the rim is added separately using Geometrical Theory of Diffraction. The results based on the present approach are compared with solutions based on the hybrid asymptotic modal method. The agreement is found to be very good for cavities made up of planar surfaces, and also for cavities with curved surfaces which are not too long with respect to their width.

Burkholder, R. J.↗

Cloud 3D Effects Evidenced in Landsat Power Spectra and Autocorrelation Functions

An analysis of nadir reflectivity Fourier spatial power spectra and autocorrelation functions at solar wavelengths and for cloudy conditions has been carried out. The data come from Landsat Thematic Mapper (TM) observations, while Monte Carlo (MC) simulations are used to aid the interpretation of the Landsat results. We show that radiative processes produce consistent signatures on power spectra and autocorrelation functions. The former take a variety of forms not shown or explained in previous observational studies. We demonstrate that the TM spectra can potentially be affected by both radiative "roughening" at intermediate scales (approx. 1 -5 km), being more prevalent at large solar zenith angles, and the already documented radiative "smoothing" at small scales (less than 1 km). These processes are wavelength dependent, as shown by systematic differences between conservative (for cloud droplets) TM band 4 (approx. 0.8 microns) and absorbing band 7 (approx. 2.2 microns): band 7 exhibits more roughening and less smoothing. This is confirmed quantitatively by comparing least-squared fitted power spectral slopes for the two bands. It is also corroborated by a slower decrease with distance of autocorrelation function values for band 4 compared to band 7. The appearance of roughening at large solar zenith angles is a result of side illumination and shadowing and adds an additional complexity to the power spectra. MC spectra are useful in illustrating that scale invariant optical depth fields can produce complex power spectra that take a variety of shapes under different conditions. We show that radiative roughening increases with the decrease of single scattering albedo and with the increase of solar zenith angle (as in the observations). For high Sun there is also a clear shift of the radiative smoothing scale to smaller values as droplet absorption increases. The shape of the power spectrum is sensitive to the magnitude and type of cloud top height variability, with the spectral signatures of decorrelation between reflectance and optical depth at large scales becoming stronger as the magnitude of cloud top variations increase. Finally, the usefulness of power spectral analysis in evaluating the skill of novel optical depth retrieval techniques in removing 3D radiative effects is demonstrated. New techniques using inverse Non-local Independent Pixel Approximation (NIPA) and Normalized Difference of Nadir Reflectivity (NDNR) yield optical depth fields which better match the scale-by-scale variability of the true optical depth field.

Oreopoulos, Lazaros↗

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗

Vestibular and Non-vestibular Contributions to Eye Movements that Compensate for Head Rotations during Viewing of Near Targets

We studied horizontal eye movements induced by en-bloc yaw rotation, over a frequency range 0.2 - 2.8 Hz, in 10 normal human subjects as th ey monocularly viewed a target located at their near point of focus. We measured gain and phase relationships between eye-in-head velocity and head velocity when the near target was either earth-fixed or head-fixed. During viewing of the earth-fixed near target,median gain was 1.49 (range 1.24 - 1.87) at 0.2 Hz for the group of subjects, but decl ined at higher frequencies, so that at 2.8 Hz median gain was 1.08 (r ange 0.68 - 1.67). During viewing of the head-fixed near target, median gain was 0.03 (range 0.01 - 0.10) at 0.2 Hz for the group of subjec ts, but increased at higher frequencies, so that at 2.8 Hz median gai n was 0.71 (range 0.28 - 0.94). We estimated the vestibular contribution to these responses (vestibulo-ocular reflex gain, Gvor) by applyin g transient head perturbations (peak acceleration> 1,000 deg's(exp 2) ) during sinusoidal rotation under the two viewing conditions. Median Gvor, estimated < 70m after the onset of head perturbation, was 0.98 (range 0.39 - 1.42) while viewing the earth-fixed near target, and 0. 97 (range 0.37 - 1.33) while viewing the head-fixed near target. For the group of subjects, 9 out of 10 subjects showed no sigificant diff erence of Gvor between the two viewing conditions ( p > 0.053 ) at all test frequencies. Since Gvor accounted for only approximately 73% of the overall response gain during viewing of the earth-fixed target, we investigated the relative contributions of non-vestibular factors. When subjects viewed the earth-fixed target under strobe illumination , to eliminate retinal image slip information, the gain of compensato ry eye movements declined compared with viewing in ambient room light . During sum-of-sine head rotations, while viewing the earth-fixed target, to minimize contributions from predictive mechanisms, gain also declined Nonetheless, simple superposition of smooth-pursuit tracking of sinusoidal target motion could not fully account for the overall r esponse at higher frequencies, suggesting other non-vestibular contributions. During binocular viewing conditions when vergence angle was s ignificantly greater than monocular viewing (p < 0.001), this gain of compensatory eye movements did not show proportional change; indeed, gain could not be correlated with vergence angle during monocular or binocular viewing. We conclude that several separate factors contribute to generate eye rotations during sinusoidal yaw head rotations whi le viewing a near target; these include the VOR, visual-tracking eye movements that utilize retinal image motion, predictive eye movements and, possibly, other unidentified nonvestibular factors. For these experiments, vergence was not an important determinant of response gain .

Han, Yanning H.↗

Projections of two-particle correlations onto transverse rapidity in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV at STAR

Two-particle correlation measurements projected onto two-dimensional, transverse rapidity coordinates (y T1 ,y T2 ), allow access to dynamical properties of the QCD medium produced in relativistic heavy-ion collisions that angular correlation measurements are not sensitive to. Here, we report non-identified charged-particle correlations for Au + Au minimum-bias collisions at $\sqrt{s_{NN}}$ = 200 GeV taken by the STAR experiment at the Relativistic Heavy-Ion Collider (RHIC). Correlations are presented as 2D functions of transverse rapidity for like-sign, unlike-sign and all charged-particle pairs, as well as for particle pairs whose relative azimuthal angles lie on the near-side, the away-side, or at all relative azimuth. The correlations are constructed using charged particles with transverse momentum p T ≥ 0.15 GeV/c, pseudorapidity from –1 to 1, and azimuthal angles from –π to π. The significant correlation structures that are observed evolve smoothly with collision centrality. The major correlation features include a saddle shape plus a broad peak with maximum near y T ≈ 3, corresponding to p T ≈ 1.5 GeV/c. The broad peak is observed in both like- and unlike-sign charge combinations and in near- and away-side relative azimuthal angles. The all-charge, all-azimuth correlation measurements are compared with the theoretical predictions of HIJING and EPOS. The results indicate that the correlations for peripheral to mid-central collisions can be approximately described as a superposition of nucleon + nucleon collisions with minimal effects from the QCD medium. Strong medium effects are indicated in mid- to most-central collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coordinate-Based Seismic Interpolation in Irregular Land Survey: A Deep Internal Learning Approach

Physical and budget constraints often result in irregular sampling, which complicates accurate subsurface imaging. Preprocessing approaches, such as missing trace or shot interpolation, are typically employed to enhance seismic data in such cases. Recently, deep learning has been used to address the trace interpolation problem at the expense of large amounts of training data to adequately represent typical seismic events. Nonetheless, most research in this area has focused on trace reconstruction, with little attention having been devoted to shot interpolation. Furthermore, existing methods assume regularly spaced receivers/sources failing in approximating seismic data from real (irregular) surveys. This work presents a novel shot gather interpolation approach which uses a continuous coordinate-based representation of the acquired seismic wavefield parameterized by a neural network. The proposed unsupervised approach, which we call coordinate-based seismic interpolation (CoBSI), enables the prediction of specific seismic characteristics in irregular land surveys without using external data during neural network training. Importantly, experimental results on real and synthetic 3-D data validate the ability of the proposed method to estimate continuous smooth seismic events in the time-space and frequency-wavenumber domains, improving sparsity or low-rank-based interpolation methods.

58 GEOSCIENCES↗

An embedded function approach for turbulent flow prediction

In conventional prediction methods for turbulent flow influence on surface properties, very small mesh sizes and considerable computational effort is required to adequately resolve the intense velocity and temperature profile variations that occur in the wall-layer region. In this study, an approach is described wherein an outer region numerical solution is smoothly matched to a set of embedded analytic profile functions in the near-wall region; these wall-layer profile functions having been derived through consideration of the coherent structure of the time-dependent near-wall flow provide analytic expressions for the mean velocity and enthalpy profiles in the wall layer. The concept is demonstrated through example applications to turbulent boundary-layer flows in two dimensions. The technique is shown to be very efficient and it is demonstrated that a reduction of approximately half the mesh points across the 2-D layer may be realized (as compared to conventional methods) with no degradation in accuracy.

Walker, J. D. A.↗

Inverse problems in diffraction

A two-dimensional problem of diffraction of a plane electromagnetic wave on a smooth 2 pi-periodic surface is considered. A numerical algorithm solving this problem is developed. An inverse problem of determination of the shape of 2 pi-periodic surface using the performance data of reverse scattering is considered. The inverse problem was solved by means of minimization of the residual functional with the help of the gradient descent method. The initial data were calculated with the help of the numerical method. On each step of the iterative method of minimization, the residual functional was calculated approximately with the help of the small slope method. The examples of the shape determination are considered.

Mikheev, Andrew G.↗

Hydrodynamic Impact of a System with a Single Elastic Mode II : Comparison of Experimental Force and Response with Theory

Hydrodynamic impact tests were made on an elastic model approximating a two-mass spring system to determine experimentally the effects of structural flexibility on the hydrodynamic loads encountered during seaplane landing impacts and to correlate the results with theory. A flexible seaplane was represented by a two-mass spring system consisting of a rigid prismatic float connected to a rigid upper mass by an elastic structure. The model had a ratio of sprung mass to hull mass of 0.6 and a natural frequency of 3.0 cycles per second. The tests were conducted in smooth water at fixed trims and included both high and low flight-path angles and a range of velocity. Theoretical and experimental comparisons indicated that the theoretical results agreed well with the experimental results.

Miller, Robert W↗

Comparison of large-scale structures and velocities in the local universe

Comparison of the large-scale density and velocity fields in the local universe shows detailed agreement, strengthening the standard paradigm of the gravitational origin of these structures. Quantitative analysis can determine the cosmological density parameter, Omega, and biasing factor, b; there is virtually no sensitivity in any local analyses to the cosmologial constant, lambda. Comparison of the dipole anisotropy of the cosmic microwave background with the acceleration due to the Infrared Astronomy Satellite (IRAS) galaxies puts the linear growth factor in the range beta approximately equals Omega (exp 0.6)/b = 0.6(+0.7/-0.3) (95% confidence). A direct comparison of the density and velocity fields of nearby galaxies gives beta = 1.3 (+0.7/-0.6), and from nonlinear analysis the weaker limit (Omega greater than 0.45 for b greater than 0.5 (again 95% confidence). A tighter limit (Omega greater than 0.3 (4-6 sigma)), is obtained by a reconstruction of the probability distribution function of the initial fluctuations from which the structures observed today arose. The last two methods depend critically on the smooth velocity field determined from the observed velocities of nearby galaxies by the POTENT method. A new analysis of these velocities, with more than three times the data used to obtain the above quoted results, is now underway and promises to tighten the uncertainties considerably, as well as reduce systematic bias.

Yahil, Amos↗

A Numerical and Experimental Study of Coflow Laminar Diffusion Flames: Effects of Gravity and Inlet Velocity

In this work, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting behavior of laminar coflow methane-air diffusion flames was investigated both computationally and experimentally. A series of flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) was assessed numerically under microgravity and normal gravity conditions with the fuel stream CH4 mole fraction ranging from 0.4 to 1.0. Computationally, the MC-Smooth vorticity-velocity formulation of the governing equations was employed to describe the reactive gaseous mixture; the soot evolution process was considered as a classical aerosol dynamics problem and was represented by the sectional aerosol equations. Since each flame is axisymmetric, a two-dimensional computational domain was employed, where the grid on the axisymmetric domain was a nonuniform tensor product mesh. The governing equations and boundary conditions were discretized on the mesh by a nine-point finite difference stencil, with the convective terms approximated by a monotonic upwind scheme and all other derivatives approximated by centered differences. The resulting set of fully coupled, strongly nonlinear equations was solved simultaneously using a damped, modified Newton's method and a nested Bi-CGSTAB linear algebra solver. Experimentally, the flame shape, size, lift-off height, and soot temperature were determined by flame emission images recorded by a digital camera, and the soot volume fraction was quantified through an absolute light calibration using a thermocouple. For a broad spectrum of flames in microgravity and normal gravity, the computed and measured flame quantities (e.g., temperature profile, flame shape, lift-off height, and soot volume fraction) were first compared to assess the accuracy of the numerical model. After its validity was established, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting tendency of laminar coflow methane-air diffusion flames was explored further by examining quantities derived from the computational results.

microgravity↗

Algebraic grid generation using tensor product B-splines

Finite difference methods are more successful if the accompanying grid has lines which are smooth and nearly orthogonal. The development of an algorithm which produces such a grid when given the boundary description. Topological considerations in structuring the grid generation mapping are discussed. The concept of the degree of a mapping and how it can be used to determine what requirements are necessary if a mapping is to produce a suitable grid is examined. The grid generation algorithm uses a mapping composed of bicubic B-splines. Boundary coefficients are chosen so that the splines produce Schoenberg's variation diminishing spline approximation to the boundary. Interior coefficients are initially chosen to give a variation diminishing approximation to the transfinite bilinear interpolant of the function mapping the boundary of the unit square onto the boundary grid. The practicality of optimizing the grid by minimizing a functional involving the Jacobian of the grid generation mapping at each interior grid point and the dot product of vectors tangent to the grid lines is investigated. Grids generated by using the algorithm are presented.

Saunders, B. V.↗

Spectroscopy of the Perseus Cluster

We present preliminary results of a XMM-Newton 50 ks observation of the Perseus Cluster that provides an unprecedented view of the central 0.5 Mpc region. The projected gas temperature declines smoothly by a factor of 2 from a maximum value of approx. 7 keV in the outer regions to just above 3 keV at the cluster center. Over this same range, the heavy-element abundance rises slowly from 0.4 to 0.5 solar as the radius decreases from 14 ft. to 5 ft., and then it rises to a peak of almost 0.7 solar at 1&farcm;25 before declining to 0.4 at the center. Th global east-west asymmetry of the gas temperature and surface brightness distributions, approximately aligned with the chain of bright galaxies, suggests an ongoing merger, although the modest degree of the observed asymmetry certainly excludes a major merger interpretation. The chain of galaxies probably traces the filament along which accretion started some time ago and is continuing at the present time. A cold and dense (low-entropy) cluster core like Perseus is probably well "protected" against the penetration of the gas of infalling groups and poor clusters, whereas in non-cooling core clusters such as Coma and A1367, infalling subclusters can penetrate deeply into the core region. In Perseus, gas associated with infalling groups may be stripped completely at the outskirts of the main cluster and only compression waves (shocks) may reach the central regions. We argue, and show supporting simulations, that the passage of such a wave(s) can qualitatively explain the overall horseshoe shaped appearance of the gas temperature map (the hot horseshoe surrounds the colder, low-entropy core) as well as other features of the Perseus Cluster core. These simulations also show that as compression waves traverse the cluster core, they can induce oscillatory motion of the cluster gas that can generate multiple sharp "edges" on opposite sides of the central galaxy. Gas motions induced by mergers may be a natural way to explain the high frequency of "edges" seen in clusters with cooling cores.

Jones, Christine↗

Compositional patterning in irradiated alloys: Effective potentials and effective interfacial energy

Compositional patterning (CP) in binary alloys during energetic particle irradiation is studied using a kinetic model that considers two competing kinetic processes, a thermally activated one promoting macroscopic phase separation (MPS) of the concentration field c(r,t) and a forced one resulting in finite-range random atomic mixing. The forced mixing is modeled by a Gaussian relocation distribution with a characteristic distance R. A series of approximate kinetic models are introduced by expanding the mixing function into a series of n terms, thus replacing the non-local evaluations of the concentration field c(r'-r,t) by local derivatives of c(r,t). This approach makes it possible to obtain exact effective potentials and build steady-state diagrams for each order-n model. Phase field (PF) simulations using these order-n models reveal that near the onset of patterning, phase evolution is accurately described using an order-3 model, which changes smoothly from an extended Cahn-Hilliard free energy in the MPS regime to a one-mode Swift-Hohenberg functional in the CP regime. Deeper into the patterning regime, higher-order models are required to achieve convergence, yielding square-like concentration profiles characteristic of a strong segregation regime. These higher-order effective free energies are analogous to multimodal Swift-Hohenberg functionals. Here, a new definition for the effective interfacial energy is proposed in the CP regime, since the interfacial area is no longer an excess quantity in that regime, precluding the use of the standard thermodynamic definition of interfacial energy.

36 MATERIALS SCIENCE↗

Vortex equations: Singularities, numerical solution, and axisymmetric vortex breakdown

A method of weighted residuals for the computation of rotationally symmetric quasi-cylindrical viscous incompressible vortex flow is presented and used to compute a wide variety of vortex flows. The method approximates the axial velocity and circulation profiles by series of exponentials having (N + 1) and N free parameters, respectively. Formal integration results in a set of (2N + 1) ordinary differential equations for the free parameters. The governing equations are shown to have an infinite number of discrete singularities corresponding to critical values of the swirl parameters. The computations point to the controlling influence of the inner core flow on vortex behavior. They also confirm the existence of two particular critical swirl parameter values: one separates vortex flow which decays smoothly from vortex flow which eventually breaks down, and the second is the first singularity of the quasi-cylindrical system, at which point physical vortex breakdown is thought to occur.

Bossel, H. H.↗

Magnetic field studies by Voyager 2 - Preliminary results at Saturn

Results of Voyager 2 studies of the magnetosphere and planetary magnetic field of Saturn are presented. Magnetometer studies have confirmed the results obtained by Voyager 1, indicating the magnetic field to be that of a centered dipole of moment 0.21 gauss Saturn radii-cubed, tilted approximately 0.8 deg from the rotation axis and a maximum measured field intensity of 1187 nT at latitude 17.3 deg N just before periapsis. Voyager 2 observed multiple bow shock and magnetopause crossings during its inbound and outbound trajectories, which were complementary to those of Voyager 1, including magnetopause crossing at 18.5 Saturn radii on the inbound trajectory, and at 48.4-50.9 Saturn radii outbound indicative of magnetospheric expansion due to changing solar wind conditions. Throughout the outbound passage, the magnetospheric field was observed to be relatively steady and smooth, with no evidence for any azimuthal asymmetry or magnetic anomaly. Results thus are incapable of accounting for the observed periodic modulation of the Saturnian kilometric radio emissions.

Ness, N. F.↗