SEARCH · Engineering Papers

Results for “Simulation Cases”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,495 records · Page 83

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Experimental Characterization of Non-Associative Plasticity Flow Rule Coefficients and Post-Peak Stress Degradation for the LS-DYNA MAT213 Model

This project is focused on developing an experimental framework for characterizing non-associative plasticity flow rule coefficients through coupon-scale tests for the LS-DYNA MAT213 model. The main objective is to characterize these coefficients based on the multi-scale (i.e., both microscopic and macroscopic) full-field measurement of the evolution of strain and stress fields. This paper focuses on presenting the experimental work on characterizing the full-scale stress-strain curves of T700/LM-PAEK composites under tension, compression, and shear loads. The experimental data set was intended to build a deformation sub-model in the MAT213 model for the material. The strain data were collected using both microscopic and macroscopic digital image correlation techniques. The microscopic technique was particularly useful for fracture cases under small strains. A preliminary simulation result obtained from the MAT213 model is also presented in the paper. The experimental framework herein will be extended to characterize post-peak stress degradation in the composite material and to develop a damage sub-model for the material. This project will contribute to developing a simulation tool based on the MAT213 model for simulating the rate-dependent impact damage in composites under multi-axial loading.

MAT213↗

DOE OSTI2021

Identification of Life Models for Li-Ion Batteries Using Penalized Regression and Bilevel Optimization

Reduced-order physics-based life models are extremely useful for rapidly predicting battery state-of-health and for simulating battery lifetime in arbitrary aging conditions. However, identification of well-parameterized models is difficult. This is because, for maximum usefulness in predicting lifetime under a variety of conditions, aging test data exhibits many degradation mechanisms, which all need to be accurately modeled. However, because aging tests are time-consuming and expensive, especially for large-format batteries, a minimum of tests are conducted while probing many stress factors. Building a well-parameterized model is then very challenging: an under-parameterized model will neglect critical degradation modes, and an over-parameterized model will extrapolate poorly to new testing conditions. To complicate this matter, the functional form of the model for any individual degradation rate can be very difficult to identify. In this work, the statistical tools of penalized regression and bilevel optimization are used to help identify both the functional forms of and optimize the parameters of reduced-order life models, accelerating identification of robust models. Model robustness is demonstrated through traditional statistics methods of cross-validation and sensitivity analysis, uncertainty quantification through bootstrap resampling and Monte-Carlo simulation, and simulation of real-world use cases.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Binary classification with noisy quantum circuits and noisy quantum data

We study the effects of single-qubit noises in the quantum circuit and the corruption in the quantum training data to the performance of binary classification problem. We find that under the presence of errors, the measurement at a qubit is affected only by the noise in the corresponding qubit, and that the errors on other qubits do not affect this outcome. Furthermore, for the task where we fit a binary classifier using a quantum training data, we show that the noise in the data can work as a regularizer, implying that we can benefit from the noise in certain cases. We support our findings with simulations.

Lee, Yonghoon↗

DOE OSTI2023

The Performance of FBMC-SS in Challenging HF Channel Environments

Provide overview of FBMC-SS waveform Discuss parameters used for simulated results herein Present simulated results in challenging channel environments Cases: 2 path channel with 2nd path frequency offset 2 path MLD and “VP” channels Large narrowband interferer Highly congested interference case Highly dynamic interference case Clarify previously presented results Make conclusions from presented results

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2023

Using Co-Optimized Machine Learned Manifolds for Modeling Chemically Reacting Flows

