Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock physics”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2024

An exploration of anomalous electrical noise in shocked cyclotrimethylenetrinitramine (RDX)-based explosives

Gas gun shock experiments on cyclotrimethylenetrinitramine (RDX)-based explosive compositions that employ embedded gauge particle velocity tracers have noted a significant amount of electrical noise when compared to other explosive formulations. This paper reexamines previously published embedded gauge data on Cyclotols (60–80 wt. % RDX) to quantify the electromagnetic behavior of these materials. The primary observation is a fourfold increase in the electrical noise when Cyclotols are shocked above 4.22 ± 0.08 GPa. Electromagnetic gauge noise is also observed within particle velocity traces in reactive growth and off-Hugoniot shocks, although at higher pressures than the direct shock case, suggesting a temperature- or kinetically dependent transition. In all cases, the electrical noise disappears upon detonation. By comparing with the static high-pressure phase diagram of RDX, we interpret this change in electromagnetic behavior to be a change in the RDX crystal structure to a piezoelectric phase, although it is uncertain whether the γ or ε phase is responsible for the observed behavior.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A detailed chemical kinetic modeling and experimental investigation of the low- and high-temperature chemistry of n-butylcyclohexane

Chemical kinetic models of gasoline, jet, and diesel fuels and their mixtures with bioderived fuels are needed to assess fuel property effects on efficiency, emissions, and other performance metrics in internal combustion and gas turbine engines. As these real fuels have too many fuel components to be included in a chemical kinetic model, surrogate fuels containing fewer components are used to represent them. These surrogate fuels mimic the chemical classes or molecular structures contained in the real fuel. One of the important chemical classes in gasoline, jet, and diesel fuels comprises cyclohexanes. Cyclohexanes comprise about 30% or more by weight in diesel fuel. Also, Mueller et al proposed n-butylcyclohexane (nBCH) as a component in a nine-component surrogate palette to represent the ignition properties, distillation curve, density, and molecular structures of a diesel certification fuel. Herein, experimental measurements of the ignition delay times (IDTs) of nBCH in a shock tube and in a rapid compression machine are reported over a wide range of temperature, pressure, and equivalence ratio important for enabling the validation of a chemical kinetic model for nBCH for combustion in diesel engines. The range of conditions are temperatures of 630–1420 K, pressures of 10, 30, and 50 bar, and equivalence ratios of 0.3, 0.5, 1.0, and 2.0 in ‘air’. A detailed chemical model is developed for nBCH to simulate its ignition at both low and high-temperature conditions and at relevant elevated pressures. The experimentally measured IDTs are used to improve and validate the chemical kinetic model.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Synthesis and investigation into explosive sensitivity for a series of new picramide explosives

Tailoring the molecular properties that govern energetic material sensitivity is essential to improve safety and help develop new energetic materials. Despite this need, understanding the complex chemistry and physics of explosive initiation and propagation is still a challenge. Recent work by our group has reinforced the view that explosive sensitivity under sub-shock conditions is connected to the strength of the weakest covalent bond in the molecule, that is, its “trigger linkage.” These correlations have been observed with different classes of energetic molecules and indicate that “trigger linkage” bond breaking, and heat of explosion are good indicators for the sensitivity trends. Herein we report the synthesis of aliphatic energetic materials with ethane, propane and neopentane backbones. Experimental and computational studies show that the trigger linkage model, based on results from quantum molecular dynamics simulations, correctly predicts trends observed in the impact sensitivity of the molecules. However, while the model predicts the impact sensitivities of the ethane series, the neopentane series has higher impact sensitivities than predicted, which is presumably influenced by crystal packing effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1991

Physical processes in the plasma mantle of Venus

The results are presented of a study, which analyzed data from 10 Pioneer Venus orbits in order to see whether similar wave particle interaction processes also exist in the corresponding region around Venus. The first conclusion is that the apparent physical processes in the mantle are indeed similar around Venus and Mars. The planetary thermal O(+) ions outside the ionopause interact with the shocked solar wind and excite electrostatic waves close to the lower hybrid frequency. These waves propagate inwards, heating first the electron and deeper down in the ionosphere the thermal ion population. The observed superthermal ions are believed to be the product of this wave particle interaction process. It is also concluded that the wave energy transferred to the thermal electrons is of the right magnitude (about 4 x 10 exp 9 eV/sq cm s) to provide the supplemental heat source necessary to reconcile observed and calculated electron temperatures in the ionosphere.

