Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive 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 361 records · Page 20

Finite Element Simulations of Two Vertical Drop Tests of F-28 Fuselage Sections

In March 2017, a vertical drop test of a forward fuselage section of a Fokker F-28 MK4000 aircraft was conducted as part of a joint NASA/FAA project to investigate the performance of transport aircraft under realistic crash conditions. In June 2017, a vertical drop test was conducted of a wing-box fuselage section of the same aircraft. Both sections were configured with two rows of aircraft seats, in a triple-double configuration. A total of ten Anthropomorphic Test Devices (ATDs) were secured in seats using standard lap belt restraints. The forward fuselage section was also configured with luggage in the cargo hold. Both sections were outfitted with two hat racks, each with added ballast mass. The drop tests were performed at the Landing and Impact Research facility located at NASA Langley Research Center in Hampton, Virginia. The measured impact velocity for the forward fuselage section was 346.8-in/s onto soil. The wing-box section was dropped with a downward facing pitch angle onto a sloping soil surface in order to create an induced forward acceleration in the airframe. The vertical impact velocity of the wing-box section was 349.2-in/s. A second objective of this project was to assess the capabilities of finite element simulations to predict the test responses. Finite element models of both fuselage sections were developed for execution in LS-DYNA(Registered Trademark), a commercial explicit nonlinear transient dynamic code. The models contained accurate representations of the airframe structure, the hat racks and hat rack masses, the floor and seat tracks, the luggage in the cargo hold for the forward section, and the detailed under-floor structure in the wing-box section. Initially, concentrated masses were used to represent the inertial properties of the seats, restraints, and ATD occupants. However, later simulations were performed that included finite element representations of the seats, restraints, and ATD occupants. These models were developed to more accurately replicate the seat loading of the floor and to enable prediction of occupant impact responses. Models were executed to generate analytical predictions of airframe responses, which were compared with test data to validate the model. Comparisons of predicted and experimental structural deformation and failures were made. Finally, predicted and experimental soil deformation and crater depths were also compared for both drop test configurations.

Jackson, Karen E.↗

Viscous flow effects on hydrogen leaks from cracks in the Orbiter Challenger main engines

An analytical model was developed to provide additional insight and understanding of the factors that influence the simulation and prediction of leak rates from small cracks in pressurized containers. Specifically, the analysis was aimed at developing an analytical model capable of predicting the hydrogen leak rates from a crack in the combustion chamber coolant discharge manifold on main engine 1 of the Orbiter Challenger that was discovered during flight readiness firings 1 and 2. This model was based on viscous pipe flow analyses and calibrated for the crack geometry by using helium leak-rate data obtained from both low- and high-pressure tests used to simulate the flight readiness firing test conditions. In addition, this model includes the effects of crack width changes caused by different working stresses associated with the different test conditions. Because of the combination of the small crack dimensions and the wide range of pressures used for the test conditions, either laminar or turbulent viscous effects dominated the flows at all test conditions. This model was used to illustrate the sensitivity of the predicted leak rates to considerations of test conditions, viscous flow effects, and geometric features of the crack. In addition, the model was certified by comparing the hydrogen leak-rate prediction for the flight readiness firing test condition to the actual measured leak rate. The prediction was within 9 percent of the measured value.

Goodrich, W. D.↗

The Use of Steady and Pulsed Detonations for Propulsion Systems

Objectives of the ODWE concept studies are: demonstrate the feasibility of the oblique detonation wave engine (ODWE) for hypersonic propulsion; demonstrate the existance and stability of an oblique detonation wave in hypersonic wind tunnels; develop engineering codes which predict the performance characteristics of the ODWE including specific impulse and thrust coefficients for various operating conditions; develop multi-dimensional computer codes which can model all aspects of the ODWE including fuel injection, mixing, ignition, combustion and expansion with fully detailed chemical kinetics and turbulence models; and validate the codes with experimental data use the simulations to predict the ODWE performance for conditions not easily obtained in wind tunnels.

Adelman, Henry G.↗

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗

Flat-Plate Solar-Collector Performance Evaluation with a Solar Simulator as a Basis for Collector Selection and Performance Prediction

The use of a solar simulator for performance determination permits collector testing under standard conditions of wind, ambient temperature, flow rate and sun. The performance results determined with the simulator have been found to be in good agreement with outdoor performance results. The measured thermal efficiency and evaluation of 23 collectors are reported which differ according to absorber material (copper, aluminum, steel), absorber coating (nonselective black paint, selective copper oxide, selective black nickel, selective black chrome), type of glazing material (glass, Tedlar, Lexan, antireflection glass), the use of honeycomb material and the use of vacuum to prevent thermal convection losses. The collectors were given performance rankings based on noon-hour solar conditions and all-day solar conditions. The determination with the simulator of an all-day collector performance was made possible by tests at different incident angles. The solar performance rankings were made based on whether the collector is to be used for pool heating, hot water, absorption air conditioning, heating, or for a solar Rankine machine.

