Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Instability control”

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 1,387 records · Page 77

Novel Power-Hardware-in-the-Loop Interface Method for Grid-Forming Inverter Systems

Power-hardware-in-the-loop (PHIL) simulations of grid-forming (GFM) inverter systems facilitate the testing of drastic scenarios, such as on-grid to off-grid transitions and islanded microgrid operations without a stiff grid. To the authors’ best knowledge, most studies in the literature focus on PHIL simulations for grid-following inverter systems. Only a few studies focus on GFM inverters, and those are challenging and problematic, especially for high-power applications. This article proposes a novel PHIL simulation platform that enables interfacing high-power GFM inverter systems. The paper proposes the concept of a virtual GFM inverter as a part of the proposed PHIL interface. This addition of a virtual GFM inverter in the PHIL interface expands the conventional ideal transformer model (ITM) method and enables it to overcome the issues of instability of existing ITM methods. In the validation stage, a PHIL experiment is conducted on a three-phase, 480-V, 125-kVA GFM inverter system with the proposed interfacing method. The results corroborate that the proposed PHIL simulation method performs well and is stable for GFM inverter systems.

droop control↗

First indirectly driven liquid-DT filled double shell implosions at the National Ignition Facility

Double shell implosions aim to explore material mixing under fusion conditions in a volume burn geometry using high-Z metal pushers. High-Z pushers are more compressible than low-Z pushers enabling high stagnation pressure which reduces the required implosion speed while maintaining a low pusher adiabat despite strong shock heating. Additionally, the use of a small fuel mass reduces the fuel internal energy required for ignition, thus achieving stable platforms with reasonable fusion output to conduct controlled experiments. These factors make volume burn in a double shell implosion highly promising. Recently, a series of liquid-DT filled, indirectly driven, double shell implosions were conducted at the National Ignition Facility with laser drives reaching up to 1.5 MJ. These experiments achieved a maximum DT neutron yield of 1.67 × 10 14 (yield—479 J), DT ion temperature of 2.6 keV, fuel areal density (ρR) of 0.15 g/cm 2 , and stagnation pressure of 79 Gbar. Over the course of these shots, the DT neutron yield has increased by an order of magnitude largely from improved mitigation of outer shell assembly joint driven instability growth using thicker gold plating at the assembly joint. Further performance improvements are expected by enhancing outer-to-inner shell kinetic energy transfer and refined mitigation of degradations from engineering features such as the outer shell joint, fill tube, and surface roughness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flame flashback in a premixed dump combustor

A stainless steel, rectangular center-dump, premixed-prevaporized combustor was used at the NASA Lewis Research Center to study flashback in controlled conditions simulating gas turbine combustor operations under steady-state and transient conditions. The combustor allowed visual access to the combustion area, and was operated at inlet air temperatures from 600 to 850 K, premixer wall temperatures from 450 to 900 K, and average premixer velocity from 40 to 80 ft/s. A parametric study revealed a slight decrease in the flashback equivalence ratio as the inlet air temperatures increased. It also indicated that the premixer wall temperature and premixer velocity are not governing parameters of flashback. It is suggested that flashback occurs through the premixer wall boundary layer flow reversal caused by combustion instability.

Proctor, M. P.↗

Detailed Multidimensional Simulations of the Structure and Dynamics of Flames

Numerical simulations in which the various physical and chemical processes can be independently controlled can significantly advance our understanding of the structure, stability, dynamics and extinction of flames. Therefore, our approach has been to use detailed time-dependent, multidimensional, multispecies numerical models to perform carefully designed computational experiments of flames on Earth and in microgravity environments. Some of these computational experiments are complementary to physical experiments performed under the Microgravity Program while others provide a fundamental understanding that cannot be obtained from physical experiments alone. In this report, we provide a brief summary of our recent research highlighting the contributions since the previous microgravity combustion workshop. There are a number of mechanisms that can cause flame instabilities and result in the formation of dynamic multidimensional structures. In the past, we have used numerical simulations to show that it is the thermo-diffusive instability rather than an instability due to preferential diffusion that is the dominant mechanism for the formation of cellular flames in lean hydrogen-air mixtures. Other studies have explored the role of gravity on flame dynamics and extinguishment, multi-step kinetics and radiative losses on flame instabilities in rich hydrogen-air flames, and heat losses on burner-stabilized flames in microgravity. The recent emphasis of our work has been on exploring flame-vortex interactions and further investigating the structure and dynamics of lean hydrogen-air flames in microgravity. These topics are briefly discussed after a brief discussion of our computational approach for solving these problems.

