Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic structure factor”

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 181 records · Page 10

Modeling structural damping for solids having distinct shear and dilatational loss factors

For steady state time harmonic problems (rigid format 8), the NASTRAN program as currently configured treats internal structural damping through the introduction of a single structural element damping coefficient that typically is viewed as the ratio of the complex to real modulus of elasticity. For problems dealing with two or three dimensional dynamic linear viscoelasticity (e.g. a Kelvin-Voigt viscoelastic model), the present NASTRAN capability cannot directly handle this situation wherein two independent damping coefficients are required to properly model the dissipation phenomenon. A technique is presented whereby the user can adapt the standard versions of NASTRAN (without resorting to either DMAP and/or FORTRAN coding changes) for the purpose of treating this class of problem.

Kalinowski, A. J.↗

Structural-Vibration-Response Data Analysis

Computer program developed as structural-vibration-response data analysis tool for use in dynamic testing of Space Shuttle. Program provides fast and efficient time-domain least-squares curve-fitting procedure for reducing transient response data to obtain structural model frequencies and dampings from free-decay records. Procedure simultaneously identifies frequencies, damping values, and participation factors for noisy multiple-response records.

Smith, W. R.↗

Advanced Silicon Avalanche Photodiodes on NASA's Global Ecosystem Dynamics Investigation (GEDI) Mission

Silicon Avalanche Photodiodes (APDs) are used in NASA’s Global Ecosystem Dynamics Investigation (GEDI) which was launched in December 2018 and is currently measuring the Earth’s vegetation vertical structure from the International Space Station. The APDs were specially made for space lidar with a much lower hole-to-electron ionization coefficient ratio (k-factor ~0.008) than that of commercially available silicon APDs in order to reduce the APD excess noise from the randomness of the avalanche gain. A silicon heater resistor was used under the APD chip to heat the device up to 70°C and improve its quantum efficiency at 1064 nm laser wavelength while maintaining a low dark current such that the overall signal to noise ratio is improved. Special APD protection circuits were used to raise the overload damage threshold to prevent device damage from strong laser return by specular surfaces, such as still water bodies, and space radiation events. The APD and a hybrid transimpedance amplifier circuit were hermetically sealed in a package with a sufficiently low leak rate to ensure multi-year operation lifetime in space. The detector assemblies underwent a series of pre-launch tests per NASA Goddard Environmental Verification Standard for space qualification. They have performed exactly as expected with GEDI in orbit. A detailed description of the GEDI detector design, signal and noise model, and test results are presented in this paper.

Sun, Xiaoli↗

Analysis of Heat Transfer from Local Heating and Cooling Sources at Cryogenic Temperatures

Understanding the dispersion of heat around a cryogenic fluid tank, specifically the interaction between the cryogenic fluid and the tank wall is critical in the analysis of long duration cryogen storage in microgravity. The heat transfer interaction between a cryogenic storage tank and heat sources from external spacecraft structures is also one of the many factors that determine how much heat enters a tank. Recent flight experiments with two-phase fluids have indicated that local concentrations of heat input (also known as “hot spots”) can cause unwanted affects including local boiling. Computational fluid dynamic (CFD) models can provide a detailed assessment of the heat transfer occurring across a cryogenic storage system. However, CFD modeling takes time to construct and run. A simpler approach that can act as initial guidance for later CFD modeling analyzes external “hot spots” as point or finite heat sources. A radial, finite element network or a local direct solution can effectively estimate the heat spread across a cryogenic storage tank by calculating the temperature and heat load as a function of distance from the heat source. This calculation accounts for the convective heat transfer between the cryogenic fluid and storage tank surface. Similar approaches can be used to determine the effectiveness of cooling from a cryocooler as a finite, local heat sink. This approach allows for quick approximations of the thermal map across a cryogenic tank as well as sensitivity analysis under a wide range of design parameters including gravitational fields as implied through natural convection coefficients.

Cryogenic Fluid Management↗

Impact of Meteorological Factors on the Mesoscale Morphology of Cloud Streets during a Cold-Air Outbreak over the Western North Atlantic