Szego, K.↗
NASA NTRS1995

The Method of Space-time Conservation Element and Solution Element: Development of a New Implicit Solver

The method of space-time conservation element and solution element is a nontraditional numerical method designed from a physicist's perspective, i.e., its development is based more on physics than numerics. It uses only the simplest approximation techniques and yet is capable of generating nearly perfect solutions for a 2-D shock reflection problem used by Helen Yee and others. In addition to providing an overall view of the new method, we introduce a new concept in the design of implicit schemes, and use it to construct a highly accurate solver for a convection-diffusion equation. It is shown that, in the inviscid case, this new scheme becomes explicit and its amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, its principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme.

Chang, S. C.↗
NASA NTRS1994

Nonthermal radio emission from hot star winds: Its origin and physical implications

Nonthermal radio emission has been observed from some of the most luminous hot star winds. It is understood to be synchrotron radiation of the relativistic electrons in the winds. To understand how the electrons are accelerated to such high energies and to correctly explain the observed radio flux and spectra require an exhaustive investigation of all the relevant physical processes involved and possibly point to a complex wind structure. In this paper we discuss the logical path toward a comprehensive model of the nonthermal radio emission from hot star winds. Based on the available observational data and fundamental theoretical considerations, we found that the only physically viable and self-consistent scenario is: the nonthermal radio emission is synchrotron radiation of relativistic electrons the electrons are accelerated by shocks via the first-order Fermi mechanism the acceleration has to be in situ in the radio emitting region and the shocks formed at the base of the winds have to propagate to beyond the radio photosphere.

Chen, Wan↗
NASA NTRS2009

CFD Analysis of Nozzle Jet Plume Effects on Sonic Boom Signature

A computational fluid dynamics study is conducted to examine nozzle exhaust jet plume effects on the Sonic boom signature of a supersonic aircraft. A simplified axisymmetric nozzle geometry, representative of the nozzle on the NASA Dryden NF-15B Lift and Nozzle Change Effects on Tail Shock research airplane, is considered. The computational fluid dynamics code is validated using available wind-tunnel sonic boom experimental data. The effects of grid size, spatial order of accuracy. grid type, and flow viscosity on the accuracy of the predicted sonic boom pressure signature are quantified. Grid lines parallel to the Mach wave direction are found to give the best results. Second-order accurate upwind methods are required as a minimum for accurate sonic boom simulations. The highly underexpanded nozzle flow is found to provide significantly more reduction in the tail shock strength in the sonic boom N-wave pressure signature than perfectly expanded and overexpanded nozzle flows. A tail shock train in the sonic boom signature is observed for the highly underexpanded nozzle flow. Axisymmetric computational fluid dynamics simulations show the flow physics inside the F-15 nozzle to be nonisentropic and complex.

Bui, Trong T.↗
NASA NTRS2011

Time-Frequency Analysis of Rocket Nozzle Wall Pressures During Start-up Transients

Surveys of the fluctuating wall pressure were conducted on a sub-scale, thrust- optimized parabolic nozzle in order to develop a physical intuition for its Fourier-azimuthal mode behavior during fixed and transient start-up conditions. These unsteady signatures are driven by shock wave turbulent boundary layer interactions which depend on the nozzle pressure ratio and nozzle geometry. The focus however, is on the degree of similarity between the spectral footprints of these modes obtained from transient start-ups as opposed to a sequence of fixed nozzle pressure ratio conditions. For the latter, statistically converged spectra are computed using conventional Fourier analyses techniques, whereas the former are investigated by way of time-frequency analysis. The findings suggest that at low nozzle pressure ratios -- where the flow resides in a Free Shock Separation state -- strong spectral similarities occur between fixed and transient conditions. Conversely, at higher nozzle pressure ratios -- where the flow resides in Restricted Shock Separation -- stark differences are observed between the fixed and transient conditions and depends greatly on the ramping rate of the transient period. And so, it appears that an understanding of the dynamics during transient start-up conditions cannot be furnished by a way of fixed flow analysis.