Patnaik, G.↗

The dynamics and control of large flexible space structures, 8

A development of the in plane open loop rotational equations of motion for the proposed Spacecraft Control Laboratory Experiment (SCOLE) in orbit configuration is presented based on an Eulerian formulation. The mast is considered to be a flexible beam connected to the (rigid) shuttle and the reflector. Frequencies and mode shapes are obtained for the mast vibrational appendage modes (assumed to be decoupled) for different boundary conditions based on continuum approaches and also preliminary results are obtained using a finite element representation of the mast reflector system. The linearized rotational in plane equation is characterized by periodic coefficients and open loop system stability can be examined with an application of the Floquet theorem. Numerical results are presented to illustrate the potential instability associated with actuator time delays even for delays which represent only a small fraction of the natural period of oscillation of the modes contained in the open loop model of the system. When plant and measurement noise effects are added to the previously designed deterministic model of the hoop column system, it is seen that both the system transient and steady state performance are degraded. Mission requirements can be satisfied by appropriate assignment of cost function weighting elements and changes in the ratio of plant noise to measurement noise.

Bainum, P. M.↗

Satellite range delay simulator for a matrix-switched time division multiple-access network simulator

The Systems Integration, Test, and Evaluation (SITE) facility at NASA Lewis Research Center is presently configured as a satellite-switched time division multiple access (SS-TDMA) network simulator. The purpose of SITE is to demonstrate and evaluate advanced communication satellite technologies, presently embodied by POC components developed under NASA contracts in addition to other hardware, such as ground terminals, designed and built in-house at NASA Lewis. Each ground terminal in a satellite communications system will experience a different aspect of the satellite's motion due mainly to daily tidal effects and station keeping, hence a different duration and rate of variation in the range delay. As a result of this and other effects such as local oscillator instability, each ground terminal must constantly adjust its transmit burst timing so that data bursts from separate ground terminals arrive at the satellite in their assigned time slots, preventing overlap and keeping the system in synchronism. On the receiving end, ground terminals must synchronize their local clocks using reference transmissions received through the satellite link. A feature of the SITE facility is its capability to simulate the varying propagation delays and associated Doppler frequency shifts that the ground terminals in the network have to cope with. Delay is achieved by means of two NASA Lewis designed and built range delay simulator (RDS) systems, each independently controlled locally with front panel switches or remotely by an experiment control and monitor (EC/M) computer.

Nagy, Lawrence A.↗

Satellite range delay simulator for a matrix-switched time division multiple-access network simulation system

The Systems Integration, Test, and Evaluation (SITE) facility at NASA Lewis Research Center is presently configured as a satellite-switched time division multiple access (SS-TDMA) network simulator. The purpose of SITE is to demonstrate and evaluate advanced communication satellite technologies, presently embodied by POC components developed under NASA contracts in addition to other hardware, such as ground terminals, designed and built in-house at NASA Lewis. Each ground terminal in a satellite communications system will experience a different aspect of the satellite's motion due mainly to daily tidal effects and station keeping, hence a different duration and rate of variation in the range delay. As a result of this and other effects such as local oscillator instability, each ground terminal must constantly adjust its transmit burst timing so that data bursts from separate ground terminals arrive at the satellite in their assigned time slots, preventing overlap and keeping the system in synchronism. On the receiving end, ground terminals must synchronize their local clocks using reference transmissions received through the satellite link. A feature of the SITE facility is its capability to simulate the varying propagation delays and associated Doppler frequency shifts that the ground terminals in the network have to cope with. Delay is ahcieved by means of two NASA Lewis designed and built range delay simulator (RDS) systems, each independently controlled locally with front panel switches or remotely by an experiment control and monitor (EC/M) computer.

Nagy, Lawrence A.↗

New millimeter-wave diagnostics to locally probe internal density and magnetic field fluctuations in National Spherical Torus Experiment-Upgrade (invited)

