researchmap

researchmap

researchmap

researchmap

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astrophysical Journal Letters  

One key aspect of understanding the solar dynamo mechanism and the evolution of solar magnetism is to properly describe the emergence of solar active regions. In this Letter, we describe the Lagrangian photospheric flows dynamics during a simulated flux emergence that produces an active region formed by pores. We analyze the lower photospheric flow organization prior, during and following the rise of an active region, uncovering the repelling and attracting photospheric structures that act as sources and sinks for magnetic element transport. Our results show that around 10 hr before the simulated emergence, considerable global changes are taking place on mesogranular scales indicated by an increase of the number of regions acting as a source to the multiple and scattered emergences of small-scale magnetic flux. At the location of active region’s appearance, the converging flows become weaker and there is an arising of a diverging region 8 hr before the emergence time. Our study also indicates that the strong concentration of magnetic field affects the flow dynamics beyond the area of the actual simulated pores, leading to complex and strongly diverging flows in the neighboring regions. Our findings suggest that the Lagrangian analysis is a powerful tool to describe the changes in the photospheric flows due to magnetic flux emergence.

DOI： 10.3847/2041-8213/acd007

Web of Science

Scopus

researchmap

その他リンク： https://iopscience.iop.org/article/10.3847/2041-8213/acd007/pdf

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Monthly Notices of the Royal Astronomical Society  

We investigate the dependence of the angular momentum transport (AMT) on the spatial scales with numerical simulation of solar-like stars. It is thought that turbulence has an essential role in constructing solar differential rotation (DR). In a widely used method to analyse the construction mechanism of DR, the flow is divided into two components, ‘mean flow’ and ‘turbulence’, where ‘turbulence’ includes a broad spectrum of spatial scales. The features of the AMT are expected to depend on the scale. In this study, we decompose the angular momentum flux (AMF) to investigate the dependence of the AMF on the spatial scale. We compare the results with anti-solar (fast pole) and solar-type (fast equator) DR. Our conclusions are summarized as (1) Radially outward AMT is seen on a large scale (60 Mm ≤ L < 120 Mm) in rotationally constrained systems. (2) Even when the scale-integrated AMF is negative, we sometimes observe positive AMF on certain scales. (3) Small-scale turbulence tends to transport the angular momentum radially inward and causes the anti-solar DR, indicating that high-resolution simulation is a negative factor for solar-like DR. Our method to decompose the AMF provides a deep understanding of the angular momentum and construction mechanism of DR.

DOI： 10.1093/mnras/stac3804

Web of Science

Scopus

researchmap

その他リンク： https://academic.oup.com/mnras/article-pdf/519/2/3091/48522790/stac3804.pdf

掲載種別：研究論文（学術雑誌）   出版者・発行元：Monthly Notices of the Royal Astronomical Society  

Strong solar flares occur in δ-spots characterized by the opposite-polarity magnetic fluxes in a single penumbra. Sunspot formation via flux emergence from the convection zone to the photosphere can be strongly affected by convective turbulent flows. It has not yet been shown how crucial convective flows are for the formation of δ-spots. The aim of this study is to reveal the impact of convective flows in the convection zone on the formation and evolution of sunspot magnetic fields. We simulated the emergence and transport of magnetic flux tubes in the convection zone using radiative magnetohydrodynamics code R2D2. We carried out 93 simulations by allocating the twisted flux tubes to different positions in the convection zone. As a result, both δ-type and β-type magnetic distributions were reproduced only by the differences in the convective flows surrounding the flux tubes. The δ-spots were formed by the collision of positive and negative magnetic fluxes on the photosphere. The unipolar and bipolar rotations of the δ-spots were driven by magnetic twist and writhe, transporting magnetic helicity from the convection zone to the corona. We detected a strong correlation between the distribution of the non-potential magnetic field in the photosphere and the position of the downflow plume in the convection zone. The correlation could be detected 20–30 h before the flux emergence. The results suggest that high free energy regions in the photosphere can be predicted even before the magnetic flux appears in the photosphere by detecting the downflow profile in the convection zone.

DOI： 10.1093/mnras/stac2635

Web of Science

Scopus

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astrophysical Journal  

Solar magnetic fields comprise an 11 yr activity cycle, represented by the number of sunspots. The maintenance of such a solar magnetic field can be attributed to fluid motion in the convection zone, i.e., a dynamo. This study conducts the mean-field analyses of the global solar dynamo simulation presented by Hotta et al. (2016). Although the study succeeds in producing coherent large-scale magnetic fields at high Reynolds numbers, the detailed physics of the maintenance of these fields have not been fully understood. This study extracts the α tensor and the turbulent magnetic diffusivity tensor β through mean-field analyses. The turbulent magnetic diffusivity exhibits a significant decrease toward high Reynolds numbers. The decrease in the turbulent magnetic diffusivity suppresses the energy conversion of large-scale field to small-scale field. This implies that the decrease in the turbulent magnetic diffusivity contributes to the maintenance of a large-scale magnetic field at high Reynolds numbers. A significant downward turbulent pumping is observed; it is enhanced in the weak phase of the large-scale field. This study proposes a cyclic reversal process of a large-scale field, which is dominantly driven by the α effect and is possibly triggered by downward pumping.

DOI： 10.3847/1538-4357/ac7e43

Web of Science

Scopus

researchmap

その他リンク： https://iopscience.iop.org/article/10.3847/1538-4357/ac7e43/pdf

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astrophysical Journal  

