Search NASA⌕ Search

SEARCH · Search NASA

Results for “consistent 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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Blanco DECam Bulge Survey (BDBS)

The Blanco DECam Bulge Survey (BDBS) has imaged more than 200 square degrees of the southern Galactic bulge, providing photometry in the ugrizy filters for ~250 million unique stars. The presence of a strong foreground disk population, along with complex reddening and extreme image crowding, has made it difficult to constrain the presence of young and intermediate age stars in the bulge population. We employed an accurate cross-match of BDBS with the latest data release (EDR3) from the Gaia mission, matching more than 140 million sources with BDBS photometry and Gaia EDR3 photometry and astrometry. We relied on Gaia EDR3 astrometry, without any photometric selection, to produce clean BDBS bulge colour-magnitude diagrams (CMDs). Gaia parallaxes were used to filter out bright foreground sources, and a Gaussian mixture model fit to Galactic proper motions could identify stars kinematically consistent with bulge membership. We applied this method to 127 different bulge fields of 1 deg 2 each, with |$\ell$| ≤ 9.5 ° and –9.5 ° ≤ b ≤ –2.5 °. The astrometric cleaning procedure removes the majority of blue stars in each field, especially near the Galactic plane, where the ratio of blue to red stars is .10%, increasing to values ~20% at higher Galactic latitudes. We rule out the presence of a widespread population of stars younger than 2 Gyr. The vast majority of blue stars brighter than the turnoff belong to the foreground population, according to their measured astrometry. We introduce the distance between the observed red giant branch bump and the red clump as a simple age proxy for the dominant population in the field, and we confirm the picture of a predominantly old bulge. Further work is needed to apply the method to estimate ages to fields at higher latitudes, and to model the complex morphology of the Galactic bulge. We also produce transverse kinematic maps, recovering expected patterns related to the presence of the bar and of the X-shaped nature of the bulge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Particle model of full-size ITER-relevant negative ion source

This work represents the first attempt to model the full-size ITER-relevant negative ion source including the expansion, extraction, and part of the acceleration regions keeping the mesh size fine enough to resolve every single aperture. The model consists of a 2.5D particle-in-cell Monte Carlo collision representation of the plane perpendicular to the filter field lines. Magnetic filter and electron deflection field have been included and a negative ion current density of j{sub H{sup −}} = 660 A/m{sup 2} from the plasma grid (PG) is used as parameter for the neutral conversion. The driver is not yet included and a fixed ambipolar flux is emitted from the driver exit plane. Results show the strong asymmetry along the PG driven by the electron Hall (E × B and diamagnetic) drift perpendicular to the filter field. Such asymmetry creates an important dis-homogeneity in the electron current extracted from the different apertures. A steady state is not yet reached after 15 μs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

DANTE as a primary temperature diagnostic for the NIF iron opacity campaign

The Opacity Platform on the National Ignition Facility (NIF) has been developed to measure iron opacities at varying densities and temperatures relevant to the solar interior and to verify recent experimental results obtained at the Sandia Z-machine, that diverge from theory. The first set of NIF experiments collected iron opacity data at ∼150 eV to 160 eV and an electron density of ∼7 × 1021 cm−3, with a goal to study temperatures up to ∼210 eV, with electron densities of up to ∼3 × 1022 cm−3. Among several techniques used to infer the temperature of the heated Fe sample, the absolutely calibrated DANTE-2 filtered diode array routinely provides measurements of the hohlraum conditions near the sample. However, the DANTE-2 temperatures are consistently low compared to pre-shot LASNEX simulations for a range of laser drive energies. We have re-evaluated the estimated uncertainty in the reported DANTE-2 temperatures and also the error generated by varying channel participation in the data analysis. An uncertainty of ±5% or better can be achieved with appropriate spectral coverage, channel participation, and metrology of the viewing slot.

Opachich, Y. P. (ORCID:0000000314323989)↗
DOE OSTI2011

THE CHANDRA DEEP FIELD-SOUTH SURVEY: 4 Ms SOURCE CATALOGS