Chemically reacting flows play a key role in a wide range of engineered systems, from chemical and polymer processing to combustion-based energy conversion technologies. Simulations of these flows involve solving a coupled set of partial differential equations for mass, momentum, energy, and all relevant chemical species in the system. Chemical reaction pathways may be extremely complex and involve hundreds or more intermediate species, with reactions that occur over timescales varying by several orders of magnitude - presenting a significant numerical stiffness challenge. The combination of these factors makes simulation of chemically reacting flows vastly more expensive than nonreactive simulations, and often makes direct solution of the governing equations intractable. It is necessary to apply lower-fidelity models in place of the detailed governing equations in order to reduce computational cost to enable reacting flow simulation tools to be used in the engineering design process. Many of the models employed for this purpose are based on reducing the dimension of the thermochemical state, motivated by the observation that the observed thermochemical states in a system lie on a low-dimensional manifold in thermochemical state space. This behavior occurs due to the fast equilibration of certain reactive and transport processes, and physics-based manifold models rely on idealized assumptions about the balance of timescales and the way in which chemistry and transport are coupled. In this work, we apply a novel method for data-driven manifold-based modeling that can leverage data from high-fidelity reacting flow simulations to improve model accuracy in cases where the physics-based modeling assumptions break down. The approach is designed to be broadly applicable across chemically reacting flow systems but is applied here to turbulent combustion modeling.

machine learning↗

DOE OSTI2022

Multi-Physics Coupled Seismic Safety Analysis of Molten-Salt Reactors

Pool-type Molten Salt Reactors (MSRs) are a promising advanced nuclear reactor concept relying on passive physical behavior to provide increased safety characteristics. The high heat capacity, high boiling point, unpressurized liquid salt contains the nuclear fuel and as it passes through the vessel achieves criticality and produces nuclear power. The salt then carries away the heat produced along with the delayed neutron precursors towards heat exchangers and primary pumps before entering the core. This forms a highly coupled multiphysics problem between neutronics and fluid flow, which makes modeling these reactors challenging. The safety/licensing case heavily relies on modeling and simulation to prove their safety, and significant effort has been expanded by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to develop the appropriate simulation tools. Their safety/licensing will require a thorough analysis of their behavior during earthquake-based transients. Earthquakes are a transient condition commonly experienced in many regions, and are safely survived by dozens of reactors every year. Nonetheless, significant conservatisms were introduced with limited predictive modeling and simulation available. The first-of-a-kind capability developed by this research will allow for engineers to analyze ranges of design-basis and beyond-design basis accident scenarios to assess the integrity of the reactor vessel and surrounding system, possibly allowing for limiting the costly conservatisms in the design. This projects uses MASTODON, Griffin, and Pronghorn in the NEAMS ecosystem to model the coupled multiphysics problem posed by an earthquake in a molten salt reactor. It marches through a number of number of coupling schemes, from uncoupled simulations to tightly two-way coupled simulations. The appropriate level of coupling for these simulations will then be determined based on the accuracy in the quantities of interest (QoI) and computational effort involved in each scheme. The QoIs chosen for this paper include the mass flow rate across the heat exchanger (for Pronghorn simulations), power output of the reactor (for Griffin simulations), and the maximum Von Mises stress (for MASTODON simulations). This summary paper reports on the current progress of this project, with 2D tightly coupled multiphysics results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Development of an Optimal Variable-Pitch Controller for Floating Axial-Flow Marine Hydrokinetic Turbines: Preprint

This article discusses the development of an optimal variable-pitch controller for floating, axial-flow marine turbines. Recently, OpenFAST, an open-source wind turbine modeling tool, has been extended to model marine turbines. A controller is necessary to simulate marine turbines for different load cases using OpenFAST, which greatly impacts the performance of the energy system. Previous studies have designed controllers using a linearized model of the marine turbine, which can be timeconsuming and require the expertise of a control engineer. In this study, we use an automated approach that uses generic models of the marine turbine to identify the controller gains, which can expedite the process of designing a controller. Using an optimizer to identify the control system parameters can additionally improve the controller's performance. The optimal controller tuned using such an approach results in a 20% reduction in the towerbase damage equivalent loading and better tracking of the rated generator speed and power.

closed-loop control↗

DOE OSTI2025

Power Quality and Load Capacity Evaluations of an Electric Vehicle for Multi-Robot System Applications