Simon, F. F.↗

Simulation studies of proposed observing systems and their impact on numerical weather prediction

A series of realistic simulation studies is being conducted as a cooperative effort between the European Centre for Medium Range Weather Forecasts (ECMWF), the National Meteorological Center (NMC), and the Goddard Laboratory for Atmospheric Sciences (GLAS) to provide a quantitative assessment of the potential impact of proposed observation systems on large scale numerical weather prediction. A special objective of this project is to avoid the unrealistic character of earlier simulation studies.

Atlas, R.↗

WRF-Chem Simulations of Lightning-NOx Production and Transport in Oklahoma and Colorado Thunderstorms Observed During DC3

The focus of this analysis is on lightning-generated nitrogen oxides (LNOx) and their distribution for two thunderstorms observed during the Deep Convective Clouds and Chemistry (DC3) field campaign in May-June 2012. The Weather Research and Forecasting Chemistry (WRF-Chem) model is used to perform cloud-resolved simulations for the May 29-30 Oklahoma severe convection, which contained one supercell, and the June 6-7 Colorado squall line. Aircraft and ground-based observations (e.g., trace gases, lightning and radar) collected during DC3 are used in comparisons against the model-simulated lightning flashes generated by the flash rate parameterization schemes (FRPSs) incorporated into the model, as well as the model-simulated LNOx predicted in the anvil outflow. Newly generated FRPSs based on DC3 radar observations and Lightning Mapping Array data are implemented in the model, along with previously developed schemes from the literature. The results of these analyses will also be compared between storms to investigate which FRPSs were most appropriate for the two types of convection and to examine the variation in the LNOx production. The simulated LNOx results from WRF-Chem will also be compared against other previously studied mid-latitude thunderstorms.

Lightning↗

Residual Strength Analyses of Monolithic Structures