We present source catalogs for the 4 Ms Chandra Deep Field-South (CDF-S), which is the deepest Chandra survey to date and covers an area of 464.5 arcmin{sup 2}. We provide a main Chandra source catalog, which contains 740 X-ray sources that are detected with WAVDETECT at a false-positive probability threshold of 10{sup -5} in at least one of three X-ray bands (0.5-8 keV, full band; 0.5-2 keV, soft band; and 2-8 keV, hard band) and also satisfy a binomial-probability source-selection criterion of P < 0.004 (i.e., the probability of sources not being real is less than 0.004); this approach is designed to maximize the number of reliable sources detected. A total of 300 main-catalog sources are new compared to the previous 2 Ms CDF-S main-catalog sources. We determine X-ray source positions using centroid and matched-filter techniques and obtain a median positional uncertainty of {approx}0.''42. We also provide a supplementary catalog, which consists of 36 sources that are detected with WAVDETECT at a false-positive probability threshold of 10{sup -5}, satisfy the condition of 0.004 < P < 0.1, and have an optical counterpart with R < 24. Multiwavelength identifications, basic optical/infrared/radio photometry, and spectroscopic/photometric redshifts are provided for the X-ray sources in the main and supplementary catalogs. Seven hundred sixteen ({approx}97%) of the 740 main-catalog sources have multiwavelength counterparts, with 673 ({approx}94% of 716) having either spectroscopic or photometric redshifts. The 740 main-catalog sources span broad ranges of full-band flux and 0.5-8 keV luminosity; the 300 new main-catalog sources span similar ranges although they tend to be systematically lower. Basic analyses of the X-ray and multiwavelength properties of the sources indicate that >75% of the main-catalog sources are active galactic nuclei (AGNs); of the 300 new main-catalog sources, about 35% are likely normal and starburst galaxies, reflecting the rise of normal and starburst galaxies at the very faint flux levels uniquely accessible to the 4 Ms CDF-S. Near the center of the 4 Ms CDF-S (i.e., within an off-axis angle of 3'), the observed AGN and galaxy source densities have reached 9800{sup +1300}{sub -1100} deg{sup -2} and 6900{sup +1100}{sub -900} deg{sup -2}, respectively. Simulations show that our main catalog is highly reliable and is reasonably complete. The mean backgrounds (corrected for vignetting and exposure-time variations) are 0.063 and 0.178 counts Ms{sup -1} pixel{sup -1} (for a pixel size of 0.''492) for the soft and hard bands, respectively; the majority of the pixels have zero background counts. The 4 Ms CDF-S reaches on-axis flux limits of {approx}3.2 x 10{sup -17}, 9.1 x 10{sup -18}, and 5.5 x 10{sup -17} erg cm{sup -2} s{sup -1} for the full, soft, and hard bands, respectively. An increase in the CDF-S exposure time by a factor of {approx}2-2.5 would provide further significant gains and probe key unexplored discovery space.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Extraction of the non-spin- and spin-transfer isovector responses via the 12 C ⁡( 10 Be, 10 B + 𝛾)⁢ 12 B reaction

The isovector response in 12 B was investigated via the 12 C ⁡( 10 Be, 10 B + 𝛾)⁢ 12 B* reaction at 100⁢𝐴MeV. By utilizing the 𝛾-decay properties of the 1.74 MeV 0 + and 0.718 MeV 1 + states in 10 B, the separate extraction of the non-spin-transfer (Δ⁢𝑆 = 0) and spin-transfer (Δ⁢𝑆 = 1) isovector responses up to an excitation energy of 50 MeV in 12 B in a single measurement is demonstrated. The experimental setup employed the S800 spectrometer to detect and analyze the 10 B ejectiles and the Gamma-Ray Energy Tracking In-beam Nuclear Array (GRETINA) for obtaining the Doppler-reconstructed spectrum for 𝛾 rays emitted in flight by 10 B. A 12 C foil was placed at the pivot point of the spectrograph. Here, the 12 B reaction product was not detected. Contributions from transitions associated with the transfer of different units of angular momentum in the non-spin- and spin-transfer responses were analyzed using a multipole decomposition analysis. The extracted non-spin-dipole (Δ⁢𝑆 = 0, Δ⁢𝐿 = 1) and spin-dipole (Δ⁢𝑆 = 1, Δ⁢𝐿 = 1) responses were found to be consistent with available data from other charge-exchange probes, validating the non-spin- and spin-transfer filters used. While statistical uncertainties and experimental resolutions were relatively large due to the modest intensity of the 10 Be secondary beam, the results show that, with the much higher intensities that will be available at new rare-isotope beam facilities, the ( 10 Be, 10 B + 𝛾) reaction and its Δ⁢𝑇 𝑧 = −1 partner, the ( 10 C, 10 B + 𝛾) reaction, are powerful tools for elucidating the isovector non-spin- and spin-transfer responses in nuclei.

Charge-exchange reactions↗
DOE OSTI2015

MO-FG-BRA-02: Modulation of Clinical Orthovoltage X-Ray Spectrum Further Enhances Radiosensitization of Cancer Cells Targeted with Gold Nanoparticles

Purpose: To assess the potential to amplify radiosensitization of cancer cells targeted with gold nanoparticles by augmenting selective spectral components of X-ray beam. Methods: Human prostate cancer cells were treated for 24h with gold nanorods conjugated to goserelin acetate or pegylated, systematically washed and irradiated with 250 kVp X-rays (25mA, 0.25mm Cu- filter, 8x8cm{sup 2} field size, 50cm SSD) with or without an additional 0.25 mm Erbium (Er) filter. As demonstrated in a companion Monte Carlo study, Er-filter acted as an external target to feed Erbium K-shell X-ray fluorescence photons (∼50 keV) into the 250 kVp beam. After irradiation, we performed measurements of clonogenic viability with doses between 0 -6Gy, irreparable DNA damage assay to measure double-strand breaks via γH2AX-foci staining, and production of stable reactive oxygen species (ROS). Results: The clonogenic assay for the group treated with conjugated nanoparticles showed radiosensitization enhancement factor (REF), calculated at the 10% survival fraction aisle, of (1.62±0.07) vs. (1.23±0.04) with/without the Er-filter in the 250 kVp beam, respectively. The group treated with pegylated nanoparticles, albeit retained in modest amounts within the cells, also showed statistically significant REF (1.13±0.09) when the Erbium filter was added to the beam. No significant radiosensitization was observed for other groups. Measurements of ROS levels showed increments of (1.9±0.2) vs. (1.4±0.1) for combined treatment with targeted nanoparticles and Er-filtered beam. γH2AX-foci showed 50% increase for the same treatment combination, confirming the enhanced radiosensitization in a consistent fashion. Conclusion: Our study demonstrates the feasibility of enhancing radiosensitization of cancer cells by combining actively targeted gold nanoparticles and modulating the X-ray spectrum in the desired energy range. The established technique will not only help develop strategies to maximize nanoparticle-mediated radiosensitization but also offer a convenient way to acquire unprecedented insights into the role of photon energy for the observed radiosensitization effects. Supported by DOD/PCRP grant W81XWH-12-1-0198.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

KIPM Detector Characterization in QUIET

Kinetic Inductance Phonon‑Mediated (KIPM) detectors are superconducting microwave resonators on silicon that sense bursts of energy through tiny shifts in their resonance. In this work, we installed new devices in the QUIET facility's cryogenic test setup and used a network analyzer to drive them across a range of power levels. At each setting, we measured the transmitted signal and fit it to a simple notch‑filter model, allowing us to extract key performance metrics. And despite the increased noise at low power, our model consistently captured the resonance behavior. We observed that the resonant frequency stayed stable across most powers, the device s internal quality factor improved with stronger signals, and its coupling strength peaked mid‑range before leveling off. Overall, these findings confirm the detector s performance, laying groundwork for future work on energy calibration with pulsed optical sources.

Savitala, Akash [Washington U., Seattle; Fermilab]↗
DOE OSTI2010

Can 3-D models explain the observed fractions of fossil and non-fossil carbon in and near Mexico City?

Abstract. A 3-D chemistry-transport model has been applied to the Mexico City metropolitan area to investigate the origin of elevated levels of non-fossil (NF) carbonaceous aerosols observed in this highly urbanized region. High time resolution measurements of the fine aerosol concentration and composition, and 12 or 24 h integrated 14C measurements of aerosol modern carbon have been performed in and near Mexico City during the March 2006 MILAGRO field experiment. The non-fossil carbon fraction (fNF), which is lower than the measured modern fraction (fM) due to the elevated 14C in the atmosphere caused by nuclear bomb testing, is estimated from the measured fM and the source-dependent information on modern carbon enrichment. The fNF contained in PM1 total carbon analyzed by a US team (f TC NF ) ranged from 0.37 to 0.67 at the downtown location, and from 0.50 to 0.86 at the suburban site. Substantially lower values (i.e. 0.24–0.49) were found for PM10 filters downtown by an independent set of measurements (Swiss team), which are inconsistent with the modeled and known differences between the size ranges, suggesting higher than expected uncertainties in the measurement techniques of 14C. An increase in the non-fossil organic carbon (OC) fraction (f OC NF ) by 0.10–0.15 was observed for both sets of filters during periods with enhanced wildfire activity in comparison to periods when fires were suppressed by rain, which is consistent with the wildfire impacts estimated with other methods. Model results show that the relatively high fraction of nonfossil carbon found in Mexico City seems to arise from the combination in about equal proportions of regional biogenic SOA, biomass burning POA and SOA, as well as non-fossil urban POA and SOA. Predicted spatial and temporal variations for f OC NF are similar to those in the measurements between the urban vs. suburban sites, and high-fire vs. low-fire periods. The absolute modeled values of f OC NF are consistent with the Swiss dataset but lower than the US dataset. Resolving the 14C measurement discrepancies is necessary for further progress in model evaluation. The model simulations that included secondary organic aerosol (SOA) formation from semi-volatile and intermediate volatility (S/IVOC) vapors showed improved closure for the total OA mass compared to simulations which only included SOA from VOCs, providing a more realistic basis to evaluate the fNF predictions. f OC NF urban sources of modern carbon are important in reducing or removing the difference in fNF between model and measurements, even though they are often neglected on the interpretation of 14C datasets. An underprediction of biomass burning POA by the model during some mornings also explains a part of the model-measurement differences. The fNF of urban POA and SOA precursors is an important parameter that needs to be better constrained by measurements. Performing faster ( 3 h) 14C measurements in future campaigns is critical to further progress in this area. To our knowledge this is the first time that radiocarbon measurements are used together with aerosol mass spectrometer (AMS) organic components to assess the performance of a regional model for organic aerosols.

09 BIOMASS FUELS↗
NASA NTRS1973

Mariner/Jupiter/Saturn navigation in the presence of massive planetary satellites

Orbit determination accuracies attainable during a Mariner/Jupiter/Saturn encounter have been established by means of covariance analyses. Advanced earth-based multistation radiometric data, optical data consisting of star/satellite pictures provided by a narrow-angle TV camera on board the spacecraft, and batch sequential filtering methods were employed. An evaluation of sequential filter sensitivities to errors in modeling small nongravitational spacecraft accelerations is presented. Selected covariances developed in the orbit determination analysis are used in an investigation of the trajectory correction costs associated with a close Ganymede encounter. The 99 percentile Delta V contours in the Ganymede aiming plane are obtained for satellite encounters before or after Jupiter closed approach, and the sources of the dominant contributions to the Delta V costs are identified.

Hildebrand, C. E.↗
NASA NTRS1979

Flight test of a VTOL digital autoland system along complex trajectories

The objective of this research is the development and flight test of a digital multi-rate VTOL autoland system. This autoland system provided the NASA/Army Fly-by-Wire CH-47 helicopter with an automatic 4-D trajectory following and landing capability. Flight tests were conducted along complex trajectories with descending, turning, and decelerating segments. The flight software consists of a trajectory generator and guidance algorithm operating at two iterations per second and a set of sensor filters and an advanced controller operating at 10 iterations per second. This low iteration rate digital controller has a proportional-integral structure that provides an autotrim function. The control gains, scheduled as a function of flight condition, are updated every two seconds. All flight software, except input/output routines, are coded in FORTRAN using floating point arithmetic and are exercised in a ROLM 1664 flight computer. The flight software, the ground and flight hardware, the autoland system's position tracking performance, and the controller's velocity and heading tracking performance are presented.

Downing, D. R.↗
NASA NTRS1981

Digital detection and processing of laser beacon signals for aircraft collision hazard warning

A low-cost collision hazard warning system suitable for implementation in both general and commercial aviation is presented. Laser beacon systems are used as sources of accurate relative position information that are not dependent on communication between aircraft or with the ground. The beacon system consists of a rotating low-power laser beacon, detector arrays with special optics for wide angle acceptance and filtering of solar background light, microprocessors for proximity and relative trajectory computation, and pilot displays of potential hazards. The laser beacon system provides direct measurements of relative aircraft positions; using optimal nonlinear estimation theory, the measurements resulting from the current beacon sweep are combined with previous data to provide the best estimate of aircraft proximity, heading, minimium passing distance, and time to closest approach.

Sweet, L. M.↗
NASA NTRS1987

Space station resistojet system requirements and interface definition study

Preliminary resistojet design requirements were established based on initial technical requirements imposed by the results of NASA and Rocketdyne studies. The requirements are directed toward long life, simplicity, flexibility, and commonality with other space station components. The resistojet assembly is comprised of eight resistojets, fluid components downstream of the waste fluid storage system, a power controller, structure, and shielding. It consists of two identical subassemblies, one of which is redundant. Each subassembly consists of four 500-W resistojets, series redundant latch values, a power controller, a water vaporizer, two pressure regulators, filters, check valves, disconnects, fluid tubing, and electrical cables. All components are packaged at the end of the stinger aft of the JEM and Columbus modules. Different flow and power control methods were studied. A constant inlet pressure and a two-power setting controller were tentatively selected based on simplicity and reasonably high specific impulse for the range of waste gas compositions that are anticipated. The constant pressure is supplied by pressure regulators. The two set point power control includes individual power supplies to each resistojet heater and water vaporizer. An embedded data processor, a multiplexer-demultiplexer, and a network interface unit that are standard space station components are included in the power controller. The total dry weight of the resistojet assembly is approximately 172 lb. The total cost for design, development, test, evaluation, qualification, and flight hardware is estimated to be $16 million.

Heckert, Bruce J.↗
NASA NTRS1988

Searching Circuit For A Servoloop

Feedback error voltage forced into range of stability. Circuit connected in series with error-voltage line of feedback loop. When loop operates normally in stable region, searching circuit does not disturb operation, and error-voltage output of loop filter is control voltage. When servo unlocks and error-voltage output of loop filter wanders into unstable range, searching circuit forces control voltage to sweep through stable region repeatedly until lock recovered. Consists of two halves identical except polarities of some parts in each half opposite of corresponding parts in other half.

Sigman, E. H.↗
NASA NTRS1997

Passive Vibration Control of Airborne Equipment using a Circular Steel Ring

Vibration isolation is needed to protect avionics equipment from adverse aircraft vibration environments. Passive isolation is the simplest means to achieve this goal. The system used here consists of a circular steel ring with a lump mass on top and exposed to base excitation. Sinusoidal and filtered zero-mean Gaussian white noise are used to excite the structure and the acceleration response spectra at the top of the ring are computed. An experiment is performed to identify the natural frequencies and modal damping of the circular ring. Comparison is made between the analytical and experimental results and good agreement is observed. The ring response is also evaluated with a concentrated mass attached to the top of the ring. The effectiveness of the ring in isolating the equipment from base excitation is studied. The acceleration response spectra of a single degree of freedom system attached to the top of the ring are evaluated and the results are compared with those exposed directly to the base excitation. It is shown that a properly designed ring could effectively protect the avionics from possible damaging excitation levels.

Ellison, Joseph↗
NASA NTRS2001

On-Orbit Calibration of Redundant Spacecraft Gyros by Optimal Reduction to Three Axes

The Aqua spacecraft will carry four single-axis gyros configured with three orthogonal axes and one skew axis. This redundancy presents a challenge for batch methods of on-orbit gyro calibration that use a spacecraft rotation model deterministically related to gyro data, in that sensor data can respond to at most three angular velocity components. When the number of gyros, N, is greater than 3, the 3xN matrix, G, that reduces the N gyro measurements to three body-frame angular-velocity components cannot be fully determined by such methods; there are many such matrices that produce essentially the same angular velocity history. In such a case, spacecraft operators require information about the Nx3 gyro linear response matrix, R, that relates gyro outputs to the body-frame angular velocities causing them. This matrix provides sufficient information to determine multiple reduced-dimension G-matrices for use in case of failure or degradation of one or more gyros, as well as to determine an optimal 3xN G for the fully-functional configuration. A viable proposal is to apply a 3xN pre-filter matrix, F, to the N gyro outputs before carrying out a conventional gyro calibration procedure. The angular-velocity history emerging from conventional calibration may then be used as input data, together with the same gyro data that generated it, to fit the alignment, scale-factor, and bias parameters of each gyro axis in turn. A difficulty of such a proposal is the arbitrariness in the choice of F. Due to gyro noise, different pre-filter matrices produce different calibrations. This paper presents a method of choosing F that is based on optimizing gyro consistency in the limit of infinite weight on gyro data, as compared to sensor data. The choice of F is independent of a priori alignment and is based on the gyro data alone. The method is applicable to any N of three or more, but reduces to conventional batch-estimation methodologies when N = 3. Results of computational comparison among calibration simulations using various choices of F will be presented for the Aqua gyro configuration with N = 4.

Radomski, M. S.↗
NASA NTRS2011

Broadband Via-Less Microwave Crossover Using Microstrip-CPW Transitions

The front-to-back interface between microstrip and CPW (coplanar waveguide) typically requires complex fabrication or has high radiation loss. The microwave crossover typically requires a complex fabrication step. The prior art in microstrip-CPW transition requires a physical vias connection between the microstrip and CPW line on a separate layer. The via-less version of this transition was designed empirically and does not have a close form solution. The prior art of the micro wave crossover requires either additional substrate or wire bond as an air bridge to isolate two microwave lines at the crossing junction. The disadvantages are high radiation loss, no analytical solution to the problem, lengthy simulation time, and complex fabrication procedures to generate air bridges or via. The disadvantage of the prior crossover is a complex fabrication procedure, which also affects the device reliability and yield. This microstrip-CPW transition is visualized as two microstrip-slotline transitions combined in a way that the radiation from two slotlines cancels each other out. The invention is designed based on analytical methods; thus, it significantly reduces the development time. The crossover requires no extra layer to cross two microwave signals and has low radiation loss. The invention is simple to fabricate and design. It produces low radiation loss and can be designed with low insertion loss, with some tradeoff with signal isolation. The microstrip-CPW transition is used as an interface to connect between the device and the circuit outside the package. The via-less microwave crossover is used to allow two signals to cross without using an extra layer or fabrication processing step to enable this function. This design allows the solution to be determined entirely though analytical techniques. In addition, a planar via-less microwave crossover using this technique was proposed. The experimental results show that the proposed crossover at 5 GHz has a minimum isolation of 32 dB. It also has low in-band insertion loss and return loss of 1.2 dB and 18 dB, respectively, over more than 44 percent of bandwidth at room temperature. This microstrip-CPW transition requires the microstrip line to be split into two sections. Each section is connected to a microstrip quarter-wavelength openended stub. A slotline is also placed perpendicular to the microstrip section. The slot is connected to a grounded-end quarter-wavelength slotline and generates a microstrip-slotline transition. When two of these sections are placed in parallel and with the microstrip section combined at transition, a microstrip- CPW transition is formed. The slotline radiation is suppressed as two slots are excited with the electric field in an opposite direction, which cancels the radiation in far field. The invention on the crossover consists of the invented microstrip-CPW transitions combined back-to-back and a microstrip low-pass filter. One signal is crossed through to the microstrip layer, while the other signal is crossed through the CPW line located on the ground plane of the microstrip line. The microstrip low-pass filter produces a narrow line at the crossing point to enhance the system isolation. It also produces broadband response in the operating frequency band. The microstrip-CPW transition allows a microwave signal to travel from microstrip line to CPW line with low radiation loss. The crossover allows two microwave signals to cross with minimal parasitic coupling.

Stevenson, Thomas↗
NASA NTRS2005

Adaptive Numerical Dissipation Control in High Order Schemes for Multi-D Non-Ideal MHD

The required type and amount of numerical dissipation/filter to accurately resolve all relevant multiscales of complex MHD unsteady high-speed shock/shear/turbulence/combustion problems are not only physical problem dependent, but also vary from one flow region to another. In addition, proper and efficient control of the divergence of the magnetic field (Div(B)) numerical error for high order shock-capturing methods poses extra requirements for the considered type of CPU intensive computations. The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free from numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multiresolution wavelets (WAV) (for the above types of flow feature). These filters also provide a natural and efficient way for the minimization of Div(B) numerical error.

Yee, H. C.↗