Baars, Woutijn J.↗
NASA NTRS2011

Earth-Affecting Solar Causes Observatory (EASCO): A Potential International Living with a Star Mission from Sun-Earth L5

This paper describes the scientific rationale for an L5 mission and a partial list of key scientific instruments the mission should carry. The L5 vantage point provides an unprecedented view of the solar disturbances and their solar sources that can greatly advance the science behind space weather. A coronagraph and a heliospheric imager at L5 will be able to view CMEs broadsided, so space speed of the Earth-directed CMEs can be measured accurately and their radial structure discerned. In addition, an inner coronal imager and a magnetograph from L5 can give advance information on active regions and coronal holes that will soon rotate on to the solar disk. Radio remote sensing at low frequencies can provide information on shock-driving CMEs, the most dangerous of all CMEs. Coordinated helioseismic measurements from the Sun Earth line and L5 provide information on the physical conditions at the base of the convection zone, where solar magnetism originates. Finally, in situ measurements at L5 can provide information on the large-scale solar wind structures (corotating interaction regions (CIRs)) heading towards Earth that potentially result in adverse space weather.

Gopalswamy, N.↗
NASA NTRSDate not supplied

High-Order Finite-Difference Large-Eddy Simulations of a Supersonic Nozzle

Large-eddy simulations of the flow emanating from supersonic nozzles were carried out for the 6th Propulsion Aerodynamics Workshop to be held at the AIAA SciTech Conference. The purpose of the work is to detail the ability of a high-order of accuracy finite-difference based numerical scheme to accurately compute this flow. Results for both a perfectly expanded and over-expanded flow case are presented. Several complicating factors are discussed: 1) the limitations on the computational grid imposed by the numerical method, 2) the effect of additional dissipation necessary to resolve shock waves and 3) the effect of the nozzle boundary layer on the flow. Results indicate that the simulations replicate the basic flow physics, but there is room for improvement in accuracy. Areas for future study are outlined.

turbulence↗
NASA NTRS1996

Molecular Dynamics Calculations

The development of thermodynamics and statistical mechanics is very important in the history of physics, and it underlines the difficulty in dealing with systems involving many bodies, even if those bodies are identical. Macroscopic systems of atoms typically contain so many particles that it would be virtually impossible to follow the behavior of all of the particles involved. Therefore, the behavior of a complete system can only be described or predicted in statistical ways. Under a grant to the NASA Lewis Research Center, scientists at the Case Western Reserve University have been examining the use of modern computing techniques that may be able to investigate and find the behavior of complete systems that have a large number of particles by tracking each particle individually. This is the study of molecular dynamics. In contrast to Monte Carlo techniques, which incorporate uncertainty from the outset, molecular dynamics calculations are fully deterministic. Although it is still impossible to track, even on high-speed computers, each particle in a system of a trillion trillion particles, it has been found that such systems can be well simulated by calculating the trajectories of a few thousand particles. Modern computers and efficient computing strategies have been used to calculate the behavior of a few physical systems and are now being employed to study important problems such as supersonic flows in the laboratory and in space. In particular, an animated video (available in mpeg format--4.4 MB) was produced by Dr. M.J. Woo, now a National Research Council fellow at Lewis, and the G-VIS laboratory at Lewis. This video shows the behavior of supersonic shocks produced by pistons in enclosed cylinders by following exactly the behavior of thousands of particles. The major assumptions made were that the particles involved were hard spheres and that all collisions with the walls and with other particles were fully elastic. The animated video was voted one of two winning videos in a competition held at the meeting of the American Physical Society's Division of Fluid Dynamics, held in Atlanta, Georgia, in November 1994. Of great interest was the result that in every shock there were a few high-speed precursor particles racing ahead of the shock, carrying information about its impending arrival. Most recently, Dr. Woo has been applying molecular dynamics techniques to the problem of determining the drag produced by the space station truss structure as it flies through the thin residual atmosphere of low-Earth orbit. This problem is made difficult by the complex structure of the truss and by the extreme supersonic nature of the flow. A fully filled section of the truss has already been examined, and drag predictions have been made. Molecular dynamics techniques promise to make realistic drag calculations possible even for very complex partially filled truss segments flying at arbitrary angles.

Source record↗
NASA NTRS1986

