Search NASA⌕ Search

SEARCH · Search NASA

Results for “modeling and filtering”

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

Digital adaptive controllers using second order models with transport lag

Design of a discrete optimal regulator requires the a priori knowledge of a mathematical model for the system of interest. Because a second-order model with transport lag is very amenable to control computations and because this type of model has been used previously to represent certain high order single input-single output processes, an adaptive controller was designed based upon adjustment of controls computed for such a model. An extended Kalman filter was utilized for tracking the model parameters which were subsequently used to update a set of optimal control gains. Favorable results were obtained in applying this procedure to the control of several examples including a ninth order nonlinear process.

Joshi, S.↗
NASA NTRS1993

Analyzing Martian winds and tracer concentrations using Mars Observer data

During the courses of a day, the Mars Observer spacecraft will acquire globally distributed profiles of the martian atmosphere. It is highly desirable that this data be assembled into synoptic weather maps (complete specifications of the atmospheric pressure, temperature, and winds at a given time), which can in turn be used as starting points in the study of many meteorological phenomena. Unfortunately, the special nature of the Mars Observer data presents several challenges above and beyond the usual difficult problem of data initialization. Mars Observer atmospheric data will consist almost exclusively of asynoptic vertical profiles of temperatures (or radiances) and pressures, whereas winds are generally in balance with horizontal gradients of these quantities (which will not be observed). It will therefore be necessary to resort to dynamical models to analyze the wind fields. As a rule, data assimilation into atmospheric models can result in the generation of spurious gravity waves, so special steps must be taken to suppress these. In addition, the asynoptic nature of the data will require a four-dimensional (space and time) data assimilation scheme. The problem is to find a full set of meteorological fields (winds and temperatures) such that, when marched forward in time in the model, they achieve a best fit (in the weighted least-squares sense) to the data. The proposed solution is to develop a model especially for the Mars Observer data assimilation problem. Gravity waves are filtered from the model by eliminating all divergence terms from the prognostic divergence equation. This leaves a diagnostic gradient wind relation between the rotational wind and the temperature field. The divergent wind is diagnosed as the wind required to maintain the gradient wind balance in the presence of the diabatic heating. The primitive equations of atmospheric dynamics (with three principal dependent variables) are thus reduced to a simpler system with a single prognostic equation for temperature - the variable that will be best observed. (This balance system was apparently first derived by Charney as a first-order Rossby number expansion of the equations of motion). Experience with a full primitive equation model of the martian atmosphere indicates that a further simplification is possible: at least for short-term integrations, the model can be linearized about the zonally symmetric basic state.

Houben, Howard C.↗
NASA NTRS2011

Lessons Learned from Assimilating Altimeter Data into a Coupled General Circulation Model with the GMAO Augmented Ensemble Kalman Filter

Satellite altimetry measurements have provided global, evenly distributed observations of the ocean surface since 1993. However, the difficulties introduced by the presence of model biases and the requirement that data assimilation systems extrapolate the sea surface height (SSH) information to the subsurface in order to estimate the temperature, salinity and currents make it difficult to optimally exploit these measurements. This talk investigates the potential of the altimetry data assimilation once the biases are accounted for with an ad hoc bias estimation scheme. Either steady-state or state-dependent multivariate background-error covariances from an ensemble of model integrations are used to address the problem of extrapolating the information to the sub-surface. The GMAO ocean data assimilation system applied to an ensemble of coupled model instances using the GEOS-5 AGCM coupled to MOM4 is used in the investigation. To model the background error covariances, the system relies on a hybrid ensemble approach in which a small number of dynamically evolved model trajectories is augmented on the one hand with past instances of the state vector along each trajectory and, on the other, with a steady state ensemble of error estimates from a time series of short-term model forecasts. A state-dependent adaptive error-covariance localization and inflation algorithm controls how the SSH information is extrapolated to the sub-surface. A two-step predictor corrector approach is used to assimilate future information. Independent (not-assimilated) temperature and salinity observations from Argo floats are used to validate the assimilation. A two-step projection method in which the system first calculates a SSH increment and then projects this increment vertically onto the temperature, salt and current fields is found to be most effective in reconstructing the sub-surface information. The performance of the system in reconstructing the sub-surface fields is particularly impressive for temperature, but not as satisfactory for salt.

