Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state”

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 775 records · Page 43

Modeling the Role of Dislocation Substructure During Class M and Exponential Creep

The different substructures that form in the power-law and exponential creep regimes for single phase crystalline materials under various conditions of stress, temperature and strain are reviewed. The microstructure is correlated both qualitatively and quantitatively with power-law and exponential creep as well as with steady state and non-steady state deformation behavior. These observations suggest that creep is influenced by a complex interaction between several elements of the microstructure, such as dislocations, cells and subgrains. The stability of the creep substructure is examined in both of these creep regimes during stress and temperature change experiments. These observations are rationalized on the basis of a phenomenological model, where normal primary creep is interpreted as a series of constant structure exponential creep rate-stress relationships. The implications of this viewpoint on the magnitude of the stress exponent and steady state behavior are discussed. A theory is developed to predict the macroscopic creep behavior of a single phase material using quantitative microstructural data. In this technique the thermally activated deformation mechanisms proposed by dislocation physics are interlinked with a previously developed multiphase, three-dimensional. dislocation substructure creep model. This procedure leads to several coupled differential equations interrelating macroscopic creep plasticity with microstructural evolution.

Raj, S. V.↗

Silicon Carbide Diodes Characterization at High Temperature and Comparison With Silicon Devices

Commercially available silicon carbide (SiC) Schottky diodes from different manufacturers rated at 200, 300, 600, and 1200 V, were electrically tested and characterized as a function of temperature up to 300 C. Electrical tests included both steady state and dynamic tests. Steady state tests produced forward and reverse I-V characteristic curves. Transient tests evaluated the switching performance of the diodes in either a hard-switched DC to DC buck converter or a half-bridge boost converter. For evaluation and comparison purposes, the same tests were performed with current state-of-the-art ultra fast silicon (Si) pn-junction diodes of similar ratings and also a Si Schottky diode. The comparisons made were forward voltage drop at rated current, reverse current at rated voltage, and turn-off peak reverse recovery current and reverse recovery time. In addition, efficiency measurements were taken for the buck DC to DC converter using both the SiC Schottky diodes and the Si pn-junction diodes at different temperatures and frequencies. The test results showed that at high temperature, the forward voltage drop for SiC Schottky diodes is higher than the forward drop of the ultra fast Si pn-junction diodes. As the temperature increased, the forward voltage drop of the SiC Schottky increased while for the ultra fast Si pn-junction diodes, the forward voltage drop decreased as temperature increased. For the elevated temperature steady state reverse voltage tests, the SiC Schottky diodes showed low leakage current at their rated voltage. Likewise, for the transient tests, the SiC Schottky diodes displayed low reverse recovery currents over the range of temperatures tested. Conversely, the Si pn-junction diodes showed increasing peak reverse current values and reverse recovery times with increasing temperature. Efficiency measurements in the DC to DC buck converter showed the advantage of the SiC Schottky diodes over the ultra fast Si pn-junction diodes, especially at the higher temperatures and higher frequencies.

Lebron-Velilla, Ramon C.↗

Numerical Calculation of the Drag Force Acting on a Solid Particle Pushed by a Solid/Liquid Interface

