Search NASA⌕ Search

SEARCH · Search NASA

Results for “Performance Rating Method”

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.

1,608 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Metallic Powder Core Tubular Wire (PCTW) Development for Additive Manufacturing

NASA Langley Research Center (LaRC) in Hampton, VA, and the Center for Welding, Joining, and Coating Research (CWJCR) in the Department of Metallurgical and Materials Engineering at the Colorado School of Mines (CSM) in Golden, CO, have collaborated for nearly fifteen years to develop metallic powder core tubular wire (PCTW) feedstock for additive manufacturing (AM) using electron beam deposition. The motivation for this development was to offset dealloying observed during deposition of titanium and aluminum alloy solid wire feedstock in the LaRC electron beam freeform fabrication (EBF3) system and to explore novel alloy compositions by tailoring the powder fill. EBF3 is a fusion-based AM process that is performed under vacuum and is consequently prone to vaporization loss of low vapor pressure alloying elements. Research spanning four graduate student research projects succeeded in developing the basic methodology to fabricate PCTW, the mass balance equations used to define the powder fill, microstructure control through particulate inoculation, and the feasibility of producing metal matrix composite materials. Metallic PCTW is comprised of a thin foil sheath wrapped around a powder blend. Fabrication of PCTW at CSM was accomplished using a tubular wire mill system to first bend a flat strip into a u-shaped cross section, deposit a powder fill, create a lap joint to seal the wire, and draw to final diameter. Powder management systems were developed to ensure uniform mixing of powder constituents and identify feed rates for optimum fill ratios. Processing parameters were successfully developed for titanium and aluminum alloy sheaths to account for material work hardening during drawing. Mass balance equations were developed to account for compositional contributions of the sheath and powder fill. PCTW composition can be adjusted through selection of the sheath or powder materials. PCTWs were successfully developed to offset Al loss in deposits of Ti-6-4 and Mg loss in Al 6061, resulting in deposits that were within the composition limits for each material. Ti-6-4 deposits were fabricated using PCTW comprised of commercially pure (CP) titanium sheath and a fill of pre-alloyed Ti-6-4 powder with additions of elemental Al and V powders. Additions of both elements were required to account for the CP Ti from the sheath, and Al was further enhanced to offset documented losses in EBF3 deposits. For Al 6061 deposits, composition was modulated through the choice of sheath material. The PCTW was made using an Al 5052 sheath to supplement Mg content and a fill composition of pre-alloyed 6061 powder with additions of Si and Cu.The slower solidification rate of the EBF3 process compared with laser powder bed fusion processes can lead to coarse-grained microstructures and mechanical property anisotropy. EBF3deposits of Ti-6-4 exhibited large epitaxially grown β grains, which result in degraded mechanical properties in the deposition direction. Microstructure refinement was achieved through a combination of electron beam modulation and the use of a PCTW designed for Ti-6-4 deposits and modified by the addition of powder made from Fe and B. Resulting deposits exhibited reductions in β grain size, α phase colony intercept length and lath width, and an increased volume fraction of β phase. The mechanism of microstructure refinement was related to precipitation of boride particles. Hardness testing confirmed that mechanical property anisotropy was reduced and an overall increase in hardness was achieved.Metal matrix composite (MMC) materials offer strength and stiffness benefits over unreinforced materials. Two methods were investigated using PCTW technology for producing Al 6061 basedMMC materials; the direct addition of SiC reinforcing particles and the in-situ generation of reinforcements through reaction synthesis of precursor powders. Al 6061 based MMCs reinforced with SiC exhibited particle clustering and formation of aluminum carbides during deposition. The application of a Ni coating to the SiC particles prior to incorporation into the powder fill eliminated aluminum carbide formation and aided particle distribution during deposition. Modulation of beam focus and power further reduced aluminum carbide formation and aided particle distribution. In reaction synthesis, precursor additions undergo an exothermic reaction during deposition to form new ceramic products, generally of finer scale than the initial precursor components. The product phases provide reinforcement in the MMC and act as nucleation sites during solidification, resulting in grain refinement, improved mechanical properties, and reduced susceptibility to solidification cracking. Higher tensile strength and stiffness were demonstrated in Al 6061 based MMCs with 2 vol% reinforcement. These combined research efforts demonstrate the potential of PCTW technology for customizing of alloy composition and achieving microstructure control. A manufacturing process for titanium- and aluminum-based PCTWs was developed based on iterative design, theoretical calculations, and empirical data. PCTW technology successfully mitigated Al loss from Ti-6-4 and Mg loss in Al 6061, resulted in grain refinement in Ti-6-4, and demonstrated fabrication of Al based MMCs