Keppenne, Christian↗
NASA NTRS2011

Numerical Simulations of Two-Phase Reacting Flow in a Single-Element Lean Direct Injection (LDI) Combustor Using NCC

A series of numerical simulations of Jet-A spray reacting flow in a single-element lean direct injection (LDI) combustor have been conducted by using the National Combustion Code (NCC). The simulations have been carried out using the time filtered Navier-Stokes (TFNS) approach ranging from the steady Reynolds-averaged Navier-Stokes (RANS), unsteady RANS (URANS), to the dynamic flow structure simulation (DFS). The sub-grid model employed for turbulent mixing and combustion includes the well-mixed model, the linear eddy mixing (LEM) model, and the filtered mass density function (FDF/PDF) model. The starting condition of the injected liquid spray is specified via empirical droplet size correlation, and a five-species single-step global reduced mechanism is employed for fuel chemistry. All the calculations use the same grid whose resolution is of the RANS type. Comparisons of results from various models are presented.

Liu, Nan-Suey↗
NASA NTRS1991

A fast algorithm for control and estimation using a polynomial state-space structure

One of the major problems associated with the control of flexible structures is the estimation of system states. Since the parameters of the structures are not constant under varying loads and conditions, conventional fixed parameter state estimators can not be used to effectively estimate the states of the system. One alternative is to use a state estimator which adapts to the condition of the system. One such estimator is the Kalman filter. This filter is a time varying recursive digital filter which is based upon a model of the system being measured. This filter adapts the model according to the output of the system. Previously, the Kalman filter has only been used in an off-line capacity due to the computational time required for implementation. With recent advances in computer technology, it is becoming a viable tool for use in the on-line environment. A distributed Kalman filter implementation is described for fast estimation of the state of a flexible arm. A key issue, is the sensor structure and initial work on a distributed sensor that could be used with the Kalman filter is presented.

Shults, James R.↗
DOE OSTI2021

Model predictive combustion control of a Gasoline Compression Ignition engine

Gasoline Compression Ignition is a novel combustion concept that derives its superiority from the high compression ratio of a compression ignition engine as well as the properties of gasoline fuel, such as longer ignition delay and higher volatility compared to diesel. Here, this combustion concept was experimentally tested on a 12.4L Class 8 truck engine. Based on these experimental data, prior efforts by the authors focused on the development of an engine model for a heavy-duty engine operating on a low-reactivity fuel. This engine model was leveraged within this study to investigate a combustion control strategy at different engine conditions and injection methods and was augmented to incorporate cycle-to-cycle combustion variations. State estimation is performed by means of a Kalman filter which feeds into a model predictive controller. The model predictive controller chooses control actions based on a predefined cost function under consideration of bounds reflecting physical constraints. The engine model was utilized to establish a state-space model that serves the Kalman filter and model predictive controller for estimation and prediction. A comparative study investigating control actions and engine behavior was performed with and without limiting in-cylinder peak pressure as well as combustion noise, which is of particular interest for early pilot injection strategies. In addition, the proposed control architecture was investigated at two different levels of cycle-to-cycle variations and compared to the performance of a control structure with input disturbance rejection. For increased cycle-to-cycle variations, disturbance estimation reduces state fluctuations and control effort. In general, this investigation highlights control aspects specific to a compression-ignited combustion regime with low-reactivity fuel. The control algorithm is able to maintain the desired references for brake mean effective pressure and combustion phasing while controlling peak in-cylinder pressure and combustion noise.

42 ENGINEERING↗
NASA NTRS1993

Simulated orbiter lift-off acceleration - A stochastic model