A set of new millimeter-wave diagnostics will deliver unique measurement capabilities for National Spherical Torus Experiment-Upgrade to address a variety of plasma instabilities believed to be important in determining thermal and particle transport, such as micro-tearing, global Alfvén eigenmodes, kinetic ballooning, trapped electron, and electron temperature gradient modes. These diagnostics include a new integrated intermediate-k Doppler backscattering (DBS) and cross-polarization scattering (CPS) system (four channels, 82.5–87 GHz) to measure density and magnetic fluctuations, respectively. The system can access reasonably large normalized wavenumbers k θ ρ s ranging from ≤0.5 to 15 (where ion sound gyroradius ρ s = 1 cm and k θ is the binormal density turbulence wavenumber). The system addresses the challenges for making useful DBS/CPS measurements with a remote control of launch polarization (X- or O-mode), probed wavenumber, polarization match of the launch beam with the edge magnetic field pitch angle, and beam steering of the launched beam for wave-vector alignment. In addition, a low-k DBS system consisting of eight fixed frequencies (34–52 GHz) and four tunable frequencies (55–75 GHz) for low-k density turbulence and fast ion physics will be located at a nearby port location. Finally, the combined systems cover the near LCFS and pedestal regions (34–52 GHz), the pedestal or mid-radius (50–75 GHz), and core plasmas (82.5–87 GHz).

47 OTHER INSTRUMENTATION↗

Acoustic properties of supersonic helium/air jets at low Reynolds numbers

Experiments have been performed with the objective of developing a greater understanding of the physics of hot supersonic jet noise. Cold helium/air jets are used to easily and inexpensively simulate the low densities of hot air jets. The experiments are conducted at low Reynolds numbers in order to facilitate study of the large-scale turbulent structures (instability waves) that cause most of the radiated noise. Experiments have been performed on Mach 1.5 and 2.1 jets of pure air, pure helium and 10 percent helium by mass. Helium/air jets are shown to radiate more noise than pure air jets due to the increased exit velocity. Microphone spectra are usually dominated by a single spectral component at a predictable frequency. Increasing the jet's helium concentration is shown to increase the dominant frequency. The helium concentration in the test chamber is determined by calculating the speed of sound from the measured phase difference between two microphone signals. Bleeding outside air into the test chamber controls the accumulation of helium so that the hot jet simulation remains valid. The measured variation in the peak radiated noise frequency is in good agreement with the predictions of the hot jet noise theory of Tam et al.

Mclaughlin, Dennis K.↗

Aging phenomena in the two-dimensional complex Ginzburg-Landau equation

The complex Ginzburg-Landau equation with additive noise is a stochastic partial differential equation that describes a remarkably wide range of physical systems which include coupled non-linear oscillators subject to external noise near a Hopf bifurcation instability and spontaneous structure formation in non-equilibrium systems, e.g., in cyclically competing populations or oscillatory chemical reactions. Here, we employ a finite-difference method to numerically solve the noisy complex Ginzburg-Landau equation on a two-dimensional domain with the goal to investigate its non-equilibrium dynamics when the system is quenched into the “defocusing spiral quadrant”. We observe slow coarsening dynamics as oppositely charged topological defects annihilate each other, and characterize the ensuing aging scaling behavior. We conclude that the physical aging features in this system are governed by non-universal aging scaling exponents. We also investigate systems with control parameters residing in the “focusing quadrant”, and identify slow aging kinetics in that regime as well. Finally, we provide heuristic criteria for the existence of slow coarsening dynamics and physical aging behavior in the complex Ginzburg-Landau equation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Rethinking PCMI Failures Under Transient Conditions

Under PCMI conditions, the cladding strain is driven by either pellet swelling during normal operation or pellet thermal expansion in transients. The loading in PCMI is displacement controlled because friction between the pellet and the cladding results in a strain compatibility at the pellet/cladding interface. Failure is known to occur under PCMI loading conditions during operation when normal power changes result in release of iodine (and other fission products) causing stress corrosion cracking where the crack-propagation rate is dependent upon the cladding stress but slow relative to strictly mechanical failures. Failures after normal power changes are often referred to as failure by pellet-cladding interaction or PCI to distinguish it from PCMI. Such failures are considered displacement-controlled failures, as noted above, and occur with limited plastic strain. This is different to what is typically thought of as PCMI failure under transient conditions, either Anticipated Operational Occurrences (AOOs) or Reactivity Insertion Accidents (RIAs), which is generally assumed to occur once the cladding reaches a mechanical strain limit. But there are two features in the typical PCMI failure in an RIA that are not consistent with the PCMI assumptions. The first is the ductile shear failure feature in cross section and the second is a large crack opening post failure (e.g., see OI-11 failure in NSRR program). The ductile feature is indicative of a plastic instability that requires very large local strains, and the large crack opening must be the result of a release of energy, both of which are not expected if stress in the cladding was only driven by thermal expansion in the pellet. It is concluded that there must be a mechanism (likely related to rod internal pressure) that eliminates strain compatibility between the pellet and cladding – allowing strain localization to occur in what are known as PCMI transient failures. In some sense, these PCMI failures may actually be early burst failures that occur when thermal expansion causes uniform strains greater than a defined limit and fission gas pressure eliminate the strain compatibility