We analyze the simulation result shown in Hotta & Kusano (2021) in which the solar-like differential rotation is reproduced. The Sun is rotating differentially with the fast equator and the slow pole. It is widely thought that the thermal convection maintains the differential rotation, but recent high-resolution simulations tend to fail to reproduce the fast equator. This fact is an aspect of one of the biggest problems in solar physics called the convective conundrum. Hotta & Kusano succeed in reproducing the solar-like differential rotation without using any manipulation with an unprecedentedly high-resolution simulation. In this study, we analyze the simulation data to understand the maintenance mechanism of the fast equator. Our analyses lead to conclusions that are summarized as follows. (1) The superequipatition magnetic field is generated by the compression, which can indirectly convert the massive internal energy to magnetic energy. (2) The efficient small-scale energy transport suppresses large-scale convection energy. (3) Non-Taylor-Proudman differential rotation is maintained by the entropy gradient caused by the anisotropic latitudinal energy transport enhanced by the magnetic field. (4) The fast equator is maintained by the meridional flow mainly caused by the Maxwell stress. The Maxwell stress itself also has a role in the angular momentum transport for the fast near-surface equator (we call it the P unching ball effect). The fast equator in the simulation is reproduced not due to the low Rossby number regime but due to the strong magnetic field. This study newly finds the role of the magnetic field in the maintenance of differential rotation.

DOI： 10.3847/1538-4357/ac7395

Web of Science

Scopus

researchmap

その他リンク： https://iopscience.iop.org/article/10.3847/1538-4357/ac7395/pdf

掲載種別：研究論文（学術雑誌）   出版者・発行元：Geophysical Research Letters  

Carbon-14 in tree rings have suggested there had been multiple extreme solar proton events (SPEs) in the past. While the largest events such as in 774–775 CE can be significantly detected by the typical precision of accelerator mass spectrometry, smaller but possibly more frequent events have been difficult to be detected. Thus, the frequency or any characteristics of such relatively smaller events are still largely unknown. In this paper, we report that large SPEs had occurred in 1261–1262, 1268–1269, and 1279–1280 CE before the onset of the Wolf minimum based on high-precision carbon-14 analyses. It is suggested that they had occurred at the maximum and the declining phase of solar cycles, and that they had occurred during the transition time of solar activity into a deep minimum. We propose that this episode may provide a unique opportunity to elucidate a potential interaction between the solar dynamo and extreme solar flares.

DOI： 10.1029/2021GL097201

Web of Science

Scopus

researchmap

その他リンク： https://onlinelibrary.wiley.com/doi/full-xml/10.1029/2021GL097201

担当区分：筆頭著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Nature Astronomy  

The Sun rotates differentially with a fast equator and slow pole1. Convection in the solar interior is thought to maintain the differential rotation. However, although many numerical simulations have been conducted to reproduce the solar differential rotation2–7, previous high-resolution calculations with solar parameters fall into the antisolar (fast-pole) differential rotation regime. Consequently, we still do not know the true reason why the Sun has a fast-rotating equator. While the construction of the fast equator requires a strong rotational influence on the convection, the previous calculations have not been able to achieve the situation without any manipulations. The problem is called the convective conundrum8. The convection and the differential rotation in numerical simulations were different from the observations. Here, we show that a high-resolution calculation succeeds in reproducing the solar-like differential rotation. Our calculations indicate that the strong magnetic field generated by a small-scale dynamo has a significant impact on thermal convection. The successful reproduction of the differential rotation, convection and magnetic field achieved in our calculation is an essential step to understanding the cause of the most basic nature of solar activity, specifically, the 11 yr cycle of sunspot activity.

DOI： 10.1038/s41550-021-01459-0

Web of Science

Scopus

researchmap

その他リンク： https://www.nature.com/articles/s41550-021-01459-0

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Earth, Planets and Space  

Although solar activity may significantly impact the global environment and socioeconomic systems, the mechanisms for solar eruptions and the subsequent processes have not yet been fully understood. Thus, modern society supported by advanced information systems is at risk from severe space weather disturbances. Project for solar–terrestrial environment prediction (PSTEP) was launched to improve this situation through synergy between basic science research and operational forecast. The PSTEP is a nationwide research collaboration in Japan and was conducted from April 2015 to March 2020, supported by a Grant-in-Aid for Scientific Research on Innovative Areas from the Ministry of Education, Culture, Sports, Science and Technology of Japan. By this project, we sought to answer the fundamental questions concerning the solar–terrestrial environment and aimed to build a next-generation space weather forecast system to prepare for severe space weather disasters. The PSTEP consists of four research groups and proposal-based research units. It has made a significant progress in space weather research and operational forecasts, publishing over 500 refereed journal papers and organizing four international symposiums, various workshops and seminars, and summer school for graduate students at Rikubetsu in 2017. This paper is a summary report of the PSTEP and describes the major research achievements it produced.[Figure not available: see fulltext.]

DOI： 10.1186/s40623-021-01486-1

Web of Science

Scopus

researchmap

その他リンク： https://link.springer.com/article/10.1186/s40623-021-01486-1/fulltext.html

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Monthly Notices of the Royal Astronomical Society  

