Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Material Response 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 1,459 records · Page 81

New insights on carbon black suspension rheology—Anisotropic thixotropy and antithixotropy

We report a detailed experimental study of peculiar thixotropic dynamics of carbon black (CB, Vulcan XC-72) suspensions in mineral oil, specifically the observation of sequential stress increase then decrease at a fixed shear rate in a step-down test. We verify that such dynamics, though peculiar, come from a true material response rather than experimental artifacts. We also reveal how this long-time stress decay is associated with antithixotropic structural change rather than viscoelastic stress relaxation by using orthogonal superposition (OSP) rheometry to probe viscoelastic moduli during the step-down tests. The orthogonal storage and loss moduli are present, showing this two-timescale recovery then decay response, which demonstrates that this response is antithixotropic, and it involves shear-induced structuring. We further show a mechanical anisotropy in the CB suspension under shear using OSP. Based on the rheological results, a microstructural schematic is proposed, considering qualitatively thixotropic structure build-up, antithixotropic densification, and anisotropic structure evolution. Our observation for these CB suspensions is outside the standard paradigm of thixotropic structure-parameter models, and the elastic response provides us with new insight into the transient dynamics of CB suspensions.

Mechanics↗

A Bezier Curve Informed Melt Pool Geometry to Model Additive Manufacturing Microstructures Using SPPARKS

Additive manufacturing is a transformative technology with the potential to manufacture designs which traditional subtractive machining methods cannot. Additive manufacturing offers fast builds at near final desired geometry; however, material properties and variability from part to part remain a challenge for certification and qualification of metallic components. AM induced metallic microstructures are spatially heterogeneous and highly process dependent. Engineering properties such as strength and toughness are significantly affected by microstructure morphologies resulting from the manufacturing process Linking process parameters to microstructures and ultimately to the dynamic response of AM materials is critical to certifying and qualifying AM built parts and components and improving the performance of AM materials. The AM fabrication process is characterized by building parts layer by layer using a selective laser melt process guided by a computer. A laser selectively scans and melts metal according to a designated geometry. As the laser scans, metal melts, fuses, and solidifies forming the final geometry in a layerwise fashion. As the laser heat source moves away, the metal cools and solidifies forming metallic microstructures. This work describes a microstructure modeling application implemented in the SPPARKS kinetic Monte Carlo computational framework for simulating the resulting microstructures. The application uses Bzier curves and surfaces to model the melt pool surface and spatial temperature profile induced by moving the laser heat source; it simulates the melting and fusing of metal at the laser hot spot and microstructure formation and evolution when the laser moves away. The geometry of the melt pool is quite flexible and we explore effects of variances in model parameters on simulated microstructures.

36 MATERIALS SCIENCE↗

Electrochemical Modeling of PID Leakage Current of PV Modules: Steady-State Current, Transient Current, and RC-Equivalent Circuit

Potential-induced degradation (PID) remains a significant reliability concern for photovoltaic (PV) modules, arising when a voltage difference between the module frame and the solar cells drives unintended leakage current through the glass-encapsulant stack. Although PID ultimately manifests as PID-s, PID-p, or PID-c, the underlying behavior of the leakage current-its magnitude and time dependence-requires clearer electrochemical interpretation. Traditional explanations attribute the initial transient current to bulk capacitive elements of the glass, encapsulant, and antireflection coatings, and the steady-state current according to their effective ohmic resistance. More recent studies, however, indicate that electrochemical charge-transfer processes at the encapsulant-metallization interface can play a dominant role in defining the leakage-current path. This paper develops a unified electrochemical framework for modeling PID leakage current. First, an RC-equivalent circuit is formulated by combining conventional RC elements with a Randles-type interface to capture transient leakage current through double-layer capacitance and faradaic processes at ionic-electronic boundaries. Second, the steady-state current-voltage behavior is explained using a linearized Butler-Volmer relationship, showing that the measured ohmic response corresponds to the low-overpotential limit of charge-transfer kinetics. Analytical results demonstrate that, for typical module materials-3.2-mm soda-lime glass and 0.45-mm encapsulant-the dominant modulators to PID leakage current are the glass surface resistance (under dry-surface conditions), the glass bulk capacitance, and the encapsulant resistance (under wet-surface conditions), with soda lime glass surface and EVA/POE encapsulant resistances primarily governing steady-state current. The proposed electrochemical model is validated against measured leakage-current data, showing good agreement in both the magnitude and the time-dependent evolution of PID leakage current.

