Mancha3D Code: Multipurpose Advanced Nonideal MHD Code for High-Resolution Simulations in Astrophysics.

The Mancha3D code is a versatile tool for numerical simulations of magnetohydrodynamic (MHD) processes in solar/stellar atmospheres. The code includes nonideal physics derived from plasma partial ionization, a realistic equation of state and radiative transfer, which allows performing high-quality realistic simulations of magnetoconvection, as well as idealized simulations of particular processes, such as wave propagation, instabilities or energetic events. The paper summarizes the equations and methods used in the Mancha3D (Multifluid (-purpose -physics -dimensional) Advanced Non-ideal MHD Code for High resolution simulations in Astrophysics 3D) code. It also describes its numerical stability and parallel performance and efficiency. The code is based on a finite difference discretization and a memory-saving Runge-Kutta (RK) scheme. It handles nonideal effects through super-time-stepping and Hall diffusion schemes, and takes into account thermal conduction by solving an additional hyperbolic equation for the heat flux. The code is easily configurable to perform different kinds of simulations. Several examples of the code usage are given. It is demonstrated that splitting variables into equilibrium and perturbation parts is essential for simulations of wave propagation in a static background. A perfectly matched layer (PML) boundary condition built into the code greatly facilitates a nonreflective open boundary implementation. Spatial filtering is an important numerical remedy to eliminate grid-size perturbations enhancing the code stability. Parallel performance analysis reveals that the code is strongly memory bound, which is a natural consequence of the numerical techniques used, such as split variables and PML boundary conditions. Both strong and weak scalings show adequate performance up to several thousands of processors (CPUs).

Parametric resonance of Alfvén waves driven by ionization-recombination waves in the weakly ionized solar atmosphere.

Parametric coupling of waves is one of the most efficient mechanisms of energy transfer that can lead to the growth or decay of waves. This transfer occurs at frequencies close to their natural frequencies. In partially ionized solar plasma, there are a multitude of waves that can undergo this process. Here, we study the parametric coupling of Alfvén waves propagating in a partially ionized solar plasma with ionization-recombination waves identified by our study to appear in a plasma in ionization non-equilibrium. Depending on the parameters that describe the plasma (density, temperature), coupling can lead to a parametric resonance. Our study determines the occurrence conditions of parametric resonance, by finding the boundaries between stable and unstable regions in the parameter space. Our results show that collisions and non-equilibrium recombination can both contribute to the onset of unstable behaviour of parametrically resonant Alfvén waves. This article is part of the theme issue 'Partially ionized plasma of the solar atmosphere: recent advances and future pathways'.

Clinical characteristics and outcomes of maintenance hemodialysis patients with COVID-19 during the Omicron wave of the pandemic in Beijing: a single center retrospective study.

BACKGROUND: The clinical manifestations and prognosis of hemodialysis patients with severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) during the Omicron wave of the pandemic infection were still unclear. This study investigated the clinical characteristics of patients undergoing maintenance hemodialysis (MHD) infected with it. METHODS: This retrospective single-center study included 151 patients undergoing MHD. Healthcare workers were selected as control group were assessed from December 1, 2022 to March 31, 2023. Clinical data, laboratory test results, treatment protocols, and prognoses were collected and analyzed. RESULTS: The study population included 146 patients with MHD, 93 (63.7%) of whom were infected with SARS-CoV-2. The number of non-severe, severe, and critical cases was 84 (90.3%), 4 (4.3%), and 5 (5.3%), respectively. Six patients (6.5%) died during the study period. The main symptoms of SARS-CoV-2 infection, including fever, cough, and fatigue, were less common in patients with MHD than the controls. During SARS-CoV-2 infection, the C-reactive protein (2.9 vs. 11.8 mg/dl, p < 0.0001) and ferritin levels(257.7 vs. 537 ng/l, p < 0.0001) were elevated. The hemoglobin(113vs 111 g/L, p = 0.0001) and albumin levels(39.4 vs. 36.1 g/L, p < 0.0001) decreased. Generally, it took two months for the hemoglobin levels to recover. Positivity rate for SARS-COV-2 serum immunoglobin G (IgG) antibodies and IgG titers were lower in dialysis patients than the controls. Age was positively associated with disease severity, while age and hyponatremia were associated with death. CONCLUSIONS: Patients with MHD and COVID-19 were primarily classified as non-severe. SARS-CoV-2 infection would soon lead to the increase of inflammation related acute response protein in dialysis patients, and then lead to the decrease of hemoglobin and albumin. About 9.6% in HD patients were severe cases and had poor prognosis. Advanced age and hyponatremia were associated with disease severity and prognosis.