Theoretical study of the X-ray emission from astrophysical shock waves

Theoretical X-ray emission spectra are needed to interpret the X-ray emission observed by many low and moderate resolution X-ray instruments, and to provide diagnosis of physical conditions for high resolution spectra. Over the past decade, a set of model codes which compute the X-ray and XUV emission for a wide set of physical conditions, including high or low densities, photoionized gas, and time-dependent ionization balance was developed. In the past year, the atomic rate coefficients in the code was improved. Further capabilities were added, and applied to several astrophysical problems.

Raymond, J.↗
NASA NTRS2008

A Comparison of EAST Shock-Tube Radiation Measurements with a New Air Radiation Model

This paper presents a comparison between the recent EAST shock tube radiation measurements (Grinstead et al., AIAA 2008-1244) and the HARA radiation model. The equilibrium and nonequilibrium radiation measurements are studied for conditions relevant to lunar-return shock-layers; specifically shock velocities ranging from 9 to 11 kilometers per second at initial pressures of 0.1 and 0.3 Torr. The simulated shock-tube flow is assumed one-dimensional and is calculated using the LAURA code, while a detailed nonequilibrium radiation prediction is obtained in an uncoupled manner from the HARA code. The measured and predicted intensities are separated into several spectral ranges to isolate significant spectral features, mainly strong atomic line multiplets. The equations and physical data required for the prediction of these strong atomic lines are reviewed and their uncertainties identified. The 700-1020 nm wavelength range, which accounts for roughly 30% of the radiative flux to a peak-heating lunar return shock-layer, is studied in detail and the measurements and predictions are shown to agree within 15% in equilibrium. The plus or minus 1.5% uncertainty on the measured shock velocity is shown to cause up to a plus or minus 30% difference in the predicted radiation. This band of predictions contains the measured values in almost all cases. For the highly nonequilibrium 0.1 Torr cases, the nonequilibrium radiation peaks are under-predicted by about half. This under-prediction is considered acceptable when compared to the order-of-magnitude over-prediction obtained using a Boltzmann population of electronic states. The reasonable comparison in the nonequilibrium regions provides validation for both the non-Boltzmann modeling in HARA and the thermochemical nonequilibrium modeling in LAURA. The N2 (+)(1-) and N2(2+) molecular band systems are studied in the 290 480 nm wavelength range for both equilibrium and nonequilibrium regimes. The non-Boltzmann rate models for these systems, which have significant uncertainties, are tuned to improve the comparison with measurements.

Johnston, Christopher O.↗
DOE OSTI2021

A theoretical approach for transient shock strengthening in high-energy-density laser compression experiments

In high-energy-density shock compression experiments, the desired state of compression is typically achieved by passing shock waves through a sequence of different materials. Here, a theoretical approach for transiently strengthening such shocks passing through interfaces of experimentally relevant impedance ratios is examined. A semi-analytical method based on characteristics analysis is used to solve the problem of a shock passing from one material to another through an intermediate region of non-uniform impedance; by appropriately designing this region, a greater shock strength can be achieved in the second material for a finite duration than in the absence of this region. When a shock passes into a material of higher impedance, the shock strengthening increases to a maximum before decreasing to an asymptotic value as the strength of the incident shock is increased. For shocks passing into materials of decreasing impedance, the shock strengthening increases monotonically as the strength of the incident shock is increased when the impedance ratio is above a critical threshold. Incorporating multiple intermediate materials can further increase the strength of the transmitted shock, with an exponential discretization of intermediate material impedances being the most effective distribution for strengthening strong shocks. The results suggest that up to 25% and 9% increases in pressure behind the leading shock can be achieved for materials of increasing and decreasing impedance, respectively. The technique is applied to the design of laser-driven dynamic compression experiments, and the results of the analysis are verified via comparison to simulations performed with the HYADES hydrodynamics code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2017

The Effect of Shock on the Amorphous Component in Altered Basalt