Marcia S Domack↗
NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗
NASA NTRSDate not supplied

Setting the Bar for the Replacement of the Probability of Collision Metric

To date, satellite conjunction assessment (CA) risk analysis has largely embraced the probability of collision (Pc) as the omnibus metric to evaluate collision likelihood, and its use in such assessments has mostly been straightforward: at the point at which a conjunction mitigation decision is required, the calculated Pc is compared to a threshold; and if the calculated Pc exceeds that threshold, then a mitigation action is warranted. With only minor variation, this approach is employed by major CA risk assessment centers (e.g., NASA, EUSST, CNES, JAXA) and is advanced as the preferred method in the published CA best practices handbooks. Despite this near unanimity of operational practice, there is a major strain of secondary literature critical of the Pc and willing to propose alternatives. Alfano (2005) pointed out the ability of the Pc to underrepresent the risk in certain situations and counselled a maximum Pc construct. Carpenter (2017, 2019) reiterated this criticism and proposed using instead a confidence interval on the miss distance. Balch et al. (2019) identified what they argued was a defect in the entire Bayesian Pc construct and believed that the use of a more conservative methodology based on covariance ellipsoid overlap was necessary. Delande (2022) introduced the framework of collision “plausibility” to the risk assessment process and sketched out how this might be used operationally. Elkantassi (2022) published a full development of the miss distance confidence interval approach and applied it to several worked examples. While these different approaches to collision risk assessment do differ in their details, they all converge on two central points: first, the Pc’s failure to give an adequate expression of the risk in dilution region situations is a fatal flaw; and second, a conjunction should be presumed risky and in need of mitigation until the evidence of the situation can establish otherwise. These criticisms, if correct, would counsel a number of modifications to current CA operational practice; as such, they force a re-examination of fundamental aspects of the CA problem, including the following: 1. Is the CA risk assessment a probability problem, a statistics problem, or something else? 2. If it is a statistics problem, does it lend itself naturally to a hypothesis test construction? 3. If it can be construed as a hypothesis test, what form should the null hypothesis take, to wit: what constraints exist on the choice of the null hypothesis, what selections are in best alignment with all of the attendant parameters of the problem, and what is implied philosophically by different choices? 4. What are the implications of using the different proposed risk assessment parameters for CA? This question should be answered both in determining how frequently the dilution region situation cited by the critics of the Pc actually appears in an operationally significant manner and the missed detection and false alarm rates of all of the proposed risk assessment metrics, compared both to the Pc and to each other. This paper explores and offers preliminary answers to the above questions, presenting a researched treatment of the philosophical nature of the CA problem and the null hypothesis choice that achieves the greatest consistency with all of the different aspects of operational CA conduct. It then profiles all of the different proposed risk assessment metrics enumerated in the earlier paragraph against an extremely large database of conjunction events at both the 550km and 700km altitudes. The combination of the philosophical exploration of the CA problem and the results of the profiling activity articulates what a risk assessment metric will need to demonstrate, in terms of both innate construction and performance, in order to be a true competitor to the Pc.

conjunction assessment↗
NASA NTRS2018

Orion Optical Navigation Progress Toward Exploration: Mission 1