Research on an Equivalent Heat Source Model of the AC Arc in the Short Gap of a Copper-Core Cable and a Fire Risk Assessment Method.

The magnetohydrodynamics (MHD) model of the alternating current (AC) arc is complex, so a simplified equivalent heat source (EHS) model can be used to replace the complex model in studying the AC arc's thermal characteristics and cable fire risk. A 2D axisymmetric AC arc MHD simulation model in the short gap of a copper-core cable is established in this paper. The AC arc voltage and current obtained by the model are consistent with experiments. The AC arc's heat source distribution obtained by the MHD model is fitted to obtain the heat source function Q of the AC arc. Q is divided into 16 independent segmented heat sources, and a correction matrix is constructed to optimize the segmented heat sources. A neural network and a genetic algorithm give the prediction model and the optimal correction matrix of the segmented heat source. The EHS model optimized by the optimal correction matrix can obtain a minimum temperature error of 5.8/4.4/4.2% with the MHD model in different AC arc peak currents 2/4/6 A. The probability of a cable fire is calculated by using AC arc's optimized EHS model when different numbers of AC arcs are generated randomly in AC half-waves. The EHS model can replace the complex MHD model to study the thermal characteristics of AC arcs and quickly calculate the probability of a cable fire caused by random AC arcs.

An influence of temperature jump and Navier's slip-on hybrid nano fluid flow over a permeable stretching/shrinking sheet with heat transfer and inclined MHD.

This research article, explores the influence of an inclined magnetic field on the fluid flow over a permeable stretching/shrinking surface with heat transfer. The study use water as a conventional base fluid, with graphene oxide (GO) and Aluminum oxide (Al2O3) nanoparticles submerged to create a nanofluid, the system of governing nonlinear partial differential equations converted into ordinary differential equations via suitable similarity conversions. This allow for the unique solution for stretching sheet/shrinking sheets to be obtained, along with the corresponding temperature solution in terms of the hypergeometric function, several parameters are included in the investigation and their contribution is graphically explained to examine physical characteristics such as radiation, inclined magnetic field, solution domain, volume fraction parameter, and temperature jump. Increasing the volume fraction and thermal radiation increases the thermal boundary layer, increasing the magnetic field parameter and inverse Darcy number increases the temperature and decays the velocity profile. The present work has many useful applications in engineering, biological and physical sciences, as well as in cleaning engine lubricants and thrust-bearing technologies.

Numerical analysis of magnetohydrodynamics in an Eyring-Powell hybrid nanofluid flow on wall jet heat and mass transfer.

The customization of hybrid nanofluids to achieve a particular and controlled growth rate of thermal transport is done to meet the needs of applications in heating and cooling systems, aerospace and automotive industries, etc. Due to the extensive applications, the aim of the current paper is to derive a numerical solution to a wall jet flow problem through a stretching surface. To study the flow problem, authors have considered a non-Newtonian Eyring-Powell hybrid nanofluid with water and CoFe2O4and TiO2nanoparticles. Furthermore, the impact of a magnetic field and irregular heat sink/source are studied. To comply with the applications of the wall jet flow, the authors have presented the numerical solution for two cases; with and without a magnetic field. The numerical solution is derived with a similarity transformation and MATLAB-based bvp4c solver. The value of skin friction for wall jet flow at the surface decreases by more than 50% when the magnetic fieldMA=0.2is present. The stream function value is higher for the wall jet flow without the magnetic field. The temperature of the flow rises with the dominant strength of the heat source parameters. The results of this investigation will be beneficial to various applications that utilize the applications of a wall jet, such as in car defrosters, spray paint drying for vehicles or houses, cooling structures for the CPU of high-processor laptops, sluice gate flows, and cooling jets over turbo-machinery components, etc.