Investigation of the geochemical and mineralogical composition of the Martian surface provides insight into the geologic history of the predominantly basaltic crust. The Chemistry and Mineralogy (CheMin) instrument onboard the Curiosity rover has returned the first X-Ray diffraction data from the Martian surface. However, large proportions (27 +/- 14 with some estimates as high as 50 weight percentage) of an amorphous component have been reported. As a remedy to this problem, mass balance equations using geochemistry, volatile chemistry, and mineralogy have been employed to constrain the geochemistry of the amorphous component. However, "the nature and number of amorphous phases that constitute the amorphous component is not unequivocally known". Multiple hypotheses have been proposed to explain the origin of this amorphous component: Allophane (Al2O); Basaltic glass (Volcanic and impact); Palagonite (Altered basaltic glass); Hisingerite (Fe (sup 3 plus)-bearing phyllosilicate); S/Cl-rich component (sulfates and/or akaganeite); Nanophase ferric oxide component (npOx). Establishing a multi-phase amorphous component from a basaltic precursor that has undergone physical and chemical weathering within geochemical constraints is of paramount importance to better understand the composition of a large portion of the Martian surface (up to 50 weight percentage). Shocked basalts from Lonar Crater in India are valuable analogs for the Martian surface because it is a well-preserved impact crater in a basaltic target. Having undergone pre- and post-shock aqueous alteration, these rocks provide crucial data regarding the effect of shock on the amorphous component in altered basalt. By conducting mass balance equations similar to what has been performed for Gale crater materials, we attempt to calculate the geochemistry of the amorphous component in altered basalts ranging from unshocked to Class 5 (Table 1). This has the potential to reveal the nature and origin (i.e. primary igneous, shock metamorphic, and/or aqueous alteration occurring before or after the impact event) of the amorphous component in shocked basalt with the goal of unravelling the history of the Martian surface.

Eckley, S. A.↗
DOE OSTI2024

Shock wave formation in radiative plasmas

The temporal evolution of weak shocks in radiative media is theoretically investigated in this work. The structure of radiative shocks has traditionally been studied in a stationary framework. Their systematic classification is complex because layers of optically thick and thin regions alternate to form a radiatively-driven precursor and a temperature-relaxation layer, between which the hydrodynamic shock is embedded. In this work, we analyze the formation of weak shocks when two radiative plasmas with different pressures are put in contact. Applying a reductive perturbative method yields a Burgers-type equation that governs the temporal evolution of the perturbed variables including the radiation field. The conditions upon which optically thick and thin solutions exist have been derived and expressed as a function of the shock strength and Boltzmann number. Below a certain Boltzmann number threshold, weak shocks always become optically thick asymptotically in time, while thin solutions appear as transitory structures. The existence of an optically thin regime is related to the presence of an overdense layer in the compressed material. Scaling laws for the characteristic formation time and shock width are provided for each regime. The theoretical analysis is supported by FLASH simulations, and a comprehensive testcase has been designed to benchmark radiative hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2017

Modeling Solar Energetic Particle Events Using ENLIL Heliosphere Simulations

Solar energetic particle (SEP) event modeling (SEPMOD) has gained renewed attention in part because of the availability of a decade of multipoint measurements from STEREO (Solar TErrestrial RElations Observatory) and L1 (Lagrangian point 1) spacecraft at 1 AU (Astronomical Unit). These observations are coupled with improving simulations of the geometry and strength of heliospheric shocks obtained by using coronagraph images to send erupted material into realistic solar wind backgrounds. The STEREO and ACE (Aerosol, Cloud systems, ocean Ecosystems) measurements in particular have highlighted the sometimes surprisingly widespread nature of SEP events. It is thus an opportune time for testing SEP models, which typically focus on protons approximately 1-100 megaelectronvolts, toward both physical insight to these observations and potentially useful space radiation environment forecasting tools. Some approaches emphasize the concept of particle acceleration and propagation from close to the Sun, while others emphasize the local field line connection to a traveling, evolving shock source. Among the latter is the previously introduced SEPMOD treatment, based on the widely accessible and well-exercised WSA-ENLIL (Wang-Sheeley-Arge-ENLIL)-cone model. SEPMOD produces SEP proton time profiles at any location within the ENLIL domain. Here we demonstrate a SEPMOD version that accommodates multiple, concurrent shock sources occurring over periods of several weeks. The results illustrate the importance of considering longer-duration time periods and multiple CME (Coronal Mass Ejection) contributions in analyzing, modeling, and forecasting SEP events.

Luhmann, J. G.↗