Postfrontal clouds (PFC) are ubiquitous in the marine boundary layer, and their morphology is essential to estimating the radiation budget in weather and climate models. Here we examine the roles of sea surface temperature (SST) and meteorological factors in controlling the mesoscale morphology and evolution of shallow clouds associated with a cold-air outbreak that occurred on 1 March 2020 during phase I of the Aerosol Cloud Meteorology Interactions over the Western Atlantic Experiment (ACTIVATE). Our results show that the simulated PFC structure and ambient conditions by the Weather Research and Forecasting (WRF) Model are generally consistent with observations from GOES-16 and dropsonde measurements. We also examine the thermodynamical and dynamical influences in the cloud mesoscale morphology using WRF sensitivity experiments driven by two meteorological forcing datasets with different domain-mean SST and spatial gradients, which lead to dissimilar values of hydrometeor water path and cloud core fraction. The SST from ERA5 leads to weaker stability and higher inversion height than the SST from FNL does. In addition, the use of large-scale meteorological forcings from ERA5 yields a distinctive time evolution of wind direction shear in the inner domain, which favors the formation and persistence of longer cloud rolls. Both factors contribute to a change in the time evolution of domain-mean water path and cloud core fraction of cloud streets. Our study takes advantage of the simulation driven by the differences between two large-scale forcing datasets to illustrate the importance of SST and wind direction shear in the cloud street morphology in a realistic scenario

Jingyi Chen↗

A predictive model of human performance.

An attempt is made to develop a model describing the overall responses of humans to exercise and environmental stresses for prediction of exhaustion vs an individual's physical characteristics. The principal components of the model are a steady state description of circulation and a dynamic description of thermal regulation. The circulatory portion of the system accepts changes in work load and oxygen pressure, while the thermal portion is influenced by external factors of ambient temperature, humidity and air movement, affecting skin blood flow. The operation of the model is discussed and its structural details are given.

Walters, R. F.↗

Aeroelastic modeling for the FIT (Functional Integration Technology) team F/A-18 simulation

As part of Langley Research Center's commitment to developing multidisciplinary integration methods to improve aerospace systems, the Functional Integration Technology (FIT) team was established to perform dynamics integration research using an existing aircraft configuration, the F/A-18. An essential part of this effort has been the development of a comprehensive simulation modeling capability that includes structural, control, and propulsion dynamics as well as steady and unsteady aerodynamics. The structural and unsteady aerodynamics contributions come from an aeroelastic mode. Some details of the aeroelastic modeling done for the Functional Integration Technology (FIT) team research are presented. Particular attention is given to work done in the area of correction factors to unsteady aerodynamics data.

Zeiler, Thomas A.↗

Tortuosity Computations of Porous Materials using the Direct Simulation Monte Carlo

Low-density carbon fiber preforms, used as thermal protection systems (TPS) materials for planetary entry systems, have permeable, highly porous microstructures consisting of interlaced fibers. Internal gas transport in TPS is important in modeling the penetration of hot boundary-layer gases and the in-depth transport of pyrolysis and ablation products. The gas effective diffusion coefficient of a porous material must be known before the gas transport can be modeled in material response solvers; however, there are very little available data for rigid fibrous insulators used in heritage TPS.The tortuosity factor, which reflects the efficiency of the percolation paths, can be computed from the effective diffusion coefficient of a gas inside a porous material and is based on the micro-structure of the material. It is well known, that the tortuosity factor is a strong function of the Knudsen number. Due to the small characteristic scales of porous media used in TPS applications (typical pore size of the order of 50 micron), the transport of gases can occur in the rarefied and transitional regimes, at Knudsen numbers above 1. A proper way to model the gas dynamics at these conditions consists in solving the Boltzmann equation using particle-based methods that account for movement and collisions of atoms and molecules.In this work we adopt, for the first time, the Direct Simulation Monte Carlo (DSMC) method to compute the tortuosity factor of fibrous media in the rarefied regime. To enable realistic simulations of the actual transport of gases in the porous medium, digitized computational grids are obtained from X-ray micro-tomography imaging of real TPS materials. The SPARTA DSMC solver is used for simulations. Effective diffusion coefficients and tortuosity factors are obtained by computing the mean-square displacement of diffusing particles.We first apply the method to compute the tortuosity factors as a function of the Knudsen number for computationally designed materials such as random cylindrical fibers and packed bed of spheres with prescribed porosity. Results are compared to literature values obtained using random walk methods in the rarefied and transitional regime and a finite-volume method for the continuum regime. We then compute tortuosity factors for a real carbon fiber material with a transverse isotropic structure (FiberForm), quantifying differences between through-thickness and in-plain tortuosities at various Knudsen regimes.