This paper evaluates the capability of a fully electric pickup truck, using the Ford F-150 Lightning as an example, to provide power to the circuit of a multi-robot system. The case study was conducted on a simulated INL Autonomous Pit Exploration System (APES) designed for the inspection of nuclear waste tank pits. Through a series of controlled tests, the vehicle’s power delivery consistency, load-handling capability, and battery performance were assessed under various conditions. First of all, the load test demonstrated that the vehicle provided stable power with low distortion and no unexpected interruptions. Second, during the operational limit test, the 240V system sustained loads up to 7.4 kW before tripping, providing insights into its operational limits. Last but not least, during a simulated full-scale APES operation, the vehicle’s battery depleted by only 6% over an hour, indicating sufficient capacity for extended use while retaining reserve power for transportation needs. This study highlights the potential of electric vehicles as reliable power sources for field operations, contributing to the advancement of sustainable technologies by reducing reliance on traditional fossil fuel generators and promoting the integration of clean energy solutions in remote and challenging environments.

Electric vehicle↗

DOE OSTI2025

IEEE SusTech 2025 Presentation

This paper evaluates the capability of a fully electric pickup truck, using the Ford F-150 Lightning as an example, to provide power to the circuit of a multi-robot system. The case study was conducted on a simulated INL Autonomous Pit Exploration System (APES) designed for the inspection of nuclear waste tank pits. Through a series of controlled tests, the vehicle’s power delivery consistency, load-handling capability, and battery performance were assessed under various conditions. First of all, the load test demonstrated that the vehicle provided stable power with low distortion and no unexpected interruptions. Second, during the operational limit test, the 240V system sustained loads up to 7.4 kW before tripping, providing insights into its operational limits. Last but not least, during a simulated full-scale APES operation, the vehicle’s battery depleted by only 6\% over an hour, indicating sufficient capacity for extended use while retaining reserve power for transportation needs. This study highlights the potential of electric vehicles as reliable power sources for field operations, contributing to the advancement of sustainable technologies by reducing reliance on traditional fossil fuel generators and promoting the integration of clean energy solutions in remote and challenging environments.

42 - ENGINEERING↗

DOE OSTI2025

Recent Improvements in Pronghorn for Advanced Reactor Modeling

Pronghorn is a thermal-hydraulics computational tool developed using the Idaho National Laboratory's Multiphysics Object-Oriented Simulation Environment (MOOSE). It is designed to support Computational Fluid Dynamics (CFD) modeling, ranging from subchannel and porous media analysis to Reynolds Averaged Navier-Stokes (RANS) turbulence modeling. As an integral part of the MOOSE-based suite of tools, Pronghorn seamlessly couples with other MOOSE-based applications to simulate a variety of physical phenomena. This article highlights recent significant enhancements to Pronghorn's CFD modeling capabilities and demonstrates their application to advanced nuclear reactor designs. The recent improvements in Pronghorn primarily focus on modifications to its turbulence modeling capabilities, near-wall corrections and numerical schemes. In terms of turbulence modeling, the two-equation $k-\epsilon$ and $k-\omega$ SST models have been implemented and validated with both equilibrium and non-equilibrium wall treatments. Additionally, corrections for wall roughness, and curvature, and wall-channeling in pebble beds have been introduced in the near-wall modeling. These developments enable more accurate simulations of advanced nuclear reactors. Two case studies are presented in this work: a pool-type Molten Chloride Reactor and a salt-cooled Pebble-Bed High Temperature Reactor. In both cases, the previous models in Pronghorn are compared with the new implementations, demonstrating the improved accuracy achieved with the updated models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

High-speed impact of micron-sized diesel drop trains—Splashing dynamics, secondary droplet formation, and effects of pre-existing film thickness

This work investigates surface impingement of monodispersed diesel drop trains using computational simulations. The three-dimensional simulations include high-speed impact of micrometer-sized fuel drops onto initially dry and wetted stainless steel substrates. The drop size and impact velocity are representative of fuel injection conditions in internal combustion engines (ICEs). The drop trains serve as a simplified representation of fuel spray. To accurately capture the film that naturally develops on initially dry substrates, a dynamic diesel-stainless steel contact angle model was developed using experiments of single drop impact under ICE representative temperature. Simulations are presented for a highly splashing case and the effects of a pre-existing film on splashing dynamics were investigated, including the temporal evolution of splashed mass and film thickness. It was concluded that for thin films, the effect of pre-existing film thickness is less significant than in single drop impingement. Secondary droplet characterization was performed on simulation results, leading to instantaneous and time-averaged distributions of secondary droplet size, velocity magnitude, and trajectory angle. It was found that for each drop impingement, approximately 58% of the splashed mass is from the impinging drop itself, while the remainder of the splashed mass is composed of film liquid. From a high resolution case, details of secondary droplet formation are observed and three distinct phases of secondary droplet formation are identified. The detailed analysis of drop train impingement under engine-relevant conditions serves as a first step toward a robust understanding of fuel impingement and the development of cleaner and more efficient ICEs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Measurement Uncertainty of Highly Asymmetrically Curved Elliptical Mirrors Using Multi-Pitch Slope Stitching Technique