Optical navigation of human spacecraft was proposed on Gemini and implemented successfully on Apollo as a means of autonomously operating the vehicle in the event of lost communication with controllers on Earth. It shares a history with the "method of lunar distances" that was used in the 18th century and gained some notoriety after its use by Captain James Cook during his 1768 Pacific voyage of the HMS Endeavor. The Orion emergency return system utilizing optical navigation has matured in design over the last several years, and is currently undergoing the final implementation and test phase in preparation for Exploration Mission 1 (EM-1) in 2019. The software development is being worked as a Government Furnished Equipment (GFE) project delivered as an application within the Core Flight Software of the Orion camera controller module. The mathematical formulation behind the initial ellipse fit in the image processing is detailed in Christian. The non-linear least squares refinement then follows the technique of Mortari as an estimation process of the planetary limb using the sigmoid function. The Orion optical navigation system uses a body fixed camera, a decision that was driven by mass and mechanism constraints. The general concept of operations involves a 2-hour pass once every 24 hours, with passes specifically placed before all maneuvers to supply accurate navigation information to guidance and targeting. The pass lengths are limited by thermal constraints on the vehicle since the OpNav attitude generally deviates from the thermally stable tail-to-sun attitude maintained during the rest of the orbit coast phase. Calibration is scheduled prior to every pass due to the unknown nature of thermal effects on the lens distortion and the mounting platform deformations between the camera and star trackers. The calibration technique is described in detail by Christian, et al. and simultaneously estimates the Brown-Conrady coefficients and the Star Tracker/Camera interlock angles. Accurate attitude information is provided by the star trackers during each pass. Figure 1 shows the various phases of lunar return navigation when the vehicle is in autonomous operation with lost ground communication. The midcourse maneuvers are placed to control the entry interface conditions to the desired corridor for safe landing. The general form of optical navigation on Orion is where still images of the Moon or Earth are processed to find the apparent angular diameter and centroid in the camera focal plane. This raw data is transformed into range and bearing angle measurements using planetary data and precise star tracker inertial attitude. The measurements are then sent to the main flight computer's Kalman filter to update the onboard state vector. The images are, of course, collected over an arc to converge the state and estimate velocity. The same basic technique was used by Apollo to satisfy loss-of-comm, but Apollo used manual crew sightings with a vehicle-integral sextant instead of autonomously processing optical imagery. The software development is past its Critical Design Review, and is progressing through test and certification for human rating. In support of this, a hardware-in-the-loop test rig was developed in the Johnson Space Center Electro-Optics Lab to exercise the OpNav system prior to integrated testing on the Orion vehicle. Figure 2 shows the rig, which the test team has dubbed OCILOT (Orion Camera In the Loop Optical Testbed). Analysis performed to date shows a delivery that satisfies an allowable entry corridor as shown in Figure 3.

Holt, Greg N.↗
NASA NTRSDate not supplied

Oxygen Isotopic Composition of Refractory Inclusions from the Miller Range (MIL) 090019 Carbonaceous Chondrite