A simple procedure for generating synthetic orbiter liftoff acceleration is described. This stochastic method is based on a filtered white noise model which accounts for the nonstationary time evolution of amplitude and frequency content of orbiter lift-off accelerations. The example of Space Transportation System STS-41 Z-acceleration is carefully examined, and the corresponding synthetic orbiter lift-off accelerations model was developed. It is shown that the resulting simulated accelerations preserve significant properties of the original record. Statistically based acceleration response spectra curves for the simulated orbiter lift-off excitations are also evaluated and the results are compared with those for the actual flight data.

Lee-Glauser, Gina↗
NASA NTRS1993

A local dynamic model for large eddy simulation

The dynamic model is a method for computing the coefficient C in Smagorinsky's model for the subgrid-scale stress tensor as a function of position from the information already contained in the resolved velocity field rather than treating it as an adjustable parameter. A variational formulation of the dynamic model is described that removes the inconsistency associated with taking C out of the filtering operation. This model, however, is still unstable due to the negative eddy-viscosity. Next, three models are presented that are mathematically consistent as well as numerically stable. The first two are applicable to homogeneous flows and flows with at least one homogeneous direction, respectively, and are, in fact, a rigorous derivation of the ad hoc expressions used by previous authors. The third model in this set can be applied to arbitrary flows, and it is stable because the C it predicts is always positive. Finally, a model involving the subgrid-scale kinetic energy is presented which attempts to model backscatter. This last model has some desirable theoretical features. However, even though it gives results in LES that are qualitatively correct, it is outperformed by the simpler constrained variational models. It is suggested that one of the constrained variational models should be used for actual LES while theoretical investigation of the kinetic energy approach should be continued in an effort to improve its predictive power and to understand more about backscatter.

Ghosal, Sandip↗
DOE OSTI2024

Seismic Signal Detection on International Monitoring System 3-Component Stations using PhaseNet

In this report we discuss training a deep learning seismic signal detection model on 3-component stations from the International Monitoring System (IMS) using the PhaseNet architecture. Using 14 years of associated signals from the International Data Centre’s (IDC) Late Event Bulletin (LEB), we auto-curated training data consisting of signal windows containing associated arrivals, and noise windows that contain no LEB-associated signals. We trained several models using different waveform window durations (30 seconds and 100 seconds), with and without bandpass filtering. We evaluated the effectiveness of our models using associated signals from the Unconstrained Global Event Bulletin (UGEB) and found that several of our models outperformed the signal detections from the IDC’s Selected Event List 3 (SEL3) arrival table. The SEL3 bulletin evaluated on the UGEB dataset with 100-second waveform windows registered a precision and recall of .15 and .48, respectively, versus .19 and .59 for our filtered-data model. For the 30-second waveform window dataset, the SEL3 bulletin achieved a precision and recall of .31 and .47, respectively, versus .32 and .60 for our filtered-data model. Finally, our models detected signals from all source-to-receiver distances, suggesting it is feasible to use a single PhaseNet model for the IMS network.

58 GEOSCIENCES↗
NASA NTRS1991

Restoration of color images by multichannel Kalman filtering

A Kalman filter for optimal restoration of multichannel images is presented. This filter is derived using a multichannel semicausal image model that includes between-channel degradation. Both stationary and nonstationary image models are developed. This filter is implemented in the Fourier domain and computation is reduced from O(Lambda3N3M4) to O(Lambda3N3M2) for an M x M N-channel image with degradation length Lambda. Color (red, green, and blue (RGB)) images are used as examples of multichannel images, and restoration in the RGB and YIQ domains is investigated. Simulations are presented in which the effectiveness of this filter is tested for different types of degradation and different image model estimates.

Galatsanos, Nikolas P.↗
DOE OSTI2020

Large Volume Airborne Contamination Monitoring To Support Nuclear Processes' Deactivation and Decommissioning

