Passive Macromodeling: a Review of Cross-Disciplinary Applications

Invited lectur

Passive Macromodeling: a Review of Cross-Disciplinary Applications

Invited lectur

Interpolatory methods for H∞\mathcal{H}_\infty model reduction of multi-input/multi-output systems

We develop here a computationally effective approach for producing high-quality $\mathcal{H}_\infty$-approximations to large scale linear dynamical systems having multiple inputs and multiple outputs (MIMO). We extend an approach for $\mathcal{H}_\infty$ model reduction introduced by Flagg, Beattie, and Gugercin for the single-input/single-output (SISO) setting, which combined ideas originating in interpolatory $\mathcal{H}_2$-optimal model reduction with complex Chebyshev approximation. Retaining this framework, our approach to the MIMO problem has its principal computational cost dominated by (sparse) linear solves, and so it can remain an effective strategy in many large-scale settings. We are able to avoid computationally demanding $\mathcal{H}_\infty$ norm calculations that are normally required to monitor progress within each optimization cycle through the use of "data-driven" rational approximations that are built upon previously computed function samples. Numerical examples are included that illustrate our approach. We produce high fidelity reduced models having consistently better $\mathcal{H}_\infty$ performance than models produced via balanced truncation; these models often are as good as (and occasionally better than) models produced using optimal Hankel norm approximation as well. In all cases considered, the method described here produces reduced models at far lower cost than is possible with either balanced truncation or optimal Hankel norm approximation

On tuning passive black-box macromodels of LTI systems via adaptive weighting

This paper discusses various approaches for tuning the accuracy of rational macromodels obtained via black-box identification or approximation of sampled frequency responses of some unknown Linear and Time-Invariant system. Main emphasis is on embedding into the model extraction process some information on the nominal terminations that will be connected to the model during normal operation, so that the corresponding accuracy is optimized. This goal is achieved through an optimization based on a suitably defined cost function, which embeds frequency-dependent weights that are adaptively refined during the model construction. A similar procedure is applied in a postprocessing step for enforcing model passivity. The advantages of proposed algorithm are illustrated on a few application examples related to power distribution networks in electronic system

The association of selected multiple sclerosis symptoms with disability and quality of life:a large Danish self-report survey

Abstract Background People with multiple sclerosis (MS) experience a wide range of unpredictable and variable symptoms. The symptomatology of MS has previously been reported in large sample registry studies; however, some symptoms may be underreported in registries based on clinician-reported outcomes and how the symptoms are associated with quality of life (QoL) are often not addressed. The aim of this study was to comprehensively evaluate the frequency of selected MS related symptoms and their associations with disability and QoL in a large self-report study. Methods We conducted a cross-sectional questionnaire survey among all patients at the Danish Multiple Sclerosis Center, Copenhagen University Hospital, Denmark. The questionnaire included information on clinical and sociodemographic characteristics, descriptors of QoL and disability, as well as prevalence and severity of the following MS symptoms: impaired ambulation, spasticity, chronic pain, fatigue, bowel and bladder dysfunction, and sleep disturbances. Results Questionnaires were returned by 2244/3606 (62%). Participants without MS diagnosis or incomplete questionnaires were excluded, n = 235. A total of 2009 questionnaires were included for analysis (mean age 49.4 years; mean disease duration 11.7 years; and 69% were women). The most frequently reported symptoms were bowel and bladder dysfunction (74%), fatigue (66%), sleep disturbances (59%), spasticity (51%) and impaired ambulation (38%). With exception of fatigue and sleep disturbances, all other symptoms increased in severity with higher disability level. Invisible symptoms (also referred to as hidden symptoms) such as fatigue, pain and sleep disturbances had the strongest associations with the overall QoL. Conclusion We found invisible symptoms highly prevalent, even at mild disability levels. Fatigue, pain and sleep disturbances had the strongest associations with the overall QoL and were more frequently reported in our study compared with previous registry-based studies. These symptoms may be underreported in registries based on clinician reported outcomes, which emphasizes the importance of including standardized patient reported outcomes in nationwide registries to better understand the impact of the symptom burden in MS