We perform radiative magnetohydrodynamic calculations for the solar-quiet region to investigate the dependence of statistical flow on magnetic properties and the three-dimensional structure of magnetic patches in the presence of large-scale flow that mimics differential rotation. It has been confirmed that strong magnetic field patches move faster in the longitudinal direction at the solar surface. Consequently, strong magnetic patches penetrate deeper into the solar interior. The motion of the deep-rooted magnetic patches is influenced by the faster differential rotation in the deeper layer. In this study, we perform realistic radiative magnetohydrodynamic calculations using R2D2 code to validate that stronger patches have deeper roots. We also add large-scale flow to mimic the differential rotation. The magnetic patches are automatically detected and tracked, and we evaluate the depth of 30 000 magnetic patches. The velocities of 2.9 million magnetic patches are then measured at the photosphere. We obtain the dependence of these values on the magnetic properties, such as field strength and flux. Our results confirm that strong magnetic patches tend to show deeper roots and faster movement, and we compare our results with observations using the point spread function of instruments at the Hinode and Solar Dynamics Observatory (SDO). Our result is quantitatively consistent with previous observational results of the SDO.

DOI： 10.1093/mnras/stab710

Web of Science

Scopus

researchmap

その他リンク： http://academic.oup.com/mnras/article-pdf/503/3/3610/36855183/stab710.pdf

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Scientific Reports  

The Sun exhibits centennial-scale activity variations and sometimes encounters grand solar minimum when solar activity becomes extremely weak and sunspots disappear for several decades. Such an extreme weakening of solar activity could cause severe climate, causing massive reductions in crop yields in some regions. During the past decade, the Sun’s activity has tended to decline, raising concerns that the Sun might be heading for the next grand minimum. However, we still have an underdeveloped understanding of solar dynamo mechanisms and hence precise prediction of near-future solar activity is not attained. Here we show that the 11-year solar cycles were significantly lengthened before the onset of the Maunder Minimum (1645–1715 CE) based on unprecedentedly high-precision data of carbon-14 content in tree rings. It implies that flow speed in the convection zone is an essential parameter to determine long-term solar activity variations. We find that a 16 year-long cycle had occurred three solar cycles before the onset of prolonged sunspot disappearance, suggesting a longer-than-expected preparatory period for the grand minimum. As the Sun has shown a tendency of cycle lengthening since Solar Cycle 23 (1996–2008 CE), the behavior of Solar Cycle 25 can be critically important to the later solar activity.

DOI： 10.1038/s41598-021-84830-5

Web of Science

Scopus

researchmap

その他リンク： http://www.nature.com/articles/s41598-021-84830-5

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Monthly Notices of the Royal Astronomical Society  

We perform a series of radiative magnetohydrodynamic simulations to understand the amplification mechanism of the exceptionally strong horizontal magnetic field in delta-type sunspots. In the simulations, we succeed in reproducing the delta-type sunspot and resulting strong magnetic field exceeding 6000 G in a light bridge between the positive and negative polarities. Our conclusions in this study are summarized as follows: (1) The essential amplification mechanism of the strong horizontal magnetic field is the shear motion caused by the rotation of two spots. (2) The strong horizontal magnetic field remains the force-free state. (3) The peak strength of the magnetic fields does not depend on the spatial resolution, top boundary condition, or Alfvén speed limit. The origin of the rotating motion is rooted in the deep convection zone. Therefore, the magnetic field in the delta-spot light bridge can be amplified to the superequipartition values in the photosphere.

DOI： 10.1093/mnras/staa2529

Web of Science

Scopus

researchmap

その他リンク： http://academic.oup.com/mnras/article-pdf/498/2/2925/33781176/staa2529.pdf

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Science Advances  

Convection in the Sun's outer envelope generates turbulence and drives differential rotation, meridional circulation, and the global magnetic cycle. We develop a greater understanding of these processes by contrasting observations with simulations of global convection. These comparisons also enhance our comprehension of the physics of distant Sun-like stars. Here, we infer toroidal flow power as a function of wave number, frequency, and depth in the solar interior through helioseismic analyses of space-based observations. The inferred flows grow with spatial wave number and temporal frequency and are confined to low latitudes, supporting the argument that rotation induces systematic differences between the poles and equator. In contrast, the simulations used here show the opposite trends-power diminishing with increasing wave number and frequency while flow amplitudes become weakest at low latitudes. These differences highlight gaps in our understanding of solar convection and point to challenges ahead.

DOI： 10.1126/sciadv.aba9639

Web of Science

Scopus

researchmap

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Monthly Notices of the Royal Astronomical Society  

We investigate the rising flux tube and the formation of sunspots in an unprecedentedly deep computational domain that covers the whole convection zone with a radiative magnetohydrodynamics simulation. Previous calculations had shallow computational boxes (<30 Mm) and convection zones at a depth of 200 Mm. By using our new numerical code Radition and RSST for Deep Dynamics (R2D2), we succeed in covering the whole convection zone and reproduce the formation of the sunspot from a simple horizontal flux tube because of the turbulent thermal convection. The main findings are as follows. (1) The rising speed of the flux tube is larger than the upward convection velocity because of the low density caused by the magnetic pressure and the suppression of the mixing. (2) The rising speed of the flux tube exceeds 250 ms-1 at a depth of 18Mm, while we do not see any clear evidence of the divergent flow 3 h before the emergence at the solar surface. (3) Initially, the root of the flux tube is filled with the downflows, and then the upflow fills the centre of the flux tube during the formation of the sunspot. (4) The essential mechanisms for the formation of the sunspot are the coherent inflow and the turbulent transport. (5) The low-temperature region is extended to a depth of at least 40 Mm in the matured sunspot, with the high-temperature region in the centre of the flux tube. Some of the findings indicate the importance of the deep computational domain for the flux emergence simulations.