Current D and D operations at Hanford have demonstrated a flaw in the current state of the art capability of defining airborne contamination boundaries - Airborne particulate emissions of Pu-239 from CM2H operations on the Hanford Pu Finishing Plant were detected well beyond areas controlled for airborne Pu, putting numerous workers at risk for radiological assimilations - Current air samplers and CAM systems were surveying insufficient volumes of air to accurately predict where respiratory protection was required - Hanford is located in a high radon environment, making air monitoring for alpha emitting actinides challenging in high alpha radon induced backgrounds. - Inexpensive HEPA based home and industrial air purifiers filter significantly higher volumes of air than commercially available continuous-air monitoring (CAM) systems. - Inexpensive models capable of filtered air volumes exceeding 1000 times that of a CAM. - Detector system can be built to detect x-rays generated from actinide decay. - Branching ratios of x-rays are 4 orders of magnitude more intense from Plutonium decay than its gamma emissions. - Detector system can be a low resolution systems, as x-ray region of interest is not terribly congested. - Detection systems evaluated will be designed based on slabs of NaI or pixelated NaI or CsI panels. - Simple graphical user interface software will be developed to operate the detection system. - Concept is to deploy some of these air-purifying systems in nonradiological areas around SRNL to confirm radon can be rejected with confidence. - Generate some filters contaminated with plutonium via electroplated Pu or lab generated simulated particles. - then deploy units in known airborne radiological environments to establish systems ability to accurately measure actinide based hot particles Lab analyses will follow up the system analyses to ensure hot particles were correctly identified. Analysis by Scintillation: - Application of scintillation media to the surface of the filters is being explored. - Evaluating slabs of ZnS(Ag), application of powdered ZnS(Cu) and spray on Perkin Elmer Enhance. - Scintillation events would then be digitized with a digital camera and quantified. - Comparing against sensitivity of a PMT or SiPD readout. - Currently evaluating a Thorlabs 8 Megapixel Monochrome Scientific CCD Camera, hermetically sealed cooled package with a wide angle lens. - Wide angle lens allows complete view of HEPA filter from 7 inches away. - Images taken of glow in the dark paint, as well as plutonium-induced fluorescence. - Plutonium was flamed mounted on a 1 inch diameter stainless steel planchet, covered with a layer of mylar and a section of scintillating ZnS obtained from Eljon. - Currently working on reducing signal-to-noise levels to boost sensitivity. - Dark box to hold contaminated or electroplated source covered filters under fabrication. Analysis by x-ray Spectroscopy: - Branching ratios for x-ray emissions from Plutonium isotopes are orders of magnitude more intense than gamma-ray emissions. - The nuclear databases are incomplete on the branching ratios of some of the isotopes. - Evaluating some custom built SrI2 x-ray spectrometers vs a conventional windowed NaI detector. - SrI2 detectors are carbon-composite-windowed 51 mm x 51 mm. - MCNP calculations to establish self absorbance of filter media on 17 keV x-ray. Radon Rejection: - Activated charcoal pre-filter acts as a Radon trap. - Adding a time lapse feature to camera to aid in radon rejection. SRS Plutonium Fuel Form D and D Operations: Currently deploying air sampler to D and D operation of Pu-238 facility at SRS to generate some field samples to analyze.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1992

Dynamic rain fade compensation techniques for the advanced communications technology satellite

The dynamic and composite nature of propagation impairments that are incurred on earth-space communications links at frequencies in and above the 30/20 GHz Ka band necessitate the use of dynamic statistical identification and prediction processing of the fading signal in order to optimally estimate and predict the levels of each of the deleterious attenuation components. Such requirements are being met in NASA's Advanced Communications Technology Satellite (ACTS) project by the implementation of optimal processing schemes derived through the use of the ACTS Rain Attenuation Prediction Model and nonlinear Markov filtering theory. The ACTS Rain Attenuation Prediction Model discerns climatological variations on the order of 0.5 deg in latitude and longitude in the continental U.S. The time-dependent portion of the model gives precise availability predictions for the 'spot beam' links of ACTS. However, the structure of the dynamic portion of the model, which yields performance parameters such as fade duration probabilities, is isomorphic to the state-variable approach of stochastic control theory and is amenable to the design of such statistical fade processing schemes which can be made specific to the particular climatological location at which they are employed.