14 SOLAR ENERGY↗

Response of Graphite to Dynamic Loading and Hypervelocity Jet Impacts

The compressive strengths of three varieties of high purity graphite, PCEA, NBG-18, and NBG-25, as well as the depth of penetration of small-scale charges into these materials was experimentally determined. These grades are similar in density, ranging from 1.80 – 1.85 g/cc, and nominal apparent porosity, ranging from 18% to 20%, but provide a wide range in maximum grain or particle size from 10s to 1000s of µm. Two very different manufacturing methods are also represented; PCEA is extruded while NBG-18 and NBG-25 are iso-molded. The quasistatic and dynamic strengths of each grade were determined on a load frame and split-Hopkinson pressure bar, respectively. The depth of penetration (DOP) of two small-scale shaped charges, the Teledyne RP-1 and RP-4, was determined against graphite. The global response of the RP-4 impacts was markedly different as the PCEA samples remained intact while all the NBG-25 samples split into 2 or 3 pieces after the jet penetration had completed. However, for all tests, the trusted DOPs fell within 2 cm. Preliminary hydrocode modeling of the penetration events used existing models that were not designed for graphite. The results can be tuned to reasonably reproduce the DOP, but the wound channel geometry is not reproduced well. A model designed for graphite would need to represent graphite’s non-linear and energy dissipation characteristics.

36 MATERIALS SCIENCE↗

Fostering a Guiding Multiscale Model for the Development of Advanced MgB 2 Hydrogen Storage Materials (Final Technical Report)