The distribution of insoluble particles in a metal casting depends primarily on the interaction of the particles with the solid/liquid interface (SLI) during the solidification process. The balance of the forces acting on the particle essentially determines whether a particle will be engulfed or pushed by the SLI. An important component of this force balance is the drag force generated by the particle motion in front of the SLI. Previously developed mathematical models for particle/SLI interaction made use of steady-state solutions of this force provided by the lubrication theory. However, our numerical model based on the SLI tracking approach shows that the steady-state approach is inappropriate to model the interaction process and that at steady-state the theoretical solution underestimates the drag force. It was found that regression analysis of steady-state numerical solutions for cylindrical particles moving normal to a flat SLI gives a relationship of the form: Abstract The distribution of insoluble particles in a metal casting depends primarily on the interaction of the particles with the solid/liquid interface (SLI) during the solidification process. The balance of the forces acting on the particle essentially determines whether a particle will be engulfed or pushed by the SLI. An important component of this force balance is the drag force generated by the particle motion in front of the SLI. Previously developed mathematical models for particle/SLI interaction made use of steady-state solutions of this force provided by the lubrication theory. However, our numerical model based on the SLI tracking approach shows that the steady-state approach is inappropriate to model the interaction process and that at steady-state the theoretical solution underestimates the drag force. It was found that regression analysis of steady-state numerical solutions for cylindrical particles moving normal to a flat SLI gives a relationship of the form: F(sub D, sup num) =sqoare root of 3(pi)(eta)V(sub p)(R(sub p)/d)(sup 10(gamma)/3). This is to compared to the solution provided by the lubrication theory: F(sub D), sup theor) = 3 square root of 2(pi)(eta)V(sub p)(R(sub p)/d)(sub 3/2), where F(sub D) is the drag force, eta is the dynamic viscosity of the fluid, V(sub p) is the particle velocity, R(sub p) is the particle radius, gamma is Euler's constant (gamma approx. = 0.577), and d is the width of the gap between the particle and the SLI. The domain on which the equations proposed by the lubrication theory are relevant is clearly identified in the paper. The numerical model was then validated against the classical lubrication theory within the domain of its validity.

Catalina, Adrian V.↗

NASA/FAA Tailplane Icing Program: Flight Test Report

This report presents results from research flights that explored the characteristics of an ice-contaminated tailplane using various simulated ice shapes attached to the leading edge of the horizontal tailplane. A clean leading edge provided the baseline case, then three ice shapes were flown in order of increasing severity. Flight tests included both steady state and dynamic maneuvers. The steady state points were 1G wings level and steady heading sideslips. The primary dynamic maneuvers were pushovers to various G-levels; elevator doublets; and thrust transitions. These maneuvers were conducted for a full range of flap positions and aircraft angle of attack where possible. The analysis of this data set has clearly demonstrated the detrimental effects of ice contamination on aircraft stability and controllability. Paths to tailplane stall were revealed through parameter isolation and transition studies. These paths are (1) increasing ice shape severity, (2) increasing flap deflection, (3) high or low speeds, depending on whether the aircraft is in a steady state (high speed) or pushover maneuver (low speed), and (4) increasing thrust. The flight research effort was very comprehensive, but did not examine effects of tailplane design and location, or other aircraft geometry configuration effects. However, this effort provided the role of some of the parameters in promoting tailplane stall. The lessons learned will provide guidance to regulatory agencies, aircraft manufacturers, and operators on ice-contaminated tailplane stall in the effort to increase aviation safety and reduce the fatal accident rate.

Ratvasky, Thomas P.↗

Small Engine Technology. Task 4: Advanced Small Turboshaft Compressor (ASTC) Performance and Range Investigation

This contact had two main objectives involving both numerical and experimental investigations of a small highly loaded two-stage axial compressor designated Advanced Small Turboshaft Compressor (ASTC) winch had a design pressure ratio goal of 5:1 at a flowrate of 10.53 lbm/s. The first objective was to conduct 3-D Navier Stokes multistage analyses of the ASTC using several different flow modelling schemes. The second main objective was to complete a numerical/experimental investigation into stall range enhancement of the ASTC. This compressor was designed wider a cooperative Space Act Agreement and all testing was completed at NASA Lewis Research Center. For the multistage analyses, four different flow model schemes were used, namely: (1) steady-state ADPAC analysis, (2) unsteady ADPAC analysis, (3) steady-state APNASA analysis, and (4) steady state OCOM3D analysis. The results of all the predictions were compared to the experimental data. The steady-state ADPAC and APNASA codes predicted similar overall performance and produced good agreement with data, however the blade row performance and flowfield details were quite different. In general, it can be concluded that the APNASA average-passage code does a better job of predicting the performance and flowfield details of the highly loaded ASTC compressor.

Hansen, Jeff L.↗

The modeling of gas phase permeation through iron and nickel membranes

The gas phase permeation of hydrogen through metal membranes encompasses many kinetic processes. This paper reviews a permeation model which incorporates second order gas-surface reaction kinetics with simple bulk diffusion. The model is used to investigate the effect of this particular surface reaction of steady-state permeation. The dependence of the steady-state permeation flux on temperature, pressure, and thickness of the membrane has been calculated. The model predicts that the bulk controlled steady-state flux will change to a surface limited steady-state flux as either the temperature or thickness of the membrane is reduced. Finally, using independently derived parameters, the model is compared with permeation measurements on iron and nickel membranes.

