Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase field simulation”

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 559 records · Page 31

Subspace-based Background Subtraction Applied to Aeroacoustic Wind Tunnel Testing

A subspace-based form of background subtraction is presented and applied to aeroacoustic wind tunnel data. A variant of this method has seen use in other fields such as climatology and medical imaging. The technique is based on an eigenvalue decomposition of the background noise cross-spectral matrix. Simulated results indicate similar performance to conventional background subtraction when the subtracted spectra are weaker than the true contaminating background levels. Superior performance is observed when the subtracted spectra are stronger than the true contaminating background levels, and when background data do not match between measurements. Experimental results show limited success in recovering signal behavior for data in which conventional background subtraction fails. The results also demonstrate the subspace subtraction technique's ability to maintain a physical coherence relationship in the modified cross-spectral matrix. Deconvolution results from microphone phased array data indicate that array integration methods are largely insensitive to subtraction type, and that background subtraction with appropriate background data is an effective alternative to diagonal removal.

deconvolution↗

Final Report for The Creation of a Physics-based Ground-effect Model, Phase 2 - Inclusion of the Effects of Wind, Stratification, and Shear into the New Ground Effect Model

The reduction of the separation of the leading and following aircrafts is desirable to enhance the airport capacity provided that there is a physics-based operational model applicable to all regions of the flight domain (out of ground effect, OGE; near ground effect, NGE; and in ground effect, IGE) and that the quality of the quantitative input from the measurements of the prevailing atmospheric conditions and the quality of the total airport operations regarding the safety and the sound interpretation of the prevailing conditions match the quality of the analysis and numerical simulations. In the absence of an analytical solution, the physics of the flow is best expressed by a mathematical model based on numerical simulations, field and laboratory experiments, and heuristic reasoning. This report deals with the creation of a sound physics-based real-time IGE model of the aircraft wake vortices subjected to crosswind, stratification and shear.

Sarpkaya, Turgut↗

The Life Cycle Application of Intelligent Software Modeling for the First Materials Science Research Rack

