Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Warm Start”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Enhancing Active Distribution Systems Resilience by Fully Distributed Self-Healing Strategy

Distributed restoration can exploit smart grid technologies to enhance the resilience of active distribution networks toward a self-healing smart grid. However, the large number of decision variables, especially the binary ones for reconfiguration, bring challenges to developing scalable distributed distribution service restoration (DDSR) strategies. This paper proposes a fully distributed solution procedure based on the alternating direction method of multipliers (ADMM) for mixed-integer programming problems and applies to develop the DDSR framework. The method consists of relax-drive-polish phases, 1) relaxing binary variables, and applying the convex ADMM as a warm start; 2) driving the solutions toward Boolean values through a proximal operator; 3) fixing the obtained binding binary variables and solving the rest of the problem to polish results and achieve a high-quality suboptimal solution. Then, an autonomous clustering strategy and consensus ADMM are integrated with the proposed method to realize the fully distributed cluster-based framework of DDSR. This framework can first determine DER scheduling and switch status for reconfiguration to energize the out-of-service areas from local faults, and then provide the load restoration solution in a distributed manner for total blackouts in large-scale distribution networks. Furthermore, the effectiveness and scalability of the proposed DDSR framework are demonstrated through testing on the IEEE 123-node, IEEE 8500-node, and synthetic 100k-node test feeders.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Inverse Design of Two-Dimensional Airfoils Using Conditional Generative Models and Surrogate Log-Likelihoods

Abstract This paper shows how to use conditional generative models in two-dimensional (2D) airfoil optimization to probabilistically predict good initialization points within the vicinity of the optima given the input boundary conditions, thus warm starting and accelerating further optimization. We accommodate the possibility of multiple optimal designs corresponding to the same input boundary condition and take this inversion ambiguity into account when designing our prediction framework. To this end, we first employ the conditional formulation of our previous work BézierGAN–Conditional BézierGAN (CBGAN)—as a baseline, then introduce its sibling conditional entropic BézierGAN (CEBGAN), which is based on optimal transport regularized with entropy. Compared with CBGAN, CEBGAN overcomes mode collapse plaguing conventional GANs, improves the average lift-drag (Cl/Cd) efficiency of airfoil predictions from 80.8% of the optimal value to 95.8%, and meanwhile accelerates the training process by 30.7%. Furthermore, we investigate the unique ability of CEBGAN to produce a log-likelihood lower bound that may help select generated samples of higher performance (e.g., aerodynamic performance). In addition, we provide insights into the performance differences between these two models with low-dimensional toy problems and visualizations. These results and the probabilistic formulation of this inverse problem justify the extension of our GAN-based inverse design paradigm to other inverse design problems or broader inverse problems.

Engineering↗
DOE OSTI2026

matsim-agents v1.0

matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗
DOE OSTI2023

A Surrogate-Based Asynchronous Decomposition Technique for Realistic Security-Constrained Optimal Power Flow Problems

Here we present a decomposition approach for obtaining good feasible solutions for the security-constrained, alternating-current, optimal power flow (SC-AC-OPF) problem at an industrial scale and under real-world time and computational limits. The approach was designed while preparing and participating in ARPA-E’s Grid Optimization Competition (GOC) Challenge 1. The challenge focused on a near-real-time version of the SC-AC-OPF problem, where a base operating point is optimized, taking into account possible single-element contingencies, after which the system adapts its operating point following the response of automatic frequency droop controllers and voltage regulators. Our solution approach for this problem relies on state-of-the-art nonlinear programming algorithms, and it employs nonconvex relaxations for complementarity constraints, a specialized two-stage decomposition technique with sparse approximations of recourse terms and contingency ranking and prescreening. The paper describes and justifies our approach and outlines the features of its implementation, including functions and derivatives evaluation, warm-starting strategies, and asynchronous parallelism. We discuss the results of the independent benchmark of our approach by ARPA-E’s GOC team in Challenge 1, where it was found to consistently produce high-quality solutions across a wide range of network sizes and difficulty, and conclude by outlining future extensions of the approach.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1983

The first winter of MAP-dynamics, 1982-1983: A winter with three different warming periods

A synoptic description is given for the first winter of the Middle Atmosphere Program's atmospheric dynamics efforts for 1982-83. After a cold early winter period three warming periods were observed at the end of December, the end of January, and the end of February. In March a pronounced late winter cooling occurred in the upper stratosphere, whereas in the lower levels the final warming started slowly, but was not accomplished before mid-April. The momentum budget, calculated from the daily height and temperature charts, is discussed in terms of the divergence of the Eliassen-Palm-Vector.