Kuhn, David K.↗

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic reconnection in the presence of sheared flow and density asymmetry: Applications to the Earth's magnetopause

Classical models of magnetic reconnection consist of a small diffusion region bounded by two symmetric slow shocks, across which the plasma is accelerated. Asymmetries often present in space plasmas are sheared plasma flow and dissimilar plasma densities on the two sides of current sheets. In this paper, we investigate magnetic reconnection in the presence of a shear flow and an asymmetric density across the current sheet using two-dimensional magnetohydrodynamic (MHD) simulations. The results demonstrate that magnetic reconnection can occur only for a plasma flow velocity (in the frame of the X line) which is below the Alfven speed in each inflow region. This limits the velocity of the X line to a certain range for a given flow shear and provides an upper limit to the total velocity shear at which reconnection ceases to operate. Depending on the direction of the flow in the adjacent inflow region, the effects from the sheared flow and from the density asymmetry will compete with or enhance each other in respect to the magnitude and location of the currents which bound the outflow regions. The results are applied to the dayside and flank regions of the magnetosphere. For the dayside region where the magnetosheath flow is slow, the magnetic field transition region is thin and the accelerated flow is earthward of the sharp current layer (magnetopause). At the flanks tailward of the X line, shear flow and density asymmetry effects compete making the magnetic field transition layer broad with the high-speed flow contained within the transition region which explains corresponding observations. At the flanks sunward of the X line, shear flow and density asymmetry effects enhance each other and lead to a strong current sheet on the magnetosheath side of the accelerated flow. The total volume affected by magnetic reconnection is much larger than the steady state region. A large bulge region precedes the steady state region. Qualitatively, the bulge and the steady state region have similar signatures and both can explain observations. We provide criteria in order to distinguish between the bulge and the steady state region in observations.

La Belle-Hamer, A. L.↗

Solving quasi-free and quadratic Lindblad master equations for open fermionic and bosonic systems

The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), here we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k-point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stabilizing Non-Abelian Topological Order Against Heralded Noise via Local Lindbladian Dynamics

An important open question for the current generation of highly controllable quantum devices is understanding which phases can be realized as stable steady states under local quantum dynamics. In this work, we show how robust steady-state phases with both Abelian and non-Abelian mixed-state topological order can be stabilized, in two spatial dimensions, against generic “heralded” noise using active dynamics that incorporate measurement and feedback, modeled as a fully local Lindblad master equation. These topologically ordered steady states are two-way connected to pure topologically ordered ground states using local quantum channels, and preserve quantum information for a time that is exponentially large in the system size. Specifically, we present explicit constructions of families of local Lindbladians for both Abelian (ℤ 2 ) and non-Abelian (𝐷 4 ) topological order whose steady states host mixed-state topological order when the noise is below a threshold strength. As the noise strength is increased, these models exhibit first-order transitions to intermediate mixed-state phases where they encode robust classical memories, followed by (first-order) transitions to a trivial steady state at high noise rates. When the noise is imperfectly heralded, steady-state order disappears but our active dynamics significantly enhances the lifetime of the encoded logical information. To carry out the numerical simulations for the non-Abelian 𝐷 4 case, we introduce a generalized stabilizer tableau formalism that permits efficient simulation of the non-Abelian Lindbladian dynamics.

Monte Carlo methods↗

Thermo-hydraulic steam pipe models for district heating simulations: Simplifications to balance accuracy and simulation speed