Marshall Space Flight Center (MSFC) has been funding development of intelligent software models to benefit payload ground operations for nearly a decade. Experience gained from simulator development and real-time monitoring and control is being applied to engineering design, testing, and operation of the First Material Science Research Rack (MSRR-1). MSRR-1 is the first rack in a suite of three racks comprising the Materials Science Research Facility (MSRF) which will operate on the International Space Station (ISS). The MSRF will accommodate advanced microgravity investigations in areas such as the fields of solidification of metals and alloys, thermo-physical properties of polymers, crystal growth studies of semiconductor materials, and research in ceramics and glasses. The MSRR-1 is a joint venture between NASA and the European Space Agency (ESA) to study the behavior of different materials during high temperature processing in a low gravity environment. The planned MSRR-1 mission duration is five (5) years on-orbit and the total design life is ten (IO) years. The MSRR-1 launch is scheduled on the third Utilization Flight (UF-3) to ISS, currently in February of 2003). The objective of MSRR-1 is to provide an early capability on the ISS to conduct material science, materials technology, and space product research investigations in microgravity. It will provide a modular, multi-user facility for microgravity research in materials crystal growth and solidification. An intelligent software model of MSRR-1 is under development and will serve multiple purposes to support the engineering analysis, testing, training, and operational phases of the MSRR-1 life cycle development. The G2 real-time expert system software environment developed by Gensym Corporation was selected as the intelligent system shell for this development work based on past experience gained and the effectiveness of the programming environment. Our approach of multi- uses of the simulation model and its intuitive graphics capabilities is providing a concurrent engineering environment for rapid prototyping and development. Operational schematics of the MSRR-1 electrical, thermal control, vacuum access, and gas supply systems, and furnace inserts are represented graphically in the environment. Logic to represent first order engineering calculations is coded into the knowledge base to simulate the operational behavior of the MSRR-1 systems. An example of engineering data provided includes electrical currents, voltages, operational power, temperatures, thermal fluid flow rates. pressures, and component status indications. These type of data are calculated and displayed at appropriate instrumentation points, and the schematics are animated to reflect the simulated operational status of the MSRR-1. The software control functions are also simulated to represent appropriate operational behavior based on automated control and response to commands received by the crew or ground controllers. The first benefit of this simulation environment is being realized in the high fidelity engineering analysis results from the electrical power system G2 model. Secondly, the MSRR-1 simulation model will be embedded with a hardware mock-up of the MSRR-1 to provide crew training on MSRR-1 integrated payload operations. G2 gateway code will output the simulated instrumentation values, termed as telemetry, in a flight-like data stream so that the crew has realistic and accurate simulated MSRR-1 data on the flight displays which will be designed for crew use. The simulation will also respond appropriately to crew or ground initiated commands, which will be part of normal facility operations. A third use of the G2 model is being planned; the MSRR-1 simulation will be integrated with additional software code as part of the test configuration of the primary onboard computer, or Master Controller, for MSRR-1. We will take advantage of the G2 capability to simulate the flight like data stream to test flight software responses and behavior. A fourth use of the G2 model will be to train the Ground Support Personnel that will monitor the MSRR-1 systems and payloads while they are operating aboard the ISS. The intuitive, schematic based environment will provide an excellent foundation for personnel to understand the integrated configuration and operation of the MSRR-1, and the anticipated telemetry feedback based on operational modes of the equipment. Expert monitoring features will be enhanced to provide a smart monitoring environment for the operators. These features include: (1) Animated, intuitive schematic-based displays which reflect telemetry values, (1) Real-time plotting of simulated or incoming sensor values, (3) High/Low exception monitoring for analog data, (4) Expected state monitoring for discrete data, (5) Data trending, (6) Automated malfunction procedure execution to diagnose problems, (7) Look ahead capability to planned MSRR-1 activities in the onboard timeline. And finally, the logic to calculate telemetry values will be deactivated, and the same environment will interface to the incoming data for the real-time telemetry stream to schematically represent the onboard hardware configuration. G2 will be the foundation for the real-time monitoring and control environment. In summary, our MSRR-1 simulation model spans many elements of the life cycle development of this project: Engineering Analysis, Test and Checkout, Training of Crew and Ground Personnel, and Real-time monitoring and control. By utilizing the unique features afforded by an expert system development environment, we have been able to synergize a powerful tool capable of addressing our project needs at every phase of project development.

Rice, Amanda↗

Computer simulations of two-pickup-ion instabilities in a cometary environment

This paper presents a one-dimensional electromagnetic hybrid computer simulation of distant cometary environments, in which protons and oxygen ions are injected continuously at the same time. It is demonstrated that, when, instead of a single component, both protons and oxygen ions are injected, the faster growing shorter-wavelength proton-proton right-hand resonant instability dominates the early time physics, leading to linear temporal growth of the fluctuating magnetic field energy density and pitch angle scattering of the protons (but causing no significant response by the nonresonant oxygen ions). At later times, as the injected oxygen ions accumulate sufficient density, the oxygen-proton right-hand resonant instability emerges, exhibiting exponential temporal growth to a relatively large peak amplitude, phase bunching, and magnetically trapping the oxygen ions. At times after peak amplitude, the oxygen-proton fluctuations stochastically accelerate some of the injected protons.

Gary, S. Peter↗

Thermal Vacuum Test Performance of the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) Variable Conductance Heat Pipe Assembly

The Variable Conductance Heat Pipe (VCHP) Assembly of the HST Wide Field Camera 3 was subjected to thermal vacuum (TN) environmental testing. The test program included both maximum and minimum environments as well as simulated on-orbit cycling. Elements of the VCHP assembly included a VCHP, an optical bench cold plate with an imbedded constant conductance heat pipe, and a VCHP reservoir radiator with a proportionally controlled heater. The purpose of the test was to characterize and demonstrate the assembly s ability to control the temperature of the cold plate, which provides a stable thermal environment for the instrument s optical bench. This paper discusses the VCHP Assembly control performance and control authority during the dynamic hot and cold 90-minute orbit cycling test phases.