Naujokat, B.↗
NASA NTRS1993

Mars: The initial emplacement of ground ice in response to the thermal evolution of its early crust

Given the geomorphic evidence for the widespread occurrence of water and ice in the early martian crust, and the difficulty involved in accounting for this distribution given the present climate, it has been suggested that the planet's early climate was originally more Earth-like, permitting the global emplacement of crustal H2O by direct precipitation as snow or rain. The resemblance of the martian valley networks to terrestrial runoff channels, and their almost exclusive occurrence in the planet's ancient (approximately 4 billion year old) heavily cratered terrain, is often cited as evidence of just such a period. An alternative school of thought suggests that the early climate did not differ substantially from that of today. Advocates of this view find no compelling reason to invoke a warmer, wetter period to explain the origin of the valley networks. Rather, they cite evidence that the primary mechanism of valley formation was ground water sapping, a process that does not require that surface water exist in equilibrium with the atmosphere. However, while sapping may successfully explain the origin of the small valleys, it fails to address how the crust was initially charged with ice as the climate evolved towards its present state. Therefore, given the uncertainty regarding the environmental conditions that prevailed on early Mars, the initial emplacement of ground ice is considered from two perspectives: (1) that the early climate started warm and wet, but gradually cooled with time; and (2) that it never differed substantially from that of today.

Clifford, Stephen M.↗
NASA NTRS1993

The evolution of the early Martian climate and the initial emplacement of crustal H2O

Given the geomorphic evidence for the widespread occurrence of water and ice in the early Martian crust, and the difficulty involved in accounting for this distribution given the present climate, it has been suggested that the planet's early climate was originally more Earth-like, permitting the global emplacement of crustal H2O by direct precipitation as snow or rain. The resemblance of the Martian valley networks to terrestrial runoff channels and their almost exclusive occurrence in the planet's ancient (approximately 4-b.y.-old) heavily cratered terrain are often cited as evidence of just such a period. An alternative school of thought suggests that the early climate did not differ substantially from that of today. Advocates of this view find no compelling reason to invoke a warmer, wetter period to explain the origin of the valley networks. Rather, they cite evidence that the primary mechanism of valley formation was groundwater sapping, a process that does not require that surface water exists in equilibrium with the atmosphere. However, while sapping may successfully explain the origin of the small valleys, it fails to address how the crust was initially charged with ice as the climate evolved towards its present state. Therefore, given the uncertainty regarding the environmental conditions that prevailed on early Mars, the initial emplacement of ground ice is considered here from two perspectives: (1) the early climate started warm and wet, but gradually cooled with time, and (2) the early climate never differed substantially from that of today.

Clifford, S. M.↗
NASA NTRS2003

Applications of Space-Filling-Curves to Cartesian Methods for CFD

The proposed paper presents a variety novel uses of Space-Filling-Curves (SFCs) for Cartesian mesh methods in 0. While these techniques will be demonstrated using non-body-fitted Cartesian meshes, most are applicable on general body-fitted meshes -both structured and unstructured. We demonstrate the use of single O(N log N) SFC-based reordering to produce single-pass (O(N)) algorithms for mesh partitioning, multigrid coarsening, and inter-mesh interpolation. The intermesh interpolation operator has many practical applications including warm starts on modified geometry, or as an inter-grid transfer operator on remeshed regions in moving-body simulations. Exploiting the compact construction of these operators, we further show that these algorithms are highly amenable to parallelization. Examples using the SFC-based mesh partitioner show nearly linear speedup to 512 CPUs even when using multigrid as a smoother. Partition statistics are presented showing that the SFC partitions are, on-average, within 10% of ideal even with only around 50,000 cells in each subdomain. The inter-mesh interpolation operator also has linear asymptotic complexity and can be used to map a solution with N unknowns to another mesh with M unknowns with O(max(M,N)) operations. This capability is demonstrated both on moving-body simulations and in mapping solutions to perturbed meshes for finite-difference-based gradient design methods.

Aftosmis, Michael J.↗
NASA NTRS2003

Applications of Space-Filling-Curves to Cartesian Methods for CFD

