Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dissipative particle 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 199 records · Page 11

Emerging supersolidity in photonic-crystal polariton condensates

A supersolid is a counter-intuitive phase of matter in which its constituent particles are arranged into a crystalline structure, yet they are free to flow without friction. This requires the particles to share a global macroscopic phase while being able to reduce their total energy by spontaneous, spatial self-organization. The existence of the supersolid phase of matter was speculated more than 50 years ago1, 2, 3–4. However, only recently has there been convincing experimental evidence, mainly using ultracold atomic Bose–Einstein condensates (BECs) coupled to electromagnetic fields. There, various guises of the supersolid were created using atoms coupled to high-finesse cavities5,6, with large magnetic dipole moments7, 8, 9, 10, 11, 12–13, and spin–orbit-coupled, two-component systems showing stripe phases14, 15–16. Here we provide experimental evidence of a new implementation of the supersolid phase in a driven-dissipative, non-equilibrium context based on exciton–polaritons condensed in a topologically non-trivial, bound state in the continuum (BiC) with exceptionally low losses, realized in a photonic-crystal waveguide. We measure the density modulation of the polaritonic state indicating the breaking of translational symmetry with a precision of several parts in a thousand. Direct access to the phase of the wavefunction allows us to also measure the local coherence of the supersolid. We demonstrate the potential of our synthetic photonic material to host phonon dynamics and a multimode excitation spectrum.

Trypogeorgos, Dimitrios↗

The Supercritical Pile GRB Model: The Prompt to Afterglow Evolution

The 'Supercritical Pile' is a very economical gamma ray burst (GRB) model that provides for the efficient conversion of the energy stored in the protons of a Relativistic Blast Wave (RBW) into radiation and at the same time produces - in the prompt GRB phase, even in the absence of any particle acceleration - a spectral peak at an energy sim 1 MeV. We extend this model to include also the evolution of the RBW Lorentz factor Gamma and thus follow the spectral and temporal features of this model into the early GRB afterglow stage. One of the novel features of the present treatment is the inclusion of the feedback of the GRB produced radiation on the evolution of Gamma with radius. This feedback and the presence of kinematic and dynamic thresholds in the model can be the sources of rich time evolution which we have begun to explore. In particular, one can this way obtain afterglow light curves with steep decays followed by the more conventional flatter afterglow slopes, while at the same time preserving the desirable features of the model, i.e. the well defined relativistic electron source and radiative processes that produce the proper peak in the nu F spectra. Furthermore, the existence of a kinematic threshold in this model provides for a operational distinction of the prompt and afterglow GRB stages; in fact, the afterglow stage sets in when the RBW Lorentz factor cannot anymore fulfill the kinematic condition for pair formation in the photon - proton pair production reactions that constitute the fundamental process for the dissipation of the blast wave kinetic energy. We present the results of a specific set of parameters of this model with emphasis on the multiwavelength prompt emission and transition to the early afterglow.

Kazanas, Demosthenes↗

High performance power handling in the absence of an H-mode edge in negative triangularity DIII-D plasmas

Abstract Experiments performed during strongly-shaped high-power diverted negative triangularity (NT) experiments in DIII-D achieved detached divertor conditions and a transient-free edge, showcasing the potential for application of NT to a core-edge integrated reactor-like scenario and providing the first characterization of the parametric dependencies for detachment onset. Detached divertor conditions will be required in future devices to mitigate divertor heat fluxes. Access to dissipative divertor conditions was investigated via an increase in upstream density. Detachment onset at the outer strike point was achieved with H-mode level confinement H 98 − y 2 ∼ 1 and reactor-relevant normalized pressures β N ∼ 2 . Confinement degradation was observed with deeper detachment, associated with the loss of an electron temperature pedestal. Differences in geometry, radial transport, impact of cross field drifts are discussed to explain differences in access to detachment in NT discharges. Higher normalized densities, with respect to equivalent discharges in positive triangularity, were necessary to achieve detachment, partially explained by the shorter parallel connection length to the targets. The effect of cross-field particle drifts ( E × B , B × ∇ B ) on access to detachment was demonstrated by the lower upstream density needed to access detachment with ion B × ∇ B drift directed outside of the active divertor (Greenwald fraction f Gw ∼ 0.9–1.0 vs f Gw ∼ 1.3). The upstream density at detachment onset was observed to increase linearly with plasma current with ion B × ∇ B drift into the divertor, consistent with the observed narrowing of the scrape-off layer heat flux width λ q . Edge fluid simulations capture separatrix densities needed to achieve detachment in NT plasma and their dependence on drift direction. The ability to reproduce detachment dynamics in NT plasma increases the confidence in future design studies for NT divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bifurcation and Hysteresis of the Magnetospheric Structure with a varying Southward IMF: Field Topology and Global Three-dimensional Full Particle Simulations