Project Goal and Objective. The demand for energy and for an upgraded energy infrastructure has steadily grown, as have the needs for energy independence and alternatives to our reliance on petroleum. Hydrogen is considered the most viable fuels for wide-scale implementation in the near future as it is less-polluting, non-toxic, and has more stored energy than petroleum. It is envisioned that hydrogen can eventually become the prime energy carrier, integrating the transportation, grid, and chemical sectors in a way that improves resiliency, diversifies feedstocks, and affords new economic opportunities. A key remaining challenge is the development materials with enhanced gravimetric and volumetric hydrogen storage capacities that offer a higher performance than compressed gas. These materials would eliminate the need for large-scale compression, thereby dramatically reducing the footprint and cost of gas storage. The high gravimetric and volumetric hydrogen capacities of complex hydrides has prompted an intensive investigation of the potential of this class of materials as hydrogen storage media over the past 25 years. Among the many complex hydrides that have been explored, magnesium borohydride, Mg(BH 4 ) 2 , has been found to possess the best combination of practical thermodynamic properties. These include a gravimetric H 2 density of 14.9 wt% H 2 and thermodynamics for the dehydrogenation of Mg(BH 4 ) 2 to MgB 2 (equation 1) (ΔH° = 39 kJ/mol H 2 , ΔS = 112 J/K mol H 2 ) which lie in the narrow window required Mg(BH 4 ) 2 $\Leftrightarrow$ MgB 2 + 4 H 2 (1) for reversibility under moderate pressure and temperature. However, overcoming the extremely slow kinetics of the reversible release of hydrogen by this material in the solid state is a daunting challenge. At temperatures greater than 400 °C, the borohydride releases up to 14 wt% hydrogen giving MgB 2 . We discovered that the direct re-hydrogenation of MgB 2 to Mg(BH 4 ) 2 can be accomplished under 950 bar H 2 at 400 °C. While this demonstrated that complete reversibility can be achieved, the conditions employed are far too extreme for commercial hydrogen storage applications. More recently, we found through US DOE funded research projects (EERE HyMARC and HySCOR), that hydrogen cycling, can be accomplish at much milder conditions upon modification of the borohydride or boride. Guided by these discoveries these discoveries, the objective of this research project was to obtain key information that will enable the development of a model of reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 . The ultimate goal of our efforts is to attain a model of this transformation that can be utilized to accelerate development further advanced materials. This project directly follows on discoveries that were made over the course of a US DOE, EERE HyMARC project that was focused on improvement of the hydrogen cycling kinetics of modified MgB 2 . We found that that mechanical milling with graphene results the desired, pronounced kinetic enhancement. The dramatic lowering of the conditions required for the hydrogenation of MgB 2 is a significant step towards overcoming its chemical inertness allowing its development as a practical onboard hydrogen storage material. However, the exact nature of the modification(s) of MgB 2 that is responsible for its activation towards hydrogenation is completely unknown. This situation is not unique, as efforts to develop hydrogen storage materials typically have a narrow focus rather than a comprehensive approach that takes atomic level bonding and structure; molecular dynamics; long range, nano- and mesoscale-structure and their interconnection all into account. The goal of this project was the development of a comprehensive, multi-scale computational model of reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 that can be utilized for development of higher performance versions of the modified material. Development of the model requires determination of: 1) the bulk, nano-scale, and meso-scale structural changes occurring at elevated pressure following mechano-chemical modification of MgB 2 ; 2) the reaction pathway of the reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 ; 3) the effect of elevated pressure and mechano-chemical modification on the chemical reaction pathways; 4) the interactions at solid-gas interfaces; and particle surfaces; and 5) the kinetics and thermodynamic parameters associated with each step of the hydrogenation reaction pathway. This investigation required advanced techniques as preliminary, standard XRD, 11 B NMR, and FTIR analysis showed no signs of material modification. In order to gain this level of understanding of modified MgB 2 , required the teaming of a diverse group of experts and state-of-the art experimental capabilities at the University of Hawaii at Manoa (UHM) and collaborating National Laboratories: Craig Jensen , Department of Chemistry (PI and Project Director), solid state, solution, and high pressure NMR spectroscopy; solid-state synthesis; and high pressure hydrogenation (collaboration with SNL); Godwin Severa , Hawaii Natural Energy Institute (co-PI) calorimetry; infrared and Raman spectroscopy (collaboration with NREL); Dera , high pressure X-ray diffraction including in situ experiments (collaboration with ANL); Hope Ishii , Hawaii Institute of Geophysics electron microscopy investigations (collaboration with LBNL); and Joe Brown , Mechanical Engineering , material electronic structure and electric field effects.

08 HYDROGEN↗

Continued Verification of MOOSE Structural Mechanics Tools for Modeling Core Bowing Phenomena in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY22, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21. Several benchmark verification examples were selected from the IAEA verification and validation report. These examples involve clusters of ducts representative of a sector of a hexagonal reactor core which bow into each other and cause contact and load pad elevations, as well as single ducts subjected to irradiation fields undergoing swelling and subsequent bowing. The MOOSE-based results were compared to both IAEA benchmark participants’ results, analytic equations as available, and NUBOW-3D, a beam model code developed by Argonne National Laboratory. In every case, the MOOSE results agreed with other simulations results, providing additional verification basis of the tools for this particular physics application.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

FY23 Status Report on MOOSE-Based Approaches to Modeling Core Bowing in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY23, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21 and FY22. Several benchmark verification examples were selected from the IAEA verification and validation report with increased number of ducts and more complex contact interaction behavior. These examples involve a full symmetric sector with restraint rings at multiple load pad locations to simulate a limited free-bow restraint system concept, as well as irradiation induced swelling and creep effects in a sector. In addition, improvements to the contact module sideset assignment were assessed and compared with previous MOOSE results to verify the contact behavior. The MOOSE-based results were compared to IAEA benchmark participants’ results. In every case, the MOOSE results agreed with other simulations results for estimating bowing behavior, providing additional verification basis of the tools for this particular physics application. Estimation of contact forces were mostly in agreement, with a few outlier results. A plan was suggested for dealing with the discrepancies with estimating contact force values. In addition, a 1-way coupling demonstration was performed using subchannel analysis code Pronghorn-SC and MOOSE on an ABR-1000-design sodium-cooled fast reactor assembly to evaluate coolant and duct wall temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Disruption modelling for engineering and physics design of ST-E1 fusion power plant