Tortuosity↗

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.↗

Development of Buffet Functions using Frequency-Dependent Coherence Factors

Launch vehicle transonic buffet environments produce large fluctuating pressures on the surface of the vehicle. Transonic buffet occurs as a launch vehicle passes through the transonic regime, typically Mach 0.8 to 1.2, where phenomena such as flow separation and shock dynamics produce large fluctuating pressures. The unsteady aerodynamic forces produced by these phenomena can excite both global and local vehicle structural response. The consequences of incorporating poorly characterized buffet forces into the design process can range from inefficient structures that can reduce the vehicle capabilities to undersized structures that, in the worst case scenario, may result in structural failure (Refs. 1, 2, 3). It is thus imperative to appropriately characterize the transonic buffet environment.

buffet↗

Multidisciplinary Dynamic Testing Challenges in Validating the NASA Artemis Architecture

NASA is in the midst of bold and exciting next steps in human exploration and spaceflight. The designs of the new Space Launch System (SLS), the Orion spacecraft and the Exploration Ground Systems (EGS) for vehicle processing and launch are essentially complete and there has been significant progress in manufacturing and assembly of specific hardware for the Artemis I and Artemis II missions. Equally as important, the program level and integrated system level testing and analyses are also well underway to support integrated verification, validation, and certificate of flight readiness (CoFR) for the first Artemis mission. Testing and analysis are key to addressing technical challenges that the Artemis missions offer. Building block approaches are required that provide the right balance between component, system, and/or element level testing that satisfies verification and validation objectives and where, uncertainties are quantified and minimized. Artemis I is a system of systems that requires a fusion of test and analysis that adeptly characterizes critical interfaces between major program elements. NASA is implementing new in-situ testing that fuse traditional aerospace structures with civil structures, such as the Integrated Modal Test for the Artemis I vehicle where the Mobile Launcher and Crawler Transporter serve as a support structure whose dynamics couple with that of the Artemis I vehicle. This new paradigm requires a closer inspection of structural behavior of the Crawler Transporter and the Mobile Launcher as they now serve multiple purposes. This requires a paradigm shift to look beyond experimental modal techniques and incorporates operational modal analysis techniques to validate dynamic models from data collected during rollout to the launch pad. A further complicating factor is the Crawler Transporter generated ground forces have numerous harmonics making extracting dynamic responses of the Artemis I, Mobile Launcher, and Crawler Transporter coupled system challenging. This discussion explores all these challenges with and attempts to understand how we best build confidence in systems and system-of-systems performance capabilities and margins and understand uncertainties.

Joel W Sills↗

Integrated Simulation Design Challenges to Support TPS Repair Operations

