Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spreading models”

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

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

At least 559 records · Page 31

The Impact of Parameterized Lateral Mixing on the Antarctic Circumpolar Current in a Coupled Climate Model

This study examines the impact of changing the lateral diffusion coefficient ARedi on the transport of the Antarctic Circumpolar Current (ACC). The lateral diffusion coefficient ARedi is poorly constrained, with values ranging across an order of magnitude in climate models. The ACC is difficult to accurately simulate, and there is a large spread in eastward transport in the Southern Ocean (SO) in these models. This paper examines how much of that spread can be attributed to different eddy parameterization coefficients. A coarse-resolution, fully coupled model suite was run with A Redi = 400, 800, 1200, and 2400 m 2 s -1 . Additionally, two simulations were run with two-dimensional representations of the mixing coefficient based on satellite altimetry. Relative to the 400 m 2 s -1 case, the 2400 m 2 s -1 case exhibits 1) an 11% decrease in average wind stress from 50° to 65°S, 2) a 20% decrease in zonally averaged eastward transport in the SO, and 3) a 14% weaker transport through the Drake Passage. The decrease in transport is well explained by changes in the thermal current shear, largely due to increases in ocean density occurring on the northern side of the ACC. In intermediate waters these increases are associated with changes in the formation of intermediate waters in the North Pacific. We hypothesize that the deep increases are associated with changes in the wind stress curl allowing Antarctic Bottom Water to escape and flow northward.

54 ENVIRONMENTAL SCIENCES↗

3D Material Response of the MSL Heatshield Using NuSil-Coated PICA

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield that used NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performance for protecting probes during planetary entry. Data from the Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. A space-grade silicone-based coating commercially known as NuSil CV-1144-0 [3] was applied to the entire MSL heatshield, including the MEDLI plugs, to mitigate the spread of dust from PICA. Modeling the thermal response of PICA-NuSil (PICA-N) system is still an open challenge. Ground testing of PICA-N models exhibited surface temperature jumps of the order of 150 K due to oxide scale formation and sub-sequent NuSil burn-off. It is therefore critical to include a validated model for the material response of the coating in engineering codes. A test campaign has been conducted at the NASA’s Langley HyMETS [4] facility to screen the response of PICA-N and gather detailed data on its behavior [5]. A first model of PICA-N thermal response has been developed using the Hy-METS experiments [6]. The objective of this work is to analyze the material response of the latest PICA-N model compared to the engineering model used to simulate the entry of MSL. The environment and material response around the MSL aeroshell during Mars atmospheric entry is simulated using a collection of tools. The Direct Simulation Monte Carlo SPARTA code [7] is used in the rarefied regime, the Data Parallel Line Relaxation (DPLR) code [8] is used in the continuum regime and radiative heating conditions are provided by the Nonequilibrium air radiation (NEQAIR) code [9] to estimate the environmental conditions. The thermal response inside the material is computed using the Porous material Analysis Toolbox based on Open-FOAM (PATO) [10,11,12]. Thermodynamic and chemistry properties are estimated using the Mutation++ library [13]. The approach implemented in PATO as a first cut PICA-N thermal response model is outlined in Figure 1. While the recession is less than the coating thickness, the Surface mass and energy balance Boundary Condition (SBC) uses the NuSil B’ tables. Once the recession removes the coating, the usual PICA B’ tables are used for the SBC. The B’ tables are computed using an equilibrium solver implemented in Mutation++, given the temperature, pressure, blowing rate, composition of the pyrolysis and environment gases, and the condensed species at the surface. Preliminary results of the 3D material response of the MSL heat-shield at the peak heating (80 sec after Entry Interface) are shown in Figure 2. Current NASA’s mission to Mars, Mars 2020, used the spare heatshield of MSL for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the PATO high-fidelity material response capability was benchmarked against flight data from MEDLI. This effort represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes. This bench-marking is awaiting the final release of the MEDLI-2 data.