This paper presents a variety of novel uses of space-filling-curves (SFCs) for Cartesian mesh methods in CFD. While these techniques will be demonstrated using non-body-fitted Cartesian meshes, many are applicable on general body-fitted meshes-both structured and unstructured. We demonstrate the use of single theta(N log N) SFC-based reordering to produce single-pass (theta(N)) algorithms for mesh partitioning, multigrid coarsening, and inter-mesh interpolation. The intermesh interpolation operator has many practical applications including warm starts on modified geometry, or as an inter-grid transfer operator on remeshed regions in moving-body simulations Exploiting the compact construction of these operators, we further show that these algorithms are highly amenable to parallelization. Examples using the SFC-based mesh partitioner show nearly linear speedup to 640 CPUs even when using multigrid as a smoother. Partition statistics are presented showing that the SFC partitions are, on-average, within 15% of ideal even with only around 50,000 cells in each sub-domain. The inter-mesh interpolation operator also has linear asymptotic complexity and can be used to map a solution with N unknowns to another mesh with M unknowns with theta(M + N) operations. This capability is demonstrated both on moving-body simulations and in mapping solutions to perturbed meshes for control surface deflection or finite-difference-based gradient design methods.

Aftosmis, M. J.↗
NASA NTRS2019

Partition-based Feasible Integer Solution Pre-computation for Hybrid Model Predictive Control

For multiparametric mixed-integer convex programming problems such as those encountered in hybrid model predictive control, we propose an algorithm for generating a feasible partition of a subset of the parameter space. The result is a static map from the current parameter to a suboptimal integer solution such that the remaining convex program is feasible. Convergence is proved with a new insight that the overlap among the feasible parameter sets of each integer solution governs the partition complexity. The partition is stored as a tree which makes querying the feasible solution efficient. The algorithm can be used to warm start a mixed integer solver with a real-time guarantee or to provide a reference integer solution in several suboptimal MPC schemes. The algorithm is tested on randomly generated systems with up to six states, demonstrating the effectiveness of the approach.

Bayard, David S.↗
NASA NTRSDate not supplied

A Reynolds-Averaged Navier-Stokes Perspective for the High Lift Common Research Model Using the LAVA Framework

An in-depth assessment of the High Lift-Common Research Model (HL-CRM) aerodynamic performance is made using a Reynolds-Averaged Navier-Stokes (RANS) methodology. Results presented here were produced in support of the 4th AIAA High Lift Prediction Workshop(HLPW-4), with additional simulation results and analysis used to support findings from the workshop. Mesh resolution, RANS closure models, and tunnel effects were among the sources of simulation sensitivity addressed, all in the pursuit of evaluating best-practice simulation methods for complex high-lift configurations. In general, free-air simulations using RANS yield accurate predictions of integrated loads in the pre-stall regime for the HL-CRM, while strong sensitivity to simulation method is observed near CLmax. Additionally, a study assessing changes in aerodynamic performance due to flap deflection angle indicated strong sensitivity to turbulence model closure and mesh resolution, particularly at the highest deflection setting. Several available corrections for the Spalart-Allmaras (SA) turbulence model were tested and generally introduced more error relative to the experiment in comparison to the baseline SA closure model. The Standard Menter Baseline Two-Equation (k-ωBSL) turbulence model was also explored, which arguably yielded results highly consistent with experiment, but at a substantially higher computational cost. Therefore, as a compromise between computational cost and accuracy, the baseline SA closure model was selected for additional HL-CRM studies that assess alternative solution methods. These strategies include “warm start” approaches, which initialize the simulation from a previously converged solution at a nearby angle of attack, unsteady calculations, and simulations that model the experimental wind tunnel. Results from these studies would reinforce the shortcomings of RANS in predicting aircraft performance at high-lift conditions, namely at high-αand highly loaded flap conditions, while also answering questions about appropriate modeling techniques to minimize the errors associated with RANS where possible

TTT↗
NASA NTRSDate not supplied

Combined Bernstein Polynomial, Optimal Reciprocal Collision Avoidance, Differential Dynamic Programming for Trajectory Replanning and Collision Avoidance for UAM Vehicles