Using a three-dimensional full electromagnetic particle model (EMPM), we have performed global simulations of the interaction between the solar wind and the terrestrial magnetosphere, and have investigated its asymptotic stability. The distance between the dayside magnetopause subsolar point and the Earth center, R(sub mp) is measured, as the intensity of southward IMF |B(sub z)| is slowly varying. Based on the field topology theory, one analyzes the variation of R(sub mp) as a reference index of the dynamics of this interaction, when IMF |B(sub z)| successively increases and decreases to its original value. Two striking results are observed. First, as the IMF |B(sub z)| increases above a critical value, the variation of R(sub mp) suddenly changes (so called 'bifurcation' process in field topology). Above this critical value, the overall magnetic field topology changes drastically and is identified as being the signature of magnetic reconnection at the subsolar point on the magnetopause. Second, this subsolar point recovers its original location R(sub mp) by following different paths as the IMF |B(sub z)| value increases (from zero to a maximum fixed value) and decreases (from this maximum to zero) passing through some critical values. These different paths are the signature of 'hysteresis' effect, and are characteristic of the so-called 'subcritical-type' bifurcation. This hysteresis signature indicates that dissipation processes take place via an energy transfer from the solar wind to the magnetosphere by some irreversible way, which leads to a drastic change in the magnetospheric field topology. This hysteresis is interpreted herein as a consequence of the magnetic reconnection taking place at the dayside magnetopause. The field topology reveals to be a very powerful tool to analyze the signatures of three-dimensional magnetic reconnection without the obligation for determining the mechanisms responsible for, and the consequences of the reconnection on the overall magnetospheric dynamics.

Cai, DongSheng↗

MARCI-observed clouds in the Hellas Basin during northern hemisphere summer on Mars: Interpretation with the NASA/Ames Legacy Mars Global Climate Model

We present a study that is motivated by a population of water ice clouds in the Hellas Basin that has been observed by the MARs Color Imager (MARCI) instrument on the Mars Reconnaissance Orbiter (MRO) to persist throughout the majority of Northern Hemisphere (NH) summer. Although water ice clouds are present in Hellas at very low opacities throughout NH spring, they noticeably thicken after L(sub s) 60° and continue to do so until their peak optical thickness is attained at L(sub s) ~120°. They dissipate rapidly from L(sub s) 120° to 150°, and by L(sub s) 150°, the only clouds in Hellas are on the southern side of the basin and are indistinguishable from the polar hood clouds. We use the NASA/Ames Legacy Mars Global Climate Model (GCM), which is supported by the Agency’s Mars Climate Modeling Center, to investigate the dynamical and microphysical mechanisms that control the formation and evolution of Hellas water ice clouds. We show that water is transported from the North Polar Residual Cap (NPRC) southward across the equator and into the Hellas region. Water vapor is confined down low by cloud formation in the aphelion cloud belt, thus limiting the effectiveness of transport by the zonal mean overturning circulation (i.e., the Hadley cell). Thus, contrary to the commonly held conceptual understanding that the Hadley cell controls cross-equatorial water transport during this season, we show that the southward transport of water is done primarily in the vapor phase by stationary eddies at nearly all latitudes, including across the equator. Clouds form near the surface in the basin as moist air mixes with cold polar air on the western side of Hellas and are subsequently transported clockwise around the north side of the basin by the low-level cyclonic circulation. The simulated Hellas clouds have small particle sizes (~3–5 μm), are very low, and exist as long as the NPRC is exposed.

M. A. Kahre↗

Onset of Reconnection in the near Magnetotail: PIC Simulations