Aerospace↗

Quantifying the Effect of Mixing on the Mean Age of Air in CCMVal-2 and CCMI-1 Models

The stratospheric age of air (AoA) is a useful measure of the overall capabilities of a general circulation model (GCM) to simulate stratospheric transport. Previous studies have reported a large spread in the simulation of AoA by GCMs and coupled chemistry-climate models (CCMs). Compared to observational estimates, simulated AoA is mostly too low. Here we attempt to untangle the processes that lead to the AoA differences between the models and between models and observations. AoA is influenced by both mean transport by the residual circulation and two-way mixing; we quantify the effects of these processes using data from the CCM inter-comparison projects CCMVal-2 (Chemistry-Climate Model Validation Activity 2) and CCMI-1 (Chemistry-Climate Model Initiative, phase 1). Transport along the residual circulation is measured by the residual circulation transit time (RCTT). We interpret the difference between AoA and RCTT as additional aging by mixing. Aging by mixing thus includes mixing on both the resolved and subgrid scale. We find that the spread in AoA between the models is primarily caused by differences in the effects of mixing and only to some extent by differences in residual circulation strength. These effects are quantified by the mixing efficiency, a measure of the relative increase in AoA by mixing. The mixing efficiency varies strongly between the models from 0.24 to 1.02. We show that the mixing efficiency is not only controlled by horizontal mixing, but by vertical mixing and vertical diffusion as well. Possible causes for the differences in the models' mixing efficiencies are discussed. Differences in subgrid-scale mixing (including differences in advection schemes and model resolutions) likely contribute to the differences in mixing efficiency. However, differences in the relative contribution of resolved versus parameterized wave forcing do not appear to be related to differences in mixing efficiency or AoA.

CCMI↗

Identifying mitigation strategies for $\mathrm{COVID}$-$\mathrm{19}$ superspreading on flights using models that account for passenger movement

Despite commercial airlines mandating masks, there have been multiple documented events of COVID-19 superspreading on flights. Conventional models do not adequately explain superspreading patterns on flights, with infection spread wider than expected from proximity based on passenger seating. An important reason for this is that models typically do not consider the movement of passengers during the flight, boarding, or deplaning. Understanding the risks for each of these aspects could provide insight into effective mitigation measures. We modeled infection risk from seating and fine-grained movement patterns – boarding, deplaning, and inflight movement. We estimated infection model parameters from a prior superspreading event. We validated the model and the impact of interventions using available data from three flights, including cabin layout and seat locations of infected and uninfected passengers, to suggest interventions to mitigate COVID-19 superspreading events during air travel. Specifically, we studied: 1) London to Hanoi with 201 passengers, including 13 secondary infections among passengers; 2) Singapore to Hangzhou with 321 passengers, including 12 to 14 secondary infections; 3) a non-superspreading event on a private jet in Japan with 9 passengers and no secondary infections. Here our results show that the inclusion of passenger movement better explains the infection spread patterns than conventional models do. We also found that FFP2/N95 mask usage would have reduced infection by 95–100%, while cloth masks would have reduced it by only 40–80%. Results indicate that leaving the middle seat vacant is effective in reducing infection, and the effectiveness increases when combined with good quality masks. However, with a good mask, the risk is quite low even without the middle seats being empty. Our results suggest the need for more stringent guidelines to reduce aviation-related superspreading events of COVID-19.

54 ENVIRONMENTAL SCIENCES↗

Implementation and Evaluation of Emission‐Driven Land‐Atmosphere Coupled Simulation in E3SMv2.1