Finite-element fracture simulation methodology predicts the residual strength of damaged aircraft structures. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic- plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intention of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels.In this regard, a series of fracture tests were conducted on both flat and curved aluminum alloy integrally-stiffened panels. These flat panels were subjected to uniaxial tension and during the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (wc) using three-dimensional code (ZIP3D) and the plane-strain core height (hJ using two-dimensional code (STAGS). These values were then used in the STAGS analysis to predict the fracture behavior of the integrally-stiffened panels. The analyses modeled stable tearing, buckling, and crack branching at the integral stiffener using different values of critical CTOA for different material thicknesses and orientation. Comparisons were made between measured and predicted load-crack extension, out-of-plane displacements and local deformations around the crack tip region. Simultaneously, three-dimensional capabilities to model crack branching and to monitor stable crack growth of multiple cracks in a large thick integrally-stiffened flat panels were implemented in three-dimensional finite element code (ZIP3D) and tested by analyzing the integrally-stiffened panels tested at Alcoa. The residual strength of the panels predicted from STAGS and ZP3D code compared very well with experimental data. In recent times, STAGS software has been updated with new features and now one can have combinations of solid and shell elements in the residual strength analysis of integrally-stiffened panels.

Forth, Scott↗

USM3D-ME Revised SA-neg Implementation and Assessment for Seventh Drag Prediction Workshop Solutions

Simulations performed by an unstructured-grid, cell-centered, finite-volume, mixed-element flow solver USM3D-MEwere presented at the Seventh Drag Prediction Workshop (DPW-7). Reynolds-averaged Navier-Stokes (RANS) solutions for the NASA high-speed common research model wing-body configuration were computed using a new implementation of the negative variant of the Spalart-Allmaras turbulence model (SA-neg) with hierarchical adaptive nonlinear iteration method (HANIM).ThisUSM3D-ME enhancement improved efficiency and robustness of RANS solutions resulting in fast iterative convergence to machine zero residuals on mixed-element unstructured grids provided by the DPW-7 organizing committee. This paper presents theoretical foundations for the new implementation of the SA-neg model and analyzesUSM3D-ME solutions for flow conditions corresponding to the required DPW-7 cases.

Boris Diskin↗

A comparison of the noise produced by a small jet on a moving vehicle with that in a free jet

A 2.54 cm (1.00 in.) nozzle supplied with nitrogen was mounted above an automobile and driven over an asphalt roadway past stationary microphones in an attempt to quantify the effects of the vehicle motion on jet mixing noise. The nozzle was then tested in the Langley anechoic noise facility with a large free jet simulating the relative motion. The results are compared for these two methods of investigating forward speed effects on jet mixing noise. The vehicle results indicate a noise with forward speed throughout the Doppler-shifted static spectrum. This decrease across the entire frequency range was also apparent in the free-jet results. The similarity of the results indicates that the effects of flight on jet mixing noise can be predicted by simulation of forward speed with a free jet. Overall sound pressure levels were found to decrease with forward speed at all observation angles for both methods of testing.

Norum, T. D.↗

Steady and Transient Performance Prediction of Gas Turbine Engines

This Lecture Series covers topics of aero-thermodynamic performance prediction methods for gas turbine engines with respect to steady and transient operation. This includes advanced cycle calculation methods, also taking into account variable cycle engine types. A very important objective is the consideration of installation effect, i.e., Reynolds number and inlet distortions, as well as advanced control concepts for increasing engine surge margins. In addition to these topics, individual papers include practical considerations in developing the engine cycle, dynamic simulation, inlet distortion effects in aircraft propulsion system integration, “smart” engines, performance and heal monitoring models. This Lecture Series, sponsored by the Propulsion and Energetics Panel of AGARD and the Consultant and Exchange Programme of AGARD presented on 27th -28th May 1992 in Cambridge, MA, United States, 9th-10th June 1992 in Neubiberg, Germany and 11th-12th June 1992 in Chàtillon/Bagneux (near Paris), France

Aircraft engines↗

Predicting Melt Properties Using Atomistic Simulations With A Highly Accurate Physically Informed Neural Network Interatomic Potential

The use of a recently developed machine learning (ML) interatomic potential for molecular dynamics simulations of aluminum melt properties will be presented. Such properties are critical for process modeling in additive manufacturing, including the melt pool size, solidification, and formation of solidification microstructures. Direct first-principles modeling of these processes is computationally prohibitive whereas simulations employing ML potentials combine the high accuracy of quantum-mechanical methods with high computational speeds. The physically-informed neural network (PINN) method used herein, integrates a high-dimensional regression implemented by an artificial neural network with a physics-based bond-order interatomic potential. PINN potentials can accurately reproduce many properties of aluminum in both crystalline-solid and liquid phases. We examine the accuracy of a PINN Al potential in predicting the density, self-diffusivity, viscosity, and the tension of the liquid surface and liquid-solid interfaces. Comparison with experimental data and ab initio molecular dynamics calculations shows very good agreement for all properties tested.

molecular dynamics↗

Piloted Simulation of a Model-Predictive Automated Recovery System

This presentation describes a model-predictive automatic recovery system for aircraft on the verge of a loss-of-control situation. The system determines when it must intervene to prevent an imminent accident, resulting from a poor approach. It estimates the altitude loss that would result from a go-around maneuver at the current flight condition. If the loss is projected to violate a minimum altitude threshold, the maneuver is automatically triggered. The system deactivates to allow landing once several criteria are met. Piloted flight simulator evaluation showed the system to provide effective envelope protection during extremely unsafe landing attempts. The results demonstrate how flight and propulsion control can be integrated to recover control of the vehicle automatically and prevent a potential catastrophe.

aircraft engine simulation↗

Predictive Modeling of Fast-Curing Thermosets in Nozzle-Based Extrusion

This work presents an approach to modeling the dynamic spreading and curing behavior of thermosets in nozzle-based extrusions. Thermosets cover a wide range of materials, some of which permit low-temperature processing with subsequent high-temperature and high-strength working properties. Extruding thermosets may overcome the limited working temperatures and strengths of conventional thermoplastic materials used in additive manufacturing. This project aims to produce technology for the fabrication of thermoset-based structures leveraging advances made in nozzle-based extrusion, such as fused deposition modeling (FDM), material jetting, and direct writing. Understanding the synergistic interactions between spreading and fast curing of extruded thermosetting materials will provide essential insights for applications that require accurate dimensional controls, such as additive manufacturing [1], [2] and centrifugal coating/forming [3]. Two types of thermally curing thermosets -- one being a soft silicone (Ecoflex 0050) and the other being a toughened epoxy (G/Flex) -- served as the test materials in this work to obtain models for cure kinetics and viscosity. The developed models align with extensive measurements made with differential scanning calorimetry (DSC) and rheology. DSC monitors the change in the heat of reaction, which reflects the rate and degree of cure at different crosslinking stages. Rheology measures the change in complex viscosity, shear moduli, yield stress, and other properties dictated by chemical composition. By combining DSC and rheological measurements, it is possible to establish a set of models profiling the cure kinetics and chemorheology without prior knowledge of chemical composition, which is usually necessary for sophisticated mechanistic modeling. In this work, we conducted both isothermal and dynamic measurements with both DSC and rheology. With the developed models, numerical simulations yielded predictions of diameter and height of droplets, along with width and height of extruded lines cured at varied temperatures. Experimental results carried out on a goniometric platform and a nozzle-based 3D printer showed agreement with the numerical simulations. Finally, this presentation will show how the models are adaptable to the planning of tool paths and designs in additive manufacturing.

Xie, Jingjin↗

Accounting for Uncertainties in Strengths of SiC MEMS Parts

A methodology has been devised for accounting for uncertainties in the strengths of silicon carbide structural components of microelectromechanical systems (MEMS). The methodology enables prediction of the probabilistic strengths of complexly shaped MEMS parts using data from tests of simple specimens. This methodology is intended to serve as a part of a rational basis for designing SiC MEMS, supplementing methodologies that have been borrowed from the art of designing macroscopic brittle material structures. The need for this or a similar methodology arises as a consequence of the fundamental nature of MEMS and the brittle silicon-based materials of which they are typically fabricated. When tested to fracture, MEMS and structural components thereof show wide part-to-part scatter in strength. The methodology involves the use of the Ceramics Analysis and Reliability Evaluation of Structures Life (CARES/Life) software in conjunction with the ANSYS Probabilistic Design System (PDS) software to simulate or predict the strength responses of brittle material components while simultaneously accounting for the effects of variability of geometrical features on the strength responses. As such, the methodology involves the use of an extended version of the ANSYS/CARES/PDS software system described in Probabilistic Prediction of Lifetimes of Ceramic Parts (LEW-17682-1/4-1), Software Tech Briefs supplement to NASA Tech Briefs, Vol. 30, No. 9 (September 2006), page 10. The ANSYS PDS software enables the ANSYS finite-element-analysis program to account for uncertainty in the design-and analysis process. The ANSYS PDS software accounts for uncertainty in material properties, dimensions, and loading by assigning probabilistic distributions to user-specified model parameters and performing simulations using various sampling techniques.

Nemeth, Noel↗

Parameter Uncertainty for Aircraft Aerodynamic Modeling using Recursive Least Squares

A real-time method was demonstrated for determining accurate uncertainty levels of stability and control derivatives estimated using recursive least squares and time-domain data. The method uses a recursive formulation of the residual autocorrelation to account for colored residuals, which are routinely encountered in aircraft parameter estimation and change the predicted uncertainties. Simulation data and flight test data for a subscale jet transport aircraft were used to demonstrate the approach. Results showed that the corrected uncertainties matched the observed scatter in the parameter estimates, and did so more accurately than conventional uncertainty estimates that assume white residuals. Only small differences were observed between batch estimates and recursive estimates at the end of the maneuver. It was also demonstrated that the autocorrelation could be reduced to a small number of lags to minimize computation and memory storage requirements without significantly degrading the accuracy of predicted uncertainty levels.

Grauer, Jared A.↗

Evaluation of Cloud Microphysics Simulated using a Meso-Scale Model Coupled with a Spectral Bin Microphysical Scheme through Comparison with Observation Data by Ship-Borne Doppler and Space-Borne W-Band Radars

Equivalent radar reflectivity factors (Ze) measured by W-band radars are directly compared with the corresponding values calculated from a three-dimensional non-hydrostatic meso-scale model coupled with a spectral-bin-microphysical (SBM) scheme for cloud. Three case studies are the objects of this research: one targets a part of ship-borne observation using 95 GHz Doppler radar over the Pacific Ocean near Japan in May 2001; other two are aimed at two short segments of space-borne observation by the cloud profiling radar on CloudSat in November 2006. The numerical weather prediction (NWP) simulations reproduce general features of vertical structures of Ze and Doppler velocity. A main problem in the reproducibility is an overestimation of Ze in ice cloud layers. A frequency analysis shows a strong correlation between ice water contents (IWC) and Ze in the simulation; this characteristic is similar to those shown in prior on-site studies. From comparing with the empirical correlations by the prior studies, the simulated Ze is overestimated than the corresponding values in the studies at the same IWC. Whereas the comparison of Doppler velocities suggests that large-size snowflakes are necessary for producing large velocities under the freezing level and hence rules out the possibility that an overestimation of snow size causes the overestimation of Ze. Based on the results of several sensitivity tests, we conclude that the source of the overestimation is a bias in the microphysical calculation of Ze or an overestimation of IWC. To identify the source of the problems needs further validation research with other follow-up observations.

Iguchi, T.↗