Plasma disruptions represent a critical challenge for high-performance tokamak operations, as they can compromise machine integrity and reduce operational availability. Although future fusion devices essentially need to incorporate strategies to minimise disruption occurrence, complete avoidance remains unattainable. Consequently, assessing and characterising unmitigated disruption consequences is fundamental for the design and qualification of next-generation fusion power plants. This work supports the pre-conceptual design of ST-E1, a low aspect-ratio Tokamak Fusion Power Plant developed by Tokamak Energy Ltd., by presenting a comprehensive disruption modelling approach applied across different design stages. The methodology integrates both physics and engineering considerations to evaluate the impact of disruptions on machine performance and structural integrity. From an engineering perspective, several ST-E1 layout options were analysed to investigate the electromagnetic response of key components under disruption-induced loads, enabling comparison between alternative design solutions. On the physics side, a broad set of disruption scenarios was explored, scanning operational space parameters, plasma-material interactions, and associated thermal loads. Furthermore, the study examined variations in disruption behaviour arising from different reference equilibria, focusing on a range starting from Double Null to Single Null configurations, reflecting the increasing up-down asymmetry consequences. The results reveal significant contrasts in plasma dynamics and structures electromagnetic behaviour between configurations, highlighting the importance of disruption modelling in guiding design choices. These analyses have proven instrumental in shaping ST-E1 development, offering critical insights for mitigating risks and optimising future fusion power plant designs.