Emissions-driven (prognostic CO 2 ) simulations are essential for representing two-way carbon-climate feedback in Earth System Models. We present an emissions-driven land–atmosphere coupled biogeochemistry (BGC) configuration (BGCLNDATM_progCO2) in version 2.1 of the Energy Exascale Earth System Model (E3SMv2.1). This is the first E3SM configuration that performs land-atmosphere emission-hindcasts. Here, we document its implementation, evaluate the model's performance against observations and other models, and propose a structured evaluation protocol for such emissions-driven simulations. We conducted transient historical simulations (1850–2014) with BGCLNDATM_progCO2 and compare them to reference simulations—a land-atmosphere coupled simulation without BGC and a standalone land simulation with BGC, both using prescribed CO 2 concentrations—and to observations. BGCLNDATM_progCO2 overestimates atmospheric CO 2 concentrations by 11–23 ppm yet stays within the 40-ppm spread CMIP6 emission-driven models and retains physical climate properties comparable to the reference runs. The CO 2 biases are partly attributed to underrepresented oceanic CO 2 uptake and inadequate representations of some terrestrial processes. In general, introducing prognostic CO 2 did not change physical climate metrics at the global scale but had larger regional effects, particularly over land where spatially heterogeneous CO 2 and prognostic leaf area index influenced surface energy balance. Finally, we propose a general evaluation protocol including spin-up assessment, atmospheric CO 2 benchmarking, physical climate evaluation, and land biogeochemical analysis to support scientific rigor and facilitate inter-model comparisons. The new configuration lays the groundwork for future enhancements, including improved terrestrial biogeochemical processes, integrated marine biogeochemistry, and additional human–Earth system interactions. These developments advance E3SM toward fully coupled emissions-driven simulations, enabling more accurate carbon–climate feedback projections and informing mitigation policy by providing physically consistent carbon-budget metrics for mitigation scenarios.

54 ENVIRONMENTAL SCIENCES↗

Transition Heat Transfer Modeling Based on the Characteristics of Turbulent Spots

While turbulence models are being developed which show promise for simulating the transition region on a turbine blade or vane, it is believed that the best approach with the greatest potential for practical use is the use of models which incorporate the physics of turbulent spots present in the transition region. This type of modeling results in the prediction of transition region intermittency which when incorporated in turbulence models give a good to excellent prediction of the transition region heat transfer. Some models are presented which show how turbulent spot characteristics and behavior can be employed to predict the effect of pressure gradient and Mach number on the transition region. The models predict the spot formation rate which is needed, in addition to the transition onset location, in the Narasimha concentrated breakdown intermittency equation. A simplified approach is taken for modeling turbulent spot growth and interaction in the transition region which utilizes the turbulent spot variables governing transition length and spot generation rate. The models are expressed in terms of spot spreading angle, dimensionless spot velocity, dimensionless spot area, disturbance frequency and Mach number. The models are used in conjunction with a computer code to predict the effects of pressure gradient and Mach number on the transition region and compared with VKI experimental turbine data.

Simon, Fred↗

Machine learning surrogate for charged particle beam dynamics with space charge based on a recurrent neural network with aleatoric uncertainty

In this work, we develop a machine learning (ML) model with aleatoric uncertainty for the low energy beam transport (LEBT) region of the LANSCE linear accelerator in which we model the transport of a space-charge-dominated 750 keV proton beam through a lattice of 22 quadrupole magnets. Our ML model is developed based on data generated by a Kapchinsky–Vladimirsky (KV) envelope model of beam transport. We show that a recurrent neural network can be used as a dynamical surrogate model for fast prediction of the LEBT beam envelope. Furthermore, we endow the model with the prediction of aleatoric uncertainty and compare three different approaches. We demonstrate that the ML-based uncertainty quantification models are well calibrated and produce good estimates of the regions where the model is less certain about its predictions. This ML framework is a necessary step in the development of a real-time virtual diagnostic tool with uncertainty quantification that can be integrated into more complex downstream tasks (e.g., adaptive control or learning flexible control policies via reinforcement learning) for improved efficiency in beam operations. In future work, we plan to expand on this preliminary study by considering more realistic envelope models that include longitudinal momentum spread and dispersive effects in bending magnets, as well as particle tracking codes with 3D space charge (such as and ). Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