Soft X-ray off-axis elliptical mirrors bring new challenges for X-ray mirror metrology. These highly asymmetrically curved elliptical cylindrical mirrors with a total slope range >10 mrad are extremely challenging to measure. Their total slope range exceeds the measuring range of most angular sensors used for X-ray mirror inspection. To overcome this problem, it is possible to stitch partial slope data by measuring the mirror at different pitch angles (multi-pitch angles). By revisiting the theory of the multi-pitch Nano-accuracy Surface Profiler (NSP), we derive the sampling position error on the mirror surface as a function of the mirror height profile and the measurement error of the pitch rotation center. When measuring “extreme”, highly asymmetrically curved, elliptical mirrors, the calculation of the mirror height profile with iterative reconstruction outperforms the classical “flat assumption” (i.e., assuming that the mirror sag is negligible). As demonstrated by our simulations, a proper tolerance evaluation on the measurement of pitch rotation center is needed to assess the measurement accuracy (systematic error) for these strongly aspherical mirrors using the multi-pitch NSP technique. Taking a real design of an “extreme” elliptical mirror as a case study, we conduct a Monte Carlo simulation to mimic the measurement and characterization process to analyze the impact of several error sources. With the measurement uncertainty of the pitch rotation center, the multi-pitch NSP measurement can estimate the grazing angle θ and the chief ray location x o with their uncertainties, as well as the slope residuals.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Criticality Analysis of FSV Spent Nuclear Fuel in the DOE Standard Canister

The U.S. Department of Energy (DOE) is responsible for managing over 300 types of spent nuclear fuel (SNF). To manage this large variety of fuel types, DOE plans to employ standardized canisters for the transportation, long-term storage, and eventual disposal of SNF. Idaho National Laboratory is currently supporting DOE’s SNF Packaging Demonstration Project, in which Fort Saint Vrain (FSV) fuel assemblies will be loaded into a DOE Standard Canister. This paper presents criticality calculations demonstrating that all four or five FSV fuel assemblies loaded into the DOE Standard Canister will remain subcritical in any expected normal or credible abnormal conditions. Previous criticality analyses were performed for one FSV fuel assembly and 12 Peach Bottom Core 2 fuel elements loaded into a DOE Standard Canister. This paper covers the criticality analysis performed for loading both four and five FSV fuel assemblies into a DOE Standard Canister. Various intact and degraded mode configurations were modeled in conducting the criticality calculations. This analysis encompassed three different configurations: (1) a single DOE Standard Canister loaded into a concrete storage overpack, (2) seven DOE Standard Canisters loaded into a concrete storage overpack, and (3) nine DOE Standard Canisters loaded into a concrete storage overpack. The overpack dimensions were varied for each of the three configurations, and transport, storage, and disposal scenarios were analyzed for each configuration. For transport scenarios, a pair of degradation cases was analyzed. In the first case, the fuel compacts became degraded and were removed from the fuel block, then deposited at the bottom of a horizontally placed canister, thereby simulating a drop event. The canister was considered to remain intact. In the second case, the spacing between horizontally placed canisters in a nine-canister overpack was reduced such that the canisters were piled on top of each other, simulating a drop event. For this case, no degradation of the canister internals or fuel was considered. For storage scenarios, the water moderator location in the system was varied to enable identification of the most reactive configurations. Dry and wet conditions were analyzed for the fuel materials, canister, and overpack. For disposal scenarios, two degradation cases were analyzed. In the first, the stainless-steel internals of the canister degraded to either hematite or goethite under both dry and wet conditions. In the second case, degraded FSV fuel formed a uranium-water slurry that filled the coolant/void holes. None of the cases presented exceeded the application specific upper subcritical limit.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Backlighting extended gas halos around luminous red galaxies: Kinematic Sunyaev-Zel’dovich effect from DESI Y1 and ACT data

