Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface damping”

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 577 records · Page 32

Study of Fluid Flow Control In Protein Crystallization Using Strong Magnetic Fields

An important component in biotechnology, particularly in the area of protein engineering and rational drug design is the knowledge of the precise three-dimensional molecular structure of proteins. The quality of structural information obtained from X-ray diffraction methods is directly dependent on the degree of perfection of the protein crystals. As a consequence, the growth of high quality macromolecular crystals for diffraction analyses has been the central focus for biochemists, biologists, and bioengineers. Macromolecular crystals are obtained from solutions that contain the crystallizing species in equilibrium with higher aggregates, ions, precipitants, other possible phases of the protein, foreign particles, the walls of the container, and a likely host of other impurities. By changing transport modes in general, i.e., reduction of convection and sedimentation, as is achieved in 'microgravity', researchers have been able to dramatically affect the movement and distribution of macromolecules in the fluid, and thus their transport, formation of crystal nuclei, and adsorption to the crystal surface. While a limited number of high quality crystals from space flights have been obtained, as the recent National Research Council (NRC) review of the NASA microgravity crystallization program pointed out, the scientific approach and research in crystallization of proteins has been mainly empirical yielding inconclusive results. We postulate that we can reduce convection in ground-based experiments and we can understand the different aspects of convection control through the use of strong magnetic fields and field gradients. Whether this limited convection in a magnetic field will provide the environment for the growth of high quality crystals is still a matter of conjecture that our research will address. The approach exploits the variation of fluid magnetic susceptibility with concentration for this purpose and the convective damping is realized by appropriately positioning the crystal growth cell so that the magnetic susceptibility force counteracts terrestrial gravity. The general objective is to test the hypothesis of convective control using a strong magnetic field and magnetic field gradient and to understand the nature of the various forces that come into play. Specifically we aim to delineate causative factors and to quantify them through experiments, analysis and numerical modeling. Once the basic understanding is obtained, the study will focus on testing the hypothesis on proteins of pyruvate dehydrogenase complex (PDC), proteins E1 and E3. Obtaining high crystal quality of these proteins is of great importance to structural biologists since their structures need to be determined.

Ramachandran, N.↗

Study of Fluid Flow Control in Protein Crystallization using Strong Magnetic Fields

An important component in biotechnology, particularly in the area of protein engineering and rational drug design is the knowledge of the precise three-dimensional molecular structure of proteins. The quality of structural information obtained from X-ray diffraction methods is directly dependent on the degree of perfection of the protein crystals. As a consequence, the growth of high quality macromolecular crystals for diffraction analyses has been the central focus for biochemists, biologists, and bioengineers. Macromolecular crystals are obtained from solutions that contain the crystallizing species in equilibrium with higher aggregates, ions, precipitants, other possible phases of the protein, foreign particles, the walls of the container, and a likely host of other impurities. By changing transport modes in general, i.e., reduction of convection and sedimentation, as is achieved in "microgravity", researchers have been able to dramatically affect the movement and distribution of macromolecules in the fluid, and thus their transport, formation of crystal nuclei, and adsorption to the crystal surface. While a limited number of high quality crystals from space flights have been obtained, as the recent National Research Council (NRC) review of the NASA microgravity crystallization program pointed out, the scientific approach and research in crystallization of proteins has been mainly empirical yielding inconclusive results. We postulate that we can reduce convection in ground-based experiments and we can understand the different aspects of convection control through the use of strong magnetic fields and field gradients. Whether this limited convection in a magnetic field will provide the environment for the growth of high quality crystals is still a matter of conjecture that our research will address. The approach exploits the variation of fluid magnetic susceptibility with concentration for this purpose and the convective damping is realized by appropriately positioning the crystal growth cell so that the magnetic susceptibility force counteracts terrestrial gravity. The general objective is to test the hypothesis of convective control using a strong magnetic field and magnetic field gradient and to understand the nature of the various forces that come into play. Specifically we aim to delineate causative factors and to quantify them through experiments, analysis and numerical modeling. Once the basic understanding is obtained, the study will focus on testing the hypothesis on proteins of pyruvate dehydrogenase complex (PDC), proteins E1 and E3. Obtaining high crystal quality of these proteins is of great importance to structural biologists since their structures need to be determined. Specific goals for the investigation are: 1. To develop an understanding of convection control in diamagnetic fluids with concentration gradients through experimentation and numerical modeling. Specifically solutal buoyancy driven convection due to crystal growth will be considered. 2. To develop predictive measures for successful crystallization in a magnetic field using analyses and numerical modeling for use in future protein crystal growth experiments. This will establish criteria that can be used to estimate the efficacy of magnetic field flow damping on crystallization of candidate proteins. 3. To demonstrate the understanding of convection damping by high magnetic fields to a class of proteins that is of interest and whose structure is as yet not determined. 4. To compare quantitatively, the quality of the grown crystals with and without a magnetic field. X-ray diffraction techniques will be used for the comparative studies. In a preliminary set of experiments, we studied crystal dissolution effects in a 5 Tesla magnet available at NASA Marshall Space Flight Center (MSFC). Using a Schlieren setup, a 1mm crystal of Alum (Aluminum-Potassium Sulfate) was introduced in a 75% saturated solution and the resulting dissolution plume was observed. The experiment was conducted both in the presence and absence of a magnetic field gradient. The magnet produces a gradient field of approx. 1 Tesla2/cm. Image analysis of the recorded images indicated an enhanced plume velocity that was of the order of the measurement limit. For this experiment, both the gradient and gravity fields are in the same direction resulting in an enhanced effective gravity that tends to accelerate the observed plume velocity. While the results are not conclusive, pending further tests, it clearly points out the inadequacy of the MSFC magnet for conducting protein crystallization experiments and the need for a stronger magnet. In spacebased experiments, however, where the gravitational effects are small, only a weak magnetic field will be required to control or mitigate the effects of convective contamination.