DOI： 10.1093/mnras/staa844

Web of Science

Scopus

researchmap

その他リンク： http://academic.oup.com/mnras/article-pdf/494/2/2523/33113724/staa844.pdf

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astrophysical Journal Letters  

https://ui.adsabs.harvard.edu/abs/2019ApJ...886L..21T/abstract

DOI： 10.3847/2041-8213/ab55e7

Web of Science

Scopus

researchmap

その他リンク： https://arxiv.org/abs/1911.03909

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astrophysical Journal  

The morphological asymmetry of leading and following sunspots is a well-known characteristic of the solar surface. In the context of the large-scale evolution of the surface magnetic field, the asymmetry has been assumed to have only a negligible effect. Using the surface flux transport (SFT) model, we show that the morphological asymmetry of leading and following sunspots has a significant impact on the evolution of the large-scale magnetic field on the solar surface. By evaluating the effect of the morphological asymmetry of each bipolar magnetic region (BMR), we observe that the introduction of asymmetry to the BMR model significantly reduces the contribution to the polar magnetic field, especially for large and high-latitude BMRs. Strongly asymmetric BMRs can even reverse regular polar field formation. The SFT simulations based on the observed sunspot record show that the introduction of morphological asymmetry reduces the root-mean-square difference from the observed axial dipole strength by 30%-40%. These results indicate that the morphological asymmetry of leading and following sunspots has a significant effect on the solar cycle prediction.

DOI： 10.3847/1538-4357/ab3b04

Web of Science

Scopus

researchmap

その他リンク： https://iopscience.iop.org/article/10.3847/1538-4357/ab3b04

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astronomy and Astrophysics  

Context. The reduced speed of sound technique (RSST) has been used for efficient simulation of low Mach number flows in solar and stellar convection zones. The basic RSST equations are hyperbolic and are suitable for parallel computation by domain decomposition. The application of RSST is limited to cases in which density perturbations are much smaller than the background density. In addition, nonconservative variables are required to be evolved using this method, which is not suitable in cases where discontinuities such as shock waves coexist in a single numerical domain. Aims. In this study, we suggest a new semiconservative formulation of the RSST that can be applied to low Mach number flows with large density variations. Methods. We derive the wave speed of the original and newly suggested methods to clarify that these methods can reduce the speed of sound without affecting the entropy wave. The equations are implemented using the finite volume method. Several numerical tests are carried out to verify the suggested methods. Results. The analysis and numerical results show that the original RSST is not applicable when mass density variations are large. In contrast, the newly suggested methods are found to be efficient in such cases. We also suggest variants of the RSST that conserve momentum in the machine precision. The newly suggested variants are formulated as semiconservative equations, which reduce to the conservative form of the Euler equations when the speed of sound is not reduced. This property is advantageous when both high and low Mach number regions are included in the numerical domain. Conclusions. The newly suggested forms of RSST can be applied to a wider range of low Mach number flows.

DOI： 10.1051/0004-6361/201834031

Web of Science

Scopus

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Science Advances  

The solar convection zone is filled with turbulent convection in highly stratified plasma. Several theoretical and observational studies suggest that the numerical calculations overestimate the convection velocity. Since all deep convection zone calculations exclude the solar surface due to substantial temporal and spatial scale separations, the solar surface, which drives the thermal convection with efficient radiative cooling, has been thought to be the key to solve this discrepancy. Thanks to the recent development in massive supercomputers, we are successful in performing the comprehensive calculation covering the whole solar convection zone. We compare the results with and without the solar surface in the local domain and without the surface in the full sphere. The calculations do not include the rotation and the magnetic field. The surface region has an unexpectedly weak influence on the deep convection zone. We find that just including the solar surface cannot solve the problem.

DOI： 10.1126/sciadv.aau2307

Web of Science

Scopus

PubMed

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Science  

The differentially rotating outer layers of stars are thought to play a role in driving their magnetic activity, but the underlying mechanisms that generate and sustain differential rotation are poorly understood. We report the measurement using asteroseismology of latitudinal differential rotation in the convection zones of 40 Sun-like stars. For the most significant detections, the stars’ equators rotate approximately twice as fast as their midlatitudes. The latitudinal shear inferred from asteroseismology is much larger than predictions from numerical simulations.

DOI： 10.1126/science.aao6571

Web of Science

Scopus

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Publications of the Astronomical Society of Japan  

Sunspot observations with telescopes in the 18th century were carried out in Japan as well as elsewhere. One of these sunspot observations is recorded in an account called Sansaizusetsu narabini Kansei irai Jissoku Zusetsu (Charts of Three Worlds and Diagrams of Actual Observations since Kansei Era). We have analyzed manuscripts of this account to show a total of 15 sunspot drawings during 1749-1750. These observations are considered to be carried out by contemporary official astronomers in Japan, with telescopes covered by zongurasus (< zonglas in Dutch, corresponding to “sunglass” in English). We counted their group number of sunspots to locate them in long-term solar activity and show that their observations were situated near the solar maximum in 1750. We also computed their locations and areas, while we have to admit differences of the variant manuscripts with one another. These observational records show the spread of sunspot observations not only in Europe, but also in Japan, and hence may contribute to crosscheck, or possibly to improve the known sunspot indices.