Cleveland, Paul E.↗

Implementation of Complex Signal Processing Algorithms for Position-Sensitive Microcalorimeters

We have recently reported on a theoretical digital signal-processing algorithm for improved energy and position resolution in position-sensitive, transition-edge sensor (POST) X-ray detectors [Smith et al., Nucl, lnstr and Meth. A 556 (2006) 2371. PoST's consists of one or more transition-edge sensors (TES's) on a large continuous or pixellated X-ray absorber and are under development as an alternative to arrays of single pixel TES's. PoST's provide a means to increase the field-of-view for the fewest number of read-out channels. In this contribution we extend the theoretical correlated energy position optimal filter (CEPOF) algorithm (originally developed for 2-TES continuous absorber PoST's) to investigate the practical implementation on multi-pixel single TES PoST's or Hydras. We use numerically simulated data for a nine absorber device, which includes realistic detector noise, to demonstrate an iterative scheme that enables convergence on the correct photon absorption position and energy without any a priori assumptions. The position sensitivity of the CEPOF implemented on simulated data agrees very well with the theoretically predicted resolution. We discuss practical issues such as the impact of random arrival phase of the measured data on the performance of the CEPOF. The CEPOF algorithm demonstrates that full-width-at- half-maximum energy resolution of < 8 eV coupled with position-sensitivity down to a few 100 eV should be achievable for a fully optimized device.

Smith, Stephen J.↗

Maglev Facility for Simulating Variable Gravity

An improved magnetic levitation apparatus ("Maglev Facility") has been built for use in experiments in which there are requirements to impose variable gravity (including zero gravity) in order to assess the effects of gravity or the absence thereof on physical and physiological processes. The apparatus is expected to be especially useful for experiments on the effects of gravity on convection, boiling, and heat transfer in fluids and for experiments on mice to gain understanding of bone loss induced in human astronauts by prolonged exposure to reduced gravity in space flight. The maglev principle employed by the apparatus is well established. Diamagnetic cryogenic fluids such as liquid helium have been magnetically levitated for studying their phase transitions and critical behaviors. Biological entities consist mostly of diamagnetic molecules (e.g., water molecules) and thus can be levitated by use of sufficiently strong magnetic fields having sufficiently strong vertical gradients. The heart of the present maglev apparatus is a vertically oriented superconducting solenoid electromagnet (see figure) that generates a static magnetic field of about 16 T with a vertical gradient sufficient for levitation of water in normal Earth gravity. The electromagnet is enclosed in a Dewar flask having a volume of 100 L that contains liquid helium to maintain superconductivity. The Dewar flask features a 66-mm-diameter warm bore, lying within the bore of the magnet, wherein experiments can be performed at room temperature. The warm bore is accessible from its top and bottom ends. The superconducting electromagnet is run in the persistent mode, in which the supercurrent and the magnetic field can be maintained for weeks with little decay, making this apparatus extremely cost and energy efficient to operate. In addition to water, this apparatus can levitate several common fluids: liquid hydrogen, liquid oxygen, methane, ammonia, sodium, and lithium, all of which are useful, variously, as rocket fuels or as working fluids for heat transfer devices. A drop of water 45 mm in diameter and a small laboratory mouse have been levitated in this apparatus.

Liu, Yuanming↗

Space Shuttle Main Engine Propellant Path Leak Detection Using Sequential Image Processing

Initial research in this study using theoretical radiation transport models established that the occurrence of a leak is accompanies by a sudden but sustained change in intensity in a given region of an image. In this phase, temporal processing of video images on a frame-by-frame basis was used to detect leaks within a given field of view. The leak detection algorithm developed in this study consists of a digital highpass filter cascaded with a moving average filter. The absolute value of the resulting discrete sequence is then taken and compared to a threshold value to produce the binary leak/no leak decision at each point in the image. Alternatively, averaging over the full frame of the output image produces a single time-varying mean value estimate that is indicative of the intensity and extent of a leak. Laboratory experiments were conducted in which artificially created leaks on a simulated SSME background were produced and recorded from a visible wavelength video camera. This data was processed frame-by-frame over the time interval of interest using an image processor implementation of the leak detection algorithm. In addition, a 20 second video sequence of an actual SSME failure was analyzed using this technique. The resulting output image sequences and plots of the full frame mean value versus time verify the effectiveness of the system.

Smith, L. Montgomery↗

A New Unsteady Model for Dense Cloud Cavitation in Cryogenic Fluids

A new unsteady, cavitation model is presented wherein the phase change process (bubble growth/collapse) is coupled to the acoustic field in a cryogenic fluid. It predicts the number density and radius of bubbles in vapor clouds by tracking both the aggregate surface area and volume fraction of the cloud. Hence, formulations for the dynamics of individual bubbles (e.g. Rayleigh-Plesset equation) may be integrated within the macroscopic context of a dense vapor cloud i.e. a cloud that occupies a significant fraction of available volume and contains numerous bubbles. This formulation has been implemented within the CRUNCH CFD, which has a compressible real fluid formulation, a multi-element, unstructured grid framework, and has been validated extensively for liquid rocket turbopump inducers. Detailed unsteady simulations of a cavitating ogive in liquid nitrogen are presented where time-averaged mean cavity pressure and temperature depressions due to cavitation are compared with experimental data. The model also provides the spatial and temporal history of the bubble size distribution in the vapor clouds that are shed, an important physical parameter that is difficult to measure experimentally and is a significant advancement in the modeling of dense cloud cavitation.

Hosangadi, A.↗

A wall interference assessment/correction system

A Wall Signature method, the Hackett method, has been selected to be adapted for the 12-ft Wind Tunnel wall interference assessment/correction (WIAC) system in the present phase. This method uses limited measurements of the static pressure at the wall, in conjunction with the solid wall boundary condition, to determine the strength and distribution of singularities representing the test article. The singularities are used in turn for estimating wall interferences at the model location. The Wall Signature method will be formulated for application to the unique geometry of the 12-ft Tunnel. The development and implementation of a working prototype will be completed, delivered and documented with a software manual. The WIAC code will be validated by conducting numerically simulated experiments rather than actual wind tunnel experiments. The simulations will be used to generate both free-air and confined wind-tunnel flow fields for each of the test articles over a range of test configurations. Specifically, the pressure signature at the test section wall will be computed for the tunnel case to provide the simulated 'measured' data. These data will serve as the input for the WIAC method-Wall Signature method. The performance of the WIAC method then may be evaluated by comparing the corrected parameters with those for the free-air simulation. Each set of wind tunnel/test article numerical simulations provides data to validate the WIAC method. A numerical wind tunnel test simulation is initiated to validate the WIAC methods developed in the project. In the present reported period, the blockage correction has been developed and implemented for a rectangular tunnel as well as the 12-ft Pressure Tunnel. An improved wall interference assessment and correction method for three-dimensional wind tunnel testing is presented in the appendix.

Lo, Ching F.↗

Wind Forcing of the North Sea Pole Tide

The Chandler wobble of the earth's rotation has a period near 14 months and sets up the 0.5 cm amplitude pole tide in the deep oceans, However, the pole tide is anomalously large in the North Sea, where the amplitude increases sharply up to 3 cm along the continental coast. It is shown here that the sea levels are well correlated with winds at the pole tide frequency. The Princeton Ocean Model is used to investigate the response of the North Sea to wind forcing. The barotropic numerical ocean model depicts realistic coastlines and bathymetry at 5 ft x 5 ft resolution, with 97 x 73 grid points. The monthly mean wind fields for the 40-year period (1958-1997) from the National Centers for Atmospheric Prediction (NCEP) reanalysis were used to force the model. The winds were converted to stress with a neutral drag coefficient that varied linearly with windspeed (instead of using the NCEP windstress). A 5-day simulation was made for each month until the resulting flow regime came into equilibrium, and model water levels at various station locations were saved for comparison with tidal-gauge observed sea levels from the Permanent Service for Mean Sea Level (PSMSL). The comparison is made for 10 North Sea stations with high quality tide gauge data. We find: (1) good agreement in annual and semi-annual phases and in the trend of amplitude w.r.t. latitude; (2) more importantly for this study, the model-predicted and observed non-seasonal sea level variations show a very significant temporal correlation as well as spectral coherence. However, a large amplitude difference exists between the two sets -- the overall amplitude variability of the observed is generally a factor of 2-3 larger than the model prediction (this same phenomenon has been reported in ocean circulation studies, although the cause is not yet clear.) Our results indicate that the wind forcing is the main cause of the observed large pole tide in the North Sea.

OConnor, W.↗

Mind-the-Gap - Part 3: Doppler Velocity Measurements From Space

Convective motions and hydrometeors microphysical properties are highly sought-after parameters for evaluating atmospheric numerical models. With most of the Earth’s surface covered by water, space-borne Doppler radars are ideal for acquiring such measurements at a global scale. While these systems have proven to be useful tools for retrieving cloud microphysical and dynamical properties from the ground, their adequacy, and specific requirements for spaceborne operation still need to be evaluated. Comprehensive forward simulations enable us to assess the advantages and drawbacks of six different Doppler radar architectures currently planned or under consideration by space agencies for the study of cloud dynamics. Radar performance is examined against state-of-the-art numerical model simulations of well-characterized shallow and deep, continental, and oceanic convective cases. Mean Doppler velocity (MDV) measurements collected at multiple frequencies (13, 35 and 94 GHz) provide complementary information in deep convective cloud systems. The high penetration capability of the 13-GHz radar enables to obtain a complete, albeit horizontally under-sampled, view of deep convective storms. The smaller instantaneous field of view (IFOV) of the 35-GHz radar captures more precise information about the location and size of convective updrafts above 5-8 km height of most systems which was determined is the portion of storms where the mass flux peak is typically located. Finally, the lower mean Doppler velocity uncertainty of displaced phase center antenna (DPCA) radars makes them an ideal system for studying microphysics in shallow convection and frontal systems, as well as ice and mixed-phase clouds. It is demonstrated that a 94-GHz DCPA system can achieve retrieval errors as low as 0.05-0.15 mm for raindrop volume-weighted mean diameter and 25% for rime fraction (for a -10 dBZ echo).

Pavlos Kollias↗

Active control of propeller induced noise fields inside a flexible cylinder

An active noise control model has been evaluated for reducing aircraft interior noise. The structural noise transmission properties of an aircraft fuselage were modelled as a flexible cylinder excited by external acoustic dipoles simulating the noise produced by twin propellers. The amplitudes of an internal distribution of monopole control sources were determined such that the area-weighted mean square acoustic pressure was minimized in the propeller plane. The noise control model was evaluated at low frequencies corresponding to the blade passage frequency and first few harmonics of a typical turbo-prop aircraft. Interior noise reductions of 20 25 dB were achieved, over a substantial region of the cylindrical cross-section, with just a few monopole control sources. The most favorable interior noise reductions were achieved when the active noise control model was used in combination with propeller source phasing.

Lester, H. C.↗

Plasma Pressure in the Topside Ionosphere

A previous three year NASA-funded project resulted in the first 2-D maps of magnetotail pressure, density and temperature. A proposal to continue the work was declined, but modest funding was provided for one year to ramp down of the work. During the phase-out year, we used a time when 5 DMSP satellites were simultaneously active to produce the first instantaneous partial image of the magnetotail. The results have been submitted to the proceedings of the 1998 Huntsville Meeting on "The New Millennium Magnetosphere: Integrating Imaging, Discrete Observations and Global Simulations". A method of inferring central plasma sheet (CPS) temperature, density, and pressure from ionospheric observations was developed under a previous 3-year grant. These particles properties are calculated from data taken by particle instruments on DMSP satellites. Ion spectra occurring in conjunction with electron acceleration events are excluded. Because of the variability of magnetotail stretching, mapping to the plasma sheet was done using a modified Tsyganenko 1989 magnetic field model adjusted to agree with the actual magnetotail stretch. On May 25, 1997, five DMSP satellites (F10-F14) passed through the southern hemisphere nightside oval within a 19 minute period. Attached is the first magnetotail image, which results from applying our technique to that data set.

Newell, Patrick T.↗

Runtime and Pressurization Analyses of Propellant Tanks

Multi-element unstructured CFD has been utilized at NASA SSC to carry out analyses of propellant tank systems in different modes of operation. The three regimes of interest at SSC include (a) tank chill down (b) tank pressurization and (c) runtime propellant draw-down and purge. While tank chill down is an important event that is best addressed with long time-scale heat transfer calculations, CFD can play a critical role in the tank pressurization and runtime modes of operation. In these situations, problems with contamination of the propellant by inclusion of the pressurant gas from the ullage causes a deterioration of the quality of the propellant delivered to the test article. CFD can be used to help quantify the mixing and propellant degradation. During tank pressurization under some circumstances, rapid mixing of relatively warm pressurant gas with cryogenic propellant can lead to rapid densification of the gas and loss of pressure in the tank. This phenomenon can cause serious problems during testing because of the resulting decrease in propellant flow rate. With proper physical models implemented, CFD can model the coupling between the propellant and pressurant including heat transfer and phase change effects and accurately capture the complex physics in the evolving flowfields. This holds the promise of allowing the specification of operational conditions and procedures that could minimize the undesirable mixing and heat transfer inherent in propellant tank operation. It should be noted that traditional CFD modeling is inadequate for such simulations because the fluids in the tank are in a range of different sub-critical and supercritical states and elaborate phase change and mixing rules have to be developed to accurately model the interaction between the ullage gas and the propellant. We show a typical run-time simulation of a spherical propellant tank, containing RP-1 in this case, being pressurized with room-temperature nitrogen at 540 R. Nitrogen, shown in blue on the right-hand side of the figures, enters the tank from the diffuser at the top of the figures and impinges on the RP-1, shown in red, while the propellant is being continuously drained at the rate of 1050 lbs/sec through a pipe at the bottom of the tank. The sequence of frames in Figure 1 shows the resultant velocity fields and mixing between nitrogen and RP-1 in a cross-section of the tank at different times. A vortex is seen to form in the incoming nitrogen stream that tends to entrain propellant, mixing it with the pressurant gas. The RP-1 mass fraction contours in Figure 1 are also indicative of the level of mixing and contamination of the propellant. The simulation is used to track the propagation of the pure propellant front as it is drawn toward the exit with the evolution of the mixing processes in the tank. The CFD simulation modeled a total of 10 seconds of run time. As is seen from Figure 1d, after 5.65 seconds the propellant front is nearing the drain pipe, especially near the center of the tank. Behind this pure propellant front is a mixed fluid of compromised quality that would require the test to end when it reaches the exit pipe. Such unsteady simulations provide an estimate of the time that a high-quality propellant supply to the test article can be guaranteed at the modeled mass flow rate. In the final paper, we will discuss simulations of the LOX and propellant tanks at NASA SSC being pressurized by an inert ullage. Detailed comparisons will be made between the CFD simulations and lower order models as well as with test data. Conditions leading to cryo collapse in the tank will also be identified.

Field, Robert E.↗

Studies of the Effects of Actuator Errors on the HCIT/PIAA Contrast Performance

The High Contrast Imaging Testbed Phase Induced Amplitude Apodization (HCIT/PIAA) coronagraph system at JPL relies on an Electric-Field Conjugation (EFC) wavefront correction algorithm to obtain the required high contrast. This algorithm works with one deformable mirror (DM) to estimate the electric-field to be controlled, and with one or multiple DM's to create a "dark-hole" in the image plane. We have investigated the effects of DM actuator errors on the efficiency of the EFC algorithm. The structural design of the optical system as well as the parameters of various optical elements used in the analysis are drawn from those of the HCIT/PIAA system that have been and will be implemented with one or two DM's. The simulation takes into account the surface errors of various optical elements. In this paper, we report our findings in the case of narrowband wavelength light.

PIAA↗

Earth-System-Model Evaluation of Cloud and Precipitation Occurrence for Supercooled and Warm Clouds Over the Southern Ocean's Macquarie Island

Over the remote Southern Ocean (SO), cloud feedbacks contribute substantially to Earth system model (ESM) radiative biases. The evolution of low Southern Ocean clouds (cloud-top heights < ∼ 3 km) is strongly modulated by precipitation and/or evaporation, which act as the primary sink of cloud condensate. Constraining precipitation processes in ESMs requires robust observations suitable for process-level evaluations. A year-long subset (April 2016–March 2017) of ground-based profiling instrumentation deployed during the Macquarie Island Cloud and Radiation Experiment (MICRE) field campaign (54.5∘ S, 158.9∘ E) combines a 95 GHz (W-band) Doppler cloud radar, two lidar ceilometers, and balloon-borne soundings to quantify the occurrence frequency of precipitation from the liquid-phase cloud base. Liquid-based clouds at Macquarie Island precipitate ∼ 70 % of the time, with deeper and colder clouds precipitating more frequently and at a higher intensity compared to thinner and warmer clouds. Supercooled cloud layers precipitate more readily than layers with cloud-top temperatures > 0 ∘C, regardless of the geometric thickness of the layer, and also evaporate more frequently. We further demonstrate an approach to employ these observational constraints for evaluation of a 9-year GISS-ModelE3 ESM simulation. Model output is processed through the Earth Model Column Collaboratory (EMC2) radar and lidar instrument simulator with the same instrument specifications as those deployed during MICRE, therefore accounting for instrument sensitivities and ensuring a coherent comparison. Relative to MICRE observations, the ESM produces a smaller cloud occurrence frequency, smaller precipitation occurrence frequency, and greater sub-cloud evaporation. The lower precipitation occurrence frequency by the ESM relative to MICRE contrasts with numerous studies that suggest a ubiquitous bias by ESMs to precipitate too frequently over the SO when compared with satellite-based observations, likely owing to sensitivity limitations of spaceborne instrumentation and different sampling methodologies for ground- versus space-based observations. Despite these deficiencies, the ESM reproduces the observed tendency for deeper and colder clouds to precipitate more frequently and at a higher intensity. The ESM also reproduces specific cloud regimes, including near-surface clouds that account for ∼ 25 % of liquid-based clouds during MICRE and optically thin, non-precipitating clouds that account for ∼ 27 % of clouds with bases higher than 250 m. We suggest that the demonstrated framework, which merges observations with appropriately constrained model output, is a valuable approach to evaluate processes responsible for cloud radiative feedbacks in ESMs.

cloud feedbacks↗

Evaluation of Mixed-Phase Microphysics Within Winter Storms Using Field Data and In Situ Observations

Snow prediction within models is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33-km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Cold-season Precipitation Experiment (GCPEx) experiment over southern Ontario, as well as a few years (12 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. The GCPex presentation will focus on verification using aircraft spirals through warm frontal snow band event on 18 February 2012. All the BMPs realistically simulated the structure of the band and the vertical distribution of snow/ice aloft, except the SBU-YLIN overpredicted slightly and Thompson (THOM) underpredicted somewhat. The Morrison (MORR) scheme produced the best slope size distribution for snow, while the Stony Brook (SBU) underpredicted and the THOM slightly overpredicted. Those schemes that have the slope intercept a function of temperature (SBU and WSM6) tended to perform better for that parameter than others, especially the fixed intercept in Goddard. Overall, the spread among BMPs was smaller than in other studies, likely because there was limited riming with the band. For the 15 cases at SBNY, which include moderate and heavy riming events, the non-spherical snow assumption (THOM and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM is approx. 0.25 m/s too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM, and MORR schemes were approx. 0.25 m/s too slow, while the SBU was 0.25 to 0.5 m/s too fast because of some excessive cloud water issues.

Colle, Brian A.↗