Cantonwine, Paul↗

Transition Core Planning and Safety Analyses in Support of LEU Fuel Conversion of the University of Missouri Research Reactor (MURR)

The University of Missouri Research Reactor (MURR®) is one of six U.S. High Performance Research Reactors (USHPRR), including one critical facility, that is working with the National Nuclear Security Administration (NNSA) Office of Material Management and Minimization (M3) Reactor Conversion Program to convert from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel. The M3 Reactor Conversion USHPRR Project objectives include the development of LEU fuel element designs that will ensure safe reactor operations and to maintain the existing experimental performance of each facility. The work is being conducted through many inter-related activities being completed by four Project Pillars: Fuel Qualification (FQ), Fuel Fabrication (FF), Reactor Conversion (RC), and Cross Cutting (CC). A new type of LEU fuel based on an alloy of uranium-10 wt% molybdenum (U-10Mo) is expected to allow the conversion of those USHPRR, like MURR, requiring higher density fuels. The very-high-density LEU U-10Mo monolithic fuel is currently undergoing irradiation testing and post-irradiation examination under a planned and documented fuel qualification effort. The FQ Pillar will document fuel property and fuel performance data and qualify the fuel for use in these reactors. The FF Pillar is fabricating fuel for ongoing and future irradiation tests, as well as conducting fabrication demonstrations to validate or update preliminary fabrication assumptions. The FF Pillar is also working to develop and install commercial manufacturing capacity with the U-10Mo monolithic fuel to produce prototypic fuel. Working with the RC Pillar at Argonne, MURR has progressed through a preliminary fuel element design using preliminary data for the proposed monolithic alloy of U-10Mo. Analyses were completed in previous work that found for typical equilibrium operations with the preliminary LEU fuel element design, in conjunction with a power uprate to 12 MW and appropriate changes to the MURR Limiting safety system settings (LSSS), MURR will have adequate margins to safety for steady-state operations and postulated transient accidents and will have experimental performance in key locations that meets or exceeds current operations with HEU fuel. The purpose of this work is to develop a sequence of transition cycles that will enable MURR to transition from operation with the reactor core loaded with fresh LEU fuel elements only to typical equilibrium operations with mixed-burnup cores following conversion while meeting operational requirements on safety and experimental performance. It is expected that the use of fresh LEU fuel at conversion and subsequent low burnup of the LEU fuel elements that will initially be available for use following conversion will result in critical control blade positions that will substantially change the axial power distribution in the core and the neutron flux available in key experimental locations relative to equilibrium LEU operations. Given the constraints of MURR safety margins, operational practices, and production and research, a novel method has been developed to identify a transition sequence that minimizes the time MURR operates atypically compared to the current prototypic cycles using HEU fuel. The proposed transition sequence moves quickly to the same sort of equilibrium cycles for the LEU fuel that have already been evaluated in documented preliminary safety analyses. Although shifting the neutron flux peak to the lower half of the core during initial cycles with LEU at 12 MW reduces the experiment performance in some key locations relative to current HEU operations at 10 MW, all LEU cores provide an average performance that meets or exceeds that of HEU. An LEU cycle is reached that meets or exceeds the level of experimental performance predicted for current HEU and equilibrium LEU operations in more than 450 key locations identified by a reactor specialist at MURR by the 23rd cycle following conversion and that afterwards will enable MURR to consistently meet its experimental performance requirements. The proposed transition sequence only requires the fabrication of 34 fresh LEU elements in the first year of operation and does not exceed the anticipated availability of fresh elements that can be produced by the fuel fabricator. By the third year after conversion, 22 fresh LEU elements will be required each year, which is the same as expected for equilibrium LEU operations and the same as current operations with HEU fuel. The proposed transition sequence thus combines a relatively short time period before equilibrium burnup is achieved, a temporary increase of fuel elements needed annually relative to typical operations that are within the production capabilities of the fuel fabricator, and demonstrates comparable experimental performance of the LEU cores relative to current HEU operations. Further measures may be taken to reduce any initial experimental performance penalty even further, where possible, by repositioning certain experiments to leverage the increased performance in the lower axial experimental positions in the initial cycles following conversion or leaving the experiments in the irradiation facilities longer in order to achieve the required neutron fluence. This analysis may require refinement depending on the experimental facilities in use at the time of conversion. Nonetheless, the results presented here, including the experimental performance, core burnup, and critical control blade positions throughout the transition cycles, show that the proposed transition cycle fuel management patterns are consistent with what is expected and desired for MURR operation with LEU U-10Mo fuel. Detailed core power distributions from the neutronics models were also used to evaluate safety margins during steady-state operations for the selected transition cycles and the equilibrium LEU core. It is shown that there are adequate safety margins for both steady-state operations and postulated accident scenarios. For the steady-state operations with the preliminary LEU fuel element design the analysis predicts at least 2.49 MW margin to the onset of flow instability at the LSSS power of 15 MW. Considering the LSSS power is 125% of full license power, the margin to OFI is sufficient. In addition, the critical heat flux ratio at LSSS power is well above the requirement of CHFR > 2.0 from NUREG-1537 for all considered cases. For postulated transient accidents, the minimum margin to the fuel temperature safety limit is at least 109 °C. In summary, the proposed sequence of core loadings for MURR operations following conversion to LEU fuel and a power uprate to 12 MW provides sufficient safety margins for both steady-state operations and postulated transient accidents during a proposed sequence of transition cycles to equilibrium operations. Analysis has shown that there are some local experimental performance penalties during the initial cycles. Although there are local shifts in the experimental performance, on average all LEU cores at 12 MW have equal or higher performance than HEU at 10 MW. Temporary adjustments are being planned that will produce suitable experimental performance during these cycles. The results indicate that for the equilibrium LEU core the experimental performance exceeds that of current HEU operations in all key locations while also demonstrating sufficient safety margins.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Vibration Challenges in the Design of NASA's Ares Launch Vehicles