DOI： 10.1093/pasj/psy066

Web of Science

Scopus

researchmap

出版者・発行元：International Conference for High Performance Computing, Networking, Storage and Analysis, SC  

Stencil computation has many applications in science and engineering, thus many optimization techniques such as temporal blocking have been developed. They are, however, rarely used in real-world applications, since a large amount of careful programming is required for even the simplest of stencils. We introduce Formura, a domain specific language that provides easy access to optimized stencil computations. Higher-order integration schemes can be defined using mathematical notations. Formura generates C code with MPI calls and performs autotuning. Hence its performance is portable to most distributed-memory computers. We show the scientific applicability of Formura by performing magnetohydrodynamics (MHD) and belowground biology simulations. Ability to reach bytes-per-flops ratio only attainable by temporal blocking is demonstrated. We also demonstrate scaling up to the full nodes of the K computer, with 1.184 Pflops, 11.62% floating-pointoperation efficiency, and 31.26% memory throughput efficiency.

DOI： 10.1109/SC.2016.2

Scopus

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Science  

The 11-year solarmagnetic cycle shows a high degree of coherence in spite of the turbulent nature of the solar convection zone. It has been found in recent high-resolutionmagnetohydrodynamics simulations that the maintenance of a large-scale coherent magnetic field is difficult with small viscosity and magnetic diffusivity (≲1012 square centimenters per second).We reproduced previous findings that indicate a reduction of the energy in the large-scalemagnetic field for lower diffusivities and demonstrate the recovery of the global-scalemagnetic field using unprecedentedly high resolution. We found an efficient small-scale dynamo that suppresses small-scale flows, which mimics the properties of large diffusivity. As a result, the global-scale magnetic field is maintained even in the regime of small diffusivities - that is, large Reynolds numbers.

DOI： 10.1126/science.aad1893

Web of Science

Scopus

researchmap

Web of Science

出版者・発行元：FHPC 2016 - Proceedings of the 5th International Workshop on Functional High-Performance Computing, co-located with ICFP 2016  

Programming in HPC is a tedious work. Therefore functional programming languages that generate HPC programs have been proposed. However, they are not widely used by application scientists, because of learning barrier, and lack of demonstrated application performance. We have designed Formura which adopts application-friendly features such as typed rational array indices. Formura users can describe mathematical concepts such as operation over derivative operators using functional programming. Formura allows intuitive expression over array elements while ensuring the program is a stencil computation, so that state-of-the-art stencil optimization techniques such as temporal blocking is always applied to Formura-generated program. We demonstrate the usefulness of Formura by implementing a preliminary below-ground biology simulation. Optimized C-code are generated from 672 bytes of Formura program. The simulation was executed on the full nodes of the K computer, with 1.184 Pflops, 11.62% floating-point-instruction efficiency, and 31.26% memory throughput efficiency.

DOI： 10.1145/2975991.2975994

Web of Science

Scopus

記述言語：英語   掲載種別：研究論文（国際会議プロシーディングス）   出版者・発行元：Space Science Reviews  

We discuss recent observational, theoretical and numerical progress made in understanding the solar global magnetism and its short and long term variability. We discuss the physical process thought to be at the origin of the solar magnetic field and its 22-yr cycle, namely dynamo action, and the nonlinear interplay between convection, rotation, radiation and magnetic field, yielding modulations of the solar constant or of the large scale flows such as the torsional oscillations. We also discuss the role of the field parity and dynamo families in explaining the complex multipolar structure of the solar global magnetic field. We then present some key MHD processes acting in the deep radiative interior and discuss the probable topology of a primordial field there. Finally we summarize how helioseismology has contributed to these recent advances and how it could contribute to resolving current unsolved problems in solar global dynamics and magnetism.

DOI： 10.1007/s11214-013-0028-0

Web of Science

Scopus

researchmap

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astrophysical Journal  

We investigate small-scale dynamo action in the solar convection zone through a series of high-resolution MHD simulations in a local Cartesian domain with 1 R⊙ (solar radius) of horizontal extent and a radial extent from 0.715 to 0.96 R⊙. The dependence of the solution on resolution and diffusivity is studied. For a grid spacing of less than 350 km, the rms magnetic field strength near the base of the convection zone reaches 95% of the equipartition field strength (i.e., magnetic and kinetic energy are comparable). For these solutions the Lorentz force feedback on the convection velocity is found to be significant. The velocity near the base of the convection zone is reduced to 50% of the hydrodynamic one. In spite of the significant decrease of the convection velocity, the reduction in the enthalpy flux is relatively small, since the magnetic field also suppresses the horizontal mixing of the entropy between up- and downflow regions. This effect increases the amplitude of the entropy perturbation and makes convective energy transport more efficient. We discuss potential implications of these results for solar global convection and dynamo simulations.

DOI： 10.1088/0004-637X/803/1/42

Web of Science

Scopus

researchmap

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astrophysical Journal  