Introduction: Primitive meteorites contain Calcium-Aluminum-rich Inclusions (CAIs)that preserve the records of the earliest times in the Solar System. CAIs record O isotopic variations with time and/or location from solar (16O-rich) to planetary (16O-poor) compositions [e.g., 1, 2].However, CAIs were also affected by parent body processing to varying degrees. The goal of this work was to identify the processes that resulted in the fine-scale spatial variations in oxygen isotopic compositions within the most ancient Solar System solids. We report oxygen isotopic composition of CAIs from MIL 090019 CO3.1 carbonaceous chondrite. MIL 090019 contains a high abundance of CAIs showing mineralogical as well as textural variations [3]. Methods: We analyzed 10 CAIs with different mineralogical assemblages, classified as corundum-bearing, grossite-bearing, hibonite-bearing, melilite-bearing, anorthite-bearing inclusions and amoeboid olivine aggregates(AOA’s). The mineralogical and petrological characterization of CAIs was performed using the JOEL Hyperprobe 8530 electron microprobe at NASA JSC. Oxygen isotopic imaging of the CAIs was done using the Cameca NanoSIMS 50L ion microprobe at NASA JSC. We followed the analytical protocol described in [4]. O-isotopic maps of the CAI were acquired by rastering a ~3 pA primary Cs+beam at 16 keVover an area of 20 ×20μm2over a period of ~7 hours. Negative secondary ions of 16O−, 17O−, 18O−, 28Si−, 24Mg16O−, 27Al16O−, and 40Ca16O−were simultaneously acquired using electron multiplier detectors at a mass resolving power sufficient to resolve the 16OH−interference from the 17O−peak, where the contribution of 16OH−was <0.1%.An electron flood gun was used to mitigate sample charging during the analyses. We used San Carlos olivine and Madagascar hibonite as isotopic standards to correct for the instrumental mass fractionation. The O-isotopic ratios were corrected for the quasi simultaneous arrival (QSA) effect and the detector dead time. All reported errors are 1 sigma. Results: The MIL 090019 carbonaceous chondrite hosts refractory inclusions varying in their mineralogies and textures[2, 3]. We found the oxygen isotopic compositions of mineral phases in MIL 090019 CAIs show large variations with Δ17O values varying from -27 to 0‰. A corundum-bearing inclusion records isotopic variations with corundum, melilite, perovskite, and anorthite showing Δ17Ovalues of -20 ±2.5, -10.2 ± 6.5, -10.7 ± 1.3,and -1.2 ± 3.2‰ respectively. A grossite-hibonite-bearing inclusion also records O-isotopic heterogeneity, where grossite, hibonite, Mg-rich spinel, and melilite show -3.7 ± 2.7, -22.9 ± 2.5, -21.4 ± 2.6, and -17.3 ± 2.8‰ respectively, whereas perovskite grains record a heterogenous isotopic composition ranging from -22.6 to -13.9‰. A hibonite-bearing inclusion shows average Δ17Ο values of -17.6 ± 5.3 ‰ and -15.1 ± 4.8 ‰ for hibonite and spinel respectively. The two perovskite grains also show relatively 16O-poor composition (Δ17Ο = -11.8 ± 6.4 ‰ and -8.6 ± 5.2 ‰). Anorthite has ~0.1 wt% Na2O and ~0.52 wt% FeO and a Δ17Ο value of -7.2 ± 5.0 ‰, whereas diopside records a Δ17Ο value of -3.5 ± 5.0 ‰. A hibonite-spinel-rich inclusion contains spinel and hibonite that are16O-rich (Δ17Οvalues -21.3 ± 2.8‰ and -19.4 ± 3.0‰, respectively). Melilite shows heterogenous oxygen isotopic composition, with an inner 16O-rich region and an outer 16O-poor region separated by a sharp boundary, with a bulk Δ17Ο= -16.2 ± 2.1‰, whereas two perovskite grains are relatively 16O-poor (Δ17Ο= -13.2 ± 3.7‰ and -7.1 ± 3.4‰).A spinel-rich inclusion contains perovskite grains with Δ17Οranging from -27.0 to -22.8 ‰and spinel with average Δ17Οvalue of -21.4 ± 2.8 ‰. In a melilite and spinel-rich CAI, perovskite, spinel, diopside, melilite, and Al-Ti-rich pyroxene record Δ17Ο values of-26.9 ± 4.6‰, -22.4 ± 4.3‰,-23.4 ± 4.9, -19.3 ± 3.4, and -19.0 ± 4.3 ‰ respectively. Olivine and spinel components of 2 AOA’s record a 16O-rich composition with Δ17O~-23‰. Discussion: The coexistence of 16O-rich and 16O-poor minerals within CAIs has been attributed to differing degrees of O isotopic exchange, reflecting a wide range of O diffusion rates[5]. However, the O isotopic heterogeneity observed in these CAIs cannot be explained by oxygen diffusion in nebular or parent body settings alone. The isotopic imaging of CAIs from the MIL 090019 meteorite shows variations in the O isotopic composition recorded in their mineral components, suggesting that these CAIs record the oxygen isotopic heterogeneity in the nebular gas from which they condensed. The coexistence of 16O-rich and 16O-poor mineral components in the CAI-forming region suggests that the gas in this region was not well-mixed and consisted of distinct 16O-rich and 16O-poor gaseous reservoirs.