The gas density profile around galaxies, shaped by feedback and affecting the galaxy lensing signal, is imprinted on the cosmic microwave background (CMB) by the kinematic Sunyaev-Zel’dovich effect (kSZ). We precisely measure this effect (𝑆/𝑁 ≈ 10) via velocity stacking with 825,283 spectroscopically confirmed luminous red galaxies (LRG) from the Dark Energy Spectroscopic Instrument Year 1 (DESI Y1) survey, which overlap with the Atacama Cosmology Telescope (ACT) Data Release 6 temperature maps over ≥ 4,000 deg 2 . We explore the kSZ dependence with various galaxy parameters and find no significant trend with redshift but clear trends with stellar mass and absolute magnitude in 𝑔, 𝑟, and 𝑧 bands. Our analysis suggests that the gas extends beyond the dark matter halo (99.5% confidence level, i.e., probability to exceed (PTE)⁢ = 0.005). We find a tentative preference for hydrodynamical simulation models with stronger feedback that drives gas further out (Illustris 𝑧 = 0.5, PTE⁢ = 0.37) over weaker-feedback cases (IllustrisTNG 𝑧 = 0.8, PTE⁢ = 0.045), though with limited statistical significance. In all cases, a free multiplicative amplitude was fit to the simulated profiles, and further modeling work is required to firm up these conclusions. We find consistency between kSZ profiles around spectroscopic and photometric LRG, with comparable statistical power, thus increasing our confidence in the photometric analysis. Additionally, we present the first kSZ measurement around the DESI Y1 bright galaxy sample (BGS) and the emission-line galaxies (ELG) whose features match qualitative expectations. Finally, we forecast 𝑆/𝑁 ∼ 50 for future stacked kSZ measurements using data from the ACT, the DESI Y3, and the Rubin Observatory. Furthermore, these measurements will serve as an input for galaxy formation models and baryonic uncertainties in galaxy lensing.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Solving the Orszag–Tang vortex magnetohydrodynamics problem with physics-constrained convolutional neural networks

We study the 2D Orszag–Tang vortex magnetohydrodynamics (MHD) problem through the use of physics-constrained convolutional neural networks (PCNNs) for forecasting the density, ρ, and the magnetic field, B, as well as the prediction of B given the velocity field v of the fluid. In addition to translation equivariance from the convolutional architecture, other physics constraints were embedded: absence of magnetic monopoles, non-negativity of ρ, use of only relevant variables, and the periodic boundary conditions of the problem. The use of only relevant variables and the hard constraint of non-negative ρ were found to facilitate learning greatly. The divergenceless condition ∇·B=0 was implemented as a hard constraint up to machine precision through the use of a magnetic potential to define B=∇×A. Residual networks and data augmentation were also used to improve performance. This allowed for some of the residual models to function as surrogate models and provide reasonably accurate simulations. For the prediction task, the PCNNs were evaluated against a physics-informed neural network, which had the ideal MHD induction equation as a soft constraint. Several models were able to generate highly accurate fields, which are visually almost indistinguishable and have low mean squared error. Only methods with built-in hard constraints produced physical fields with ∇·B=0. The use of PCNNs for MHD has the potential to produce physically consistent real-time simulations to serve as virtual diagnostics in cases where inferences must be made with limited observables.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Hanford Site Composite Analysis: Recharge and Inventory Sensitivity Analyses for the A Trenches, BC Cribs and Trenches, and PUREX Area Vadose Zone Models

The main objectives of the composite analysis (CA) compliance case vadose zone (VZ) facet are to simulate transport of radionuclide releases from the surface to the water table, and to provide radionuclide transfer rates to groundwater for inclusion into the Plateau-to-River (P2R) saturated zone (SZ) model (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). The CA compliance case approach to address the scale and distribution of contaminant sources in the VZ facet is to subdivide the Hanford Site Central Plateau into 26 models that contain contaminant sources and liquid discharges likely to commingle during migration through the VZ to the water table

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