Ramachandran, Narayanan↗

Vibration mitigation in interlocking metasurfaces

Interlocking metasurfaces (ILMs) are a new class of mechanical metamaterials with patterned surface features that enable joining of two adjacent bodies, serving as an alternative to welds, fasteners, or adhesives. While the strength of ILMs has been explored previously, it is also possible to envision ILMs that possess the ability to dampen vibration or tailor energy transmission across the interface. In this study, we explore the vibration transmission under frequency sweep and random excitation of three ILM designs. This is the first study to characterize the response of ILMs under vibration. We find that ILMs which are tightly coupled (e.g., small interference fit) behave as a linear system and depart only marginally from the behavior of a solid body. However, through strategic manipulation of the ILM topology, and in particular the fit clearances, we show that vibrations can be damped by over an order of magnitude. Design strategies to employ ILMs for vibration mitigation are developed using a first-order model.

Bolmin, O. [Carnegie Mellon Univ., Pittsburgh, PA ↗

A vortex-lattice method for calculating lifting-surface interference

The interference effect of closely coupled lifting surfaces was investigated. A nonsteady vortex-lattice method was applied to calculate the steady-state and nonsteady lift characteristics of configurations having interference effects. The configurations consist of low-aspect-ratio delta wings, at low and high angles of attack. Leading-edge separation and wake roll-up are simulated by a time-dependent wake-shedding procedure. For steady-state flow conditions, this numerical procedure saves a considerable amount of computer time, compared with iterative methods, and yields the same results. A better understanding of the interference effect can be gained by using the method to study the transient behavior. The nonsteady approach also offers the capability of calculating various nonsteady motions, as is demonstrated in the calculation of the longitudinal damping.

Levin, D.↗

Dynamic Shape Reconstruction of Three-Dimensional Frame Structures Using the Inverse Finite Element Method

A robust and efficient computational method for reconstructing the three-dimensional displacement field of truss, beam, and frame structures, using measured surface-strain data, is presented. Known as shape sensing , this inverse problem has important implications for real-time actuation and control of smart structures, and for monitoring of structural integrity. The present formulation, based on the inverse Finite Element Method (iFEM), uses a least-squares variational principle involving strain measures of Timoshenko theory for stretching, torsion, bending, and transverse shear. Two inverse-frame finite elements are derived using interdependent interpolations whose interior degrees-of-freedom are condensed out at the element level. In addition, relationships between the order of kinematic-element interpolations and the number of required strain gauges are established. As an example problem, a thin-walled, circular cross-section cantilevered beam subjected to harmonic excitations in the presence of structural damping is modeled using iFEM; where, to simulate strain-gauge values and to provide reference displacements, a high-fidelity MSC/NASTRAN shell finite element model is used. Examples of low and high-frequency dynamic motion are analyzed and the solution accuracy examined with respect to various levels of discretization and the number of strain gauges.

Gherlone, Marco↗

Active Struts With Variable Spring Stiffness and Damping

An ultrasonic rock-abrasion tool (URAT) was developed using the same principle of ultrasonic/sonic actuation as that of the tools described in two prior NASA Tech Briefs articles: Ultrasonic/ Sonic Drill/Corers With Integrated Sensors (NPO-20856), Vol. 25, No. 1 (January 2001), page 38 and Ultrasonic/ Sonic Mechanisms for Drilling and Coring (NPO-30291), Vol. 27, No. 9 (September 2003), page 65. Hence, like those tools, the URAT offers the same advantages of low power demand, mechanical simplicity, compactness, and ability to function with very small axial loading (very small contact force between tool and rock). Like a tool described in the second of the cited previous articles, a URAT includes (1) a drive mechanism that comprises a piezoelectric ultrasonic actuator, an amplification horn, and a mass that is free to move axially over a limited range and (2) an abrasion tool bit. A URAT tool bit is a disk that has been machined or otherwise formed to have a large number of teeth and an overall shape chosen to impart the desired shape (which could be flat or curved) to the rock surface to be abraded. In operation, the disk and thus the teeth are vibrated in contact with the rock surface. The concentrated stresses at the tips of the impinging teeth repeatedly induce microfractures and thereby abrade the rock. The motion of the tool induces an ultrasonic transport effect that displaces the cuttings from the abraded area. The figure shows a prototype URAT. A piezoelectric-stack/horn actuator is housed in a cylindrical container. The movement of the actuator and bit with respect to the housing is aided by use of mechanical sliders. A set of springs accommodates the motion of the actuator and bit into or out of the housing through an axial range between 5 and 7 mm. The springs impose an approximately constant force of contact between the tool bit and the rock to be abraded. A dust shield surrounds the bit, serving as a barrier to reduce the migration of rock debris to sensitive instrumentation or mechanisms in the vicinity. A bushing at the tool-bit end of the housing reduces the flow of dust into the actuator and retains the bit when no axial load is applied.

Farley, Gary L.↗

Influence of thickness and camber on the aeroelastic stability of supersonic throughflow fans: An engineering approach

An engineering approach was used to include the nonlinear effects of thickness and camber in an analytical aeroelastic analysis of cascades in supersonic acial flow (supersonic leading-edge locus). A hybrid code using Lighthill's nonlinear piston theory and Lanes's linear potential theory was developed to include these nonlinear effects. Lighthill's theory was used to calculate the unsteady pressures on the noninterference surface regions of the airfoils in cascade. Lane's theory was used to calculate the unsteady pressures on the remaining interference surface regions. Two airfoil profiles was investigated (a supersonic throughflow fan design and a NACA 66-206 airfoil with a sharp leading edge). Results show that compared with predictions of Lane's potential theory for flat plates, the inclusion of thickness (with or without camber) may increase or decrease the aeroelastic stability, depending on the airfoil geometry and operating conditions. When thickness effects are included in the aeroelastic analysis, inclusion of camber will influence the predicted stability in proportion to the magnitude of the added camber. The critical interblade phase angle, depending on the airfoil profile and operating conditions, may also be influenced by thickness and camber. Compared with predictions of Lane's linear potential theory, the inclusion of thickness and camber decreased the aerodynamic stifness and increased the aerodynamic damping at Mach 2 and 2.95 for a cascade of supersonic throughflow fan airfoils oscillating 180 degrees out of phase at a reduced frequency of 0.1.

Ramsey, John K.↗

Influence of thickness and camber on the aeroelastic stability of supersonic throughflow fans

An engineering approach was used to include the nonlinear effects of thickness and camber in an analytical aeroelastic analysis of casecades in supersonic acial flow (supersonic leading-edge locus). A hybrid code using Lighthill's nonlinear piston theory and Lane's linear potential theory was developed to include these nonlinear effects. Lighthill's theory was used to calculate the unsteady pressures on the noninterference surface regions of the airfoils in cascade. Lane's theory was used to calculate the unsteady pressures on the remaining interference surface regions. Two airfoil profiles were investigated (a supersonic throughflow fan design and a NACA 66-206 airfoil with a sharp leading edge). Results show that compared with predictions of Lane's potential theory for flat plates, the inclusion of thickness (with or without camber) may increase or decrease the aeroelastic stability, depending on the airfoil geometry and operating conditions. When thickness effects are included in the aeroelastic analysis, inclusion of camber will influence the predicted stability in proportion to the magnitude of the added camber. The critical interblade phase angle, depending on the airfoil profile and operating conditions, may also be influenced by thickness and camber. Compared with predictions of Lane's linear potential theory, the inclusion of thickness and camber decreased the aerodynamic stiffness and increased the aerodynamic damping at Mach 2 and 2.95 for a cascade of supersonic throughflow fan airfoils oscillating 180 degrees out of phase at a reduced frequency of 0.1.

Ramsey, John K.↗

How well are hazards associated with derechos reproduced in regional climate simulations?

Abstract. A 15-member ensemble of convection-permitting regional simulations of the fast-moving and destructive derecho of 29–30 June 2012 that impacted the northeastern urban corridor of the USA is presented. This event generated 1100 reports of damaging winds, generated significant wind gusts over an extensive area of up to 500 000 km2, caused several fatalities, and resulted in widespread loss of electrical power. Extreme events such as this are increasingly being used within pseudo-global-warming experiments to examine the sensitivity of historical, societally important events to global climate non-stationarity and how they may evolve as a result of changing thermodynamic and dynamic contexts. As such it is important to examine the fidelity with which such events are described in hindcast experiments. The regional simulations presented herein are performed using the Weather Research and Forecasting (WRF) model. The resulting ensemble is used to explore simulation fidelity relative to observations for wind gust magnitudes, spatial scales of convection (as is manifest in high composite reflectivity, cREF), and both rainfall and hail production as a function of model configuration (microphysics parameterization, lateral boundary conditions (LBCs), start date, use of nudging, compiler choice, damping, and number of vertical levels). We also examine the degree to which each ensemble member differs with respect to key mesoscale drivers of convective systems (e.g., convective available potential energy and vertical wind shear) and critical manifestations of deep convection, e.g., vertical velocities, cold-pool generation, and how those properties relate to the correct characterization of the associated atmospheric hazards (wind gusts and hail). Use of a double-moment, seven-class scheme with number concentrations for all species (including hail and graupel) results in the greatest fidelity of model-simulated wind gusts and convective structure to the observations of this event. All ensemble members, however, fail to capture the intensity of the event in terms of the spatial extent of convection and the production of high near-surface wind gusts. We further show very high sensitivity to the LBCs employed and specifically that simulation fidelity is higher for simulations nested within ERA-Interim compared to ERA5. Excess convective available potential energy (CAPE) in all ensemble members after the derecho passage leads to excess production of convective cells, wind gusts, cREF > 40 dBZ, and precipitation during a frontal passage on the subsequent day. This event proved very challenging to forecast in real time and to reproduce in the 15-member hindcast simulation ensemble presented here. Future work could examine if simulations with other initial and lateral boundary conditions can achieve greater fidelity.

Shepherd, Tristan (ORCID:0000000186276419)↗

The NASA/GISS Mars general circulation model: Preliminary experiments

The NASA/GISS Mars General Circulation Model (GCM) is an adapted version of the GISS Global Climate/Middle Atmosphere Model, specifically developed for the diagnostic validation and objective analysis of measured atmospheric temperatures from the Mars Observer Pressure Modulator Infrared Radiometer (PMIRR) experiment. The GISS Mars GCM has 23 vertical layers extending from the surface to approximately 80 km altitude, representing a vertical resolution of about 0.3 scale heights. The primitive (vertically hydrostatic) equations are solved in finite difference form on the Krakawa B grid, with a horizontal resolution of 8 deg x 10 deg (latitude-longitude). The model includes a diurnal solar cycle, heat transport within a two-layer ground, and a high-order 'slopes-scheme' for the advection of heat in the upper atmosphere. The radiative transfer scheme is based on the correlated k distribution method for the treatment of nongray gaseous absorption thermal emission, and multiple scattering, including options for suspended dust. A special feature of the model of particular importance for Mars is a parameterization of gravity-wave-induced drag incorporating orographic forcing, wind shear, convection, and radiative damping. The implementation of the GISS Mars model includes global maps of topography, roughness, and albedo.

Allison, Michael↗

Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications

Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind-farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, which contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind-farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches relies heavily on ad hoc processes. For instance, the widely used Rayleigh damping method requires ad hoc tuning to produce a setup that may only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct a large-eddy simulation (LES) study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics. In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This, in turn, will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths (the representative wavelengths of the entire wave spectrum) to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency. Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upward-propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow–outflow boundaries is of equal importance to damping the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.

17 WIND ENERGY↗

Characterization of Field-Exposed Photovoltaic Modules Featuring Signs of Contact Degradation

Here, this work investigates several photovoltaic (PV) modules that have shown signs of metal contact corrosion due to field exposure in a hot and humid climate. This includes two multicrystalline silicon aluminum back surface field systems with 10 and 14 years of exposure and one monocrystalline silicon passivated emitter and rear cell system with four years of exposure. A comprehensive, multiscale characterization process is used to evaluate these PV modules in great detail. Current–voltage (I−V), Suns-V OC measurements, electroluminescence imaging, infrared imaging, and ultraviolet fluorescence imaging were performed, and locations of interest were cored and analyzed using cross-sectional scanning electron microscopy (SEM). A rigorous, quantitative analysis procedure for the cross-sectional SEM images is proposed and implemented. Careful characterization does reveal that some of these PV modules do indeed exhibit the same classic signs of acetic-acid-based corrosion of the glass frit that is present at the silver/silicon interface, which have been observed previously in PV modules exposed to damp heat in an environmental chamber.

14 SOLAR ENERGY↗

Chaotic /strange/ and periodic behavior in instability saturation by the oscillating two-stream instability

The nonlinear Schroedinger equation with linear growth and damping is truncated to three waves. The resulting system of nonlinear ordinary differential equations describes the excitation of linearly damped waves by the oscillating two-stream instability driven by a linearly unstable pump wave. This system represents a simple model for the nonlinear saturation of a linearly unstable wave. The model is examined analytically and numerically as a function of the dimensionless parameters of the system. It is found that the model can exhibit a wealth of characteristic dynamical behavior including stationary equilibria, Hopf bifurcations to periodic orbits, period doubling bifurcations, chaotic solutions characteristic of a strange attractor, tangent bifurcations from chaotic to periodic solutions, transient chaos, and hysteresis. Many of these features are shown to be explainable on the basis of one-dimensional maps. In the case of chaotic solutions, evidence for the presence of a strange attractor is provided by demonstrating Cantor set-like structure (i.e., scale invariance) in the surface of section.

Russell, D. A.↗

Fermi surface and Berry phase analysis for Dirac nodal line semimetals: Cautionary tale from SrGa 2 and BaGa 2

A Berry phase of odd multiples of 𝜋 inferred from quantum oscillations (QOs) has often been treated as evidence for nontrivial reciprocal space topology. However, disentangling the Berry phase values from the Zeeman effect and the orbital magnetic moment is often challenging. In centrosymmetric compounds, the case is simpler as the orbital magnetic moment contribution is negligible. Although the Zeeman effect can be significant, it is usually overlooked in most studies of QOs in centrosymmetric compounds. Here, we present a detailed study on the nonmagnetic centrosymmetric SrGa 2 and BaGa 2 , which are predicted to be Dirac nodal line semimetals based on density functional theory (DFT) calculations. Evidence of the nontrivial topology is found in magnetotransport measurements. The Fermi surface topology and band structure are carefully studied through a combination of angle-dependent QOs, angle-resolved photoemission spectroscopy (ARPES), and DFT calculations, where the nodal line is observed in the vicinity of the Fermi level. Strong de Haas–van Alphen fundamental oscillations associated with higher harmonics are observed in both compounds, which are well fitted by the Lifshitz-Kosevich (LK) formula. However, even with the inclusion of higher harmonics in the fitting, we found that the Berry phases cannot be unambiguously determined when the Zeeman effect is included. We revisit the LK formula and analyze the phenomena and outcomes that were associated with the Zeeman effect in previous studies. Our experimental results confirm that SrGa 2 and BaGa 2 are Dirac nodal line semimetals. Additionally, we highlight the often overlooked role of spin-damping terms in Berry phase analysis.

36 MATERIALS SCIENCE↗

Effect of Baffle on Gravity-Gradient-Excited Slosh Waves and Spacecraft Moment and Angular-Momentum Fluctuations in Microgravity

The dynamical behavior of fluids affected by the asymmetric gravity gradient acceleration has been investigated. In particular, the effects of surface tension on partially filled rotating fluids applicable to a full-scale Gravity Probe-B Spacecraft dewar tank with and without baffles are studied. Results of slosh wave excitation along the liquid-vapor interface induced by gravity gradient acceleration indicate that the gravity gradient acceleration is equivalent to the combined effect of a twisting force and a torsional moment acting on the spacecraft. The results are clearly seen from one-up one-down and one-down one-up oscillations in the cross-section profiles of two bubbles in the vertical (r, z)-plane of the rotating dewar, and from the eccentric contour of the bubble rotating around the axis of the dewar in a horizontal (r, theta)-plane. As the viscous force, between liquid and solid interface, greatly contributes to the damping of slosh wave excitation, a rotating dewar with baffles provides more areas of liquid-solid interface than that of a rotating dewar without baffles. Results show that the damping effect provided by the baffles reduces the amplitude of slosh wave excitation and lowers the degree of asymmetry in liquid-vapor distribution. Fluctuations of angular momentum and fluid moment caused by the slosh wave excited by gravity gradient acceleration with and without baffle boards are also investigated. It is also shown that the damping effect provided by the baffles greatly reduces the amplitudes of angular momentum and fluid moment fluctuations.

Hung, R. J.↗

Two-dimensional simulations of internal gravity waves in a $\mathrm{5M_⊙}$ zero-age-main-sequence model

Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91 percent of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low-frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to 10 4 times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.

79 ASTRONOMY AND ASTROPHYSICS↗

Error-Field Penetration Thresholds in Ohmically Heated ITER and SPARC Plasmas

The critical n=1, n=2, and n=3 error-field amplitudes needed to trigger error-field penetration in ITER and SPARC, steady-state, ohmically heated plasmas are calculated using a standard asymptotic matching approach. The calculation incorporates plasma impurities, trapped particles, the bootstrap current, and neoclassical poloidal flow-damping. The energy confinement time is specified by the neo-Alcator scaling law in the low-density linear ohmic confinement (LOC) regime, and by the ITER-89P L-mode scaling law in the high-density saturated ohmic confinement regime (SOC) regime. The response of the plasma in the inner region is calculated using a linearized version of the four-field model. At the normal operating electron number density, diamagnetic levels of rotation are found to be sufficient to protect ITER and SPARC ohmically heated plasmas from m=2/n=1 error-field penetration. On other hand, SPARC, and especially ITER, ohmically heated plasmas may be vulnerable to n>1 error-field penetration. ITER and SPARC ohmically heated plasmas are also slightly more susceptible to error-field penetration when the electron fluid at the rational surface rotates in the ion diamagnetic direction, rather than the electron diamagnetic direction. At electron number densities that are sufficiently low that the plasma lies in the LOC confinement regime, the error-field penetration threshold increases with increasing density. However, as soon as the electron number density becomes large enough that the plasma enters the SOC regime, the increase in the error-field penetration threshold with increasing density levels off.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamics of the Trans-Neptune Region: Apsidal Waves in the Kuiper Belt

The role of apsidal density waves propagating in a primordial trans-Neptune disk (i.e., Kuiper belt) is investigated. It is shown that Neptune launches apsidal waves at its secular resonance near 40 AU that propagate radially outward, deeper into the particle disk. The wavelength of apsidal waves is considerably longer than waves that might be launched at Lindblad resonances, because the pattern speed, g(sub s), resulting from the apsis precession of Neptune is much slower than its mean motion, Omega(sub s). If the early Kuiper belt had a sufficient surface density, sigma, the disk's wave response to Neptune's secular perturbation would have spread the disturbing torque radially over a collective scale lambda(sub *) approx. = r(2(mu)(sub d)Omega/ absolute value of r dg/dr)(sup 1/2), where mu(sub d)equivalent pi(sigma)r(exp 2)/(1 solar mass) and Omega(r) and g(r) are respectively the mean motion and precession frequency of the disk particles. This results in considerably smaller eccentricities at resonance than had the disk particles been treated as noninteracting test particles. Consequently, particles are less apt to be excited into planet-crossing orbits, implying that the erosion timescales reported by earlier test-particle simulations of the Kuiper belt may be underestimated. It is also shown that the torque the disk exerts upon the planet (due to its gravitational attraction for the disk's spiral wave pattern) damps the planet's eccentricity and further inhibits the planet's ability to erode the disk. Key words: celestial mechanics, stellar dynamics - comets: general minor planets, asteroids

Ward, William R.↗