Using 2.5-dimensional particle-in-cell (PIC) simulations of magnetotail dynamics, we investigate the onset of reconnection in two-dimensional tail configurations with finite Bz. Reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. We found a clear distinction between stable and unstable cases, dependent on deformation amplitude and ion/electron mass ratio. The threshold appears consistent with electron tearing. The evolution prior to onset, as well as the evolution of stable cases, are largely independent of the mass ratio, governed by integral flux tube entropy conservation as imposed in MHD (magnetohydrodynamics). This suggests that ballooning instability in the tail should not be expected prior to the onset of tearing and reconnection. The onset time and other onset properties depend on the mass ratio, consistent with expectations for electron tearing. At onset,we found electron anisotropies T⊥⁄ T∥ (bottom tail divided by parallel tail) equals 1.1-1.3, raising growth rates and wavenumbers. Our simulations have provided a quantitative onset criterion that is easily evaluated in MHD simulations, provided the spatial resolution is sufficient. The evolution prior to onset and after the formation of a neutral line does not depend on the electron physics, which should permit an approximation by MHD simulations with appropriate dissipation terms.

PIC↗

Cavity-assisted magnetization switching in a quantum spin-phonon chain

Néel order switching in antiferromagnets has typically required intense optical driving, leading to substantial heating and limited efficiency. By placing antiferromagnets inside a terahertz-driven optical cavity, we propose a multi-particle mechanism for Néel order switching that benefits from reduced heating. Our analysis reveals that phonons are indispensable to this mechanism. A driven cavity mode couples to a spin-phonon chain, with all excitations dissipating energy through external baths. Mean-field analysis shows that cavity photons induce sublattice spin-density imbalance—an intrinsic symmetry-breaking effect absent without the cavity. Contrary to known (1–10 V/nm) laser fields required to switch the Néel order, our mechanism enables switching at remarkably low laser fields (1–5 V/μm), selectively targeting low-energy and perpendicular magnon modes. By virtue of the suppressed heating, the switching remains highly tunable through laser fluence, damping, and photon loss, establishing a low-dissipation route toward cavity-assisted opto-spintronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Forced Reconnection in the Near Magnetotail: Onset and Energy Conversion in PIC and MHD Simulations

Using two-dimensional particle-in-cell (PIC) together with magnetohydrodynamic (MHD) Q1 simulations of magnetotail dynamics, we investigate the evolution toward onset of reconnection and the subsequent energy transfer and conversion. In either case, reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. The boundary deformation leads to the formation of thin embedded current sheets, which are bifurcated in the near tail, converging to a single sheet farther out in the MHD simulations. The thin current sheets in the PIC simulation are carried by electrons and are associated with a strong perpendicular electrostatic field, which may provide a connection to parallel potentials and auroral arcs and an ionospheric signal even prior to the onset of reconnection. The PIC simulation very well satisfies integral entropy conservation (intrinsic to ideal MHD) during this phase, supporting ideal ballooning stability. Eventually, the current intensification leads to the onset of reconnection, the formation and ejection of a plasmoid, and a collapse of the inner tail. The earthward flow shows the characteristics of a dipolarization front: enhancement of Bz, associated with a thin vertical electron current sheet in the PIC simulation. Both MHD and PIC simulations show a dominance of energy conversion from incoming Poynting flux to outgoing enthalpy flux, resulting in heating of the inner tail. Localized Joule dissipation plays only a minor role.

two-dimensional↗

Fluctuation dynamo in a weakly collisional plasma

The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle–particle collisions allows departures from local thermodynamic equilibrium. These departures – pressure anisotropies – exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present an extensive numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart, particularly when the magnetic field is dynamically weak. If instead the parallel viscous stress is left unabated – a situation relevant to recent kinetic simulations of the fluctuation dynamo and, we argue, to the early stages of the dynamo in a magnetized ICM – the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics reminiscent of the saturated state of the large-Prandtl-number ( ${Pm}\gtrsim {1}$ ) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches simulated dynamo growth rates and magnetic-energy spectra. A prediction of this model, confirmed by our numerical simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. Finally, in all cases that do exhibit a viable dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗

Simultaneous velocity and density measurements of fully developed Rayleigh-Taylor mixing