During the Orbiter Repair Maneuver (ORM) operations planned for Return to Flight (RTF), the Shuttle Remote Manipulator System (SRMS) must grapple the International Space Station (ISS), undock the Orbiter, maneuver it through a long duration trajectory, and orient it to an EVA crewman poised at the end of the Space Station Remote Manipulator System (SSRMS) to facilitate the repair of the Thermal Protection System (TPS). Once repair has been completed and confirmed, then the SRMS proceeds back through the trajectory to dock the Orbiter to the Orbiter Docking System. In order to support analysis of the complex dynamic interactions of the integrated system formed by the Orbiter, ISS, SRMS, and SSMS during the ORM, simulation tools used for previous ‘nominal’ mission support required substantial enhancements. These upgrades were necessary to provide analysts with the capabilities needed to study integrated system performance. Prevalent throughout this ORM operation is a dynamically varying topology. In other words, the ORM starts with the SRMS grappled to the mated Shuttle/ISS stack (closed loop topology), moves to an open loop chain topology consisting of the Shuttle, SRMS, and ISS, and then, at the repair configuration, extends the chain topology to one consisting of the Shuttle, SMS, ISS, and SSRMS/EVA crewman. The resulting long dynamic chain of vehicles and manipulators may exhibit significant motion between the Shuttle worksite and the EVA crewman due to the system flexibility throughout the topology (particularly within the SRMS/SSRMS joints and links). Since the attachment points of both manipulators span the flexible structure of the ISS, simulation analysis may also need to take that into consideration. Moreover, due to the lengthy time duration associated with the maneuver and repair, orbital effects become a factor and require the ISS vehicle control system to maintain active attitude control. Several facets of the ORM operation make the associated analytical efforts different from previous mission support, including: (1) the magnitude of the SRMS handled payload (Le., Orbiter class), (2) the orbital effects induced on the integrated system consisting of the large Shuttle and ISS masses connected by a light flexible SRMS, (3) long duration environmental consequences due to the lengthy operational times associated with the maneuver and repair of the TPS, (4) active attitude control (as opposed to free drift) interacting with the SRMS and SSRMS manipulators (also due to the length of the maneuver and repair), (5) relative dynamics between the EVA crewman and thc worksite influenced by the extended flexible topology. In order to meet these analysis challenges, an ORM simulation architecture was developed leveraging upon numerous pre-existing simulation elements to analyze the various subsystems individually. For example, core manipulator subsystem simulations for both the SRMS and SSRMS were originally combined to provide the dual-arm dynamics topology simulation (in the absence of orbital dynamics and vehicle control). This capability was later merged with the simulation used to analyze SRMS loading with a heavy payload in the orbital environment with an active payload control system (in this case, the ISS Attitude Control System (ACS)), configured for the ORM. The resulting worksite dynamics simulation, based off of the modified ORM simulation, provided the extended topological chain of vehicles and manipulators, while taking into account the orbital effects of both the Shuttle and ISS (as well as its ACS). Verification and validation (V&V) of these integrated simulations became a challenge in itself. A systematic approach needed to be developed such that integration simulation results could be tested against previous constituent simulations upon which these simulations were built. General V&V categories included: (1) core orbital state propagation, (2), stand-alone SRMS, (3) stand-alone SSRMS, (4) stand-alone ISS ACS, (5) integrated Shuttle, SRMS, ISS (with active ACS) in the orbital environment, and (5) dual-arm SRMS/SSRMS dynamics topology. Integrated simulation V&V run suites were created and correlated to verification runs from subsystem simulations, in order to establish the validity of the results. This paper discusses the simulation design challenges encountered while developing simulation capabilities to mirror the ORM operations. The paper also describes the incremental build approach that was utilized, starting with the subsystem simulation elements and integration into increasing more complex simulations until the resulting ORM worksite dynamics simulation had been assembled. Furthermore, the paper presents an overall integrated simulation V&V methodology based upon a subsystem level testing, integrated comparisons, and phased checkout.

Thermal Protection System↗

A new state space model for the NASA/JPL 70-meter antenna servo controls

A control axis referenced model of the NASA/JPL 70-m antenna structure is combined with the dynamic equations of servo components to produce a comprehansive state variable (matrix) model of the coupled system. An interactive Fortran program for generating the linear system model and computing its salient parameters is described. Results are produced in a state variable, block diagram, and in factored transfer function forms to facilitate design and analysis by classical as well as modern control methods.

Hill, R. E.↗

Molecular Dynamics study of Pb overlayer on Cu(100)

Isothermal-isobaric Molecular Dynamics (MD) simulation of a submonolayer Pb film in c(2x2) ordered structure adsorbed on a Cu(100) substrate showed retention of order to high T. The Embedded Atom Method (EAM) calculated the energy of atoms of overlayer and substrate. The time-averaged squared modulus of the two dimensional structure factor for the Pb overlayer measured the order of the overlayer. The results are for increasing T only, and require verification by simulated cooling.