Steam piping networks are essential for optimizing performance in industrial processes and district heating systems. However, dynamic models that balance thermo-hydraulic accuracy with computational efficiency remain limited. In response, this paper presents a new discretized steam pipe model based on the plug flow approach, capturing key thermo-hydraulic behaviors while simplifying steam phase change processes. Implemented in Modelica, the model accurately calculates temperature and pressure distributions along steam pipelines. To improve computational efficiency for district-scale simulations, five model simplifications are introduced: lumped thermo-hydraulic functions, empirical correlations, fluid state approximations, steady-state dynamics and inclusion of flow derivatives. These simplified models achieve 85%-98% accuracy in predicting pressure drop and condensation losses, including dynamic condensate behavior during pipe warm-up—a factor often overlooked in existing models. The models support diverse network configurations, scaling effectively to systems with multiple distribution pipes and connected building loads. Discrete models provide detailed insights but exhibit a cubic increase in simulation time as the network scales by N connected building O(N 2.42 ). In contrast, lumped models simulate 10–28 times faster than discrete, offering quadratic scaling of simulation time O(N 1.73 ). However, they still require 6 times more computation time than a lossless network, highlighting the inherent computational challenges of modeling compressible fluid flow. In conclusion, the steady-state lumped variant, with its near-linear scalability in computational time O(N 1.01 ), emerges as an efficient solution for preliminary design evaluations and extensive parametric studies.

15 GEOTHERMAL ENERGY↗

Electric propulsion interactive effects with spacecraft science payloads.

Interactive effects between electric thruster units and a Solar Electric Pioneer (SEP) spacecraft science payload are investigated. Studied effects include magnetic contamination (steady state), electric contamination (steady state and time varying) and particle contamination of regions in and around the spacecraft. Conducted electromagnetic interference is measured for thrusters powered with prototype power processing units and radio frequency interference effects between the thruster and the omnidirectional low gain antenna are determined. Material deposition and transport of thruster propellant ions and atoms and charge exchange metal ions and atoms are also experimentally determined. Large scale environmental effects between the plasma thrust beam and the particles and fields in interplanetary space are discussed. Contaminant reduction techniques are proposed, and both laboratory and in-space interactive effects tests are suggested.

Sellen, J. M., Jr.↗

A computer simulation of the midlatitude plasmasphere and ionosphere

A computer model has been developed to simulate species density, temperature, and plasma flow in the ionosphere and plasmasphere. A new approach, the flux preserving method involving a Newton iteration, is used to solve a system of governing equations comprising four second-order partial differentials. Tests are performed to demonstrate that the simulation converges to a stable steady state solution, then steady state ion fluxes are analyzed. Finally, simulations of the collapse of the sunset ionosphere are presented. Comparisons with satellite and radar data show good agreement in a number of cases.

Young, E. R.↗

Modeling and characterization of the Earth Radiation Budget Experiment (ERBE) nonscanner and scanner sensors

The Earth Radiation Budget Experiment (ERBE) is making high-absolute-accuracy measurements of the reflected solar and Earth-emitted radiation as well as the incoming solar radiation from three satellites: ERBS, NOAA-9, and NOAA-10. Each satellite has four Earth-looking nonscanning radiometers and three scanning radiometers. A fifth nonscanner, the solar monitor, measures the incoming solar radiation. The development of the ERBE sensor characterization procedures are described using the calibration data for each of the Earth-looking nonscanners and scanners. Sensor models for the ERBE radiometers are developed including the radiative exchange, conductive heat flow, and electronics processing for transient and steady state conditions. The steady state models are used to interpret the sensor outputs, resulting in the data reduction algorithms for the ERBE instruments. Both ground calibration and flight calibration procedures are treated and analyzed. The ground and flight calibration coefficients for the data reduction algorithms are presented.

Halyo, Nesim↗

Integrated flow and structural modeling for rocket engine component test facility propellant systems

A set of PC-based computational Dynamic Fluid Flow Simulation models is presented for modeling facility gas and cryogenic systems. Data obtained provide important information regarding performance envelope parameters for the facility using different engine components; time-dependent valve setting for controlling steady-state, quasi-steady state, and transient profiles; optimum facility pipe and pipe component sizes and parameters; momentum transfer loads; and fluid conditions at critical points. A set of COSMIC NASTRAN-based finite element models is also presented to evaluate the loads and stresses on test facility piping systems from fluid and gaseous effects, thermal chill down, and occasional wind loads. The models are based on Apple Macintosh software which makes it possible to change numerous parameters.

Dequay, L.↗

Micromechanical studies of composites by BEM