Here, the dynamics of molecular mixing and the energy transfer process in the Rayleigh-Taylor instability (RTI) are studied through the collection of simultaneous density-velocity measurements. These experiments provide simultaneous density-velocity field measurements, in contrast to previous point measurements. Statistically stationary experiments are performed in a “convective-type” gas tunnel facility, with density contrast achieved through the injection of helium into the bottom stream. Three experiments at Atwood number ≈ 0.1 are captured at three outer scale Reynolds numbers Re = 520, 2260, and 4050. Particle image velocimetry and laser induced fluorescence are employed simultaneously. Statistics of the density and velocity show self-similar collapse of RTI profiles at large Reynolds number Re >2000. Flat velocity profiles indicate homogeneous turbulence characteristics in the core of the mixing region. Significant anisotropy develops in the flow, with horizontal velocity fluctuations being only 60% of the vertical velocity fluctuations. The turbulent mass flux is found to be asymmetric about the centerline, with increased peak towards the spike. Measurements of the molecular mixing show that mixing is maximized at the core of the flow and increases with increased Reynolds number. The transport equation of density-specific-volume correlation b shows that it is mostly produced in the core of the mixing region, and that the spatial evolution of its profile is the result of transport by bulk motion of the bubble and spike. Energy transfer from gravitational potential energy to turbulent kinetic energy and viscous dissipation is observed to occur in the experiment with a ratio of dissipated energy to potential energy released of 38%. The analysis of the turbulent kinetic energy transport equation budget reveals that production is the primary mechanism towards the growth of turbulent kinetic energy in the core of the flow, and is asymmetrically slightly skewed towards the spike. However, it is through the transport that the strong advection at the edges of the mixing region is maintained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamics of multilayer Rayleigh–Taylor instability at moderately high Atwood numbers

This paper investigates the multilayer Rayleigh–Taylor instability (RTI) using statistically stationary experiments conducted in a gas tunnel. Employing diagnostics such as particle image velocimetry (PIV) and planar laser induced fluorescence (PLIF), we make simultaneous velocity–density measurements to study how dynamics and mixing are linked in this variable density flow. Experiments are conducted in a newly built, blow-down three-layer gas tunnel facility. Mixing between three gas streams is studied, where the top and bottom streams are comprised of air, and the middle stream is an air–helium mixture. Shear is minimized between these streams by matching their inlet velocities. The four experimental conditions investigated here consist of two different density ratios (Atwood numbers 0.3 and 0.6), each investigated at two instability development times (or equivalently, two streamwise locations), and all experiments are with the same middle stream thickness of 3 cm. The growth of the middle layer is measured using laser-based planar Mie scattering visualization. The mixing width is found to grow linearly with time at late times. Various quantitative measures of molecular mixing indicate a very high degree of molecular mixing at late times in the multilayer RTI flow. The vertical turbulent mass flux $a_y$ is calculated. In addition to mostly negative values of $a_y$ , typical of buoyancy-dominated flows due to negative correlation between velocity and density fluctuations, positive regions are also observed in profiles of $a_y$ due to entrainment and erosion at the lower edge of the mixing region. Global energy budgets are calculated for the multilayer RTI flow at late times and it is found that the majority of potential energy released has been dissipated due to viscous effects, and a large value of mixing efficiency ( $\sim$ 60 %) is observed.

Mechanics↗

Computational fluid dynamic modeling to determine the indoor environment of an electron-ion collider service building

The design for the Electron-Ion Collider (EIC) calls for several service buildings that house various power supplies and control electronics for the collider ring itself. In order to operate within specified conditions, the ambient air entering the power supplies needs to be within a certain temperature range while dissipating the heat from losses. Proper cooling is therefore a necessity in the service buildings to ensure that every aspect of the EIC works as intended. Since the EIC is in the design stage, we are evaluating the indoor environment of the service building using the current design specifications. We have researched multiple forms of literature and performed the necessary calculations to compile a list of boundary conditions that accurately represent the situation at hand. We are using computational fluid dynamics modeling to solve the conservation equations for mass, momentum, and energy (Navier-Stokes). This allows us to perform a finite element analysis which will give us the flow distribution in the room as well as temperature profiles throughout the building. We have obtained a simulation result giving us the temperature profiles for the building and it shows that the placement of the racks and supply vents are essential to obtaining an even temperature distribution. This model will provide a basis for design decisions which will affect the overall cooling of the service building without extending the schedule and avoiding a costly reworking of the cooling system.