This paper focuses on the vibration challenges inherent in the design of NASA s Ares launch vehicles. A brief overview of the launch system architecture is provided to establish the context for the discussion. Following this is a general discussion of the design considerations and analytical disciplines that are affected by vibration. The first challenge discussed is that of coupling between the vehicle flight control system and fundamental vibrational modes of the vehicle. The potential destabilizing influence of the vibrational dynamics is described along with discussion of the typical methods employed to overcome this issue. Next is a general discussion of the process for developing the design loads for the primary structure. This includes quasi-steady loads and dynamic loads induced by the structural dynamic response. The two principal parts of this response are the gust induced responses of the lower frequency modes and the buffet induced responses of the higher frequency modes. Structural dynamic model validation will also be addressed. Following this, discussions of three somewhat unique topics of Pogo Instability, Solid Booster Thrust Oscillation, and Liquid Rocket Engine Turbopump Rotordynamic Stability and Response are presented.

Ryan, Stephen G.↗

Preliminary Results of Stability and Control Investigation of the Bell X-5 Research Airplane

During the acceptance tests of the Bell X-5 airplane, measurements of the static stability and control characteristics and horizontal-tail loads were obtained by the NACA High-Speed Flight Research Station. The results of the stability and control measurements are presented in this paper. A change in sweep angle between 20 deg and 59 deg had a minor effect on the longitudinal trim, with a maximum change of about 2.5 deg in elevator deflection being required at a Mach number near 0.85; however, sweeping the wings produced a total stick-force change of about 40 pounds. At low Mach numbers there was a rapid increase in stability at high normal-force coefficients for both 20 0 and 1100 sweepback, whereas a condition of neutral stability existed for 58 0 sweepback at high normal-force coefficients. At Mach numbers near 0.8 there was an instability at normal-force coefficients above 0.5 for all sweep angles tested. In the low normal-force-coefficient range a high degree of stability resulted in high stick forces which limited the maximum load factors attainable in the demonstration flights to values under 5g for all sweep angles at a Mach number near 0.8 and an altitude of 12,000 feet. The aileron effectiveness at 200 sweepback was found to be low over the Mach number range tested.

Finch, Thomas W↗

