Search NASA⌕ Search

SEARCH · Search NASA

Results for “2D Gravity”

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 37 records · Page 2

Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets

Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.

79 ASTRONOMY AND ASTROPHYSICS↗

MHD Modelling of Coronal Loops: Injection of High-Speed Chromospheric Flows

Context. Observations reveal a correspondence between chromospheric type II spicules and bright upward-moving fronts in the corona observed in the extreme-ultraviolet (EUV) band. However, theoretical considerations suggest that these flows are probably not the main source of heating in coronal magnetic loops. Aims. We investigate the propagation of high-speed chromospheric flows into coronal magnetic flux tubes and the possible production of emission in the EUV band. Methods. We simulated the propagation of a dense 104 K chromospheric jet upward along a coronal loop by means of a 2D cylindrical MHD model that includes gravity, radiative losses, thermal conduction, and magnetic induction. The jet propagates in a complete atmosphere including the chromosphere and a tenuous cool (approximately 0.8 MK) corona, linked through a steep transition region. In our reference model, the jet initial speed is 70 km per second, its initial density is 10(exp 11) per cubic centimeter, and the ambient uniform magnetic field is 10 G. We also explored other values of jet speed and density in 1D and different magnetic field values in 2D, as well as the jet propagation in a hotter (approximately 1.5 MK) background loop. Results. While the initial speed of the jet does not allow it to reach the loop apex, a hot shock-front develops ahead of it and travels to the other extreme of the loop. The shock front compresses the coronal plasma and heats it to about 10(exp 6) K. As a result, a bright moving front becomes visible in the 171 Angstrom channel of the SDO/AIA mission. This result generally applies to all the other explored cases, except for the propagation in the hotter loop. Conclusions. For a cool, low-density initial coronal loop, the post-shock plasma ahead of upward chromospheric flows might explain at least part of the observed correspondence between type II spicules and EUV emission excess.

Sun: chromosphere↗

Evidence for an Early Formation of Serenitatis Basin at 4.25 Ga Shifts Lunar Chronology

Troctolite sample 76535, collected in Serenitatis basin during Apollo 17, formed at least 50 km deep, experienced maximum shock pressures of 6 GPa, and has a 40 Ar/ 39 Ar excavation age of 4.25 Ga. Previous work attributed 76535 to the South Pole-Aitken (SPA) basin, presumably dating the SPA-impact and constraining lunar bombardment history. Here we use the iSALE-2D shock-physics code and gravity inversion modeling to determine if instead the Serenitatis impact event excavated 76535. We find nearly 140,000 km 3 of material (∼2% of near-surface ejecta) matching the depth and pressure constraints of 76535 is displaced to the surface during crater collapse of a Serenitatis-like impact event. We conclude that the Serenitatis impact event possibly excavated 76535, redefining its age to 4.25 Ga, 300 My older than the consensus age based on Apollo 17 samples. This finding would provide an important anchor point where lunar chronology where bombardment flux is especially uncertain.