We present a high-resolution, highly stratified numerical simulation of rotating thermal convection in a spherical shell. Our aim is to study in detail the processes that can maintain a near surface shear layer (NSSL) as inferred from helioseismology. Using the reduced speed of sound technique, we can extend our global convection simulation to 0.99 R⊙ and include, near the top of our domain, small-scale convection with short timescales that is only weakly influenced by rotation. We find the formation of an NSSL preferentially in high latitudes in the depth range of r = 0.95-0.975 R⊙. The maintenance mechanisms are summarized as follows. Convection under the weak influence of rotation leads to Reynolds stresses that transport angular momentum radially inward in all latitudes. This leads to the formation of a strong poleward-directed meridional flow and an NSSL, which is balanced in the meridional plane by forces resulting from the <υ′rυ′θ> correlation of turbulent velocities. The origin of the required correlations depends to some degree on latitude. In high latitudes, a positive correlation <υ′rυ′θ> is induced in the NSSL by the poleward meridional flow whose amplitude increases with the radius, while a negative correlation is generated by the Coriolis force in bulk of the convection zone. In low latitudes, a positive correlation <υ′rυ′θ> results from rotationally aligned convection cells ("banana cells"). The force caused by these Reynolds stresses is in balance with the Coriolis force in the NSSL.

DOI： 10.1088/0004-637X/798/1/51

Web of Science

Scopus

researchmap

DOI： 10.1007/978-4-431-55399-1_3

Web of Science

DOI： 10.1007/978-4-431-55399-1_4

Web of Science

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（国際会議プロシーディングス）   出版者・発行元：ASTRONOMICAL SOC PACIFIC  

We investigate the solar differential rotation with special interest for the near surface shear layer (NSSL) in a high-resolution hydrodynamic numerical calculation. The sun is rotating differentially. Helioseismology has revealed the detailed structure of the solar differential rotation. One of the most important features is the NSSL. It is thought that the solar differential rotation is maintained by the turbulent thermal convection. In the NSSL convection time scales are short, leading to a regime with weak influence of rotation on convection. In order to reproduce the NSSL by the numerical calculations, we must use a large number of grids and integrate a large number of time steps for covering the broad spatial and temporal scales. This requirements for the NSSL is achieved using our recent efficient numerical method. In the calculation, the global scale and the 10 Mm-scale convection is established simultaneously. Then the solar like NSSL is partially reproduced. Around the NSSL, the convection transports the angular momentum radially inward and generates the poleward meridional flow. The small scale convection acts as the turbulent viscosity on the meridional flow. The turbulent viscous stress balances with the Coriolis force in the NSSL.

Web of Science

researchmap

DOI： 10.1007/978-4-431-55399-1_2

Web of Science

DOI： 10.1007/978-4-431-55399-1_1

Web of Science

DOI： 10.1007/978-4-431-55399-1_5

Web of Science

出版者・発行元：Astrophysical Journal  

We carry out non-rotating high-resolution calculations of the solar global convection, which resolve convective scales of less than 10 Mm. To cope with the low Mach number conditions in the lower convection zone, we use the reduced speed of sound technique (RSST), which is simple to implement and requires only local communication in the parallel computation. In addition, the RSST allows us to expand the computational domain upward to about 0.99 R ⊙, as it can also handle compressible flows. Using this approach, we study the solar convection zone on the global scale, including small-scale near-surface convection. In particular, we investigate the influence of the top boundary condition on the convective structure throughout the convection zone as well as on small-scale dynamo action. Our main conclusions are as follows. (1) The small-scale downflows generated in the near-surface layer penetrate into deeper layers to some extent and excite small-scale turbulence in the region >0.9 R ⊙, where R ⊙is the solar radius. (2) In the deeper convection zone (<0.9 R ⊙), the convection is not influenced by the location of the upper boundary. (3) Using a large eddy simulation approach, we can achieve small-scale dynamo action and maintain a field of about 0.15B eq-0.25B eq throughout the convection zone, where B eq is the equipartition magnetic field to the kinetic energy. (4) The overall dynamo efficiency varies significantly in the convection zone as a consequence of the downward directed Poynting flux and the depth variation of the intrinsic convective scales. © 2014. The American Astronomical Society. All rights reserved.

DOI： 10.1088/0004-637X/786/1/24

Scopus

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）  

researchmap

記述言語：日本語   出版者・発行元：一般社団法人 日本物理学会  

DOI： 10.11316/jpsgaiyo.68.1.2.0_257_2

CiNii Research

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astronomy and Astrophysics  

Using two-dimensional magnetohydrodynamics calculations, we investigate a twist generation mechanism on a magnetic flux tube at the base of the solar convection zone based on the idea of Choudhuri (2003, Sol. Phys., 215, 31) in which a toroidal magnetic field is wrapped by a surrounding mean poloidal field. During generation of the twist, the flux tube follows four phases. (1) It quickly splits into two parts with vortex motions rolling up the poloidal magnetic field. (2) Owing to the physical mechanism similar to that of the magneto-rotational instability, the rolled-up poloidal field is bent and amplified. (3) The magnetic tension of the disturbed poloidal magnetic field reduces the vorticity, and the lifting force caused by vortical motion decreases. (4) The flux tube gets twisted and begins to rise again without splitting. Investigation of these processes is significant because it shows that a flux tube without any initial twist can rise to the surface in relatively weak poloidal fields. © 2012 ESO.

DOI： 10.1051/0004-6361/201220108

Web of Science

Scopus

researchmap

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astrophysical Journal Letters  