Borowiec, Katarzyna [ORNL] (ORCID:0000000335911739↗

New self-consistent homogenization for thermo-elastic polycrystals with imperfect interfaces

Here we present a new thermo-elastic self-consistent formulation to calculate the thermo-mechanical response of polycrystalline aggregates with imperfect interfaces. The new model is based on an improved treatment of the effect of imperfect interfaces on the micromechanical response of the grains compared to pre-existing models. New expressions/algorithms are derived for the effective elastic and thermal properties, localization relations, and calculation of intragranular field fluctuations. In addition, derivation of a novel integral equation for the strain field in a heterogeneous elastic medium subjected to arbitrary thermal strain field and displacement jumps is presented. The new model is applied to study the effective thermo-elastic properties and average field fluctuations as a function of the interface properties for imperfectly-bonded triaminotrinitrobenzene (TATB) polycrystalline granular aggregates. The proposed formulation is general in the sense that it can be extended to other material systems and deformation regimes that require a proper consideration of the effect of imperfect interface properties on the material's thermo-mechanical response.

36 MATERIALS SCIENCE↗

Engineered Polymer Architectures for Thermo–Responsive Desiccants in Dehumidification Applications

Humidity control is relevant to many processes in buildings and manufacturing and is responsible for over 600 million tons of CO 2 annually. As such, desiccants that can efficiently remove moisture from air can greatly reduce emissions. A new class of desiccant materials, thermo–responsive desiccants, has been of increased interest in the literature. These materials are generally based on thermo–responsive polymers and exhibit a temperature–dependent isotherm, an attribute that is small in traditional desiccants. In this work, seven key parameters needed to create an effective thermo–responsive polymer desiccant are identified using a combination of modeling and experimental data. Further, it is found that more thermo–responsiveness is not necessarily better for energy efficiency, and that an optimal, “medium” value is desired. Additionally, polymer composition and architecture play key roles in adjusting the seven key parameters. Finally, based on findings from both modeling and experimental work, guidelines on synthesizing future thermo–responsive polymer desiccants are presented.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Magnesium abundances in cool metal-polluted white dwarfs

ABSTRACT The accretion of rocky material is responsible for the presence of heavy elements in the atmospheres of a large fraction of white dwarf stars. Those objects represent a unique opportunity to infer the bulk composition of exoplanetesimals. This chemical characterization requires the use of detailed atmosphere models to determine the elemental abundances at the photospheres of white dwarfs. In this work, we use a state-of-the-art model atmosphere code to reanalyse the first large survey of metal-polluted white dwarfs for which abundances are found for multiple elements. We show that the improved constitutive physics of our models lead to systematically higher Mg abundances than previous analyses. We find an average $\log \, {\rm Mg/Ca}$ number abundance ratio of 1.5. This value is significantly above the reference abundance for chondrites, which is expected as current diffusion models predict that for the cool helium-atmosphere white dwarfs of our sample, Mg should remain in the atmosphere longer than Ca. This helps resolve a recently identified Mg depletion problem, where the planetesimals accreted by white dwarfs were reported to be Mg-deficient compared to the expected composition of their planetary systems.

79 ASTRONOMY AND ASTROPHYSICS↗

Group Conductivity and Nonadiabatic Born Effective Charges of Disordered Metals, Warm Dense Matter, and Hot Dense Plasma

The average ionization state is a critical parameter in plasma models for charged particle transport, equations of state, and optical response. The dynamical or nonadiabatic Born effective charge (NBEC), calculated via first principles time-dependent density-functional theory, provides exact ionic partitioning of bulk electron response for both metallic and insulating materials. The NBEC can be transformed into a “group conductivity,” i.e., the electron conductivity ascribed to a subset of ions. We show that for disordered metallic systems, such as warm dense matter (WDM) and hot dense plasma, the static limit of the NBEC is different from the average ionization states, but that the ionization state can be extracted from the group conductivity even in mixed systems. Here we demonstrate this approach using a set of archetypical examples, including cold and warm aluminium, low- and high-density WDM carbon, and a WDM carbon-beryllium-hydrogen mixture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design of a Single Layer Metamaterial for Pressure Offloading of Transtibial Amputees

While using a prosthesis, transtibial amputees may experience pain and discomfort brought on by large pressure gradients at the interface between the residual limb and prosthetic socket. Current prosthetic interface solutions attempt to alleviate these pressure gradients by using soft homogenous liners to reduce and distribute pressures. This research investigates an additively manufactured metamaterial inlay with a tailored mechanical response to reduce peak pressure gradients around the limb. The inlay uses a hyperelastic behaving metamaterial (US10244818) comprised of triangular pattern unit cells which can be 3D printed with walls of various thicknesses controlled by draft angles. The hyperelastic material properties are modeled using a Yeoh 3rd order model. The 3rd order coefficients can be adjusted and optimized, which corresponds to a change in the unit cell wall thickness to create an inlay that can meet the unique offloading needs of an amputee. Finite element analysis simulations evaluated the pressure gradient reduction from: 1) a common homogenous silicone liner, 2) a prosthetist's inlay prescription that utilizes three variations of the metamaterial, and 3) a metamaterial solution with optimized Yeoh 3rd order coefficients. When compared to a traditional homogenous silicone liner for two unique limb loading scenarios, the prosthetist prescribed inlay and the optimized material inlay can achieve equal or greater pressure gradient reduction capabilities. These results show the potential feasibility of implementing this metamaterial as a method of personalized medicine for transtibial amputees by creating a customizable interface solution to meet the unique performance needs of an individual patient.

60 APPLIED LIFE SCIENCES↗

Detailed Design and Cost Estimation of a 300 MWe Oxy-Fuel sCO2 Turbine

The detailed design of a 300 MWe, utility scale oxy-fuel turbine has been completed for purposed operation in the sCO2 direct fired Allam-Fetvedt cycle, targeting near-zero emissions and a 50% LHV system efficiency. The turbine and its supporting plant aim to offer a lower levelized cost of energy than a natural gas combined cycle plant employing carbon capture. The oxy-fuel turbine conditions include an inlet temperature of 1150°C and inlet pressure of 305 bar, representing temperatures near that of a gas turbine simultaneously with pressures near an ultra-supercritical steam turbine. The combustor housing and turbine designs were completed according to the ASME BPVC; the turbine case specifically incorporates a multi-body design with inner high-pressure barrel case and low-pressure (30 bar) horizontally split outer case of low-chromium steel material. Lateral rotordynamic evaluation demonstrated acceptable vibration response for a range of imbalance conditions per API standards. The cooling flow required in the six-stage turbine flowpath for 30,000 hr. blade and stator lifetime is predicted through thermal and structural modeling of the first stage. The provided cost estimate of the turbine is formed through a combination of scaled up-costs from procured 10 MWe scale sCO2 turbomachinery hardware, and vendor provided budgetary quotes of larger components including the turbine case requiring casting, welding, and final machining processes. The performance and cost estimation of the oxy-fuel turbine predicted for the completed detailed design provides important information towards future development needs for market penetration of utility scale direct fired sCO2 power cycles.