Bjonnes, Evan [Lawrence Livermore National Laborat↗

Formation, levitation, and stability of prominences in the magnetized solar atmosphere

The dynamic formation of prominences in the initial magnetothermal equilibrium and their stability to sideward displacements are investigated focusing on the structure of the 2D solar atmosphere in the presence of coronal arcades or loops. A model based on 2D magnetohydrodynamic equations takes into account gravity, compressible flows, heating, radiation, anisotropic thermal conduction, and coupling to a deep chromosphere. It is found that prominences in simple arcades characterized by magnetic field with significant curvature at the apex are unstable to a lateral displacement.

Drake, J. F.↗

Monsoonal MCS Initiation, Rainfall, and Diurnal Gravity Waves over the Bay of Bengal: Observation and a Linear Model

Abstract Previous observational studies have indicated that mesoscale convective systems (MCSs) contribute the majority of precipitation over the Bay of Bengal (BoB) during the summer monsoon season, yet their initiation and propagation remain incompletely understood. To fill this knowledge gap, we conducted a comprehensive study using a combination of 20-yr satellite observations, MCS tracking, reanalysis data, and a theoretical linear model. Satellite observations reveal clear diurnal propagation signals of MCS initiation frequency and rainfall from the west coast of the BoB toward the central BoB, with the MCS rainfall propagating slightly slower than the MCS initiation frequency. Global reanalysis data indicate a strong association between the offshore-propagating MCS initiation frequency/rainfall and diurnal low-level wind perturbations, implying the potential role of gravity waves. To verify the hypothesis, we developed a 2D linear model that can be driven by realistic meteorological fields from reanalysis. The linear model realistically reproduces the characteristics of offshore-propagating diurnal wind perturbations. The wind perturbations, as well as the offshore propagation signals of MCS initiation frequency and rainfall, are associated with diurnal gravity waves emitted from the coastal regions, which in turn are caused by the diurnal land–sea thermal contrast. The ambient wind speed and vertical wind shear play crucial roles in modulating the timing, propagation, and amplitude of diurnal gravity waves. Using the linear model and satellite observations, we further show that the stronger monsoonal flows lead to faster offshore propagation of diurnal gravity waves, which subsequently control the offshore propagation signals of MCS initiation and rainfall. Significance Statement Rainfall over the Bay of Bengal (BoB) is primarily contributed by large and organized rainfall systems in the summer monsoon season. During this season, these systems are commonly observed over the east coast of India around midnight, the western BoB in the morning, and the central BoB in the afternoon. This eastward rainfall propagation is confirmed by observations, reanalysis data, and a theoretical model to have a strong association with atmospheric diurnal gravity waves. These waves are caused by land–sea thermal contrast and can trigger rainfall systems over the offshore regions. We also found that the diurnal gravity waves, as well as their triggered rainfall systems, can be greatly modulated by the large-scale monsoonal flows. All the above findings improve our understanding of diurnal rainfall cycle over the BoB and may contribute to the future improvement of rainfall forecast over the region.

Meteorology & Atmospheric Sciences↗

Effect of gravity on the stability and structure of lean hydrogen-air flames

Detailed, time-dependent, 2D numerical simulations with full hydrogen-oxygen chemistry are used to investigate the effects of gravity on the stability and structure of laminar flames in lean, premixed hydrogen-air mixtures. The calculations show that the effects of gravity becomes more important as the lean flammability limit is approached. In a 12 percent hydrogen-air mixture, gravity plays only a secondary role in determining the multidimensional structure of the flame with the stability and structure of the flame controlled primarily by the thermo-diffusive instability mechanism. However, in leaner hydrogen-air mixtures gravity becomes more important. Upward-propagating flames are highly curved and evolve into a bubble rising upwards in the tube. Downward-propagating flames are flat or even oscillate between structures with concave and convex curvatures. The zero-gravity flame shows only cellular structures. Cellular structures which are present in zero gravity can be suppressed by the effect of buoyancy for mixtures leaner than 11 percent hydrogen. These observations are explained on the basis of an interaction between the processes leading to buoyancy-induced Rayleigh-Taylor instability and the thermo-diffusive instability.

Patnaik, G.↗

The QBO as Potential Amplifier and Conduit to Lower Altitudes of Solar Cycle Influence?

The solar cycle (SC) effect in the lower atmosphere has been linked observationally to the Quasi-biennial Oscillation (QBO), which is generated primarily by small-scale gravity waves (GW). Following up on a 2D study with our Numerical Spectral Model (NSM), we discuss here a 3D study in which we simulated the QBO under the influence of the SC. For a SC period of 10 years, the relative amplitude of radiative forcing is taken to vary exponentially with height, i.e., 0.2% at the surface, 2% at 50 km, 20% at 100 km and above. Applying spectral analysis to filter out and identify the SC signature, this model produces a relatively large modulation in the QBO of the lower stratospheric circulation, which is in qualitative agreement with the results obtained by Salby and Callaghan (2000) who analyzed zonal wind observations covering more than 40 years. The modeled SC modulation of the QBO extends to high latitudes where it produces temperature variations of about 1 K in the troposphere below 10 km. We report that the SC also generates in the model a hemispherically symmetric Equatorial Annual Oscillation (EAO, with 12-month period) that is largely confined to low latitudes. Under the influence of the GWs, the SC modulated EA0 propagates down into the lower stratosphere like the QBO. As is the case for the QBO, the energy of this EA0 is partially redistributed by the meridional circulation and planetary waves, presumably, to generate measurable SC signatures in the tropospheric temperature of the polar regions, which may be related to the so called Arctic Oscillation (Thompson and Wallace, 1998). The larger SC influence at higher altitudes is apparently transferred to the lower and denser regions of the atmosphere by tapping the momentum from the upward propagating GWs that drive the oscillations. The SC modulation of the QBO period could prove to be very effective in this process, as our earlier 2D study indicated. Further studies are needed, (1) to make sure that the SC effects are real and the numerical results are robust, and (2) to explore more fully the mechanism(s) that may conspire to amplify the SC effect. Quasi-decadal oscillations, generated internally by the QBO interacting with the seasonal cycles, may interfere with, or aid, the SC influence.

Mayr, Hans G.↗

Combustion Characteristics in a Non-Premixed Cool-Flame Regime of n-Heptane in Microgravity

A series of distinct phenomena have recently been observed in single-fuel-droplet combustion tests performed on the International Space Station (ISS). This study attempts to simulate the observed flame behavior numerically using a gaseous n-heptane fuel source in zero gravity and a time-dependent axisymmetric (2D) code, which includes a detailed reaction mechanism (127 species and 1130 reactions), diffusive transport, and a radiation model (for CH4, CO, CO2, H2O, and soot). The calculated combustion characteristics depend strongly on the air velocity around the fuel source. In a near-quiescent air environment (< or = 2 mm/s), with a sufficiently large fuel injection velocity (1 cm/s), a growing spherical diffusion flame extinguishes at ≈1200 K due to radiative heat losses. This is typically followed by a transition to the low-temperature (cool-flame) regime with a reaction zone (at ≈700 K) in close proximity to the fuel source. The 'cool flame' regime is formed due to the negative temperature coefficient in the low-temperature chemistry. After a relatively long period (≈18 s) of the cool flame regime, a flash re-ignition occurs, associated with flame-edge propagation and subsequent extinction of the re-ignited flame. In a low-speed (≈3 mm/s) airstream (which simulates the slight droplet movement), the diffusion flame is enhanced upstream and experiences a local extinction downstream at ≈1200 K, followed by steady flame pulsations (≈0.4 Hz). At higher air velocities (4-10 mm/s), the locally extinguished flame becomes steady state. The present axisymmetric computational approach helps in revealing the non-premixed 'cool flame' structure and 2D flame-flow interactions observed in recent microgravity droplet combustion experiments.

diffusion flame↗

Morphological and developmental variation of red mizuna grown under blue light in ground-based microgravity analogues through the utilization of specialized clinostats

Microgreens, such as red (purple) mizuna, have been considered as spaceflight candidate crop plants to be integrated in bioregenerative life support systems for long-duration spaceflight. This is due to their high nutrient density, relatively fast developmental cycles, low horticultural energy, water needs and minimal space requirements. It is therefore important to develop space plant life-support systems capable of mitigating any potentially detrimental effects associated with crop responses to spaceflight related stressors; such as microgravity and light variation. This study focuses on identifying phenotypic variation in the development of red (purple) mizuna grown under blue light in a 2D clinostat to simulate the lack of a gravity vector found within true microgravity.

Gilbert Cauthorn↗

Paramyxovirus Infection Mimics In Vivo Cellular Dynamics in Three-Demensional Human Bronchio-Epithelial Tissue-Like Assemblies

Respiratory syncytial virus and parainfluenza virus cause severe respiratory disease, especially in infants, children and the elderly. An in vitro model that accurately mimics infection of the human respiratory epithelium (HRE) would facilitate vaccine development greatly. Monolayer cultures traditionally used to study these viruses do not accurately and precisely differentiate the replication efficiencies of wild type and attenuated viruses. Therefore, we engineered novel three-dimensional (3D) tissue-like assemblies (TLAs) of human broncho-epithelial (HBE) cells to produce a more physiologically relevant in vitro model of the HRE. TLAs resemble HRE structurally and by expression of differentiated epithelial cell markers. Most significantly, wild type viruses exhibited a clear growth advantage over attenuated strains in TLAs unlike monolayer cultures. In addition, the TLAs responded to virus infection by secreting pro-inflammatory mediators similar to the respiratory epithelia of infected children. These characteristics make the TLA model a valuable platform technology to develop and evaluate live, attenuated respiratory virus vaccine candidates for human use. Respiratory virus diseases, the most frequent and least preventable of all infectious diseases, range in severity from the common cold to severe bronchiolitis and pneumonia . Two paramyxoviruses, respiratory syncytial virus (RSV) and parainfluenza virus type 3 (PIV3), are responsible for a majority of the most severe respiratory diseases of infants and young children. RSV causes 70% of all bronchiolitis cases and is a major cause of morbidity and mortality worldwide, especially in infants. PIV3 causes 10-15% of bronchiolitis and pneumonia during infancy, second only to RSV, and 40% of croup in infants To date, licensed vaccines are not available to prevent these respiratory diseases. At present, traditional monkey kidney (Vero and LLC-MK2) and human (HEp-2) tissue culture cells and small animal models (mouse, cotton rat, guinea pig, ferret, and hamster) fail to accurately imitate viral replication and human disease states (8). Lacking an authentic model has impeded the development and evaluation of live, attenuated vaccine candidates. Development of a physiologically relevant in vitro tissue culture model that reproduces characteristics of the HRE, the primary target of RSV and PIV3, would aid in predicting clinical attenuation and safety of vaccine candidates. Successful tissue engineering of a 3D human intestinal model using novel NASA technology inspired the development of a tri-culture 3D model for the HRE. Sequential layering of primary mesenchymal cells (comprised of normal human fibroblasts and endothelial cells) followed by BEAS-2B epithelial cells derived from human bronchi and tracheae were recapitulated on Cultisphere and/or cytodex3 microcarriers in cylindrical vessels that rotate horizontally creating an organized epithelial structure. Horizontal rotation randomizes the gravity vector modeling aspects of microgravity. Mesenchymal and epithelial cells grown under these conditions reproduce the structural organization, multi-cellular complexity, and differentiation state of the HRE. The opportunity to study respiratory viruses in a nasal epithelium model is invaluable because the most promising respiratory virus vaccine candidates are live attenuated viruses for intranasal administration. Here we characterize the interactions of respiratory viruses and epithelial cells grown under modeled microgravity in comparison to gravity-ladened monolayers. 3D HBE TLAs and traditional monolayers (2D) are infected at 35 C, the upper temperature of the upper HRE, to simulate in vivo infection conditions. Growth kinetics of wild type (wt) RSV and PIV3 viruses were compared in 2D and 3D cells to that of strains attenuated in humans or rhesus macaques. This novel 3D HBE model also offers an opportunity to study whether the epithelial cell function, especially in host defenses recapitulated by mimicking the structural organization of the HRE. In vivo, airway epithelial cells play a significant and dynamic role in host defense by blocking paracellular permeability and modulating airway function through cellular interactions or tight junctions. As regulators of the innate immune response, epithelial cells constitutively express cytokines, chemokines, and colony stimulating factors including RANTES, IL-8, IL-6, GM-CSF, and G-CSF for proactive host defense. In response to viral infection, epithelial cells induce potent immuno-modulatory and pro-inflammatory cytokines that recruit phagocytic and inflammatory cells to clear the virus and enhance protection. Although disease pathogenesis is classically attributed to the cytopathic effects of the pathogen, severe disease states associated with RSV and PIV3 are attributed to the inflammatory response, especially in infants. RSV is a potent inducer of cytokines and pro-inflammatory mediators in epithelial cells in vivo. A differentiated human epithelial model independent of the complete functional immune system will help elucidate the role of epithelial cells in respiratory disease. We reported here, virus and host cell interactions in 3D HBE TLAs are similar to that in vivo. Because the epithelial cell organization of the TLAs impacts not only the expression of airway epithelial characteristics, but also cellular communication, the TLAs represent a more physiologically relevant model of the HRE than BEAS-2B or other non-tumour monolayer models of respiratory disease. As a result, wild type respiratory viruses have a clear growth advantage over attenuated viruses in TLAs unlike traditional monolayers. In addition, the TLAs respond to wild type virus infection by secreting pro-inflammatory mediators characteristic of infected HRE. TLAs expressing microbial defense mechanisms provide an excellent model to study the interactions of respiratory pathogens with their host and to identify the innate immunity mediators. Therefore, 3D HBE TLAs offer advantages for the study of respiratory viruses and the development of viral vaccine candidates.

Deatly, Anne M.↗

Modeling the Dynamics of the Middle Atmosphere and Lower Thermosphere Under the Influence of Gravity Waves

Our Numerical Spectral Model (NSM), which extends from the ground up into the thermosphere, is non-linear, time-dependent and has been employed for 2D and 3D applications. The standard version of the NSM incorporates Hines' Doppler Spread Parameterization for small scale gravity waves (GW), but planetary waves generated in the troposphere have also been incorporated. The NSM has been applied to describe: (1) the anomalous seasonal variations of the zonal circulation and temperature in the upper mesosphere, (2) the equatorial oscillations (quasi-biennial and semi-annual oscillations (QBO and SAO)) extending from the stratosphere into the upper mesosphere, (3) the diurnal and semi-diurnal tides, and (4) the planetary waves that are excited in the mesosphere. With the emphasis to provide understanding, we present here results from numerical experiments with the NSM that shed light on the GW processes that are of central importance in the mesosphere and lower thermosphere. These are our conclusions: (1) The large semiannual variations in the diurnal tide (DT), with peak amplitudes observed around equinox, are produced primarily by GW interactions that involve, in part, planetary waves. The DT, like planetary waves, tends to be amplified by GW momentum deposition, which reduces also the vertical wavelength, but variations in eddy viscosity associated with GW interactions are also important. (2) The semidiurnal tide (SDT) and its phase in particular, is strongly influenced by the mean zonal circulation. The SDT, individually, is also amplified by GW. But the DT filters out GW such that the GW interaction effectively reduces the amplitude of the SDT, producing a strong nonlinear interaction between the DT and SDT. (3) Without external time dependent energy or momentum sources, planetary waves (PW) are generated in the model for zonal wavenumbers 1 to 4, which have amplitudes in the mesosphere above 50 km as large as 40 m/s and periods between 50 and 2 days. The waves are generated primarily during solstice conditions, which indicates that the baroclinic instability (associated with the GW induced reversal in the latitudinal temperature gradient) is playing an important role. Numerical experiment show that GW, directly, also greatly amplify the PW. A common feature of the PW generated in summer and winter is that their vertical wavelengths throughout the mesosphere are large, which indicates that the waves are not propagating freely but are generated throughout the region. Another common feature is that the PW propagate preferentially westward in summer and eastward in winter, being launched from the westward and eastward zonal winds that prevail respectively in summer and winter at altitudes below 80 km. (4) Planetary waves generated internally by baroclinic instability and GW interaction produce large amplitude modulations of the DT and SDT. In summary we conclude that GW play major roles in generating and amplifying the dynamical components in the MLT region and, acting principally through wave filtering, produce important non-linear interactions between the components.

Mayr, H. G.↗

Vibration control of a manipulator tip on a flexible body

Vibration control of a rigid manipulator tip on a main flexible uniform beam is examined. It is proposed to add a compensator between the manipulator and the beam to rotate and extend/retrieve the manipulator during the control period. The 2D station-keeping maneuvers within the linear range without gravity and damping are considered. The compensatory open-loop control law, which depends on the amplitudes of the beam's flexible deformations at the connection joint, is synthesized using linear quadratic regulator techniques. After introducing the compensatory control into the system, system control is still stable, and the tip coordinates of the manipulator can be made to closely follow the rigid beam motion, which is assumed to be a desired motion.

Xu, J.↗

NASA’s Ground-Based Microgravity Simulation Facility

Since opportunities to conduct experiments in space are scarce, various microgravity simulators and analogs have been widely used in space biology ground studies. Even though microgravity simulators do not produce all of the biological effects observed in the true microgravity environment, they provide alternative test platforms that are effective, affordable, and readily available to facilitate microgravity research. The Microgravity Simulation Support Facility (MSSF) at the National Aeronautics and Space Administration (NASA) John F. Kennedy Space Center (KSC) has been established for conducting short duration experiments, typically less than 1 month, utilizing a variety of microgravity simulation devices for research at different gravity levels. The simulators include, but are not limited to, 2D Clinostats, 3D Clinostats, Random Positioning Machines, and Rotating Wall Vessels. In this presentation, we will provide an overview of current MSSF capabilities, research conducted in MSSF, new technology developments, customized sample holder designs, and the physical characteristics of these microgravity simulators.

Ye Zhang↗

NASA's Ground-Based Microgravity Simulation Facility

Since opportunities to conduct experiments in space are scarce, various microgravity simulators and analogs have been widely used in space biology ground studies. Even though microgravity simulators do not produce all of the biological effects observed in the true microgravity environment, they provide alternative test platforms that are effective, affordable, and readily available to facilitate microgravity research. The Microgravity Simulation Support Facility (MSSF) at the National Aeronautics and Space Administration (NASA) John F. Kennedy Space Center (KSC) has been established for conducting short duration experiments, typically less than 1 month, utilizing a variety of microgravity simulation devices for research at different gravity levels. The simulators include, but are not limited to, 2D Clinostats, 3D Clinostats, Random Positioning Machines, and Rotating Wall Vessels. In this chapter, we will discuss current MSSF capabilities, development concepts, customized sample holder designs, and the physical characteristics of these microgravity simulators.

Ye Zhang↗

Comparison of marine gravity from shipboard and high-density satellite altimetry along the Mid-Atlantic Ridge, 30.5-35.5 deg S

We compare new marine gravity fields derived from satellite altimetry with shipboard measurements over a region of more than 120,000 sq km in the central South Atlantic. Newly declassified satellite data were employed to construct free-air anomaly maps on 0.05 degree grids. An extensive gravity and bathymetry data set from four cruises along the Mid-Atlantic Ridge from 30.5-35.5 deg S provides a benchmark for testing the 2D resolution and accuracy of the satellite measurements where their crosstrack spacing is near their widest. The satellite gravity signal is coherent with bathymetry in this region down to wavelengths of 26 km, compared to 12.5 km for shipboard gravity. Residuals between the shipboard and satellite data sets have a roughly normal distribution. The standard deviation of satellite gravity with respect to shipboard measurements is nearly 7 mGal in a region of 140 mGal total variation, whereas the internal standard deviation at crossovers for GPS-navigated shipboard data is 1.8 mGal. The differences between shipboard and satellite data are too large to use satellite gravity to determine crustal thickness variations within a typical ridge segment.

Neumann, Gregory A.↗

Effects of buoyancy-driven flow and thermal boundary conditions on physical vapor transport

A 2D numerical model was developed in order to ascertain if reduced gravity conditions are beneficial to physical vapor transport (PVT) and to determine its tolerance limits to residual accelerations. This was solved using the PHOENICS finite-volume code. Reduction of gravitational accelerations to less than 0.1 g0 was found to be sufficient to suppress buoyancy-driven convection to an extent that diffusion was the dominant transport mode, whence a greater uniformity in the growth rate could be obtained. It is shown that a uniform temperature gradient on the ampoule walls causes the vapor to be supersaturated throughout the ampoule, potentially resulting in undesirable nucleation at the walls. A 'hump' in the wall temperature profile can be used to avoid this. The prevailing transport conditions determine the size of the hump needed.

Nadarajah, Arunan↗

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↗

A novel transformation of the ice sheet Stokes equations and some of its properties and applications

We introduce a novel transformation of the Stokes equations into a form closely resembling the shallow Blatter–Pattyn equations. The two forms differ by only a few additional terms, while their variational formulations differ only by a single term in each horizontal direction. Specifically, the variational formulation of the Blatter–Pattyn model drops the vertical velocity in the second invariant of the strain rate tensor. Here we make use of the new transformation in two ways. First, we consider incorporating the transformed equations into a code that can be very easily converted from a Stokes to a Blatter–Pattyn model, and vice versa, by switching these terms on or off. This may be generalized so that the Stokes model is switched on adaptively only where the Blatter–Pattyn model loses accuracy. Second, the key role played by the vertical velocity in the Blatter–Pattyn approximation motivates new approximations. Two examples are presented. These require a mesh that enables the discrete continuity equation to be invertible for the vertical velocity in terms of the horizontal velocity components. Examples of such meshes, such as the first-order P1–E0 mesh and the second-order P2–E1 mesh, are given in both 2D and 3D. However, the transformed Stokes model has the same type of gravity forcing as the Blatter–Pattyn model, determined by the ice surface slope, thereby forgoing some of the mesh generality of the traditional formulation of the Stokes model.

58 GEOSCIENCES↗