An Improved Approach to the Predictability & Reliability of the Onset of Turbulence With Shocks

The construction of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and for (b) obtaining correct propagation speed of discontinuities in the presence of stiff source terms share one important ingredient – minimization of numerical dissipation while maintaining numerical stability. The dual requirements to achieve both numerical stability and minimal numerical dissipation are often conflicting since existing shock capturing schemes were designed mainly to be robust for rapidly developed turbulence-free flows and for shock waves without stiff source term. For the past two decades, Yee and collaborators have focused on an improved understanding of the nonlinear behavior of different high order shock-capturing methods. It was found that even very high order methods without proper nonlinear stability and numerical dissipation control can either numerically smear the onset of turbulence due to excess numerical dissipation, or induce (onset) numerical turbulence that is not physical turbulence due to lack of proper numerical dissipation to improve nonlinear stability for long time integration. Our approach is to combine (I) and (II) below for obtaining the physically correct onset of turbulence with shocks, including problems with stiff source terms: (I) Nonlinear dynamics is utilized to complement the traditional linearized stability theory (Yee & Sweby, Yee et al., Griffiths et al., Lafon & Yee, Yee, Wang et al., Kotov et al. 1990- 2015) in order to (i) Minimize numerically induced false transition to turbulence, (ii) Minimize numerical instability due to long time integration of turbulent flows, (iii) Minimize numerically induced standing wave solutions, and (iv) Minimize wrong propagation of speed of discontinuities due to the presence of stiff source terms. (II) Our recently developed physical preserving (structural preserving) high order methods with improved nonlinear stability & accuracy that are essential in minimizing spurious numerics are used.

HECC↗

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

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

Bainum, P. M.↗

Impact of pedestal density gradient and collisionality on ELM dynamics

BOUT++ turbulence simulations are conducted to capture the underlying physics of small ELM characteristics achieved by increasing separatrix density via controlling strike points from vertical to horizontal divertor plates for three EAST discharges. BOUT++ linear simulations show that the most unstable modes change from high-n ideal ballooning modes to intermediate-n peeling–ballooning modes and eventually to peeling–ballooning stable plasmas in the pedestal. Nonlinear simulations show that the fluctuation is saturated at a high level for the lowest separatrix density. The ELM size decreases with increasing separatrix density, until the fraction of this energy lost during the ELM crash becomes less than 1% of the pedestal stored energy, leading to small ELMs. Simulations indicate that small ELMs can be triggered either by the marginally peeling–ballooning instability near the peak pressure gradient position inside the pedestal or by a local instability in the pedestal foot with a larger separatrix density gradient. The pedestal collisionality scan for type-I ELMs with steep pedestal density gradient shows that both linear growth rate and ELM size decrease with increasing collisionality. On the contrary, the pedestal collisionality and pedestal density width scan with a weak pedestal density gradient indicate small ELMs can either be triggered by a high-n ballooning mode or by a low-n peeling mode in a low collisionality region 0.04–0.1. We report the simulations indicate the weaker the linear unstable modes near marginal stability with small linear growth rate, the lower nonlinearly saturated fluctuation intensity and the smaller turbulence spreading from the linear unstable zone to stable zone in the nonlinear saturation phase, leading to small ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Studying marine stratus with large eddy simulation

Data sets from field experiments over the stratocumulus regime may include complications from larger scale variations, decoupled cloud layers, diurnal cycle, or entrainment instability, etc. On top of the already complicated turbulence-radiation-condensation processes within the cloud-topped boundary layer (CTBL), these complexities may sometimes make interpretation of the data sets difficult. To study these processes, a better understanding is needed of the basic processes involved in the prototype CTBL. For example, is cloud top radiative cooling the primary source of the turbulent kinetic energy (TKE) within the CTBL. Historically, laboratory measurements have played an important role in addressing the turbulence problems. The CTBL is a turbulent field which is probably impossible to generate in laboratories. Large eddy simulation (LES) is an alternative way of 'measuring' the turbulent structure under controlled environments, which allows the systematic examination of the basic physical processes involved. However, there are problems with the LES approach for the CTBL. The LES data need to be consistent with the observed data. The LES approach is discussed, and results are given which provide some insights into the simulated turbulent flow field. Problems with this approach for the CTBL and information from the FIRE experiment needed to justify the LES results are discussed.

Moeng, Chin-Hoh↗