Taylor-Couette flow for astrophysical purposes.

A concise review is given of astrophysically motivated experimental and theoretical research on Taylor-Couette flow. The flows of interest rotate differentially with the inner cylinder faster than the outer, but are linearly stable against Rayleigh's inviscid centrifugal instability. At shear Reynolds numbers as large as [Formula: see text], hydrodynamic flows of this type (quasi-Keplerian) appear to be nonlinearly stable: no turbulence is seen that cannot be attributed to interaction with the axial boundaries, rather than the radial shear itself. Direct numerical simulations agree, although they cannot yet reach such high Reynolds numbers. This result indicates that accretion-disc turbulence is not purely hydrodynamic in origin, at least insofar as it is driven by radial shear. Theory, however, predicts linear magnetohydrodynamic (MHD) instabilities in astrophysical discs: in particular, the standard magnetorotational instability (SMRI). MHD Taylor-Couette experiments aimed at SMRI are challenged by the low magnetic Prandtl numbers of liquid metals. High fluid Reynolds numbers and careful control of the axial boundaries are required. The quest for laboratory SMRI has been rewarded with the discovery of some interesting inductionless cousins of SMRI, and with the recently reported success in demonstrating SMRI itself using conducting axial boundaries. Some outstanding questions and near-future prospects are discussed, especially in connection with astrophysics. This article is part of the theme issue 'Taylor-Couette and related flows on the centennial of Taylor's seminal Philosophical Transactions paper (Part 2)'.

Development and validation of a multivariate model for predicting heart failure hospitalization and mortality in patients receiving maintenance hemodialysis.

BACKGROUND: Heart failure (HF) in patients undergoing maintenance hemodialysis (MHD) increases their hospitalization rates, mortality, and economic burden significantly. We aimed to develop and validate a predictive model utilizing contemporary deep phenotyping for individual risk assessment of all-cause mortality or HF hospitalization in patients on MHD. MATERIALS AND METHODS: A retrospective review was conducted from January 2017 to October 2022, including 348 patients receiving MHD from four centers. The variables were adjusted by Cox regression analysis, and the clinical prediction model was constructed and verified. RESULTS: The median follow-up durations were 14 months (interquartile range [IQR] 9-21) for the modeling set and 14 months (9-20) for the validation set. The composite outcome occurred in 72 (29.63%) of 243 patients in the modeling set and 39 (37.14%) of 105 patients in the validation set. The model predictors included age, albumin, history of cerebral hemorrhage, use of angiotensin-converting enzyme inhibitors/angiotensin II receptor blockers/"sacubitril/valsartan", left ventricular ejection fraction, urea reduction ratio, N-terminal prohormone of brain natriuretic peptide, and right atrial size. The C-index was 0.834 (95% CI 0.784-0.883) for the modeling set and 0.853 (0.798, 0.908) for the validation set. The model exhibited excellent calibration across the complete risk profile, and the decision curve analysis (DCA) suggested its ability to maximize patient benefits. CONCLUSION: The developed prediction model offered an accurate and personalized assessment of HF hospitalization risk and all-cause mortality in patients with MHD. It can be employed to identify high-risk patients and guide treatment and follow-up.

Mathematical modelling and analysis of thermoregulation effects on blood viscosity under magnetic effects and thermal radiation in a permeable stretching capillary.