43 PARTICLE ACCELERATORS↗

Transient Mixing Driven by Buoyancy Flows

Mixing driven by buoyancy-induced flows is of particular interest to microgravity processes, as the body force that governs the intensity of flow fields can be directly controlled. We consider a model experimental system to explore the dynamics of mixing which employs two miscible liquids inside a cavity separated initially by a divider. The two liquids are oriented vertically inside a rectangular cavity with constant width and height, and varying depths to span the range of a Hele-Shaw cell to a 3-[) configuration. The two miscible liquids can be sufficiently diluted and died, for example water and deuterium oxide, such that a distinct interface exists across the divider. The transient mixing characteristic of the two fluids is addressed by following the Lagrangian history of the interface for various aspect ratios in the z-plane (depth variation) as well as a range of pulling velocities of the divider. The mixing characteristics of the two fluids are quantified from measurement of the length stretch of the interface and its flow field using respectively image processing techniques and Particle Imaging Velocimetry. Scaling analysis shows that the length stretch depends on four governing parameters, namely the Grashof number (Gr), Schmidt number (Sc), aspect ratio (Ar), and Reynolds number (Re). Variation of the Schmidt number is taken into account through thermophysical property variation. Thus our problem reduces to a codimension three bifurcation in parametric space for Gr, Ar, and Re. Our experimental results show that for Gr on the order of 106 and a nominal cavity aspect ratio Ar = 0.2, the net effect of removal of the divider and the overwhelming buoyancy force causes an overturning motion which stretches and folds the interface to produce ml internal breakwave. The structure of the breakwave is similar to the ubiquitous Rayleigh-Taylor instability morphology. The breakwave is dissipated either through internal or wall collision depending on the impulsive velocity of the divider as prescribed by the Reynolds number. The decay of the collision event occurs through sloshing oscillations over a short time scale. The two fluids then become stably stratified with a diffusive band at the interface indicating local mass transport. The local bifurcation of the internal breakwave is investigated as a function of aspect ratio. Results show that for narrow cavities on the order of 2mm (Ar = 0.04) folding does not occur, the interface only stretches. As the cavity size increases folding occurs through a supercritical bifurcation. Insight into the mechanism of folding is obtained from measurement of the flow field which shows that in the neighborhood of the folding event, there exists hyperbolic points caused by multiple vortex interactions. The global length stretch of the interface as a function of time is nearly Gaussian; calculations of finite-time Liapunov exponents as well as construction of horseshoe maps indicate the likelihood of a chaotic transient.

W.M.B. Duval↗

Two-fluid simulations of galaxy formation

We investigate the formation of galaxies and larger structure with a simulation modeling two gravitationally coupled fluids representing dark matter and baryons. The baryon gas dynamics are calculated with a smoothed particle hydrodynamics (SPH) method, and the physics modeled includes thermal pressure, shock heating, and radiative cooling. We simulate a 16 Mpc periodic cube with 64(exp 3) particles in each fluid and 10% baryon mass fraction. We confirm, for the first time experimentally, disk formation as a natural consequence of hierarchical clustering in a large-scale cosmological environment. The majority of isolated galaxies exhibit centrifugally supported disks. A power-law relation between cold baryonic mass and maximum rotation velocity is found, M varies as nu(sub rot)(exp alpha) with alpha = 2.5 after correcting for differential numerical resolution. Both the spatial and velocity distributions of the simulated galaxies are biased with respect to the dark matter. A counts-in-cells analysis indicates that an unphysical degree of merging in the central cluster is likely responsible for the antibias signal in the correlation function. A robust, scale-dependent velocity bias is measured. The ratio of galaxy to dark matter pairwise velocity dispersions on a scale of 1 Mpc is 0.7. The amplitude is only mildly dependent on redshift or mass cutoff and scales with separation as r(exp 0.2). The degree to which these results depend on numerical parameters is discussed. Mass resolution plays a key role in controlling the resulting fraction of cold, dense baryons. The mass fraction associated with galaxies decreases by a factor of approximately greater than 3 when the mass per particle is increased by a factor 8. Photoionization and energy input from supernova will have to be included to determine more carefully the fraction of highly dissipated material and the characteristics of the stellar component of galaxies.