Marshall, Michael [Southwest Research Institute, S↗

Analysis methods for Kevlar shield response to rotor fragments

Several empirical and analytical approaches to rotor burst shield sizing are compared and principal differences in metal and fabric dynamic behavior are discussed. The application of transient structural response computer programs to predict Kevlar containment limits is described. For preliminary shield sizing, present analytical methods are useful if insufficient test data for empirical modeling are available. To provide other information useful for engineering design, analytical methods require further developments in material characterization, failure criteria, loads definition, and post-impact fragment trajectory prediction.

Gerstle, J. H.↗

Magnetic bearings for a spaceflight optical disk recorder

The development and testing of a magnetic bearing system for the translator of the read/write head in a magneto-optic disk drive are discussed. The asymmetrical three-pole actuators with permanent magnet bias support the optical head, and its tracking and focusing servos, through their radial excursion above the disk. The specifications for the magnetic bearing are presented, along with the configuration of the magnetic hardware. Development of a five degree of freedom collision model is examined which allowed assessment of the system response during large scale transients. Experimental findings and the results of performance testing are presented, including the roll-off of current-to-force due to eddy current loss in the magnetic materials.

Hockney, Richard↗

Sensitivity analysis for the dynamic response of thermoviscoplastic shells of revolution

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric, fully-coupled, thermoviscoplastic response of shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of large-strain, transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be viscoplasticity with strain hardening and thermal hardening, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used for the shell equations with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. The energy-balance equation is solved using a Galerkin procedure, with the temperature as the fundamental unknown. Spatial discretization is performed in one dimension (meridional direction) for the momentum and constitutive equations of the shell, and in two dimensions (meridional and thickness directions) for the energy-balance equation. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) for the momentum equation; a predictor-corrector version of the trapezoidal rule is used for the energy-balance equation; and an explicit scheme consistent with the central difference method is used to integrate the constitutive equations. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading.

Kulkarni, Makarand↗

Global Aspects of Heliosphere-Geosphere Coupling

The magnetosphere serves both to hold off the solar wind and to couple it selectively to the Earth through the auroral zones. It has long been understood that the plasmasphere consists of geogenic plasmas that expand out of the sunlit low latitude ionosphere, but it was initially assumed that the hot low density plasma beyond the plasmasphere is largely of solar wind origin with a minor admixture of ionospheric plasmas, discovered via mass spectrometric observations in the early 70's. Since then, elaborate simulation models have been developed that have taken us from an era of cartoon physics to a new era of quantitative global comparisons between observations and theory. In most current global circulation models of the magnetosphere, the ionospheric load on the system is taken to lie exclusively in the thin F layer of the ionosphere. This layer is coupled with solar wind and magnetospheric plasmas via Maxwell stresses communicated by field aligned current systems, and with the thermosphere via ion-neutral charge exchange and Coulomb collisions. However, recent observations have shown us that ionospheric plasmas flow sporadically in various forms into the plasmasphere and up into the high latitude circulation cells of the outer magnetosphere, as driven by the solar wind and its variable intensity and magnetic field. Under some conditions, ionospheric material is observed to be the principal component of plasmas at the dayside magnetopause and in the plasma sheet and ring current regions. Given a global model of magnetospheric circulation, it is relatively straightforward to investigate the behaviors of ionospheric plasmas in response to solar wind drivers, and we report the results of such efforts here. We find that ionospheric plasmas dominate the plasma pressure in the magnetosphere in some regions, particularly when the solar wind is especially intense and-or southward directed. This result violates the assumption that the ionospheric load is confined to the F layer, and shows that the ionosphere is often an important dynamic element of the solar wind-ionosphere system throughout the magnetosphere, especially during larger geospace storms. This means that future global circulation models must account for ionospheric plasma inertia, heat capacity, and pressure, to be quantitatively credible.

Moore, Thomas E.↗