Discovery of GeV emission from the Circinus galaxy with the Fermi large area telescope

We report the discovery of γ-ray emission from the Circinus galaxy using the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. Circinus is a nearby (∼4 Mpc) starburst with a heavily obscured Seyfert-type active nucleus, bipolar radio lobes perpendicular to the spiral disk, and kpc-scale jet-like structures. Our analysis of 0.1-100 GeV events collected during 4 yr of LAT observations reveals a significant (≅ 7.3σ) excess above the background. We find no indications of variability or spatial extension beyond the LAT point-spread function. A power-law model used to describe the 0.1-100 GeV γ-ray spectrum yields a flux of (18.8 ± 5.8) × 10{sup –9} photon cm{sup –2} s{sup –1} and photon index 2.19 ± 0.12, corresponding to an isotropic γ-ray luminosity of 3 × 10{sup 40} erg s{sup –1}. This observed γ-ray luminosity exceeds the luminosity expected from cosmic-ray interactions in the interstellar medium and inverse Compton radiation from the radio lobes. Thus, the origin of the GeV excess requires further investigation.

79 ASTRONOMY AND ASTROPHYSICS↗

Data-driven efficient network and surveillance-based immunization

Given a contact network and coarse-grained diagnostic information such as electronic Healthcare Reimbursement Claims (eHRC) data, can we develop efficient intervention policies from data to control an epidemic? Immunization is an important problem in multiple areas, especially epidemiology and public health. However, most existing studies rely on assuming prior epidemiological models to develop pre-emptive strategies, which may fail to adapt to the change in new epidemiological patterns and the availability of rich data such as eHRC. In practice, disease spread is usually complicated, hence assuming an underlying model may deviate from true spreading patterns, leading to possibly inaccurate interventions. Additionally, the abundance of health care surveillance data (such as eHRC) makes it possible to study data-driven strategies without too many restrictive assumptions. Therefore, such a data-driven intervention approach can help public-health experts take more practical decisions. In this work, we take into account propagation log and contact networks for controlling propagation. Furthermore, different from previous model-based approaches, our solutions are solely data driven in a sense that we develop immunization strategies directly from the network and eHRC without assuming classical epidemiological models. In particular, we formulate the novel and challenging data-driven immunization problem. To solve it, we first propose an efficient sampling approach to align surveillance data with contact networks, then develop an efficient algorithm with the provably approximate guarantee for immunization. Finally, we show the effectiveness and scalability of our methods via extensive experiments on multiple datasets, and conduct case studies on nation-wide real medical surveillance data.

60 APPLIED LIFE SCIENCES↗

Are Solar Activity Variations Amplified by the QBO: A Modeling Study

Solar cycle activity effects (SCAE) in the lower and middle atmosphere, reported in several studies, are difficult to explain on the basis of the small changes in solar radiation that accompany the 11-year cycle. It is therefore natural to speculate that dynamical processes may come into play to produce a leverage. Such a leverage may be provided by the Quasi-Biennial Oscillation (QBO) in the zonal circulation of the stratosphere, which has been linked to solar activity variations. Driven primarily by wave mean flow interaction, the QBO period and its amplitude are variable but are also strongly influenced by the seasonal cycle in the solar radiation. This influence extends to low altitudes and is referred to as 'downward control'. Small changes in the solar radiative forcing may produce small changes in the period and phase of the QBO, but these in turn may produce measurable differences in the wind field. Thus, the QBO may be an amplifier of solar activity variations and a natural conduit of these variations to lower altitudes. To test this hypothesis, we conducted experiments with a 2D version of our Numerical Spectral Model that incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GW). Solar cycle radiance variations (SCRV) are accounted for by changing the radiative heating rate on a logarithmic scale from 0.1% at the surface to 1% at 50 km to 10% at 100 km. With and without SCRV, but with the same GW flux, we then conduct numerical experiments to evaluate the magnitude of the SCAE in the zonal circulation. The numerical results indicate that, under certain conditions, the SCAE is significant and can extend to lower altitudes where the SCRV is small. For a modeled QBO period of 30 months, we find that the seasonal cycle in the solar forcing acts as a strong pacemaker to lock up the phase and period of the QBO. The SCAE then shows up primarily as a distinct but relatively weak amplitude modulation. But with a different QBO period between 30 and 34 (or less than 30, presumably) months, the seasonal phase lock is weak. Solar flux variations in the seasonal cycle then cause variations in the QBO period and phase. These amplify the SCAE to produce relatively large variations in the wind field. The SCAE in this case extends to mid-latitudes.