We investigate an intensification mechanism for the magnetic field near the base of the solar convection zone that does not rely on differential rotation. Such mechanism in addition to differential rotation has been suggested by studies of flux emergence, which typically require field strength in excess of those provided by differential rotation alone. We study here a process in which potential energy of the superadiabatically stratified convection zone is converted into magnetic energy. This mechanism, known as the "explosion of magnetic flux tubes," has been previously studied in thin flux tube approximation as well as two-dimensional magnetohydrodynamic (MHD) simulations; here we expand the investigation to three-dimensional MHD simulations. Our main result is that enough intensification can be achieved in a three-dimensional magnetic flux sheet as long as the spatial scale of the imposed perturbation normal to the magnetic field is sufficiently large. When this spatial scale is small, the flux sheet tends to rise toward the surface, resulting in a significant decrease of the magnetic field amplification. © 2012. The American Astronomical Society. All rights reserved.

DOI： 10.1088/2041-8205/759/1/L24

Web of Science

Scopus

researchmap

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astrophysical Journal Letters  

We investigate the value of horizontal turbulent diffusivity η by numerical calculation of thermal convection. In this study, we introduce a new method whereby the turbulent diffusivity is estimated by monitoring the time development of the passive scalar, which is initially distributed in a given Gaussian function with a spatial scale d 0. Our conclusions are as follows: (1) assuming the relation η = L c v rms/3, where v rms is the root-mean-square (rms) velocity, the characteristic length L c is restricted by the shortest one among the pressure (density) scale height and the region depth. (2) The value of turbulent diffusivity becomes greater with the larger initial distribution scale d 0. (3) The approximation of turbulent diffusion holds better when the ratio of the initial distribution scale d 0 to the characteristic length L c is larger. © 2012. The American Astronomical Society. All rights reserved..

DOI： 10.1088/2041-8205/751/1/L9

Web of Science

Scopus

researchmap

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astronomy and Astrophysics  

Context. The anelastic approximation is often adopted in numerical calculations with low Mach numbers, such as those including stellar internal convection. This approximation requires so-called frequent global communication, because of an elliptic partial differential equation. Frequent global communication is, however, negative factor for the parallel computing performed with a large number of CPUs. Aims. We test the validity of a method that artificially reduces the speed of sound for the compressible fluid equations in the context of stellar internal convection. This reduction in the speed of sound leads to longer time steps despite the low Mach number, while the numerical scheme remains fully explicit and the mathematical system is hyperbolic, thus does not require frequent global communication. Methods. Two-and three-dimensional compressible hydrodynamic equations are solved numerically. Some statistical quantities of solutions computed with different effective Mach numbers (owing to the reduction in the speed of sound) are compared to test the validity of our approach. Results. Numerical simulations with artificially reduced speed of sound are a valid approach as long as the effective Mach number (based on the lower speed of sound) remains less than 0.7. © 2012 ESO.

DOI： 10.1051/0004-6361/201118268

Web of Science

Scopus

researchmap

記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Nature Communications  

The accretion of small bodies in the Solar System is a fundamental process that was followed by planet formation. Chronological information of meteorites can constrain when asteroids formed. Secondary carbonates show extremely old 53Mn-53Cr radiometric ages, indicating that some hydrous asteroids accreted rapidly. However, previous studies have failed to define accurate Mn/Cr ratios; hence, these old ages could be artefacts. Here we develop a new method for accurate Mn/Cr determination, and report a reliable age of 4,563.4+0.4/-0.5 million years ago for carbonates in carbonaceous chondrites. We find that these carbonates have identical ages, which are younger than those previously estimated. This result suggests the late onset of aqueous activities in the Solar System. The young carbonate age cannot be explained if the parent asteroid accreted within 3 million years after the birth of the Solar System. Thus, we conclude that hydrous asteroids accreted later than differentiated and metamorphosed asteroids. © 2012 Macmillan Publishers Limited. All rights reserved.

DOI： 10.1038/ncomms1635

Web of Science

Scopus

researchmap

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（国際会議プロシーディングス）   出版者・発行元：ASTRONOMICAL SOC PACIFIC  

The flux-transport dynamo model for the solar sunspot cycle is revised and is demonstrated by using the axisymmetric kinematic simulations. The flux-transport dynamo has succeeded to explain the general cyclic behaviors of the sunspots especially in the gradual shift of the sunspot toward the equator and the poleward migration of the surface magnetic field. It has been known, however, that previous models failed to avoid the strong polar surface field and the strong toroidal field at the base in the high latitude, both of which are not consistent with observations. With an additional intense diffusivity profile near the surface two problematic features can be avoided. The surface poloidal field generated by the a effect is transported down to the base of the convection zone not by the meridional flow but by the surface diffusion mainly in the mid-latitude. This prevents the concentration of the polar surface field and the amplification of the toroidal field at the high latitude. The condition to obtain the proper magnetic field strength near the pole is eta(surf)/u(0) &gt; 2 x 10(9)cm, where eta(surf) and u(0) are the surface diffusivity and the meridional flow speed, respectively. We also do some parameter studies to ensure the importance of the surface strong diffusivity.

Web of Science

researchmap

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astrophysical Journal  

We investigate differential rotation in rapidly rotating solar-type stars by means of an axisymmetric mean field model that was previously applied to the Sun. This allows us to calculate the latitudinal entropy gradient with a reasonable physical basis. Our conclusions are as follows. (1) Differential rotation approaches the Taylor-Proudman state when stellar rotation is faster than solar rotation. (2) Entropy gradient generated by the attached subadiabatic layer beneath the convection zone becomes relatively small with a large stellar angular velocity. (3) Turbulent viscosity and turbulent angular momentum transport determine the spatial difference of angular velocity ΔΩ. (4) The results of our mean field model can explain observations of stellar differential rotation. © 2011. The American Astronomical Society. All rights reserved.