Manning, Robert M.↗
NASA NTRS1994

The Rb 780-nanometer Faraday anomalous dispersion optical filter: Theory and experiment

The Faraday anomalous dispersion optical filter may provide ultra-high background noise rejection for free-space laser communications systems. The theoretical model for the filter is reported. The experimental measurements and their comparison with theoretical results are discussed. The results show that the filter can provide a 56-dB solar background noise rejection with about a 2-GHz transmission bandwidth and no image degradation. To further increase the background noise rejection, a composite Zeeman and Faraday anomalous dispersion optical filter is designed and experimentally demonstrated.

Yin, B.↗
NASA NTRS2009

Satellite-Scale Snow Water Equivalent Assimilation into a High-Resolution Land Surface Model

An ensemble Kalman filter (EnKF) is used in a suite of synthetic experiments to assimilate coarse-scale (25 km) snow water equivalent (SWE) observations (typical of satellite retrievals) into fine-scale (1 km) model simulations. Coarse-scale observations are assimilated directly using an observation operator for mapping between the coarse and fine scales or, alternatively, after disaggregation (re-gridding) to the fine-scale model resolution prior to data assimilation. In either case observations are assimilated either simultaneously or independently for each location. Results indicate that assimilating disaggregated fine-scale observations independently (method 1D-F1) is less efficient than assimilating a collection of neighboring disaggregated observations (method 3D-Fm). Direct assimilation of coarse-scale observations is superior to a priori disaggregation. Independent assimilation of individual coarse-scale observations (method 3D-C1) can bring the overall mean analyzed field close to the truth, but does not necessarily improve estimates of the fine-scale structure. There is a clear benefit to simultaneously assimilating multiple coarse-scale observations (method 3D-Cm) even as the entire domain is observed, indicating that underlying spatial error correlations can be exploited to improve SWE estimates. Method 3D-Cm avoids artificial transitions at the coarse observation pixel boundaries and can reduce the RMSE by 60% when compared to the open loop in this study.

De Lannoy, Gabrielle J.M.↗
NASA NTRSDate not supplied

Examination of L-Band Brightness Temperature Forecasts in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission [4] has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model [5] using a spatially distributed ensemble Kalman filter, the NASA Global Modeling and Assimilation Office generates the SMAP Level-4 Soil Moisture (L4_SM) product, which provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications [6]. The L4_SM product also includes estimates of soil temperature, land surface fluxes, and assimilation diagnostics such as the model forecast and observed Tb values [7]. The output from the L4_SM system is routinely monitored by the L4_SM team. Such monitoring provides valuable information; instances of unusually large Tb observation-minus-forecast (O-F) residuals can indicate events during which soil moisture conditions are poorly described in the land modeling system [7,8]. For example, repeated occurrences of very large Tb O-F values in central Australia were traced back to deficiencies in the gauge-based precipitation product used in the land modeling system through L4_SM Version 5; this discovery prompted the use, beginning in Version 6, of satellite- and gauge-based precipitation observations outside of North America [9]. Recently, a statistical analysis of the Tb O-F residuals revealed systematic errors in the “tau-omega” L-band radiative transfer model that converts the land-model simulated soil moisture and temperature into the forecast Tb prior to the L4_SM analysis. The L4_SM analysis is built on the ensemble Kalman filter and assumes unbiased forecast errors. The presence of systematic Tb forecast errors could thus adversely impact the quality of the analyzed soil moisture. In this paper, we examine the Tb O-F residuals of the latest Version 7 L4_SM data (Science Version ID Vv7030 and Vv7032) [10,11,12]. In the L4_SM Version 7 algorithm, key parameters of the L-band radiative transfer model, including soil roughness, scattering albedo, and a (seasonally varying) climatology of vegetation opacity, are obtained from the SMAP Level-2 Radiometer retrieval product (Version 5) [1]. Additionally, we also use ground measurements of surface soil moisture and soil temperature from the SMAP core validation sites as in situ reference of soil conditions [2,3]. As part of the L4_SM system calibration, the seasonally-varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations [7]. That is, the L4_SM system is designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (potentially erroneous) climatology. The Tb O-F residuals examined here are thus computed using SMAP Tb observations after they are rescaled to the mean seasonal cycle of Tb from the modeling system without SMAP data assimilation. Consequently, the long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. Despite the small time-average values of the Tb O-F residuals, the model forecast Tb was nevertheless found to exhibit undesirable systematic errors. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to further improve the L4_SM algorithm.