The current study deals with the mathematical modelling and analysis of the effects of thermoregulation on blood viscosity under magnetic and thermal radiation effects on a permeable, stretching blood capillary. The model comprises the governing equations of the resulting boundary layer problem, which is a set of nonlinear partial differential equations, and this is transformed into a coupled system of nonlinear differential equations using similarity transformation. The numerical solution of the problem is attained by the fourth-order Runge Kutta method with the shooting method. The current work examines the velocity and temperature profile of blood flow with the impact of varying temperature-dependent blood viscosity, involving the variation of several physical parameters such as magnetic field, radiation parameter, permeability parameter, Prandtl number, and surface temperature, as well as their physical interpretation. During certain therapies, the thermoregulation mechanism in humans causes a change in blood flow behaviour. Also, blood flow must be regulated during pathological conditions with the help of radiation, heat, or magnetic effects. It is found that, due to the varying viscosity parameter, the magnetic field effect inhibits heat transport in the body. It has also been found that increasing the permeability parameter enhances mass and heat transfer. The results of this paper could be important in analysing and regulating blood flow and body temperature during hypothermia and hyperthermia therapies.

Fluid-Kinetic Variations in the Storm-Time Inner Magnetosphere.

Storm-time broadband electromagnetic field variations along the interface between the dipolar field of the Earth's inner-magnetosphere and the stretched fields of the plasma-sheet are decomposed as a superposition of fluid-kinetic modes. Using model eigen-vectors operating on the full set of Van Allen Probes fields measurements it is shown how these variations are composed of a broad spectrum of dispersive Alfvén waves with significant spectral energy densities in the fast and slow modes over scales extending into the kinetic range. These modes occupy volumes in k -space that define the field variations observed at each spacecraft frame frequency ( f s c ). They are in aggregate not necessarily planar and often comprise filamentary structures with no distinct propagation direction in the perpendicular plane. Within these volumes the characteristic parallel phase speeds of the fast and Alfvénic modes coincide over a broad range of f s c suggestive of coupling/conversion between modes.

Mental health disorders in newly diagnosed non-functional pitutary adenoma under initial observation: an observational cohort study using the nationwide MarketScan database.

PURPOSE: Nonfunctioning pituitary adenomas account for 15-30% of pituitary tumors. Studies exploring the role of an intracranial tumor diagnosis, specifically nonfunctioning pituitary adenomas, on mental health disorders (MHDs) in patients have been limited. We characterize the incidence and factors affecting the development of MHDs in untreated pituitary adenomas. METHODS: Utilizing a large-scale private payor database, MarketScan, we performed a retrospective study of patients with an untreated pituitary adenomas and corresponding MHD. RESULTS: We found that in patients diagnosed with an untreated pituitary adenomas, approximately 15% were newly diagnosed with a MHD within 1 year of the pituitary adenoma diagnosis. Independent risk factors included female gender and substance abuse. Headaches, visual symptoms, and higher Charlson Co-morbidity indexes were also independently associated with a subsequent diagnosis of MHD. On multivariable analysis, patients in the pituitary tumor cohort were more likely to be diagnosed with a MHD than those in the matched cohort (aOR: 1.31, CI: 1.19-1.44). CONCLUSION: By identifying risk factors, advanced screening can focus on non-operative pituitary adenoma patients at high-risk for the development of MHD.

Numerical investigation of blood flow and red blood cell rheology: the magnetic field effect.

In this study, the motion and deformation of a red blood cell in a Poiseuille flow through microvessels under the effect of a uniform transverse magnetic field is comprehensively investigated to get a better insight into blood hemorheology. The rheology of the RBC and the surrounding blood flow are examined numerically in two dimensions using a Finite Element Method. It is essential to know that the flow patterns of blood change in the presence of an RBC. The simulation results demonstrate that the magnetic field has significant influence on the flow stream and the behavior of the RBC, including the motion and the cells deformation.

A Modular Grad-Div Stabilization Method for Time-Dependent Thermally Coupled MHD Equations.

In this paper, we consider a fully discrete modular grad-div stabilization algorithm for time-dependent thermally coupled magnetohydrodynamic (MHD) equations. The main idea of the proposed algorithm is to add an extra minimally intrusive module to penalize the divergence errors of velocity and improve the computational efficiency for increasing values of the Reynolds number and grad-div stabilization parameters. In addition, we provide the unconditional stability and optimal convergence analysis of this algorithm. Finally, several numerical experiments are performed and further indicated these advantages over the algorithm without grad-div stabilization.