Mengel, J. G.↗

Solar Cycle Variations and Equatorial Oscillations: Modeling Study

Solar cycle activity effects (SCAE) in the lower and middle atmosphere, reported in several studies, are difficult to explain on the basis of the small changes in solar radiation that accompany the 11-year cycle, It is therefore natural to speculate that dynamical processes may come into play to produce a leverage. Such a leverage may be provided by the Quasi-Biennial Oscillation (QBO) in the zonal circulation of the stratosphere, which has been linked to solar activity variations. Driven primarily by wave mean flow interaction, the QBO period and its amplitude are variable but are also strongly influenced by the seasonal cycle in the solar radiation. This influence extends to low altitudes referred to as "downward control". Relatively small changes in solar radiative forcing can produce small changes in the period and phase of the QBO, but this in turn can produce measurable differences in the wind field. Thus, the QBO may be an amplifier of solar activity variations and a natural conduit of these variations to lower altitudes. To test this hypothesis, we conducted experiments with a 2D (two-dimensional) version of our Numerical Spectral Model that incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GW). Solar cycle radiance variations (SCRV) are accounted for by changing the radiative heating rate on a logarithmic scale from 0.1 % at the surface to 1 % at 50 km to 10% at 100 km. With and without SCRV, but with the same GW flux, we then conduct numerical experiments to evaluate the magnitude of the SCAE in the zonal circulation. The numerical results indicate that, under certain conditions, the SCAE is significant and can extend to lower altitudes where the SCRV is inconsequential. At 20-km the differences in the modeled wind velocities are as large as 5 m/s. For a modeled QBO period of 30 months, we find that the seasonal cycle in the solar forcing (through the Semi-annual Oscillation (SAO)) acts as a strong pacemaker to lockup the phase and period of the QBO. The SCAE then shows up primarily as a distinct but relatively weak amplitude modulation. But with the QBO period between 30 and 34 (or less than 30, presumably) months, the seasonal phase lock is weak. Solar flux radiance variations in the seasonal cycle then cause variations in the QBO period and phase that amplify the SCAE to produce relatively large variations in the wind field. These variations also extend to mid latitudes.

Mayr, H. G.↗

Modelling of equatorial phase and amplitude scintillations from OGO-6 and AE irregularity data

The technique of modelling equatorial phase and amplitude scintillations from in-situ measurements of irregularity amplitude, ambient ionization density and irregularity power spectrum is outlined. The in-situ measurements do not directly provide information on irregularity layer thickness for which radar backscatter observations are utilized. The equatorial scintillation models developed from OGO-6 and AE-C in-situ measurements during the two solstices are presented. Considerable differences in the longitude variation of scintillations are predicted during the two solstices. The model estimates are in good agreement with the available ground based phase and amplitude measurements. Problems associated with the use of bottomside spread-F data for transionospheric propagation modelling at VHF/UHF are also discussed.

Basu, S.↗

Comprehensive compartmental model and calibration algorithm for the study of clinical implications of the population-level spread of COVID-19: a study protocol