Evrard, August E.↗

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗

Comparison of Classical and Charge Storage Methods for Determining Conductivity of Thin Film Insulators

Conductivity of insulating materials is a key parameter to determine how accumulated charge will distribute across the spacecraft and how rapidly charge imbalance will dissipate. Classical ASTM and IEC methods to measure thin film insulator conductivity apply a constant voltage to two electrodes around the sample and measure the resulting current for tens of minutes. However, conductivity is more appropriately measured for spacecraft charging applications as the "decay" of charge deposited on the surface of an insulator. Charge decay methods expose one side of the insulator in vacuum to sequences of charged particles, light, and plasma, with a metal electrode attached to the other side of the insulator. Data are obtained by capacitive coupling to measure both the resulting voltage on the open surface and emission of electrons from the exposed surface, as well monitoring currents to the electrode. Instrumentation for both classical and charge storage decay methods has been developed and tested at Jet Propulsion Laboratory (JPL) and at Utah State University (USU). Details of the apparatus, test methods and data analysis are given here. The JPL charge storage decay chamber is a first-generation instrument, designed to make detailed measurements on only three to five samples at a time. Because samples must typically be tested for over a month, a second-generation high sample throughput charge storage decay chamber was developed at USU with the capability of testing up to 32 samples simultaneously. Details are provided about the instrumentation to measure surface charge and current; for charge deposition apparatus and control; the sample holders to properly isolate the mounted samples; the sample carousel to rotate samples into place; the control of the sample environment including sample vacuum, ambient gas, and sample temperature; and the computer control and data acquisition systems. Measurements are compared here for a number of thin film insulators using both methods at both facilities. We have found that conductivity determined from charge storage decay methods is 102 to 104 larger than values obtained from classical methods. Another Spacecraft Charging Conference presentation describes more extensive measurements made with these apparatus. This work is supported through funding from the NASA Space Environments and Effects Program and the USU Space Dynamics Laboratory Enabling Technologies Program.

Swaminathan, Prasanna↗

Nonlinear dynamics and dissipation of a curvilinear vortex driven by a strong time-dependent Meissner current

Here, we report numerical simulations of large-amplitude oscillations of a trapped vortex line under a strong ac magnetic field $\textit{H(t) = H}$ sin$\textit{ωt}$ parallel to the surface. The power dissipated by an oscillating vortex segment driven by the surface ac Meissner currents was calculated by taking into account the nonlinear vortex line tension, vortex mass, and a nonlinear Larkin-Ovchinnikov (LO) viscous drag coefficient $\textit{η(v)}$. We show that the LO decrease of $\textit{η(v)}$ with the vortex velocity $\textit{v}$ can radically change the field dependence of the surface resistance R i ( H ) caused by trapped vortices. At low frequencies, $R_i(H)$ exhibits a conventional increase with $\textit{H}$, but as $\textit{ω}$ increases, the surface resistance becomes a nonmonotonic function of $\textit{H}$ which decreases with $\textit{H}$ at higher fields. The effects of frequency, pin spacing, and the mean-free path $l_i$ on the field dependence of $R_i(H)$ were calculated. It is shown that, as the surface gets dirtier and $l_i$ decreases, the anomalous drop of $R_i(H)$ with $\textit{H}$ shifts to lower fields which can be much smaller than the lower critical magnetic field. Our numerical simulations also show that the LO decrease of $\textit{η(v)}$ with $\textit{v}$ can cause a vortex bending instability at high field amplitudes and frequencies, giving rise to the formation of dynamic kinks along the vortex. Measurements of $R_i(H)$ caused by sparse vortices trapped perpendicular to the surface can offer opportunities to investigate an extreme nonlinear dynamics of vortices driven by strong current densities up to the depairing limit at low temperatures. The behavior of $R_i(H)$, which can be tuned by varying the rf frequency or concentration of nonmagnetic impurities, is not masked by strong heating effects characteristic of dc or pulse transport measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