Karimi, M.↗

3D Representation of UAV-obstacle Collision Risk Under Off-nominal Conditions

Safe operations of autonomous unmanned aerial vehicles (UAVs) in low-altitude airspace with beyond visual line-of-sight (BVLOS) flights demand robust risk monitoring of airspace as well as of people and property on ground. One of the safety critical factors for UAV flights is the risk of collision with static and dynamic obstacles in proximity to its flight path. This paper presents a detailed formulation of risk of obstacle collision incorporating the effects of off-nominal conditions introduced by component failures, degraded controllability and environmental disturbances such as wind gusts. The risk is represented in terms of a matrix with rows corresponding to the likelihood of occurrence of collision and columns representing severity of collision to the vehicle and surrounding structures. Risk likelihood is generated using a Bayesian Belief Network (BBN) that compiles knowledge from related Failure Modes and Effects Analysis (FMEAs) and Subject Matter Experts (SMEs) to determine the probability of collision based on on-board sensor measurements indicative of vehicle health and controllability. Risk severity is computed utilizing a point-mass 3D kinematic model of the vehicle in presence of wind. The proposed risk factor is demonstrated on real flight data from experimental flights of an octocopter at NASA Langley Research Center in presence of simulated obstacles and wind conditions. Effect of varying wind conditions, level of controllability and obstacle measurement noise on the risk factor is demonstrated. The proposed approach enables risk-informed decision making for timely mitigation of current and future unsafe events in autonomous systems.

risk analysis↗

Gas dynamics of semidetached binaries. II - The vertical structure of the stream

The letter considers the three-dimensional dynamics of the stream in semidetached binaries undergoing Roche-lobe overflow. Results of an earlier communication are generalized to include the dynamic effects in the direction perpendicular to the orbital plane. It is found that the scale height of the stream characteristically exceeds its corresponding hydrostatic value by a significant factor because the inertia of the gas prevents it from responding instantaneously to the changing gravitational field. This effect is important for the interpretation of observations relating to stream-disk impacts in cataclysmic variables and in binary X-ray sources of low total mass.

Lubow, S. H.↗

Environmental and high strain rate effects on composites for engine applications

The complex environmental and loading conditions experienced by many turbine engine components impose severe durability and damage tolerance requirements for component materials. Programs are being conducted to establish the structural performance of composite materials under anticipated engine operating environments. A description is presented of the results obtained in connection with several of these programs. A comparison of predicted and measured hygrothermal effects is considered along with hygrothermal effects on defect growth, high strain rate effects on composite mechanical properties, dynamic stress intensity factors for composites, and the indentation laws for composite impact analysis.

Chamis, C. C.↗

Plume mechanics and stratocumulus convection

The First ISCCP Regional Experiment (FIRE) marine stratocumulus Intensive Field Observation (IFO) held in July, 1987 produced a data set that is far more comprehensive than data sets from previous stratocumulus experiments. One exciting new development was the use of the 10.6 micron lidar system for cloud-top mapping that was available on the National Center for Atmospheric Research (NCAR) Electra. This system provided a unique look at the small scales of the turbulence in the clouds, images of the turbulent structures that are quantitatively revealed by conditional sampling (e.g., Khalsa and Greenhut, 1987). The behavior of these updrafts and downdrafts is central to the dynamics of the stratocumulus-topped marine boundary layer. FIRE's objectives of understanding cloud dynamics and how they affect the cloud optical depth - which, in turn, is the crucial factor in determining the clouds' albedo - therefore require the investigation of these drafts. Given here are initial results from a simple model capable of simulating moist, entraining plumes that are subject to water phase changes and radiative heating and cooling. These results correspond, therefore, to the 'dry cloud' case discussed by Lilly and Schubert (1980). The model's simplicity limits the realism of the results - the plumes are assumed not to interact, for example - but the role of radiative processes in influencing the plume dynamics is clear. Also revealed is the role the plumes play is maintaining the cloud-top inversion. The model equations and methodology used, and qualitative results are given.

Hanson, Howard P.↗