The complex dynamics of the coronavirus disease 2019 (COVID-19) pandemic has made obtaining reliable long-term forecasts of the disease progression difficult. Simple mechanistic models with deterministic parameters are useful for short-term predictions but have ultimately been unsuccessful in extrapolating the trajectory of the pandemic because of unmodelled dynamics and the unrealistic level of certainty that is assumed in the predictions. We propose a 22-compartment epidemiological model that includes compartments not previously considered concurrently, to account for the effects of vaccination, asymptomatic individuals, inadequate access to hospital care, post-acute COVID-19 and recovery with long-term health complications. Additionally, new connections between compartments introduce new dynamics to the system and provide a framework to study the sensitivity of model outputs to several concurrent effects, including temporary immunity, vaccination rate and vaccine effectiveness. Subject to data availability for a given region, we discuss a means by which population demographics (age, comorbidity, socioeconomic status, sex and geographical location) and clinically relevant information (different variants, different vaccines) can be incorporated within the 22-compartment framework. Considering a probabilistic interpretation of the parameters allows the model's predictions to reflect the current state of uncertainty about the model parameters and model states. We propose the use of a sparse Bayesian learning algorithm for parameter calibration and model selection. This methodology considers a combination of prescribed parameter prior distributions for parameters that are known to be essential to the modelled dynamics and automatic relevance determination priors for parameters whose relevance is questionable. This is useful as it helps prevent overfitting the available epidemiological data when calibrating the parameters of the proposed model. Population-level administrative health data will serve as partial observations of the model states.

59 BASIC BIOLOGICAL SCIENCES↗

Measures of model performance based on the log accuracy ratio

Quantitative assessment of modeling and forecasting of continuous quantities uses a variety of approaches. We review existing literature describing metrics for forecast accuracy and bias, concentrating on those based on relative errors and percentage errors. Of these accuracy metrics, the mean absolute percentage error (MAPE) is one of the most common across many fields and has been widely applied in recent space science literature and we highlight the benefits and drawbacks of MAPE and proposed alternatives. We then introduce the log accuracy ratio, and derive from it two metrics: the median symmetric accuracy; and the symmetric signed percentage bias. Robust methods for estimating the spread of a multiplicative linear model using the log accuracy ratio are also presented. The developed metrics are shown to be easy to interpret, robust, and to mitigate the key drawbacks of their more widely-used counterparts based on relative errors and percentage errors. Their use is illustrated with radiation belt electron flux modeling examples.

54 ENVIRONMENTAL SCIENCES↗

Fluid Dynamics Effects on Microstructure Prediction in Single Laser Tracks for Additive Manufacturing

The Laser Powder Bed Fusion Additive Manufacturing (LPBFAM) is one of the most important processes for the production of lightweight, cost-effective, complex, and high-performance enduse parts. At present, the cost and time associated with LPBFAM process development is very high due to a lack of fundamental process understanding. In this project, a multi-physics model was developed on a highly parallel open-source code, Truchas, with the ultimate goal of providing experimentally validated process maps for tailoring microstructure to achieve desired performance for LPBFAM. As a critical step towards fully LPBFAM modeling, modeling of single-track laser fusion (STLF) were conducted. Multi-physics simulations were conducted using Truchas for STLF by considering heat transfer, phase-change, fluid dynamics, surface tension phenomena, and evaporation. In order to assess the effect of the fluid dynamics on the solidification and ensuing microstructure, a heat transfer-and-solidification-only (HTS) model and a fully coupled heat-transfer-solidification and fluid-dynamics (HTS-FD) model were considered. In the HTS_FD model, the fluid flow was considered to be laminar while the molten alloy surface is assumed to be flat and non-deformable. This study is one the first attempt to understand the effect of the fluid flow on microstructure in STLF and LPBFAM. The results show that the fluid flow affects the solidification and ensuing microstructure in STLF. In order to validate Truchas for STLF and LPBFAM modeling, experimental data, which was obtained at GE Global Research (GEGR) for liquid pool shape and microstructure, were compared with those from numerical simulation results. For the keyhole regime case, numerical simulation results indicate that the melt-pool shape was typical to that of the conduction case with very large deviations from the measured melt-pool depths. Two sensitivity studies were conducted by varying the evaporation flux and the surface tension coefficient at fixed laser absorptivities of 0.5 and 0.3, respectively. The minimum value for the overall combined error, between the calculated values with the HTS_FD model and measured values for both the width and height of the melt-pool, was attained for a surface tension coefficient of -0.75e-4 N/m. An analytical solidification model for the columnar-to-equiaxed transition (CET) in rapid solidification was used to assess the microstructure variation within the melt-pool. From numerical simulation results, thermal gradient, G, and solidification velocity, V, were obtained in order to predict the microstructure type (e.g., dendritic, cellular). For the fluid-flow HTS_FD model, G was found to exhibit a maximum at the extremity of the solidified pool (i.e., at the free surface). By contrast, for HTS simulations, G was found to exhibit a maximum around the entire edge of the solidified pool. For the HTS_FD simulations, the minimum values of the cooling rate, GV, were found to be approximately half than their corresponding values for HTS simulations. By contrast with the min GV values, the maximum values for GV were found to occur for the HTS_FD simulations. Thus, HTS_FD simulations were found to exhibit a wider range of cooling rates than the HTS simulations. Concerning the solidification maps, the variation of the solidification velocity, V, as function of the thermal gradient, G, was obtained. It was found that the fluid flow model results (HTS_FD) exhibited an increased spread in the V(G) variation within the melt pool with respect to the HTS model results (without the fluid flow). A preliminary correlation model for the primary dendrite arm spacing (PDAS) based on a power law dependence on thermal gradient and solidification velocity was used to estimate PDAS from HTS and HTS_FD simulations. For low laser powers and low laser speeds, the PDAS obtained 6 with the fluid dynamics model (HTS_FD) was larger by more than 30% with respect to the PDAS calculated with the simple HTS model. In the second part of the project, fifty-seven simulations were conducted in order to obtain data on microstructure variables that can be used for process map development. In order to cover the entire processing space, thirty-three cases were selected STLF process simulations at six power levels and six laser scanning speeds. The 57 STLF process simulations were conducted as follows: 33 simulations with the heat-transfer-only (HTS) model and 24 simulations with the fluid flow model (HTS_FD). The solidification map data shows that for all simulations, a columnar dendritic microstructure would be expected. It was found that the minimum, average, and maximum thermal gradient exhibit exponential variations with respect to a process variable defined as the ratio between the power and square root of scan speed.

36 MATERIALS SCIENCE↗

Diurnal Forcing of Planetary Atmospheres

Much progress has been made on calculations of the Martian seasonal water cycle using the Mars Climate Model developed for this purpose. Two papers, documenting the model and the water transport results obtained with it have been published in the Journal of Geophysical Research - Planets. An additional paper describing results related to the evolution of the seasonal water cycle as a result of orbital changes was published in Advances in Space Research. Since that time, further studies have concentrated on the consequences of the soil adsorption required to match the observed water cycle and its relation to the stability of ground ice and other potential water reservoirs. Earth-related studies have concentrated on incorporating an efficient and realistic microphysical model into the Ames Stratospheric General Circulation Model used to simulate the spread of the ML Pinatubo and other volcanic clouds in the stratosphere. In addition, visualizations of the simulations are being incorporated into a video describing the UARS mission. A paper describing the new stratospheric aerosol microphysics package (and its consequences for volcanic cloud evolution) will be submitted in the near future. The paper will discuss the relative importance of condensation and coagulation to early particle growth and the separation of the cloud by sedimentation of the larger particles. A more general paper which highlights the observation that particle number densities did not increase dramatically after the ML Pinatubo eruption is planned. Simulations of atmospheric transport will be extended to include studies of terrestrial tropospheric tracers using the Fifth-Generation Penn State/NCAR Mesoscale Model.

Houben, Howard C.↗