Error Analysis of a PFEM Based on the Euler Semi-Implicit Scheme for the Unsteady MHD Equations.

In this article, we mainly consider a first order penalty finite element method (PFEM) for the 2D/3D unsteady incompressible magnetohydrodynamic (MHD) equations. The penalty method applies a penalty term to relax the constraint "∇·u=0", which allows us to transform the saddle point problem into two smaller problems to solve. The Euler semi-implicit scheme is based on a first order backward difference formula for time discretization and semi-implicit treatments for nonlinear terms. It is worth mentioning that the error estimates of the fully discrete PFEM are rigorously derived, which depend on the penalty parameter ϵ, the time-step size τ, and the mesh size h. Finally, two numerical tests show that our scheme is effective.

Two-Level Finite Element Iterative Algorithm Based on Stabilized Method for the Stationary Incompressible Magnetohydrodynamics.

In this paper, based on the stabilization technique, the Oseen iterative method and the two-level finite element algorithm are combined to numerically solve the stationary incompressible magnetohydrodynamic (MHD) equations. For the low regularity of the magnetic field, when dealing with the magnetic field sub-problem, the Lagrange multiplier technique is used. The stabilized method is applied to approximate the flow field sub-problem to circumvent the inf-sup condition restrictions. One- and two-level stabilized finite element algorithms are presented, and their stability and convergence analysis is given. The two-level method uses the Oseen iteration to solve the nonlinear MHD equations on a coarse grid of size H, and then employs the linearized correction on a fine grid with grid size h. The error analysis shows that when the grid sizes satisfy h=O(H2), the two-level stabilization method has the same convergence order as the one-level one. However, the former saves more computational cost than the latter one. Finally, through some numerical experiments, it has been verified that our proposed method is effective. The two-level stabilized method takes less than half the time of the one-level one when using the second class Nédélec element to approximate magnetic field, and even takes almost a third of the computing time of the one-level one when adopting the first class Nédélec element.

Optimal Convergence Analysis of Two-Level Nonconforming Finite Element Iterative Methods for 2D/3D MHD Equations.

Several two-level iterative methods based on nonconforming finite element methods are applied for solving numerically the 2D/3D stationary incompressible MHD equations under different uniqueness conditions. These two-level algorithms are motivated by applying the m iterations on a coarse grid and correction once on a fine grid. A one-level Oseen iterative method on a fine mesh is further studied under a weak uniqueness condition. Moreover, the stability and error estimate are rigorously carried out, which prove that the proposed methods are stable and effective. Finally, some numerical examples corroborate the effectiveness of our theoretical analysis and the proposed methods.

Four-Objective Optimization of an Irreversible Magnetohydrodynamic Cycle.

Based on the existing model of an irreversible magnetohydrodynamic cycle, this paper uses finite time thermodynamic theory and multi-objective genetic algorithm (NSGA-II), introduces heat exchanger thermal conductance distribution and isentropic temperature ratio of working fluid as optimization variables, and takes power output, efficiency, ecological function, and power density as objective functions to carry out multi-objective optimization with different objective function combinations, and contrast optimization results with three decision-making approaches of LINMAP, TOPSIS, and Shannon Entropy. The results indicate that in the condition of constant gas velocity, deviation indexes are 0.1764 acquired by LINMAP and TOPSIS approaches when four-objective optimization is performed, which is less than that (0.1940) of the Shannon Entropy approach and those (0.3560, 0.7693, 0.2599, 0.1940) for four single-objective optimizations of maximum power output, efficiency, ecological function, and power density, respectively. In the condition of constant Mach number, deviation indexes are 0.1767 acquired by LINMAP and TOPSIS when four-objective optimization is performed, which is less than that (0.1950) of the Shannon Entropy approach and those (0.3600, 0.7630, 0.2637, 0.1949) for four single-objective optimizations, respectively. This indicates that the multi-objective optimization result is preferable to any single-objective optimization result.

Dietary Plant Protein and Mortality Among Patients Receiving Maintenance Hemodialysis: A Cohort Study.