Rolf Reichle↗
NASA NTRS1997

An Extended Kalman Filter to Assimilate Altimetric Data into a Non-Linear Model of the Tropical Pacific

A new implementation of the extended Kaman filter is developed for the purpose of assimilating altimetric observations into a primitive equation model of the tropical Pacific. Its specificity consists in defining the errors into a reduced basis that evolves in time with the model dynamic. Validation by twin experiments is conducted and the method is shown to be efficient in quasi real conditions. Data from the first 2 years of the Topex/Poseidon mission are assimilated into the Gent & Cane [1989] model. Assimilation results are evaluated against independent in situ data, namely TAO mooring observations.

Gourdeau, L.↗
DOE OSTI2020

DPM: A deep learning PDE augmentation method with application to large-eddy simulation

A framework is introduced that leverages known physics to reduce overfitting in machine learning for scientific applications. The partial differential equation (PDE) that expresses the physics is augmented with a neural network that uses available data to learn a description of the corresponding unknown or unrepresented physics. Training within this combined system corrects for missing, unknown, or erroneously represented physics, including discretization errors associated with the PDE's numerical solution. For optimization of the network within the PDE, an adjoint PDE is solved to provide high-dimensional gradients, and a stochastic adjoint method (SAM) further accelerates training. Additionally, the approach is demonstrated for large-eddy simulation (LES) of turbulence. High-fidelity direct numerical simulations (DNS) of decaying isotropic turbulence provide the training data used to learn sub-filter-scale closures for the filtered Navier–Stokes equations. Out-of-sample comparisons show that the deep learning PDE method outperforms widely-used models, even for filter sizes so large that they become qualitatively incorrect. It also significantly outperforms the same neural network when a priori trained based on simple data mismatch, not accounting for the full PDE. Measures of discretization errors, which are well-known to be consequential in LES, point to the importance of the unified training formulation's design, which without modification corrects for them. For comparable accuracy, simulation runtime is significantly reduced. A relaxation of the typical discrete enforcement of the divergence-free constraint in the solver is also successful, instead allowing the DPM to approximately enforce incompressibility physics. Since the training loss function is not restricted to correspond directly to the closure to be learned, training can incorporate diverse data, including experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Earth Filter Mathematica Toolbox V4 User’s Guide: Simplifying, Visualizing and Quick Modeling of All Fission Products from 60 Nuclear Sources

The Earth Filter Mathematica toolbox is a set of user-friendly routines that accurately calculate the ingrowth and decay of any of the 908 fission products for an infinite combination of sources over a large dynamic range and large time scales. This toolbox was developed by Dr. W. P. Brug at Los Alamos National Laboratory and eliminates the need for the user to define and solve differential equations by automatically scanning the ENDF-349 dataset and creating a complete set of differential equations (including delayed neutron branches to both the ground and metastable states) for every fission product. Additionally, this code allows the user to choose the data source for half-lives from several key references (i.e., TORI, IAEA, ENDF or NUDAT.) The complete set of ENDF-349 independent and cumulative fission yields for all 60 isotopic targets and 4 different neutron spectrums are included. This allows nearly any nuclear fission source to be accurately modeled as a linear combination of endpoint sources. The Earth Filter Toolbox has no limit as to the length of split list, and split lists of 5 or more endpoint sources can be quite common.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