DOI： 10.1088/0004-637X/740/1/12

Web of Science

Scopus

researchmap

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astrophysical Journal Letters  

We investigated the dependence of the solar magnetic parity between the hemispheres on two important parameters, the turbulent diffusivity and the meridional flow, by means of axisymmetric kinematic dynamo simulations based on the flux-transport dynamo model. It is known that the coupling of the magnetic field between hemispheres due to turbulent diffusivity is an important factor for the solar parity issue, but the detailed criterion for the generation of the dipole field has not been investigated. Our conclusions are as follows. (1) The stronger diffusivity near the surface is more likely to cause the magnetic field to be a dipole. (2) The thinner layer of the strong diffusivity near the surface is also more apt to generate a dipolar magnetic field. (3) The faster meridional flow is more prone to cause the magnetic field to be a quadrupole, i.e., symmetric about the equator. These results show that turbulent diffusivity and meridional flow are crucial for the configuration of the solar global magnetic field. © 2010. The American Astronomical Society. All rights reserved.

DOI： 10.1088/2041-8205/714/2/L308

Web of Science

Scopus

researchmap

担当区分：筆頭著者,　責任著者   記述言語：英語   掲載種別：研究論文（学術雑誌）   出版者・発行元：Astrophysical Journal  

A revision to the flux-transport dynamo model for the solar sunspot cycle is proposed and is demonstrated by using the axisymmetric kinematic simulations. The flux-transport dynamo has succeeded to explain the general cyclic behaviors of the sunspots. It has been known, however, that previous models failed to avoid the strong polar surface field and the strong toroidal field at the base in the high latitude, both of which are not consistent with observations. We propose a new regime of the flux-transport dynamo model by assuming an additional intense diffusivity profile near the surface. The surface poloidal field generated by the α effect is transported down to the base of the convection zone not by the meridional flow but by the surface diffusion mainly in the mid-latitude. With a moderate α quenching, this prevents the concentration of the polar surface field and the amplification of the toroidal field at the high latitude. The condition to obtain the proper magnetic field strength near the pole is ηsurf /u0 > 2 × 109 cm, where ηsurf and u0 are the surface diffusivity and the meridional flow speed, respectively. We also do some parameter studies to ensure the importance of the surface strong diffusivity. In addition, the dependence of the cycle period on free parameters, the speed of meridional flow and the surface diffusivity, is investigated.

DOI： 10.1088/0004-637X/709/2/1009

Web of Science

Scopus

researchmap

開催年月日： 2024年2月

記述言語：日本語   会議種別：ポスター発表  

開催年月日： 2023年12月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

開催年月日： 2023年12月

記述言語：日本語   会議種別：口頭発表（一般）  

開催年月日： 2023年12月

記述言語：日本語   会議種別：公開講演，セミナー，チュートリアル，講習，講義等  

開催年月日： 2023年12月

記述言語：日本語   会議種別：公開講演，セミナー，チュートリアル，講習，講義等  

開催年月日： 2023年11月

記述言語：英語   会議種別：口頭発表（招待・特別）  

開催年月日： 2023年9月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

開催年月日： 2023年9月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

開催年月日： 2023年9月

記述言語：日本語   会議種別：口頭発表（一般）  

開催年月日： 2023年8月

記述言語：英語   会議種別：口頭発表（招待・特別）  

開催年月日： 2023年6月

記述言語：日本語   会議種別：公開講演，セミナー，チュートリアル，講習，講義等  

開催年月日： 2023年6月

記述言語：英語   会議種別：口頭発表（一般）  

開催年月日： 2022年11月

記述言語：英語   会議種別：口頭発表（招待・特別）  

開催年月日： 2022年9月

記述言語：英語   会議種別：口頭発表（招待・特別）  

開催年月日： 2022年8月

記述言語：英語   会議種別：口頭発表（招待・特別）  

開催年月日： 2022年8月

記述言語：英語   会議種別：口頭発表（招待・特別）  

開催年月日： 2022年2月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

開催年月日： 2022年1月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

開催年月日： 2021年12月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

開催年月日： 2021年11月

記述言語：英語   会議種別：口頭発表（基調）  

開催年月日： 2021年11月

記述言語：英語   会議種別：口頭発表（招待・特別）  

開催年月日： 2021年9月

researchmap

開催年月日： 2021年8月

記述言語：日本語  

researchmap

開催年月日： 2021年6月

researchmap

開催年月日： 2021年3月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2021年3月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2020年8月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2020年2月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2019年12月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2019年9月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2019年8月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2019年7月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2018年10月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2018年10月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2018年8月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2018年3月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2018年2月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2017年12月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2017年11月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2017年11月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2017年10月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2017年9月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2017年3月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2017年3月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2017年2月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2016年11月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2016年10月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2016年9月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2016年8月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2016年7月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2016年6月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2016年3月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2016年3月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2016年2月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2015年12月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2015年12月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2015年11月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2015年8月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2015年6月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2015年3月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2014年11月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2014年6月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2013年12月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2013年11月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2013年6月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2013年2月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2011年12月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2011年7月

記述言語：英語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2010年9月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap

開催年月日： 2010年1月

記述言語：日本語   会議種別：口頭発表（招待・特別）  

researchmap