1 Introduction
Steady-state plasma flow and magnetic field can be regarded as a realization of potential field in the planetary magnetosheath region when the vorticity and the electric current are ignored. In such a case, the potential is obtained by solving the Laplace equation, which was elegantly and analytically solved by for a parabolic shape of magnetosheath (hereafter KF). The KF potential was further extended to the stream function in the magnetosheath by . The KF solution made a series of breakthroughs in magnetosheath research. One of the most successful applications is the ability to track the plasma parcel along the streamline in the modeled magnetosheath. The tracking method was extensively used to observationally study the mirror mode growth
The KF potential is obtained using the assumption that the planetary bow shock and magnetopause have a parabolic shape sharing the same focal point. Empirical models of the bow shock and magnetopause (fitted to the spacecraft data), on the other hand, are not necessarily parabolically or co-focally shaped. For example, the empirical Earth bow shock model by and has a parabolic shape, but the focal point differs from that of the KF solution. The empirical magnetopause model by applies a power-law scaling to the parabolic shape such that the magnetic field lines appear stretched in the tail region. The gap between the KF parabolic magnetosheath and the empirical magnetosheath needs to be filled when applying the KF potential in the empirical magnetosheath.
Naively speaking, one wishes to find a conformal mapping (angle-preserving mapping) from the KF parabolic magnetosheath onto a non-parabolic empirical magnetosheath shape such as the analytic extension of magnetopause shape . However, no general mathematical algorithm is known so far to obtain the conformal mapping when the spatial domain is not properly bounded. The problem lies in the fact that the magnetosheath is bounded only by two sides, i.e., the standing shock and the magnetopause in the radial direction to the planet, but it is not bounded along the flow in the tail region. The algorithms of numerical conformal mapping are so far proposed for spatially bounded domains or domains with a closed shape of internal boundaries .
Here we address the question, how can we utilize the KF magnetosheath model by transforming into a more general, empirical, non-parabolic magnetosheath geometry? To achieve the goal, we develop a mapping method which provides a semi-analytic estimate of steady-state flow velocity and magnetic field in the empirical magnetosheath domain. Our scalar-potential mapping method is computationally inexpensive by using the analytic expression as much as possible, and it is applicable to various planetary magnetosheath domains.
This work is organized in the following fashion. After reviewing the magnetosheath model constructed by (Sect. ) and discussing different mapping methods (Sect. ), the detailed procedure of the magnetopause-normal mapping is presented (Sect. ) and discussed (Sect. ), with concluding remarks (Sect. ).
2 Revisiting the magnetosheath scalar potential
2.1 Parabolic coordinates
In the KF parabolic coordinates, the shell variable (iso-contour lines enveloping the magnetosphere) and the connector (iso-contour lines connecting from the bow shock to the magnetopause) play an important role in computing the flow velocity and magnetic field in the magnetosheath. These variables are explicitly evaluated using Cartesian coordinates and the radial distance from the focal point as where is the distance to the focus at :
3 The focus is along the axis and is defined as 4 , , and are the Cartesian representation of the KF magnetosheath model (i.e., with the pre-fixed bow shock and magnetopause shapes) obtained by projecting the position vector onto the unit vectors , , and : To complete the variable set for computing the potentials and the stream function, the azimuthal angle is introduced as 8
2.2 Velocity potentialIn the frame of potential field theory, the flow velocity is obtained either from the velocity potential (a scalar potential) or from the stream function (also a scalar potential) as
9
The symbol is the unit vector in the azimuthal direction around the symmetry axis (Sun-to-planet direction). and obtained the analytic expression of the velocity potential using the shell variable and the connector variable . 10 where is the upstream flow velocity, is the shell variable at the magnetopause, is the shell variable at the bow shock, is the shell variable, is the connector variable, and is a free parameter (integration constant) which is set to zero without loss of generality. The boundary shell values and contain the information on the stand-off distances ( and ) in the subsolar region, and they are defined by as
2.3 Stream functiontransformed the KF potential and obtained analytically the stream function as a function of the shell variable and the connector variable:
13
Hereafter, one may set so that the velocity potential is normalized to the upstream velocity. Iso-contour lines of the stream function represent the streamline.
2.4 Magnetic scalar potentialThe magnetic field in the magnetosheath is derived from the scalar potential in the same fashion as the flow velocity; that is,
14
The magnetic potential is a function of the shell variable and the connector : 15 where is the sunward component of the upstream magnetic field (corresponding to the GSE-X in near-Earth space), and and are two components of the upstream magnetic field perpendicular to the x direction. is the azimuthal angle of the position around the symmetry axis (the direction is given by the angle ). The integration constant is chosen as . The magnetic potential cannot be further transformed into the form of a stream function since the magnetic field distribution is essentially three-dimensional in the magnetosheath.
3 Mapping method comparison3.1 Mapping problem
Our task is to find the shell variable and the connector in the empirical magnetosheath by finding a suitable mapping of the position vector from the empirical magnetosheath (denoted by ) onto the KF parabolic system (denoted by ). The evaluated and are then readily used to obtain the scalar potentials and the stream function. The flow velocity and the magnetic field in the empirical magnetosheath are obtained by computing the gradient of the respective potential.
A practically useful mapping procedure to utilize the KF potential is proposed by by using the radial direction as a reference. While the radial mapping can reasonably (i.e., with a relatively high accuracy) transform the KF potential into the empirical magnetosheath domain on the dayside, the mapping quality becomes degraded in the flank region due to the conversion effect associated with the non-orthogonal grid construction. Our approach differs from the radial mapping by using the magnetopause-normal direction as a reference. The azimuthal coordinate is still orthogonal to the and coordinates. We briefly compare between the two mapping methods here.
3.2 Radial direction as reference
presented in their pioneering work an algorithm of radial mapping by transforming the KF magnetosheath model into a general, empirical magnetosheath shape by referring to the radial direction from the planet and scaling the radial position in the magnetosheath to the KF model. While the radial mapping can reasonably (i.e., with a relatively high accuracy) transform the KF potential into the empirical magnetosheath domain on the dayside, the mapping quality becomes degraded in the flank region due to the strongly non-orthogonal grids. Figure displays a comparison of the radial grids between the KF magnetosheath model and the empirical magnetosheath model. The grids span the radial direction to the planet and transfinite interpolation between the bow shock and the magnetopause. The radial mapping has a drawback in a stronger grid non-orthogonality effect, which causes an artificial converging flow pattern in the flank region (velocity potential shown in Fig. ) when the scalar potential is directly transformed. In the Soucek–Escoubet method, the problem of the flow conversion effect was avoided by solving the magnetohydrodynamics (MHD) Rankine–Hugoniot relation and tracking the streamline iteratively between the KF parabolic magnetosheath model and the empirical magnetosheath model.
Although the method introduced by and later adapted by is computationally less expensive than global magnetosheath simulations, the density and velocity fields from the bow shock to a given point in the magnetosheath still need to be computed along the streamline in an incremental way. Moreover, the Rankine–Hugoniot relations need to be solved before calculating iteratively the streamline, the flow velocity vector, to track the plasma density flow velocity along the streamline. Naturally, the uncertainty in this calculation depends on the step size (larger uncertainties for larger step sizes), and the errors accumulate along the streamline. The method introduced in this work overcomes this issue by constructing a magnetopause-normal grid system such that computational efforts are improved (no need to solve the Rankine–Hugoniot relations, and the error does not accumulate in the flank region).
Figure 1
Grid pattern generated by the radial mapping for the Kobel–Flückiger parabolic magnetosheath (a) and the non-parabolic, empirical magnetosheath (b).
[Figure omitted. See PDF]
Figure 2
Velocity potential in the Kobel–Flückiger model (a) and its radial mapping onto the non-parabolic empirical boundaries (b). Note that the same function is plotted here for different mapping methods. Contours represent the velocity potential normalized to the solar wind, , which is in units of the planetary radii. The gray range is chosen for the visual demonstration purpose from to (a) and from to (b).
[Figure omitted. See PDF]
3.3 Magnetopause-normal direction as referenceOur mapping method differs from the radial mapping method in that the magnetopause-normal direction is used as a reference to the magnetopause. Our method guarantees the grid orthogonality around the magnetopause both on the dayside and in the flank region. The magnetopause-normal grids are shown in Fig. for the KF magnetosheath model (with parabolic boundaries) and the empirical magnetosheath model (with non-parabolic boundaries). Even though the exact conformal mapping is not available, the magnetopause-normal mapping method retains the grid orthogonality around the magnetopause. This feature (orthogonality around the magnetopause) plays a crucial role in mapping the scalar potentials. An example of the scalar-potential mapping by referring to the magnetopause-normal direction (our final results) is shown in Sect. .
Figure 3
Mesh pattern used in the magnetopause-normal mapping in this work for the parabolic boundaries (a) and the non-parabolic, empirical boundaries (b).
[Figure omitted. See PDF]
4 Magnetopause-normal mapping4.1 Overview of the procedure
The magnetopause-normal mapping is performed with two transformations. In the first transformation, the position vector is mapped from the empirical magnetosheath onto the KF magnetosheath model . This is based on the assumption that the distance to the magnetopause along the magnetopause-normal direction when scaled to the magnetosheath thickness (defined as the distance from the magnetopause to the bow shock along the magnetopause-normal direction) remains constant. The azimuthal angle is the same between the empirical magnetosheath and the KF model. The first transformation is divided into computing the distance to the magnetopause (step 1), the thickness of the empirical magnetosheath (step 2), the thickness of the KF magnetosheath (step 3), and the mapping of the position vector onto the KF model (step 4).
In the second transformation, the mapped position vector is used to compute the shell variable and the connector variable (step 5) and to obtain the potentials and the stream function in the empirical magnetosheath using Eqs. (), (), and () (step 6). Here again, the azimuthal angle is treated as the same.
Figure illustrates the mapping procedure and graphically explains the variables that need to be determined to perform the mapping such as the zenith angle associated with the minimum distance to the magnetopause the distance from the planet to the bow shock , the distance from the planet to the magnetopause , the relative distance to the magnetopause , and the magnetosheath thickness . The position vector , the bow shock stand-off distance , the bow shock shape, the magnetopause stand-off distance , and the magnetopause shape are assumed to be known in our mapping.
Figure 4
Variables used in the magnetopause-normal mapping with the zenith angle along the direction nearest to the magnetopause, the radial distance to the bow shock and magnetopause along the magnetosheath-normal direction ( and , respectively), the distance from the magnetosheath to the magnetopause , and the magnetopause thickness . The position vector is denoted by . The bow shock and magnetopause stand-off distances are denoted by and , respectively.
[Figure omitted. See PDF]
4.2 SetupWe begin with a position vector in the empirical magnetosheath domain, and we express the position vector as . Hereafter, we present the mapping procedure in the two-dimensional plane spanning the and directions for simplicity, but the computation in three dimensions is straightforward by representing the component of position vector in the cylindrical fashion as and the component as using the distance to the axis. The boundaries (bow shock and magnetopause) are specified by the users and do not need to be parabolic. In this paper, we use the following bow shock and magnetopause models.
-
The empirical bow shock position expressed in GSE (Geocentric Solar Ecliptic) coordinates proposed and discussed by and :
16where is the bow shock stand-off distance, and is the empirical flaring parameter. We note here that the original Farris empirical bow shock model is not a paraboloid model; it is an ellipsoid model (with an eccentricity of 0.81), describing the bow shock on the dayside. It is not a proper representation of the far-flank bow shock. Also, the Cairns paraboloid bow shock model does not properly represent the far-flank bow shock. The distant bow shock shape approaches that of a hyperboloid.
-
The empirical magnetopause position by :17in the Cartesian representation and18in the polar representation.
We use a specific exponent for the Shue model (with an alpha exponent of 0.5) in an effort to show that the analytic model is “simple”. The solar wind conditions for which this exponent is applicable is not often encountered (e.g., interplanetary magnetic field has the component larger than nT, with specific values of solar wind dynamic pressure).
In our setup, the radial distance from the planet to the bow shock is expressed as 19
By introducing the zenith angle and inserting and in Eq. (), we obtain the equation for the radial distance to the empirical bow shock: 20
Equation () can be algebraically solved, and we take the positive value of the solution as presented in Eq. ().
The radial distance to the magnetopause is given conveniently by Eq. (). The Shue model reproduces the magnetopause stand-off distance in the subsolar direction (), and the cylindrical distance asymptotes to in the tail. It is worth noting here that one needs to compute the radial distance from the planet to the bow shock or magnetopause as a function of the zenith angle when using different shapes.
4.3 Step 1: measuring the distance to magnetopauseIn the first step, the distance from the given position in the magnetosheath to the nearest magnetopause is computed (see Fig. ). We express the position vector along the magnetopause-normal direction as
21 where is the magnetopause position nearest to the position vector, and is the unit vector in the magnetopause-normal direction. The unit vector points away from the planet and satisfies the condition 22 The symbol is the distance to the magnetopause along the magnetopause-normal direction in the empirical magnetosheath.
The nearest magnetopause position is obtained by searching for the zenith angle for the minimum distance from the sample position to the magnetopause. The distance is defined as 23 The search for the minimum distance is implemented in a brute-force fashion as a function of in our study.
Figure 5
Measuring the distance to the empirical magnetopause (step 1).
[Figure omitted. See PDF]
Using the minimum distance to the magnetopause and the zenith angle , we are ready to compute the magnetopause-normal direction and the distance . To obtain the magnetopause-normal direction, we define the magnetopause shape function as 24 and we compute the normal direction by the gradient of as
The magnetopause-normal direction is obtained by normalizing the gradient vector and representing it with the basis vectors ( and ) as 27 evaluated at the magnetopause ( and ). The magnetopause-normal vector has a unit length, and the sign () is chosen such that the normal vector is pointing outward (Eq. ). The distance to the magnetopause along the normal direction is obtained from Eq. () as 28 Equation () is constructed to be robust against the singular behavior on the dayside () and in distant tail ().
4.4 Step 2: computing the thickness of empirical magnetosheathIn the second step, the magnetosheath thickness is computed using the position vector and the magnetopause normal direction (Fig. ). For our mapping purpose, the distance is normalized to the magnetosheath thickness such that the relative distance serves as an invariant of the mapping from the empirical magnetosheath onto the KF magnetosheath. To achieve this, we combine Eq. () with Eq. () and analytically determine the thickness from the bow shock to the magnetopause in the empirical magnetosheath. That is, the thickness is obtained by rewriting the bow shock quadratic equation (Eq. ) for the position vector using the variable (Eq. ) extended to the bow shock location. The equation is again quadratic, and the solution is algebraically obtained as:
29 where is an auxiliary variable defined as 30
In the subsolar direction (), the thickness is simply given as 31
Equation () becomes singular in the subsolar direction, and Eq. () needs to be set separately.
Figure 6
Computing the magnetosheath thickness in the empirical model (step 2).
[Figure omitted. See PDF]
4.5 Step 3: computing the magnetosheath thickness in the KF systemIn the third step, the magnetosheath thickness is computed in the KF model (Fig. ). We repeat the procedures of steps 1 and 2 for the KF system and determine the KF magnetosheath thickness as reference. We treat the zenith angle and the relative distance as invariants of the mapping between the empirical magnetosheath and the KF system. The KF bow shock location is given as
32 where the KF bow-shock flaring parameter is pre-fixed as 33
The radial distance from the planet to the KF bow shock is 34
The KF magnetopause is defined in as 35
From Eq. (), the radial distance from the planet to the KF magnetopause is computed as 36
To obtain the magnetopause-normal direction in the KF system, we compute the gradient of the magnetopause shape function: 37
The gradient is analytically given as
The magnetopause-normal direction is then obtained by applying Eqs. () and () to Eq. (), which reads as 40
Figure 7
Computing the magnetosheath thickness in the KF model (step 3). The same zenith angle as that in step 2 is used.
[Figure omitted. See PDF]
The thickness in the KF system is determined by combining the bow shock shape (Eq. ) with the position vector at the bow shock: 41
Equation () becomes again a quadratic equation with respect to the thickness , and the solution reads 42 where the auxiliary variable is defined as 43
4.6 Step 4: mapping the position vector onto the KF systemIn the fourth step, the mapping of the position vector is performed from the empirical magnetosheath onto the KF system (Fig. ). An assumption is made such that the relative distance to the magnetopause along the magnetopause-normal direction is the same between the two systems. The distance from the magnetosheath position vector to the magnetopause along the magnetopause-normal direction in the KF system is then determined by the relative distance in the empirical magnetosheath , the thickness of the empirical magnetosheath , and magnetosheath thickness in the KF system as
44 The mapped position vector is then computed as 45 using the nearest magnetopause position (Eq. ), the magnetosheath-to-magnetopause distance (Eq. ), and the magnetopause-normal direction (Eq. ).
Figure 8
Mapping the position vector onto the KF magnetosheath model (step 4).
[Figure omitted. See PDF]
4.7 Step 5: evaluating the shell and connector variablesIn the fifth step, the shell variable and the connector variable are computed from the mapped position vector using Eqs. () and (), respectively. The variables and are the same as the parabolic coordinates used in the KF potential with a focus at . In our algorithm, the focus is explicitly given in Eqs. (), (), and (). The azimuthal angle around the symmetry axis is treated in the same way as in the KF model.
Figure 9
Evaluating the shell variable and the connector variable in the KF magnetosheath model (step 5).
[Figure omitted. See PDF]
Figure compares the iso-contours of the shell and the connector represented in the KF system (left panel) and the empirical magnetosheath (right panel) for a bow shock stand-off distance of 12.8 , a bow shock flaring of , and a magnetopause stand-off distance 9.8 . The shell variable is characterized by the lines with the curvature center on the right side in the panel and contains the parabolic bow shock (at ) and magnetopause (at ) marked by thick lines. The connector variable has the curvature center on the left side in the panel, and the iso-contour lines are orthogonal to the bow shock and magnetopause. The computation of the and variables and their gradient and curl is performed in Cartesian coordinates so that the connection represented by the Christoffel symbol vanishes in the computation. Computation in the Cartesian domain is also beneficial to the practical application because spacecraft trajectories are often represented in Cartesian coordinates.
Figure 10
Iso-contour lines with (center of curvature on the left side) and that with (center of curvature on the right side) in the KF magnetosheath model (a) and the empirical magnetosheath model (b). The bow shock stand-off distance is 12.8 Earth radii and the magnetopause stand-off distance is 9.8 Earth radii.
[Figure omitted. See PDF]
4.8 Step 6: computing the potentials and stream functionThe scalar potentials (velocity potential and magnetic potential) and the stream function are obtained from the shell and the connector using Eqs. (), (), and (). The velocity potential (normalized to the upstream flow) is displayed in Fig. left panel, and the stream function is displayed in the right panel. The iso-contours of the velocity potential represent the lines of the same flow velocity. The iso-contours of the stream function represent the streamlines in the magnetosheath. The flow is deflected around the nose of the magnetopause (the subsolar point), and the streamlines are tangential to the magnetopause. Due to the grid orthogonality around the magnetopause, the streamlines are constructed as tangential to the magnetopause shape, which qualifies the magnetopause-normal mapping method as a useful tool for the magnetosheath model.
Figure 11
Velocity potential (a) and stream function (b) in the empirical magnetosheath domain obtained by mapping onto the shell variable and the connector variable . Note that two different functions are plotted here for the same grid and mapping method. The gray range is chosen for the visual demonstration from to (a) and from to (b).
[Figure omitted. See PDF]
The magnetic potential and the derived magnetic field are displayed in Fig. . The magnetic potential and the magnetic field (the gradient of the potential multiplied by the minus sign) depend on the upstream field. Figure shows an example with an upstream field angle of 135° to the axis (i.e., 45° to the upstream flow direction). The magnetic field is computed using the central difference scheme. Near the boundaries (bow shock or magnetopause), the mesh resolution is enhanced so that the mesh points do not cross the boundary when computing with the central difference scheme. The upstream field is deflected on the positive side (right panel, lower half plane) and is draping the magnetopause on the negative side (right panel, upper half plane).
Figure 12
Magnetic potential for the upstream magnetic field with an angle of 135° to the axis (45° to the upstream flow direction (a) and sampled magnetic field vectors obtained by the negative gradient of the magnetic potential (b). The gray range of the magnetic potential is from to .
[Figure omitted. See PDF]
5 Discussion5.1 Extension of the mapping approach
It is straightforward to extend our method to different shapes of the Shue magnetopause model. The following form of gradient can be used for a general value of the Shue exponent : where in Eqs. ()–() is defined as
48
Mesh pattern applied to different values of the Shue exponent for an open-type magnetopause (Fig. a) and a closed-type magnetopause (Fig. b).
Figure 13
Mesh pattern applied to different values of the Shue exponent for an open-type magnetopause (a) and a closed-type magnetopause (b).
[Figure omitted. See PDF]
Also, our method can be extended to a three-dimensional, non-axisymmetric geometry of magnetosheath
It is possible to obtain the steady-state magnetosheath potential in different ways.
-
First, one may numerically solve the Laplace equation for a given set of boundary shapes (bow shock and magnetopause). Various numerical solvers are known for solving the Laplace equation such as the Jacobi method, the Gauss–Seidel method, and the successive over-relaxation (SOR) method. These Laplace solvers are numerically more expensive than the mapping method, but the computation in 3-D is feasible with the contemporary computational resources. On the other hand, the magnetosheath is not bounded but extends in the tail direction. The challenge here is thus to construct a properly bounded area for the Laplace equation.
-
Second, one may expand the magnetosheath magnetic field in different orthogonal functions. The KF solution makes use of the Bessel functions . For flexible magnetopause and bow shock boundary models, a magnetosheath magnetic field model is constructed by making use of Legendre polynomials .
-
Third, one may introduce a suitable conformal mapping by limiting the magnetosheath modeling to a complex plane (two-dimensional domain). The harmonic functions such as the KF solution are transformed from the parabolic magnetosheath shape into the non-axisymmetric magnetosheath shape. The problem here is that finding the conformal mapping is not an easy task, because the magnetosheath is not a spatially-closed domain and one has to set the boundary in the magnetosheath to complete the domain bounded by the bow shock, the subsolar axis, and the magnetopause.
-
Fourth, one may solve the Rankine–Hugoniot relations and track the streamline stepwise by referring to the KF solution, as is done in . This method computes the magnetosheath flow and magnetic field along the streamline, and the computing for the entire magnetosheath domain is numerically expensive.
6 Concluding remarks
Our potential mapping method may be regarded as an updated version of the radial mapping method by retaining the orthogonality near the magnetopause in the flank to tail region and also by computing the field through the potential mapping. Velocity potential, stream function, and magnetic potential are evaluated in the empirical magnetosheath.
The advantages of our methods are as follows:
-
The method makes extensive use of the exact solution of the Laplace equation (the Kobel–Flückiger potential and the Guicking stream function). The plasma flow and magnetic field can be determined semi-analytically in a wide range of zenith angles in the magnetosheath when the solar wind conditions and the boundary shapes are given.
-
The method is applicable to a non-parabolic shape of magnetosheath domain, opening the door to develop a tool to assist numerical simulations and spacecraft observations of not only the Earth but also the planetary magnetosheath domain.
-
The method is computationally inexpensive. In particular, if the shapes of bow shock and magnetopause are analytically given, most of the computational steps in the potential mapping method have an analytic expression.
As stated in Sect. , one ideally needs to find a conformal mapping from the KF magnetosheath model onto the empirical magnetosheath. While the conformal mapping is known both for the empirical bow shock and the empirical magnetopause, the conformal mapping of the entire magnetosheath domain still remains a challenge. There are two problems with this approach. First, the closing boundary (the contours) connecting between the bow shock and the magnetopause is not known, and moreover, the uniqueness of finding such a boundary is not guaranteed. Second, the gradients along are not the same between the empirical bow shock and the empirical magnetopause such that a naive transfinite interpolation ends up with highly non-orthogonal grids in the magnetosheath.
Our method of computing the plasma flow and magnetic field should be compared against the observations and simulations. For example, THEMIS and ARTEMIS spacecraft and MMS spacecraft are providing a huge amount of data on both sides of the bow shock in the equatorial plane; Cluster spacecraft are collecting data in polar orbit; ACE spacecraft data may be used as an upstream monitor; and Earth Flyby data of planetary missions (such as Cassini, BepiColombo) cover the far-distance tail region. In reality, non-axisymmetric structure arises in the magnetosheath. There is no restriction regarding the choice of the magnetopause model. The magnetopause-normal direction needs to be computed either analytically using the gradient of the magnetopause function as or numerically for a user-defined magnetopause shape.
Code and data availability
The IDL (Interactive Data Language) codes are available upon request to the authors.
Author contributions
YN, DS, and ST developed the idea of the potential mapping method; checked mathematics; and wrote the manuscript. YN prepared the figures. All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
Financial support
This open-access publication was funded by Technische Universität Braunschweig.
Review statement
This paper was edited by Minna Palmroth and reviewed by Steven Petrinec, Octav Marghitu, and two anonymous referees.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2024. This work is published under https://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
The steady-state magnetosheath model has various applications for studying the plasma and magnetic field profile around the planetary magnetospheres. In particular, the magnetosheath model is analytically obtained by solving the Laplace equation for parabolic boundaries (bow shock and magnetopause). We address the question, how can we utilize the magnetosheath model by transforming into a more general, empirical, non-parabolic magnetosheath geometry? To achieve the goal, we develop the scalar-potential mapping method which provides a semi-analytic estimate of steady-state flow velocity and magnetic field in the empirical magnetosheath domain. The method makes use of a coordinate transformation from the empirical magnetosheath domain into the parabolic magnetosheath domain and evaluates a set of variables (shell variable and connector variable) to utilize the solutions of the Laplace equation obtained for the parabolic magnetosheath domain. Our model uses two invariants of transformation: the zenith angle within the magnetosheath with respect to the direction to the Sun and the ratio of the distance to the magnetopause and the thickness of magnetosheath along the magnetopause-normal direction. The use of magnetopause-normal direction makes a marked difference from the earlier model construction using the radial direction as reference. The plasma flow and magnetic field can be determined as a function of the upstream condition (flow velocity or magnetic field) in a wide range of zenith angles. The scalar-potential mapping method is computationally inexpensive, using analytic expressions as much as possible, and is applicable to various planetary magnetosheath domains.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer