Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “neutral model”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

An investigation of the solar zenith angle variation of D-region ionization

Model calculations are carried out with a view to interpreting the solar zenith angle variation of D-region ionization. A model is developed for the neutral chemistry including the transport terms relating to molecular and eddy diffusion. The diurnal behavior is described of the minor neutral constituents formed in an oxygen-hydrogen-nitrogen atmosphere, in the height interval between 30 and 120 km. Computations carried out for two cases of the eddy diffusion coefficients models indicate that the constituents which are important for the D-region positive-ion chemistry do not show a significant variation with zenith angle for values up to 75 deg over the D-region heights. In the ion chemistry model, ion-pair production rates are calculated for solar X-rays between 1 A and 100 A, EUV radiations from 100 A up to the Lyman-alpha line, precipitating electrons, and galactic cosmic rays. The solar zenith angle variation of the positive-ion composition, negative-ion composition, and the electron densities are described up to 75 deg zenith angle, in the height interval between 60 and 100 km.

Ratnasiri, P. A. J.↗
NASA NTRS2023

Numerical Investigation of Occupant Injury Risks During A Realistic Transport Aircraft Crash Conditions

Researchers at the National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC) have conducted a full-scale crash test of a Fokker F28 MK1000 aircraft to investigate the performance of transport aircraft under realistic crash conditions. This crash test was computationally recreated using finite element (FE) human body models (HBMs) to further explore potential injury risks to occupants and analyze the utilization of HBMs in the aerospace crash environment. The Global Human Body Model Consortium (GHBMC) male 50th percentile occupant detailed model (v6.0) and the Toyota Human Model for Safety (THUMS) male 50th percentile occupant model (v6.1) were selected to be used in the crash simulations. The HBMs were simulated in conditions matching those of anthropomorphic test device (ATD) experiments included within the aircraft cabin during the crash test. Seven occupant locations within the cabin were simulated utilizing each of the models. The models were positioned in a neutral upright posture with hands resting on the legs and the feet contacting the floor. Head, brain, neck, and lumbar vertebra injury metrics were calculated for all trials. Both HBM models required minor modifications to stabilize these simulations. The GHBMC model required added erosion for six parts while the THUMS model only required one in order to complete the full simulation. The THUMS model, however, required a much smaller timestep for stability and therefore took significantly more computational time. In addition the GHBMC model includes integrated instrumentation while the THUMS model requires development and implementation of instrumentation. Both models predicted 100% injury risk for lumbar vertebra fracture in all test conditions. This prediction was in family with high lumbar load values measured by the ATDs during the crash test. The THUMS model consistently predicted lower injury risks than the GHBMC model in all three other metrics varying depending on the crash pulse. Overall, the THUMS model required less modifications to allow for this study. However, the GHBMC models significantly faster run time and integrated instrumentation make it a more intuitive model for this research.

Crashworthiness↗
DOE OSTI2025

Guiding Principles for Geochemical/Thermodynamic Model Development and Validation in Nuclear Waste Disposal: A Close Examination of Recent Thermodynamic Models for H + —Nd 3+ —NO 3 - (—Oxalate) Systems

Development of a defensible source-term model (STM), usually a thermodynamical model for radionuclide solubility calculations, is critical to a performance assessment (PA) of a geologic repository for nuclear waste disposal. Such a model is generally subjected to rigorous regulatory scrutiny. In this article, we highlight key guiding principles for STM model development and validation in nuclear waste management. We illustrate these principles by closely examining three recently developed thermodynamic models with the Pitzer formulism for aqueous H + —Nd 3+ —NO 3 - (—oxalate) systems in a reverse alphabetical order of the authors: the XW model developed by Xiong and Wang, the OWC model developed by Oakes et al., and the GLC model developed by Guignot et al., among which the XW model deals with trace activity coefficients for Nd(III), while the OWC and GLC models are for concentrated Nd(NO 3 ) 3 electrolyte solutions. The principles highlighted include the following: (1) Principle 1. Validation against independent experimental data: A model should be validated against experimental data or field observations that have not been used in the original model parameterization. We tested the XW model against multiple independent experimental data sets including electromotive force (EMF), solubility, water vapor, and water activity measurements. The results show that the XW model is accurate and valid for its intended use for predicting trace activity coefficients and therefore Nd solubility in repository environments. (2) Principle 2. Testing for relevant and sensitive variables: Solution pH is such a variable for an STM and easily acquirable. All three models are checked for their ability to predict pH conditions in Nd(NO 3 ) 3 electrolyte solutions. The OWC model fails to provide a reasonable estimate for solution pH conditions, thus casting serious doubt on its validity for a source-term calculation. In contrast, both the XW and GLC models predict close-to-neutral pH values, in agreement with experimental measurements. (3) Principle 3. Honoring physical constraints: Upon close examination, it is found that the Nd(III)-NO 3 association schema in the OWC model suffers from two shortcomings. Firstly, its second stepwise stability constant for Nd(NO 3 ) 2+ (log K 2 ) is much higher than the first stepwise stability constant for NdNO 3 2+ (log K 1 ), thus violating the general rule of (log K 2 –log K 1 ) < 0, or $\frac{K1}{K2}$>1. Secondly, the OWC model predicts abnormally high activity coefficients for Nd(NO 3 ) 2 + (up to ~900) as the concentration increases. (4) Principle 4. Minimizing degrees of freedom for model fitting: The OWC model with nine fitted parameters is compared with the GLC model with five fitted parameters, as both models apply to the concentrated region for Nd(NO 3 ) 3 electrolyte solutions. The latter appears superior to the former because the latter can fit osmotic coefficient data equally well with fewer model parameters. The work presented here thus illustrates the salient points of geochemical model development, selection, and validation in nuclear waste management.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

An mRNA vaccine encoding the Ebola virus glycoprotein induces high neutralizing antibody titers and provides strong protection against lethal infections in mouse models

Ebola virus (EBOV) is the causative agent of Ebola disease (EBOD), a viral hemorrhagic fever with a notably high case fatality rate. Current treatments for EBOD are limited to monoclonal antibodies or two licensed viral vector vaccines, a recombinant vesicular stomatitis virus (rVSV)-vectored vaccine or an adenovirus and modified vaccinia Ankara regimen. However, comparisons of protection, efficacy, and durability with alternative nucleotide platforms remain understudied. Here, we evaluated the immunogenicity of an mRNA vaccine expressing the EBOV glycoprotein (GP) in parallel with rVSV- and DNA-based vaccine platforms. The mRNA EBOV-GP vaccine, formulated in lipid nanoparticles, elicited significantly higher levels of total IgG and neutralizing antibody titers compared to the rVSV-EBOV-GP vaccine. Linear antibody epitope analysis indicated a preference for targeting the mucin-like domain in EBOV-GP1 following rVSV-based vaccination, while the mRNA platform distinctly targeted the internal fusion loop of EBOV-GP2. After characterizing the immunogenicity of the mRNA vaccine, two models of EBOD were used to demonstrate its protective efficacy: a surrogate rVSV-based challenge model of EBOD using type-I interferon deficient C57BL/6 mice and infection of BALB/c mice with authentic mouse-adapted EBOV. In both studies, the EBOV mRNA vaccine fully protected the mouse cohorts against morbidity and mortality. Additionally, the EBOV mRNA vaccine produced greater neutralizing antibody titers compared to the DNA EBOV-GP vaccine. These results suggest that an mRNA vaccine expressing EBOV-GP can induce robust, functional humoral responses that are protective against EBOD, warranting further development as an alternative to, or as part of a vaccine strategy including, viral vectored vaccines.

DNA vaccines↗
DOE OSTI2026

Role of Neutrals Versus Transport in Determining the Pedestal Density Structure: Final Technical Report

In fusion devices the plasma density plays a crucial role in determining the fusion reaction rate and has a direct impact on the fusion gain of a given device. This density is in general regulated by the particle sources and transport near the plasma edge, which give rise to an edge density pedestal. When predicting the performance of future devices, this density pedestal is often prescribed, rather than predicted, due to a lack of models which allow confident extrapolation. This project aims to advance these models through the focused validation of theoretical models related to the transport of fueling neutral particles, and through interpretive transport modeling in present day fusion plasmas, in which the penetration of neutrals is altered to better simulate future reactor-like conditions. Achievements in theory and model validation under this project have advanced our understanding of how much of the edge density profile is set by transport versus direct ionization, enabling interesting projections to future burning plasma devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Signals with six bottom quarks for charged and neutral Higgs bosons

In extensions of two Higgs doublet models with vectorlike quarks, the decays of vectorlike quarks may easily be dominated by cascade decays through charged or neutral Higgs bosons leading to signatures with 6 top or bottom quarks. Since top quark decays also contain bottom quarks, a common signature for many possible decay chains is 6 bottom quarks in the final state. We present a search strategy focusing on this final state and find the mass ranges of vectorlike quarks and Higgs bosons that can be explored at the Large Hadron Collider. Among other results, the sensitivity to the charged and neutral Higgs bosons, extending to about 2 TeV, stands out when compared to models without vectorlike matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Isolation and characterization of cross-neutralizing coronavirus antibodies from COVID-19+ subjects

SARS-CoV-2 is one of three coronaviruses that have crossed the animal-to-human barrier and caused widespread disease in the past two decades. The development of a universal human coronavirus vaccine could prevent future pandemics. We characterize 198 antibodies isolated from four COVID-19+ subjects and identify 14 SARS-CoV-2 neutralizing antibodies. One targets the N-terminal domain (NTD), one recognizes an epitope in S2, and 11 bind the receptor-binding domain (RBD). Three anti-RBD neutralizing antibodies cross-neutralize SARS-CoV-1 by effectively blocking binding of both the SARS-CoV-1 and SARS-CoV-2 RBDs to the ACE2 receptor. Using the K18-hACE transgenic mouse model, we demonstrate that the neutralization potency and antibody epitope specificity regulates the in vivo protective potential of anti-SARS-CoV-2 antibodies. All four cross-neutralizing antibodies neutralize the B.1.351 mutant strain. Thus, our study reveals that epitopes in S2 can serve as blueprints for the design of immunogens capable of eliciting cross-neutralizing coronavirus antibodies.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1993

A quantitative model of plasma in Neptune's magnetosphere

A model encompassing plasma transport and energy processes is applied to Neptune's magnetosphere. Starting with profiles of the neutral densities and the electron temperature, the model calculates the plasma density and ion temperature profiles. Good agreement between model results and observations is obtained for a neutral source of 5 x 10 exp 25/s if the diffusion coefficient is 10 exp -8 L3R(N)/2s, plasma is lost at a rate 1/3 that of the strong diffusion rate, and plasma subcorotates in the region outside Triton.

Richardson, John D.↗
DOE OSTI2021

Quantifying Polaron Mole Fractions and Interpreting Spectral Changes in Molecularly Doped Conjugated Polymers

Molecular doping of conjugated polymers causes bleaching of the neutral absorbance and results in new polaron absorbance transitions in the mid and near infrared. In this report the concentration dependent changes in the spectra for a series of molecularly doped diketopyrrolopyrrole (DPP) co-polymers with a series of ultra-high electron affinity cyanotrimethylenecyclopropane-based dopants is analyzed. With these strong dopants the polaron mole fraction (Θ) reaches saturation. Analysis of the full spectrum enables separation of neutral and polaron signals and quantification of the polaron mole fraction using a simple noninteracting site model. The peak ratios for both neutral and polaron peaks change systematically with increasing polaron mole fraction for all measured polymers. Analysis of the spectral changes indicates that the polaron mole fraction can be quantified to within 5%. While the total change in the absorbance spectrum with increasing polaron mole fraction is linear, the lowest energy polaron peak (P1) grows nonlinearly, which indicates increased polarization/delocalization. Molecular doping of polymers that form either H- or J-aggregates shows systematically different spectral changes in the vibronic peak ratios of the neutral spectra and provides insights into the polymer configuration at undoped sites in the film.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Cascade decays of heavy Higgs bosons through vectorlike quarks in two Higgs doublet models

We study cascade decays of heavy neutral Higgs bosons through vectorlike quarks. We focus on scenarios where decay modes into pairs of vectorlike quarks are not kinematically open which extends the sensitivity of the LHC to larger masses. Assuming only mixing with the third family of standard model quarks the new decay modes of heavy Higgs bosons are: H → t 4 t → Wbt, Ztt, htt and H → b 4 b → Wtb, Zbb, hbb, where t 4 (b 4 ) is the new up-type (down-type) quark mass eigenstate. In the numerical analysis we assume the CP even Higgs boson in the two Higgs doublet model type-II but the signatures are relevant for many other scenarios. We identify the region of the parameter space where these decay modes are significant or can even dominate, and thus they provide the best opportunities for the simultaneous discovery of a new Higgs boson and vectorlike quarks. We further explore the reach of the High Luminosity LHC for two representative decay modes,t 4 → Zt → ℓℓt and b 4 → Zb → ℓℓb, and found that cross sections at a 0.1 fb level can be probed with simple cut based analyses. We also find that the rates for Higgs cascade decays can be much larger than the rates for a single production of vectorlike quarks. Furthermore, the reach for vectorlike quarks in Higgs cascade decays and pair production extends to comparable masses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2000

Using Model-Based Reasoning for Autonomous Instrument Operation

Multiprobe missions are an important part of NASA's future: Cluster, Magnetospheric Multi Scale, Global Electrodynamics and Magnetospheric Constellation are representatives from the Sun-Earth Connections Theme. To make such missions robust, reliable, and affordable, ideally the many spacecraft of a constellation must be at least as easy to operate as one spacecraft is today. To support this need for scalability, science instrumentation must become increasingly easy to operate, even as this same instrumentation becomes more capable and advanced. Communication and control resources will be at a premium for future instruments. Many missions will be out of contact with ground operators for extended periods either to reduce operations cost or because of orbits that limit communication to weekly perigee transits. Autonomous capability is necessary if such missions are to effectively achieve their operational objectives. An autonomous system is one that acts given its situation in a mission appropriate manner without external direction to achieve mission goals. To achieve this capability autonomy must be built into the system through judicious design or through a built-in intelligence that recognizes system state and manages system response. To recognize desired or undesired system states, the system must have an implicit or explicit understanding of its expected states given its history and self observations. The systems we are concerned with, science instruments, can have stringent requirements for system state knowledge in addition to requirements driven by health and safety concerns. Without accurate knowledge of the system state, the usefulness of the science instrument may be severely limited. At the same time, health and safety concerns often lead to overly conservative instrument operations further reducing the effectiveness of the instrument. These requirements, coupled with overall mission requirements including lack of communication opportunities and tolerance of environmental hazards, frame the problem of constructing autonomous science instruments. we are developing a model of the Low Energy Neutral Atom instrument (LENA) that is currently flying on board the Imager for Magnetosphere-to-Aurora Global Exploration (IMAGE) spacecraft. LENA is a particle detector that uses high voltage electrostatic optics and time-of-flight mass spectrometry to image neutral atom emissions from the denser regions of the Earth's magnetosphere. As with most spacecraft borne science instruments, phenomena in addition to neutral atoms are detected by LENA. Solar radiation and energetic particles from Earth's radiation belts are of particular concern because they may help generate currents that may compromise LENA's long term performance. An explicit model of the instrument response has been constructed and is currently in use on board IMAGE to dynamically adapt LENA to the presence or absence of energetic background radiations. The components of LENA are common in space science instrumentation, and lessons learned by modelling this system may be applied to other instruments. This work demonstrates that a model-based approach can be used to enhance science instrument effectiveness. Our future work involves the extension of these methods to cover more aspects of LENA operation and the generalization to other space science instrumentation.

Johnson, Mike↗
NASA NTRS2003

Modeling Variability of Plasma Conditions in the Io Torus

Telescopic observations an in situ measurements of the Io plasma torus show the density, temperature and composition to vary over time, sometimes up to a factor of 2. While previous models of the physical and chemical processes in the Io plasma torus have reasonably modeled the conditions of the Voyager 1 era, their authors have not addressed the observed variability nor explored the sensitivity of torus conditions to input parameters. In this paper we present a homogeneous torus model parameterized by five variables (transport timescale, neutral source strength, ratio of oxygen sulfur to atoms in the source, fraction of superthermal electrons, temperature of these hot electrons). The model incorporates the most recent data for ionization, recombination, charge exchange and radiative energy losses for the major torus species (S, S(sup +), S(sup ++), S(sup +++), O, O(sup +), O(sup ++). We solve equations of conservation of mass and energy to find equilibrium conditions for a set of input parameters. We compare model plasma conditions with those observed by Voyager 1 Voyager 2, and Cassini. Furthermore, we explore the sensitivity of torus conditions to each parameter. We find that (1) torus conditions are distinctly different for the Voyager 1, Voyager 2 and Cassini eras, (2) unique torus input parameters for any given era are poorly constrained given the wide range of solution space that is consistent with the range of observed torus conditions, (3) ion composition is highly sensitive to the specification of a non-thermal electron distribution, (4) neutral O/S source ratio is highly variable with model values ranging between 1.7 for Cassini to 4.0 for Voyager 1 conditions, (5) transport times range between 23 days for Voyager 2 to 50 days for Voyager 1 and Cassini, (6) neutral source strengths range between 7 to 30 x 10(sup -4) cm (sup -3) s(sup -1) which corresponds to a net production of 0.4 to 1.3 tons/s for a torus volume of 1.4 x 10(sup 31) cm(sup 3), or 38 R(sub j)(sup 3).

Delamere, P. A.↗
DOE OSTI2022

A reduced model for the ITER divertor based on SOLPS solutions for ITER Q = 10 baseline conditions: B. A reduced model based on reversed-direction two point modeling

Edge codes such as SOLPS coupled to neutral codes such as EIRENE have become so comprehensive and sophisticated that they now constitute, in effect, 'code-experiments' that, as for actual experiments, can benefit from interpretation using simple models and conceptual frameworks, i.e. reduced models. The first task is the identification of options for the reduced model control parameters that are best suited for control of the action of the divertor, i.e. for control of target power loading and sputter–erosion, primarily. A strong correlation between the electron temperature at the divertor target, T e,t , and the neutral deuterium D 2 density at the target, n D2,t , flux-tube resolved, has recently been reported for a number of code studies including SOLPS-4.3 modeling of a set of ~50 ITER baseline cases: Q DT = 10, q 95 = 3, P SOL = 100 MW, metallic walls, and Ne seeding (Pitts et al 2019). Part A of the present study reports new results for largely the same ITER cases, confirming the strong correlation reported earlier between local values of T e,t , and (i) n D2,t , and (ii) normalized volumetric losses of power and pressure in the divertor. Strong correlations have now also been found, and are reported here for the first time, between T e,t and all of the divertor target quantities of practical interest. A physical explanation for this surprising result has not as yet been fully identified; nevertheless it has encouraging implications for reduced modeling of the ITER divertor. For such ITER conditions, (i) the global Ne injection rate, Inj Ne (Ne s –1 ), and (ii) the electron temperature at the location on the target where the peak power deposition occurs, T e,t @q ⊥,pk (eV), are found to be promising reduced model control parameters. In this part B, a reduced model for the ITER divertor is developed and described in detail, based on reversed-direction two point modeling, Rev2PM. The input to the reduced model is a value of the variable pair (T e,t @q ⊥,pk , Inj Ne )for a chosen case and the output are values of the various target as well as divertor-entrance quantities of practical interest, e.g. q ⊥,pk , the electron density at the X-point, n e,Xpt , etc. The reduced model was quantitatively characterized using one half of the code cases; it was then used to successfully predict (replicate) the code values of e.g. n e,Xpt for the other half.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

The ionization of the diffuse ionized gas

The diffuse ionized gas (DIG), responsible for producing the faint but pervasive H(sub alpha), requires 15% of the power of all Galactic O stars for its ionization. No other source of ionization seems practical. The spectrum is strong in low stages of ionization (N II, SII), weak in (O III), and very weak in (O I) lambda 6300, at least in one well-observed direction, all significantly different from denser H II regions. We give low-density, low-excitation photoionization models that explain the observed spectrum. The lambda 6300 observation eliminates the simplest models for the DIG (Mathis 1986), in which neutral H extends beyond the edge of the very dilute stellar radiation field. Our present models include two components. One, representing the edges of interstellar H I clouds, extends to the point where H becomes neutral. In the second, the fraction of H(sup o) is not allowed to exceed 0.05 to 0.1. Both have very low values of the ionization parameter, or ratio of the number densities of ionizing photons to electrons. The ionization parameter required by our models is shown to be compatible with observed values of photon and electron densities in the diffuse interstellar medium. Interstellar dust is not important in either model. Predictions of the model are that the Galaxy is leaking about 4% of the ionizing radiation of at least the local O stars into the intergalactic medium, that (O II) lambda 3727 approximately equals 1.1 H(sub alpha), and that He(sup +)/H (sup +)) approximately equals 0.6 He/H. A major difficul ty of this picture is having ionizing radiation propagate from the ionizing stars in the plane of the Galaxy to the DIG. We suggest, as do Miller & Cox (1993), that the ionizing radiation from O stars produces holes or tunnels in the distribution of neutral H as seen from their positions, while we see a sky covered with H I sheets or filaments because we are not near a source of ionizing radiation. With our models, several interfaces between ionized and neutral gas are allowed before the observational limit on lambda 6300 is exceeded. Another difficulty is the gas-phase N(S(sup +))/N(H(sup o)) ratio toward the high-latitude star HD 93521 (Spitzer & Fitzpatrick 1993). We require (S/H) to be somewhat lower than solar, either because of depletion onto grains or a Galactic abundance gradient, so that the S(sup +) can be produced in each of separate H(sup o) and H(sup +) regions.

Domgorgen, H.↗
DOE OSTI2024

Verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations

We present the development and verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations. We apply the corrected binary collision model originally developed for charged-particles collisions to all considered scattering channels, including Coulomb collisions, elastic neutral–neutral and charged–neutral collisions, ionization, excitation, and fusion. The model's implementation is described and verified through a series of simulations, including charged- and neutral-particle thermal equilibration, slowing of electrons in warm solid-density aluminum, collisional damping of a Langmuir wave, helium gas breakdown in an applied electric field, and thermonuclear and beam–target fusion. Then, we demonstrate the model within simulations of hydrogen plasma formation in the Princeton Field-Reversed Configuration as well as of the burning of aneutronic fusion fuel p-11B. The latter includes measurement of the fusion power density in a low-density plasma and fusion production due to the stopping of a proton ignitor beam in a compressed boron target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Profile Evolution Simulation in Etching Systems Using Level Set Methods

Semiconductor device profiles are determined by the characteristics of both etching and deposition processes. In particular, a highly anisotropic etch is required to achieve vertical sidewalls. However, etching is comprised of both anisotropic and isotropic components, due to ion and neutral fluxes, respectively. In Ar/Cl2 plasmas, for example, neutral chlorine reacts with the Si surfaces to form silicon chlorides. These compounds are then removed by the impinging ion fluxes. Hence the directionality of the ions (and thus the ion angular distribution function, or IAD), as well as the relative fluxes of neutrals and ions determines the amount of undercutting. One method of modeling device profile evolution is to simulate the moving solid-gas interface between the semiconductor and the plasma as a string of nodes. The velocity of each node is calculated and then the nodes are advanced accordingly. Although this technique appears to be relatively straightforward, extensive looping schemes are required at the profile corners. An alternate method is to use level set theory, which involves embedding the location of the interface in a field variable. The normal speed is calculated at each mesh point, and the field variable is updated. The profile comers are more accurately modeled as the need for looping algorithms is eliminated. The model we have developed is a 2-D Level Set Profile Evolution Simulation (LSPES). The LSPES calculates etch rates of a substrate in low pressure plasmas due to the incident ion and neutral fluxes. For a Si substrate in an Ar/C12 gas mixture, for example, the predictions of the LSPES are identical to those from a string evolution model for high neutral fluxes and two different ion angular distributions.(2) In the figure shown, the relative neutral to ion flux in the bulk plasma is 100 to 1. For a moderately isotropic ion angular distribution function as shown in the cases in the left hand column, both the LSPES (top row) and rude's string method (bottom row) predict tapered profiles. The LSPES uses an AD with a FWHM = 13.5 degrees, and rude's model uses a ratio of sheath voltage to ion temperature of 50. The more anisotropic IADs produce profiles with more vertical sidewalls and a wider bottom surface, as shown in the right hand column. Here, the LSPES has an AD with a FWHM = 2 degrees, and rude's model uses a ratio of 500 of sheath voltage to ion temperature. The agreement between the LSPES and rude's model is excellent in both cases. We will present etching profiles generated by the LSPES, including calculations of the re-emitted fluxes of both neutrals and ions off of the profile walls. In addition, we will show the effect of geometric structures (overhangs, etc.) on the etching profiles. Other physical aspects, such as surface diffusion, will also be included in the model.

Hwang, Helen H.↗
NASA NTRS1991

The Marshall Engineering Thermosphere model atmosphere Statistical Analysis Mode (MET-SAM)

The minimum, mean, and maximum exospheric temperature on the globe were calculated for every three hour period from 1947 through 1989 using the algorithms in the Marshall Engineering Thermosphere (MET) model and the appropriate solar activity input parameters. Cumulative percent frequency (CPF) distributions were then calculated for each of these temperatures at five levels of solar activity as defined by the 13-month smoothed values of the 10.7-cm solar radio noise flux. Next, the 50, 95, 97.7, and 100 percentile temperature values in each of these five levels of solar activity were curve fit as a function of the 13-month smoothed 10.7-cm flux. The resulting algorithms are used to compute the exospheric temperature in the MET model instead of the technique developed by Jacchia in his 1970 model. These temperatures are then used to enter tables to determine the total mass density and/or the atomic oxygen number density for application to engineering problems. Users can specify the risk level they are willing to accept in the results of analyses that require neutral atmosphere parameters inputs. The model eliminates the guess work in how to combine the solar activity input parameters to insure that the results provide answers at the proper risk levels.

Smith, Robert E.↗