P Mane↗
NASA NTRS2018

Towards an Enhanced Droplet Activation Scheme for Multi-Moment Bulk Microphysics Schemes

Initial droplet spectra produced upon activation impact the ensuing chain of microphysical processes andtherefore play a crucial role in cloud evolution. This work re-examines dependencies of newly formed clouddroplet size distribution (CDSD) characteristics on environmental and aerosol properties via parcel model simulationsthat serve as the basis for a multi-moment bulk microphysics droplet activation scheme suitable for acloud-resolving model (CRM). It is found that applying a fixed size threshold to define activated droplets versusemploying physical considerations can lead to erroneous activation and overly broad CDSDs for high aerosolconcentration and weak updraft conditions. Aerosol distributions characterized by larger median sizes and/orincreased solubility can result in greater activated droplet numbers, whereas impacts of these parameters onCDSD spectral width depend on both aerosol number concentration and updraft velocity. An expansion of theactivation scheme to include CDSD spectral width is proposed to aid efforts to extend high-order momentprediction to cloud droplet categories in CRMs as well as better represent variability in the activation process onthe cloud scale.simulations to investigate the regime dependence of the relative dispersion(d)1 of newly activated CDSDs, where d is the ratio of dropletradius standard deviation (σ) to the mean radius (r ). C16 demonstratedthat increasing Na resulted in increasing (decreasing) d values via reducedcondensational narrowing (spectral broadening) rates within theAL (UL) regime, with d values peaking in the TR regime. Their findingssuggest a similar regime dependence for d as R09 noted for Nc and helpexplain reportedly conflicting relationships between Na and CDSDspectral characteristics (cf. Hudson and Noble, 2014; Liu et al., 2014),although the applicability of these results within bulk microphysicalschemes was not addressed.Simulating aerosol-cloud interactions with CRMs employing bulkmicrophysics requires that the scheme minimally predict two CDSDparameters, namely mass and number concentrations, and represent thedroplet activation process. Various activation schemes aim to determineNc from aerosol and environmental properties and include analyticalexpressions (e.g., Abdul-Razzak et al., 1998; Morrison et al., 2005) aswell as lookup tables (LUTs) based on detailed parcel model calculations(e.g., Saleeby and Cotton, 2004, hereafter SC04; Segal and Khain,2006; Thompson and Eidhammer, 2014). Expressions to diagnose CDSDspectral width from Nc (Grabowski, 1998; Liu et al., 2006; Morrison andGrabowski, 2007) or cloud water content (Geoffroy et al., 2010) havealso been developed, although more robust methods to obtain CDSDspectral width upon activation are presently lacking. This latter point isrelevant for triple-moment (3 M) bulk microphysics that aim to predictdistribution spectral width alongside number and mass concentrations(e.g., Loftus et al., 2014; Milbrandt and Yau, 2005).The current work extends the findings of C16 to the current LUTbasedaerosol activation scheme used in the Regional AtmosphericModeling System (RAMS) (Cotton et al., 2003; SC04; Saleeby and vanden Heever, 2013, hereafter SvdH13) and additionally examinesaerosol size and solubility impacts on newly activated CDSD properties.Because early cloud development processes such as condensationalgrowth, evaporation, and droplet self-collection depend on and impactCDSD spectral width (Hudson and Yum, 1997; Seifert and Beheng 2001;Lu and Seinfeld, 2006; Igel and van den Heever, 2017), an expansion ofthe activation LUTs to include CDSD spectral width is proposed as apreliminary step for extending 3M prediction to CDSDs in CRMs forimproved simulations of aerosol-cloud interactions.2. MethodologyThe current RAMS two-moment microphysics module determinesthe fractional number of aerosol particles that activate to cloud dropletsfrom five-dimensional LUTs based on model predicted air temperature(T), w, Na, and the geometric median radius (rg) and soluble fraction (ε)of the aerosol size distribution (SvdH13). These LUTs are created offlineusing a one-dimensional Lagrangian adiabatic parcel model (Feingoldand Heymsfield, 1992; Heymsfield and Sabin, 1989; SC04) to simulateexplicit droplet activation and initial CDSD growth for a range of ambientatmospheric conditions [T, w] and binned lognormal aerosol sizedistributions given by= ⎡⎣ ⎢− ⎤⎦ ⎥N r Nr π σr rσ( )2 lnexp[ln( / )]2(ln )aggg22 (1)where r is the dry aerosol particle bin radius and σg is the geometricstandard deviation of the distribution. As the parcel model simulationsfocus on the activation process, other processes such as coalescence,sedimentation, and mixing are not considered. Details of the parcelmodel can be found in SC04 and SvdH13, and only a brief description isprovided here. At the onset of parcel model calculations, the initiallydry aerosol particles in all bins first deliquesce and reach theirequilibrium diameters in a sub-saturated environment based on theKöhler equation for solution droplets. The parcel is then lifted at a fixedupward velocity w, and particle growth by vapor diffusion, along withconcurrent changes in the ambient environment, are iteratively computedusing the Variable-coefficient Ordinary Differential Equation(VODE) solver (Brown et al., 1989). The time resolution of these calculationsis determined within the VODE solver, and the frequency atwhich the solver is called is controlled by a longer model time stepbased on fixed upward parcel displacement increments (Δz) at thespecified w (Δt=Δz/w). Model calculations proceed until the parcelreaches a height 50m beyond the level of maximum saturation ratio(Smax) or total parcel displacement exceeds 2 km. Upon model termination,Smax and the fractional number of aerosols (factv) resulting innewly formed cloud droplets, defined as particles having diameters of atleast 2 μm, are cataloged in the LUTs according to the specified T, w, Na,rg, and ε parameter values.A critical point regarding the creation of these LUTs is the use of afixed minimum diameter (Dmin) to define cloud droplets in the parcelmodel, which can produce erroneous CDSD characteristics, particularlywithin the UL regime. For aerosol distributions with large rg valuesunder low SS conditions, for example, deliquesced aerosols within thelarge tail of the distribution can exceed 2 μm in diameter yet remainunactivated as ‘haze’ particles (Levin and Cotton, 2009; McFigganset al., 2006). For this study, aerosol particles activate to cloud dropletsbased on the critical diameter Dcrit as a function of parcel supersaturationratio (Sr) as in R09:D = σ MS RTρ83 ln( ) critsol wr w (2)where σsol is the surface tension of a solution droplet, Mw and ρw are themolar mass and density of liquid water, respectively, and R is theuniversal gas constant. Additionally, at relatively large w values withinthe AL regime, Nc stabilizes shortly after reaching supersaturation.However, parcel ascent and condensational growth continue beyondthe level of Smax, potentially causing additional narrowing of the CDSD.In the current work, model calculations terminate upon reaching Smaxas changes in Nc are negligible with continued ascent (Peng et al., 2007;R09).Parcel model simulations are performed to examine the sensitivitiesof CDSD characteristics to w, Na, rg, and ε, with the ranges for theseparameters listed in Table 1. Aerosols are assumed to be a mix of solubleand insoluble material of equal density, specified by ε, where fullysoluble aerosols correspond to ammonium sulfate with hygroscopicityparameter κ=0.61 (Petters and Kreidenweis, 2007). FollowingSvdH13, aerosol geometric standard deviation is fixed at σg=1.8, andaerosol distributions (Eq. 1) are partitioned into 100 logarithmicallyspacedbins spanning a size range specific to each rg value. For all simulations,Δz=1 m, and initial values of relative humidity, air temperatureand pressure are set to RH=0.99, T=10 °C and p=900 hPa,respectively.

Loftus, Adrian M.↗