A sharp interface Lagrangian-Eulerian method for rigid-body fluid-structure interaction

Herein this paper introduces a sharp interface method to simulate fluid-structure interaction (FSI) involving rigid bodies immersed in viscous incompressible fluids. The capabilities of this methodology are benchmarked using a range of test cases and demonstrated using large-scale models of biomedical FSI. The numerical approach developed herein, which we refer to as an immersed Lagrangian-Eulerian (ILE) method, integrates aspects of partitioned and immersed FSI formulations by solving separate momentum equations for the fluid and solid subdomains, as in a partitioned formulation, while also using non-conforming discretizations of the dynamic fluid and structure regions, as in an immersed formulation. A simple Dirichlet-Neumann coupling scheme is used, in which the motion of the immersed solid is driven by fluid traction forces evaluated along the fluid-structure interface, and the motion of the fluid along that interface is constrained to match the solid velocity and thereby satisfy the no-slip condition. To develop a practical numerical method, we adopt a penalty approach that approximately imposes the no-slip condition along the fluid-structure interface. In the coupling strategy, a separate discretization of the fluid-structure interface is tethered to the volumetric solid mesh via stiff spring-like penalty forces. Our fluid-structure coupling scheme relies on an immersed interface method (IIM) for discrete geometries, which enables the accurate determination of both velocities and stresses along complex internal interfaces. Numerical methods for FSI can suffer from instabilities related to the added mass effect, but computational tests indicate that the methodology introduced here remains stable for selected test cases across a broad range of solid-fluid mass density ratios, including extremely small, nearly equal, equal, and large density ratios. Biomedical FSI demonstration cases include results obtained using this method to simulate the dynamics of a bileaflet mechanical heart valve in a pulse duplicator, and to model transport of blood clots in a patient-averaged anatomical model of the inferior vena cava.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1998

Simulation of Ultra-Small MOSFETs Using a 2-D Quantum-Corrected Drift-Diffusion Model

The continued down-scaling of electronic devices, in particular the commercially dominant MOSFET, will force a fundamental change in the process of new electronics technology development in the next five to ten years. The cost of developing new technology generations is soaring along with the price of new fabrication facilities, even as competitive pressure intensifies to bring this new technology to market faster than ever before. To reduce cost and time to market, device simulation must become a more fundamental, indeed dominant, part of the technology development cycle. In order to produce these benefits, simulation accuracy must improve markedly. At the same time, device physics will become more complex, with the rapid increase in various small-geometry and quantum effects. This work describes both an approach to device simulator development and a physical model which advance the effort to meet the tremendous electronic device simulation challenge described above. The device simulation approach is to specify the physical model at a high level to a general-purpose (but highly efficient) partial differential equation solver (in this case PROPHET, developed by Lucent Technologies), which then simulates the model in 1-D, 2-D, or 3-D for a specified device and test regime. This approach allows for the rapid investigation of a wide range of device models and effects, which is certainly essential for device simulation to catch up with, and then stay ahead of, electronic device technology of the present and future. The physical device model used in this work is the density-gradient (DG) quantum correction to the drift-diffusion model [Ancona, Phys. Rev. B 35(5), 7959 (1987)]. This model adds tunneling and quantum smoothing of carrier density profiles to the drift-diffusion model. We used the DG model in 1-D and 2-D (for the first time) to simulate both bipolar and unipolar devices. Simulations of heavily-doped, short-base diodes indicated that the DG quantum corrections do not have a large effect on the IN characteristics of electronic devices without heteroj unction s. On the other hand, ultra-small MOSFETs certainly exhibit important quantum effects that the DG model will include: quantum repulsion of the inversion and gate charges from the oxide interfaces, and quantum tunneling through thin gate oxides. We present initial results of 2-D DG simulations of ultra-small MOSFETs. Subtle but important issues involving the specification of the model, boundary conditions, and interface constraints for DG simulation of MOSFETs will also be illuminated.

Biegal, Bryan A.↗