Over the last five years BEM has been developed into a very effective tool for micromechanical studies of composites. The analyses developed range from static elastic analysis, steady state heat transfer, steady state and transient thermoelastic analysis, periodic and transient dynamic analysis. Currently some of these analyses are being extended to nonlinear micromechanical analyses involving viscoplasticity, elastoplasticity, combined plasticity fracture and creep. The paper describes some of the results for elastic analysis, heat transfer and thermoelastic analyses.

Banerjee, P. K.↗

Elastic Model Transitions: A Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes a hybrid approach for determining physical displacements by augmenting the original quadratically constrained least squares (LSQI) algorithm with Direct Shape Mapping (DSM) and modifying the energy constraints. The approach presented is applicable to simulation of the elastic behavior of launch vehicles and other structures that utilize discrete LTI finite element model (FEM) derived mode sets (eigenvalues and eigenvectors) that are propagated throughout time. The time invariant nature of the elastic data presents a problem of how to properly transition elastic states from the prior to the new model while preserving motion across the transition and ensuring there is no truncation or excitation of the system. A previous approach utilizes a LSQI algorithm with an energy constraint to effect smooth transitions between eigenvector sets with no requirement that the models be of similar dimension or have any correlation. This approach assumes energy is conserved across the transition, which results in significant non-physical transients due to changing quasi-steady state energy between mode sets, a phenomenon seen when utilizing a truncated mode set. The computational burden of simulating a full mode set is significant so a subset of modes is often selected to reduce run time. As a result of this truncation, energy between mode sets may not be constant and solutions across transitions could produce non-physical transients. In an effort to abate these transients an improved methodology was developed based on the aforementioned approach, but this new approach can handle significant changes in energy across mode set transitions. It is proposed that physical velocities due to elastic behavior be solved for using the LSQI algorithm, but solve for displacements using a two-step process that independently addresses the quasi-steady-state and non-steady-state contributions to the elastic displacement. For structures subject to large external forces, such as thrust or atmospheric drag, it is imperative to capture these forces when solving for elastic displacement. To simplify the mathematical formulation, assumptions are made regarding mass matrix normalization, constant external forcing, and constant viscous damping. These simplifications allow for direct solutions to the quasi-steady-state displacements through a process titled Direct Shape Mapping. DSM solves for the displacements using the eigenvalues of the elastic modes and the external forcing and returns a set of elastic displacements dictated by the eigenvectors of the post-transition mode set. For the non-steady-state contributions to displacement we formulate a LSQI problem that is constrained by energy of the non-steady state terms. The contributions from the quasi-steady-state and non-steady state solutions are then combined to obtain the physical displacements associated with the new set of eigenvectors. Results for the LSQI-DSM approach show significant reduction/complete removal of transients across mode set transitions while maintaining elastic motion from the prior state. For time propagation applications employing discrete elastic models that need to be transitioned in time and where running with full a full mode set is not feasible, the method developed offers a practical solution to simulating vehicle elasticity.

Hannan, Mike R.↗

Early stage growth of amorphous thin film: Average kinetics, nanoscale dynamics, and pressure dependence

We used the coherent grazing incidence small angle x-ray scattering technique to study the average kinetics and nanoscale dynamics during early-stage a-WSi 2 sputter deposition. The kinetic and dynamic properties are examined as a function of pressure, which is known to be a critical factor in determining final surface roughness. Surface growth kinetics and dynamics are characterized by time parameters extracted from the height-height structure factor and correlation functions. The roughness at a given length scale reaches a maximum before relaxing down to a steady state. The lateral length scale dependence and pressure dependence are then compared among the measured kinetics and dynamics time parameters. Surfaces grown at lower pressures are smoother, leading to longer correlation times. The time parameters to reach a dynamic steady state and a kinetic steady state show contrasting pressure dependence. A dynamic steady state is reached earlier than the kinetic steady state at high pressure. A more random deposition direction and lower kinetic energy at higher pressures can explain these phenomena, along with the hypothesis that larger nanoclusters form in vapor before arriving at the surface. Here, a continuum model is applied to simulate the overall behavior with mixed success.

36 MATERIALS SCIENCE↗