RATIONALE & OBJECTIVE: Although greater dietary intake of protein has been associated with beneficial health effects among patients receiving maintenance hemodialysis (MHD), the effects of plant protein intake are less certain. We studied the association of the proportion of protein intake derived from plant sources with the risk of mortality among patients receiving MHD and explored factors that may modify these associations. STUDY DESIGN: Prospective observational cohort study. SETTING & PARTICIPANTS: 1,119 Chinese hemodialysis patients aged over 18 years receiving MHD in 2014-2015. PREDICTORS: The proportion of plant protein intake to total protein intake. OUTCOMES: All-cause mortality and cardiovascular disease (CVD) mortality. ANALYTICAL APPROACH: Segmented regression models were fit to examine the association of plant protein intake proportion with the risk of all-cause mortality and CVD mortality. Multivariable-adjusted Cox proportional and cause-specific hazards models were used to estimate the hazard ratios (HR) and 95% CI for these outcomes. RESULTS: The means of plant protein intake normalized to ideal body weight and plant protein intake proportion were 0.6±0.2 (SD) g/kg per day and 0.538±0.134, respectively. During a median follow-up period of 28.0 months, 249 deaths occurred, with 146 of these deaths resulting from CVD. Overall, there was a U-shaped association between plant protein intake proportion and the risk of all-cause mortality, with an inflection point at 45%. Among patients with a plant protein intake proportion<45%, there was a 17% lower rate of mortality with each 5% greater plant protein intake proportion (HR, 0.83 [95% CI, 0.73-0.96]). Among patients with plant protein intake proportion≥45%, there was a 9% greater rate of mortality with each 5% greater plant protein intake proportion. A similar U-shaped association was observed for CVD mortality, with an inflection point at 44%. LIMITATIONS: Observational study, potential unmeasured confounding. CONCLUSIONS: There was a U-shaped association between plant protein intake proportion and the risk of all-cause and cardiovascular mortality in MHD patients. If confirmed, these findings suggest a potential avenue to improve outcomes in this patient population.

Transverse motions in sunspot super-penumbral fibrils.

Sunspots have played a key role in aiding our understanding of magnetohydrodynamic (MHD) wave phenomena in the Sun's atmosphere, and it is well known they demonstrate a number of wave phenomena associated with slow MHD modes. Recent studies have shown that transverse wave modes are present throughout the majority of the chromosphere. Using high-resolution Ca II 8542 Å observations from the Swedish Solar Telescope, we provide the first demonstration that the chromospheric super-penumbral fibrils, which span out from the sunspot, also show ubiquitous transverse motions. We interpret these motions as transverse waves, in particular the MHD kink mode. We compile the statistical properties of over 2000 transverse motions to find distributions for periods and amplitudes, finding they are broadly consistent with previous observations of chromospheric transverse waves in quiet Sun fibrils. The very presence of the waves in super-penumbral fibrils raises important questions about how they are generated, and could have implications for our understanding of how MHD wave energy is transferred through the atmosphere of a sunspot. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

Spectropolarimetric fluctuations in a sunspot chromosphere.

The instrumental advances made in this new era of 4 m class solar telescopes with unmatched spectropolarimetric accuracy and sensitivity will enable the study of chromospheric magnetic fields and their dynamics with unprecedented detail. In this regard, spectropolarimetric diagnostics can provide invaluable insight into magneto-hydrodynamic (MHD) wave processes. MHD waves and, in particular, Alfvénic fluctuations associated with particular wave modes were recently recognized as important mechanisms not only for the heating of the outer layers of the Sun's atmosphere and the acceleration of the solar wind, but also for the elemental abundance anomaly observed in the corona of the Sun and other Sun-like stars (also known as first ionization potential) effect. Here, we take advantage of state-of-the-art and unique spectropolarimetric Interferometric BIdimensional Spectrometer observations to investigate the relation between intensity and circular polarization (CP) fluctuations in a sunspot chromosphere. Our results show a clear link between the intensity and CP fluctuations in a patch which corresponds to a narrow range of magnetic field inclinations. This suggests the presence of Alfvénic perturbations in the sunspot. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.