This paper presents an integration of Differential Dynamic Programming (DDP) with the Optimal Reciprocal Collision Avoidance (ORCA) algorithm as the basis for a new algorithm, titled Combined Bernstein Polynomial Optimal Reciprocal Collision Avoidance DDP (COBRA-DDP), for trajectory replanning and collision avoidance for Urban Air Mobility (UAM) vehicles. State-constrained variants of DDP provide the ability to plan trajectories while avoiding obstacles, but these methods require a large increase in computational time per iteration which hinders the overall speed of the algorithm. ORCA utilizes simplified dynamics to recognize potential collisions along a trajectory and provides an optimal velocity for the avoidance of multiple vehicles. These velocity commands, however, may not result in a dynamically feasible trajectory for DDP to plan around. As such, a Bernstein polynomial curve that considers the general dynamic constraints of the vehicle is generated to approximate a trajectory based on the velocity commands. COBRA-DDP optimizes this suggested trajectory via unconstrained DDP to provide a dynamically feasible trajectory that provides collision avoidance. This new trajectory can be applied to the vehicle or used to warm start the state constrained DDP algorithms to decrease computation time. Its benefits and effectiveness of the algorithm are demonstrated on a UAM Vertical Takeoff and Landing (VTOL) vehicle simulation with highly nonlinear dynamics.

Optimal Reciprocal Collision Avoidance↗
DOE OSTI2026

Distributed Quantum-Enhanced Optimization: A Topographical Preconditioning Approach for High-Dimensional Search

Optimization problems become fundamentally challenging as the number of variables increases. Because the volume of the search space grows exponentially, classical algorithms frequently fail to locate the global minimum of non-convex functions. While quantum optimization offers a potential alternative, mapping continuous problems onto near-term quantum hardware introduces severe scaling limits and barren plateaus. To bridge this gap, we propose the Distributed Quantum-Enhanced Optimization (D-QEO) framework. Instead of forcing the quantum processor to find the exact minimum, we use it simply as a topographical preconditioner. The QPU maps the landscape to locate the most promising basin of attraction, generating high-quality seed points for a classical GPU-accelerated solver to refine. To make this approach viable for utility-scale problems, we exploit the mathematical structure of separable functions. This allows us to cut a 50-qubit (i.e., $2^{50}$) global search space into independent and manageable sub-spaces using 5-qubit subcircuits. By executing these fragments concurrently with CUDA-Q, we completely bypass the overhead of cross-register entanglement and classical tensor knitting for separable functions. Benchmarks on the 10-dimensional Rastrigin and Ackley functions show that D-QEO prevents the exponential failure rates observed in purely classical algorithms. Furthermore, this quantum warm-start significantly reduces the number of classical BFGS iterations required to converge, providing a highly practical blueprint for utilizing near-term quantum resources in complex global search.

Soos, Dominik [Old Dominion U.]↗
DOE OSTI2026

Surrogates for Valve-Controlled Pipe Flow: Accelerating Nuclear Reactor Design

Neural surrogate models are developed to replace expensive steady-state RANS CFD simulations for valve-controlled pipe flow in nuclear reactor design. Using parametric CFD data generated with MOOSE Pronghorn across a range of valve geometry and flow conditions, three approaches are compared: a POD-based reduced-order model, a structured UNet on a cylindrical grid, and unstructured models (DeepONet and BiStride MeshGraphNet) on nondimensionalized point clouds. POD achieves the highest accuracy (99%) with fast inference but requires storing all solution snapshots, while the DeepONet and BSMS-GNN both achieve ~89% accuracy at sub-second inference, with the BSMS-GNN offering superior geometric generalizability. These surrogates enable rapid ranking of candidate valve designs and can warm-start CFD solvers to accelerate convergence, supporting agentic design iteration on the Prometheus platform.

42 - ENGINEERING↗
DOE OSTI2022

Improving Unit Flexibility Utilizing Plasma Ignitors

Coal fired steam generation operators face increasing market challenges, including competition from low cost generation, renewables, and regulatory pressure. These market conditions are forcing utilities to operate their coal assets in a more flexible mode, including more frequent starts and stops, faster ramp rates, frequent cycling, and extended operation at the lowest possible loads. Without upgrades to firing and control systems, pressure parts, and auxiliary systems, CAPEX and OPEX costs will increase significantly. Typically, expensive support fuels (Oil or Gas) are required to maintain safe, stable coal ignition when operating below 25% to 35% MCR (Maximum Continuous Rating). The exact minimum load without support fuel differs depending upon unit design specifications and fuel being burned. Frequent starts also increase the use of support fuels. The firing and burner control systems must be designed to allow for proper operation over the unit’s load range, ensuring that all mechanical components operate properly, and proper control of air and fuel streams can be maintained. To address the need for frequent starts and the need for stable low load operation on coal fired steam generators, this paper discusses options available, including, Firing Systems upgrades, Digital solutions, and pressure part upgrades. Several case studies are included that highlight options available allowing units to stay in operation and reduce CAPEX and OPEX costs.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Plasma Ignition and Combustion Stabilization Technology to Improve Flexible Operation, Reliability and Economics of an Existing Coal-Fired Boiler