Progress in measuring detonation wave profiles in PBX9501

The authors have measured detonation wave profiles in PBX9501 (95 wt% HMX and 5 wt% binders) using VISAR. Planar detonations were produced by impacting the explosive with projectiles launched in a 72 mm bore gas gun. Particle velocity wave profiles were measured at the explosive/window interface using two VISARs with different fringe constants. Windows with very thin vapor deposited aluminum mirrors were used for all experiments. PMMA windows provided an undermatch, and LiF (Lithium Fluoride) windows provided an overmatch to the explosive, reacted and unreacted. While the present experiments do not have adequate time resolution to adequately resolve the ZND spike condition, they do constrain it to lie between 38.7 and 53.4 Gpa or 2.4 and 3.3 km/s. Accurate knowledge of the CJ state places the reaction zone length at 35 {+-} 12 ns ({approx} 0.3 mm). The present experiments do not show any effect of the window on the reaction zone; both window materials result in the same reaction zone length

Detonation wave profiles in HMX based explosives

Detonation wave profiles have been measured in several HMX based plastic bonded explosives including PBX9404, PBX9501, and EDC-37, as well as two HMX powders (coarse and fine) pressed to 65% of crystal density. The powders had 120 and 10 {micro}m average grain sizes, respectively. Planar detonations were produced by impacting the explosive with projectiles launched in a 72-mm bore gas gun. Impactors, impact velocity, and explosive thickness were chosen so that the run distance to detonation was always less than half the explosive thickness. For the high density plastic bonded explosives, particle velocity wave profiles were measured at an explosive/window interface using two VISAR interferometers. PMMA windows with vapor deposited aluminum mirrors were used for all experiments. Wave profiles for the powdered explosives were measured using magnetic particle velocity gauges. Estimates of the reaction zone parameters were obtained from the profiles using Hugoniots of the explosive and window

Ignition of thermally sensitive explosives between a contact surface and a shock

The dynamics of ignition between a contact surface and a shock wave is investigated using a one-step reaction model with Arrhenius kinetics. Both large activation energy asymptotics and high-resolution finite activation energy numerical simulations are employed. Emphasis is on comparing and contrasting the solutions with those of the ignition process between a piston and a shock, considered previously. The large activation energy asymptotic solutions are found to be qualitatively different from the piston driven shock case, in that thermal runaway first occurs ahead of the contact surface, and both forward and backward moving reaction waves emerge. These waves take the form of quasi-steady weak detonations that may later transition into strong detonation waves. For the finite activation energies considered in the numerical simulations, the results are qualitatively different to the asymptotic predictions in that no backward weak detonation wave forms, and there is only a weak dependence of the evolutionary events on the acoustic impedance of the contact surface. The above conclusions are relevant to gas phase equation of state models. However, when a large polytropic index more representative of condensed phase explosives is used, the large activation energy asymptotic and finite activation energy numerical results are found to be in quantitative agreement

Observations of shock-induced reaction in liquid bromoform up to 11 GPA

Shock measurements on bromoform (CHBr{sub 3}) over the past 33 years at Los Alamos have led to speculation that this material undergoes a shock-induced reaction. Ramsay observed that it became opaque after a 1 to 2 {micro}s induction time when shocked to pressures above 6 GPa. McQueen and Isaak observed that it is a strong light emitter above 25 GPa. Hugoniot data start to deviate from the anticipated liquid Hugoniot at pressures above 10 GPa. The authors have used electromagnetic particle velocity gauging to measure wave profiles in shocked liquid bromoform. At pressures below 9 GPa, there is no mechanical evidence of reaction. At a pressure slightly above 10 GPa, the observed wave profiles are similar to those observed in initiating liquid explosives such as nitromethane. Their characteristics are completely different from the two-wave structures observed in shocked liquids where the products are more dense than the reactants. As with explosives, a reaction producing products which are less dense than the reactants is indicated. BKW calculations also indicate that a detonation type reaction may be possible