GE Steam Power, Inc. (GE) proposed to improve reliability, flexibility, and economics of an existing coal-fired power plant by applying a new advanced technology developed by GE, a plasma-assisted pulverized fuel firing system. The objective of this program is to demonstrate the achievement of lower load by improved flame stabilization and therefore lower operating costs in a full-scale field installation at coal-fired electric utility. GE’s Plasma Ignition and Combustion Stabilization System is designed to operate continuously to support low load operation. With the plasma on, the flame will be attached and stable, removing the firing system as a limitation to low load operation. In addition, GE’s exclusively from ABENZ company licensed AC based technology has a 90+% system efficiency compared to all other systems at which are DC and operate with ~75% efficiency. Plant operating costs are lowered by eliminating use of expensive support fuel as well as the ability to operate at lower loads. The utilities’ ability to better match the demand curve will result in significant savings. Maintenance is lower for an AC system than a DC system as it operates at lower current. This eliminates the need for a demineralized cooling water system and provides longer electrode life which translates into both material and labor savings. It is the objective of GE to not only demonstrate the additional low load achievable with a plasma system after best achievable tuning, sensor and software approach has been exhausted, but also the increased stability of the flame at all loads with plasma assistance as well as cost savings at all low loads using plasma instead of oil. Upon successful completion of this project, GE will have sufficient field experience to rapidly deploy the Plasma Technology. The project objectives were achieved through the implementation of plasma ignitor technology at PacifiCorp Hunter Station Unit 3. A plasma ignitor system was retrofitted on ten wall-fired burners, Mill 3-4 combustion system. The Plasma Ignitors installed at Hunter proved that this GE technology is a direct and complete replacement for the original oil ignitors. The Hunter Unit 3 burner management system allows the plasma system to be used in all applications that originally required oil to be burned. This includes any time the Mill 3-4 is started or stopped for any reason including boiler starts, load changes, and low load support.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

A Modeling-Based Flammable Risk Treatment of Refrigerant Leakage from a Commercial R-290 Refrigeration Machine

Because of serious concerns about global warming, manufacturers have started phasing out high global warming potential (GWP) refrigerants in commercial refrigeration equipment (e.g., R-134a). As a potential replacement, propane (R-290) is an environmentally friendly refrigerant for commercial refrigeration equipment because its GWP is only three. However, propane is flammable and is therefore classified as a Class A3 refrigerant per ASHRAE Standards, so safety is a very important consideration when propane-based equipment is designed and deployed in buildings. In the event of a refrigerant leak, flammability of the refrigerant depends on the refrigerant’s local concentration, which is highly affected by the indoor air environment, including temperature and air flow. In this study, a ventilation system attached to a commercial R-290 refrigeration device was designed to eliminate the flammability risk. Moreover, a computational fluid dynamics (CFD) model was developed to investigate the refrigerant leak, thereby evaluating effects of the ventilation system. The CFD model can visualize the flammable zones owing to the leak.

42 ENGINEERING↗
NASA NTRS2000

Cirrus Parcel Model Comparison Project

The Cirrus Parcel Model Comparison (CPMC) is a project of the GEWEX Cloud System Study Working Group on Cirrus Cloud Systems (GCSS WG2). The primary goal of this project is to identify cirrus model sensitivities to the state of our knowledge of nucleation and microphysics. Furthermore, the common ground of the findings may provide guidelines for models with simpler cirrus microphysics modules. We focus on the nucleation regimes of the warm (parcel starting at -40 C and 340 hPa) and cold (-60 C and 170 hPa) cases studied in the GCSS WG2 Idealized Cirrus Model Comparison Project. Nucleation and ice crystal growth were forced through an externally imposed rate of lift and consequent adiabatic cooling. The background haze particles are assumed to be lognormally-distributed H2SO4 particles. Only the homogeneous nucleation mode is allowed to form ice crystals in the HN-ONLY runs; all nucleation modes are switched on in the ALL-MODE runs. Participants were asked to run the HN-lambda-fixed runs by setting lambda = 2 (lambda is further discussed in section 2) or tailoring the nucleation rate calculation in agreement with lambda = 2 (exp 1). The depth of parcel lift (800 m) was set to assure that parcels underwent complete transition through the nucleation regime to a stage of approximate equilibrium between ice mass growth and vapor supplied by the specified updrafts.

Lin, Ruei-Fong↗