Electro-magneto-hydrodynamic lubrication

Krzysztof Wierzcholski; Andrzej Miszczak

doi:10.1515/phys-2018-0040

Artikel Open Access

Electro-magneto-hydrodynamic lubrication

Krzysztof Wierzcholski und Andrzej Miszczak

Veröffentlicht/Copyright: 30. Mai 2018

Veröffentlicht von

Veröffentlichen auch Sie bei De Gruyter Brill

Manuskript einreichen Informationen für Autor*innen Erkunden Sie dieses Fachgebiet

Aus der Zeitschrift Open Physics Band 16 Heft 1

Abstract

The topic of the presented paper aims to demonstrate a new principle of hydrodynamic lubrication in mechanical, thermal and electro-magnetic fields. Up till now, when dealing with the hydrodynamic theory lubrication, many authors of scientific papers have assumed the constant oil dynamic viscosity value without variations caused by temperature crosswise the film thickness. Simultaneously, due to the numerous AFM measurements, it appears that oil temperature gradients and oil viscosity changes in the bearing gap height directions cannot be omitted. Therefore, in this paper, the problem of the viscosity changes across the lubricant thin layer was resolved as the main novelty in principles of mechanical thermal lubrication. The method of solving the mentioned problem was manifested by a general model of semi-analytical solutions of isothermal electro-magneto-elastohydro-dynamic and non-Newtonian, lubrication problem formulated for two deformable rotational surfaces in curvilinear, co-ordinates.

Keywords: 3D-lubrication; viscosity variations crosswise film thickness; semi-analytical solutions; curvilinear coordinates

PACS: 47.50.-d; 47.65.-d; 47.10.ab; 47.10.ad; 47.15.gm

1 Introduction and governing equations

This paper presents a semi analytical method of solution of the asymmetrical, laminar, steady, non-Newtonian lubrication flow problem between two rotational and deformable, curvilinear orthogonal movable surfaces in conjugated elasto-hydro-electro-magnetic fields. The fluid flow between the two above mentioned solid surfaces in the electromagnetic field will be described by the 3 momentum equations of equilibrium in a vector form, a fluid continuity equation, and by equation of a conservation of energy equation in a scalar form. Hence, we obtain the following system [1, 2]:

DivS+μo(N∇)H+12μorotN×H+J×B=ρdvdt(1)

∂ρ∂t+divρv=0,(2)

divκgradT+φF=ρddtcvT+μoTΞv∇H+J2/σ.(3)

The above-mentioned system of equations is completed by magneto-thermo-elasticity equations describing problem in stresses for two solid surfaces restricting the thin fluid layer. This system consists of the three partial differential equations (4) in a vector form. Heat conductivity equation in a solid body (5) is added to this set of equations and we obtain a system in the following form [2, 3, 4]:

DivS∗+J∗×B∗+μoN∗∇H∗=ρ∗∂2u∂t2(4)

divκ∗gradT∗=ρ∗cv∗∂T∗∂t(5)

Moreover, we add the Maxwell and Ohm equations as well for the two surfaces as for the thin boundary liquid layer between two surfaces. Thus we get such equations [5]:

∇⋅B∗=0,∇⋅B=0,∇×H∗=J∗+∂D∂t∗,∇×H=J+∂D∂t(6)

∇×E∗=−∂B∗∂t,∇×E=−∂B∂t,J∗=σ∗E∗,J=σ(E+v×B)(7)

We assume the following notations: μ_o – magnetic permeability in vacuum H/m, T – fluid temperature in K, T^* – solid body temperature in K, B – magnetic induction vector in T, N - magnetization vector A/m, E – electric intensity vector V/m, H – magnetic intensity vector A/m, (– first derivative of magnetization vector with respect to temperature A/mK, σ - electrical conductivity coefficient S/m, J - electric current density in A/m², D - electric induction vector As/m, ρ - fluid density kg/m³, κ - thermal conductivity coefficient W/mK, v – fluid velocity vector in m/s, φ_F – dissipation of energy in W/m³, S - stress tensor in the fluid in Pa, u - displacement vector of the solid body in m, t - time in s, c_v - specific heat in J/kgK. The symbols with an asterisk are related to the solid body.

The relationship between stress tensor S and strain rate tensor 2T_d = A₁i.e. constitutive equations are assumed for the lubricant in the following form [3]:

S=−pδ+ηpA1(8)

whereas unit tensor δ, strain rate tensor A₁ have the following components: δ_ij, Θ_ij in s^–1. For power law of compressible fluid the apparent viscosity η_p has the form [2, 3, 4]:

ηp≡ηpr=2n−1m(n)12I12−I2n−12,φF=4ηp12I12−I2−pI1,(9a)

I1=divv,I2=Θ11Θ22−Θ122+Θ11Θ33−Θ132+Θ22Θ33−Θ232,(9b)

We assume: n – dimensionless flow index (0.5, 1.2), m = m(n) fluid consistency coefficient in Pasⁿ.

The Duhamel Neumann relations between the components τ_ij of the stress tensor S^* of the elastic body on the sleeve and the strain tensor components ϵ_ij take the following form:

τij=2Gεij+Λεkk−3αTKΔT∗δij(10)

for i,j = 1,2,3, where: δ_ij – Kronecker unit tensor component (δ_ij = 1 for i = j and δ_ij = 0 for i≠j), K = Λ + 2G/3, G = E_Y/[2(1+ ϑ)], Λ = E_Yϑ/(1+ϑ)(1–2ϑ), whereas E_Y - Young‘s Modulus in Pa, α_T – thermal coefficient of linear expansionin K^–1, ϑ – Poisson ratio.

And by virtue of known linear geometrical relations, the strain tensor components Θ_ij(v_i) in the oil, and strain components ϵ_ij(u_i) in the sleeve are mutually connected with the oil velocity components v_i and displacement vector of solid body components u_i, respectively:

Θij=0.5(vij+vji),εij=0.5(uij+uji).(11)

2 Electro-magneto-hydrodynamic equations

We assume: unsymmetrical, incompressible (for invariant I₁ = 0), steady, pseudo-plastic lubrication in electro-magnetic field for power law model without viscous-elastic properties. Oil apparent viscosity η_p varies in length, width and gap-height directions and depends on pressure, temperature flow shear ratio, and electro-magnetic field [1, 2, 3, 4, 5]. Inertia forces and terms of energy convections are neglected, pressure is constant in gap height direction. The total distance ϵ_T between the two surfaces is signifucantly smaller than other dimensions of the considered surfaces. Taking into account the layer boundary simplification i.e. neglecting in (1)-(7) the negligibly small terms of order Ψ = 0.001 (radial clearance) presenting the quotient of characteristic gap height (₀ to the radius R of journal, the problem is made in local curvilinear and orthogonal coordinates (α₁, α₂, α₃) connected with one of the movable surfaces, where α₂ denotes the direction of gap height. The parallel and longitudinal intersections of the cooperating deformed, rotational surfaces have curvilinear and non-monotone generating lines. Hence, Lame coefficients have the form: h₁ = h₁(α₃), h₂ = 1, h₃ = h₃(α₃).

We put physical and geometrical dependencies (9), (10), (11) into expanded boundary simplified equations (1)-(7). Hence, we obtain the system of equations of conservation of momentum, continuity, energy, thermo-elasticity, heat transfer, Maxwell simplified results in the following form of curvilinear orthogonal co-ordinates (α₁, α₂, α₃):

0=−1hi∂p∂αi+∂∂α2ηpα1,α2,α3∂vi∂α2+Mi,∂p∂α2=0,(12)

0=1h1∂v1∂α1+∂v2∂α2+1h1h3∂∂α3h1v3(13)

∂∂α2κ∂T∂α2+[ηpα1,α2,α3]Vv1,v3=MT.(14)

∂∂α2κ∗∂T∗∂α2=0,∂∂α22G+Λδ2i∂ui∂α2==δ2i∂∂α23KαT∗T∗+J∗×B∗i−μ0N∗∇Hi∗(15)

ηpα1,α2,α3≡mn,T,B,E⋅Vg(v1,v3),V≡∂v1∂α22+∂v3∂α22,(16)

gn≡n−12,MTα1,α3≡μ0TΞv∇H+J2/σ,Miα1,α3≡μ0N∇Hi+12μorotN×Hi+J×Bi,(17)

where the length, width and gap- height directions, are limited respectively: 0 < α₁ ≤ 2π, –b_m ≤ α₃b_s, 0 ≤ α₂ ≤ ϵ_T and i = 1,2,3. The system of equations (12)((17) contains the following unknown dimensional values: pressure p(α₁,α₃), temperature in oil T(α₁,α₂,α₃) and sleeve T^*, three oil velocity components v_i (α₁,α₂,α₃) and three layer displacement components u_i (α₁,α₂,α₃), for i = 1,2,3 in three curvilinear, orthogonal dimensional directions: α₁,α₂,α₃.

3 Boundary conditions

The lubricant flow in the bearing gap is generated by the rotation of a rotational, curvilinear journal. Hence the boundary conditions for lubricant velocity components have the form:

v1=ωh1forα2=0,v1=0(18)

for α₂ = ϵ_T, v₂ = 0 for α₂ = 0,

and for α₂ = ϵ_T, v₃ = 0 for α₂ = 0, v₃ = 0 for α₂ = ϵ_T,

where ω denotes the angular velocity of rotational journal, in the circumferential direction.

Decrements or increments that are above the characteristic environmental temperature T₀ have constant value f_c on the journal surface and variable unknown values f_p (α₁, α₃) on the sleeve surface. Heat flux is transferred from the rotational surface of bearing journal into lubricant. Hence the boundary temperature values are as follows:

T(α1,α2,α3)=T0+fcforα2=0,(19a)

T(α1,α2,α3)=T0+fp(α1,α3)forα2=ϵT(19b)

κ∂T∂α2=−qc≡−υΔfforα2=0,(19c)

where Δf – temperature difference across the film thickness, υ – heat transfer coefficient in W/m²K, q_c - heat flux density in oil on journal surface in W/m².

The elastic layer of bearing alloy is laying on the rigid ring in place α₂ = c₃, and therefore the contact surface of these bodies is not deformed by the pressure. The contact surface α₂ = c₃ is yet deformed by variable temperature, which implies the shape and the total volume change of the rigid ring. These deformations are not considered now.

The temperatures and surfaces localization is presented in Figure 1.

Figure 1

Motionless solid A₂ and moving body A₁ loaded by means of pressure and temperature and magnetic field; N, S - magnetic poles; a) cylindrical bearing, b) other bearing

The boundary conditions in the elastic layer have the following form:

τij∗(α1,α2=c2,α3)=−p(α1,α3)δi2δij(20)

T∗(α1,α2=c3,α3)=T0+e2(α1,α3)(21)

T∗(α1,α2=c2,α3)=T0+f2(α1,α3)(22)

κ∗∂T∗∂α2=−qp,(23)

where i = 1,2,3. Moreover: T₀ - average ambient temperature, q_p - heat flux density across the elastic layer located on the sleeve surface, f₂(α₁, α₃) - temperature increases above the temperature T₀ on the external surface of the elastic layer for α₂ = c₂, additionally e₂(α₁,α₃) temperature increases in the excess of temperature T₀ on the bottom surface of the elastic layer in the place of α₂ = c₃.

4 Oil consistency and viscosity variations

Temperature decreases and magnetic induction, electric intensity increases the oil consistency m and dynamic viscosity coefficient η. This phenomenon is described in the following form:

m(α1,α2,α3)=m(n,T,B,E)=η0η1(T,B,E)⋅v0/ε01−n,η1(T,B,E)=expδBB+δEE−δTT,(24)

with coefficient δ_B, δ_E of magnetic, electric influence on the oil viscosity and coefficient δ_T of temperature influence on the oil viscosity Mentioned phenomena are confirmed by virtue of authors own measurements.

5 Sketch of semi-analytical solutions

After double integration of two equations (15) in α₂ direction we apply the low heat flux continuity (q_p = q_c) in the contact region between the elastic alloy layer and lubricant, and finally we obtain the external displacement of the sleeve surface in the following dimensional form:

u02=u2α1,α2=c2,α3==εs1+ϑ1−ϑ⋅1−2ϑEp−fp−12εsqpκ∗εsαT∗1+ϑ1−ϑ++12⋅1−2ϑE⋅1+ϑ1−ϑεs2N∗∇B2∗+J∗×B∗2(25)

where ϵ_s = c₃ – c₂, 0 < α₁ < 2π, –b_m < α₃ < +b_m. Solution (25) valid if and only if E_Y - Young modulus and solid body coefficients α_T, κ^*, ϑ - Poisson ratio, are independent of α₂, The hydrodynamic pressure p, the load carrying capacity and temperature T distributions with oil velocity components v₁, v₂, v₃ in machine slide journal bearing are solved by the means of the quasi half numerical Frobenius small parameter method implemented by the finite difference method using Mathcad 15 Professional Program. A system of partial differential equations (12)-(14) defined in thin space (α₁, α₂, α₃, t) between two movable rotational surfaces and for continuous m(n), η == m(n = 1), has an analytical solution in the form of the following infinite uniform convergent power function series in following form:

vi(α1,α2,α3)=vi(0)+gvi(1)+…+gkvi(k)+…T(α1,α2,α3)=T(0)+gT(1)+…+gkT(k)+…(26a)

p(α1,α3)=p(0)+gp(1)+…+gkp(k)+…qc=qc(0)+gqc(1)+…+gkqc(k)+…(26b)

where qc(0)=qc(1)=…=qc(k)≡q, with q_c = q/(1 – g); for i = 1,2,3; k = 0,1,2,…, g = g(n)≡ (n-1)/2 where for 0 < n ≤ 3/2 the small parameter g is less than +1/4 and greater than –1/2. In order to achieve the linearization of apparent viscosity (16), we put functions (26a) into Eq. (16) and expand function in closed interval [1, n] or [n, 1] for 0 < n ≤ 3/2 using Taylor series in the neighborhood of point n = 1 with respect to the small parameter in the following form:

ηpα1,α2,α3=η0⋅η1T,B,E⋅1+gηpr1+…+gkηprk+…(27)

where η₁(T, B, E) = η₁(α₁, α₂, α₃) and functions η_prk ≡ η_prk(v₁, v₃) for k = 1, 2, … – are the dimensionless expansion coefficients of dynamic viscosity dependent on α₁, α₂, α₃, whereas for k = 0 η_pr0 = 1. Such functions after involved calculations are determined in analytical form. Putting the infinite series (26), (27) into the system (12)-(14), multiplying the series by Cauchy methods and equating the terms by the same powers k = 0,1,2,… of the small parameter g, and then using the L. Kronecker symbol i.e. δ_k0 = 1 for k = 0, δ_k0 = 0 for k = 1,2,…, we obtain the linear differential systems of equations for k = 0,1,2,… i = 1,3 in following form:

1hi∂p(k)∂α1=∂∂α2η0η1B,E,T(k)∂vi(k)∂α2++∂∂α2η0η1B,E,T(k)Si(k)+Miδk0,∂p(k)∂α2=0,(28)

1h1∂∂α1v1(k)+∂∂α2v2(k)+1h1h3∂∂α3h1v3(k)=0,(29)

∂∂α2κ∂T(k)∂α2=−η0η1B,E,T(k)VT(k)+MTδk0.(30)

with

η1B,E,T(k)=η1BEB,Ee−δTT(k)α1,α2,α3,η1BEB,E≡η1BEα1,α3,Si(0)=0,VT(0)≡Vv1(0),v3(0),Si(1)≡∂vi(0)∂α2ε0v02VT(0),(31a)

VT(1)≡VT(0)⋅lnε0v02VT(0)+2∂v1(0)∂α2∂v1(1)∂α2+2∂v3(0)∂α2∂v3(1)∂α2,(31b)

For k = 0 the abovementioned system of equations (28)-(31) determines functions: vi(0),Ti(0),p(0) for Newtonian oil taking into account the boundary conditions in steady motion (18), (19). For k = 1,2,… and i = 1,2,3 the system of equations (28)-(31) determines the first, second,… corrections: vi(k),T^(k), p^(k) caused by the non-Newtonian oil. The boundary conditions for k = 0 are as follows:

v1(0)=ωh1,v2(0)=0,v3(0)=0,T(0)=T0+fc(32a)

for α₂ = 0,

v1(0)=0,v2(0)=0,v3(0)=0,T(0)=T0+fp(0)

for α₂ = ε_T;

and for k = 1,2,… we have:

vi(k)=0,T(k)=sfpk,fors=α2/εT,α2=0,α2=εT;(32b)

The total gap height ϵ_T includes surface displacement (25). Integrating twice (28) and (29) once with respect to the variable α₂, then the oil velocity components for k = 0,1,2,… by virtue of boundary conditions (32a, b) have the following form:

vi(k)=1η0η1BE⋅1hi∂p(k)∂αiKi(k)(α1,α2,α3)++K2(k)(α1,α2,α3)⋅∫0εTSi(k)dα2−∫0α2Si(k)dα2++ωh1[1−K2(k)(α1,α2,α3)]δi1δk0+δk0ΔMi,(33)

v2(k)=−1h1∫0α2∂v1(k)∂α1dα2−1h1h3∫0α2∂∂α3h1v3(k)dα2,(34)

for i = 1, 3; 0 < α₂₁ < α₂ < ϵ_T, 0 < α₁ < 2π, –b_m < α₃ < +b_m, whereas we denote:

ΔMi=−1η0η1BE∫0α2∫0α21Midα21dα2+1η0η1BEs∫0εT∫0α21Midα21dα2=εT22η0η1BEMi(s−s2),(35a)

K1(k)(α1,α2,α3)≡detΩ(k)Ω1(k)(α1,α2=εT,α3),K2(k)(α1,α2,α3)≡Ω1(k)(α1,α2,α3)Ω1(k)(α1,α2=εT,α3),(35b)

Ω(k)≡Ω2(k)(α1,α2,α3)Ω2(k)(α1,α2=εT,α3)Ω1(k)(α1,α2,α3)Ω1(k)(α1,α2=εT,α3),Ω1(k)(α1,α2,α3)≡∫0α2e−δT⋅T(k)dα2,Ω2(k)(α1,α2,α3)≡∫0α2α2⋅e−δT⋅T(k)dα2.(35c)

It is easy to see, that:

K1(k)(α1,α2=0,α3)=0,K1(k)(α1,α2=εT,α3)=0,K2(k)(α1,α2=0,α3)=0,K2(k)(α1,α2=εT,α3)=1.(35d)

Imposing the boundary condition (32ab) in point α₂ = ϵ_T on the solution (34) and taking into account the law of integral differentiation with variable limits of integration, we obtain the following modified Reynolds equation, which determines the unknown function p^(k)(α₁, α₃):

1h1∂∂α11η0η1BE(α1,α3)∂p(k)(α1,α3)∂α1∫0εTK1(k)(α1,α2,α3)dα2+1h3∂∂α31η0η1BE(α1,α3)h1h3∂p(k)(α1,α3)∂α3∫0εTK1(k)(α1,α2,α3)dα2++∂∂α1ΠNNi(k)−δk0εT3M1η0η1BE(α1,α3)+∂h1h3∂α3ΠNN3(k)−δk0εT3M3η0η1BE(α1,α3)==ωh1∂∂α1∫0εTK2(k)(α1,α2,α3)dα2−εT,(36)

with

ΠNNi(k)(α1,α3)≡∫0εTSi(k)dα2∫0εTK2(k)(α1,α2,α3)dα2−∫0εT∫0α2Si(k)dα2dα2,(37)

for i = 1,3; k = 1,2,3.

To obtain the temperature functions T^(k)(α₁α₂, α₃) for k = 0,1,2,…, we put solutions of oil velocity components (33), (34) for i = 1,3; into the temperature equation (30). After double integration with respect to the variable α₂i.e. in film thickness direction, using boundary conditions (19ab) in the form: T^(k) = f_cδ_k0 for α₂ = 0, T^(k) = f_pδ_k0 for (α₂ = ϵ_T, and after the term ordering, we finally obtain:

T(k)=η0η1BEs∫0εT1κ∫0α2VT(k)e−δTT(k)dα2dα2−∫0α21κ∫0α2VT(k)e−δTT(k)dα2dα2+∫0α21κ∫0α2MTdα2dα2−s∫0εT1κ∫0α2MTdα2dα2dα2δk0+δk0fc1−s+fps,(38)

where s ≡ α₂/ϵ_T, 0 ≤ s ≤ 1. Non linear differential equation (38) determines temperature functions T^(k) for k = 0,1,2,… Imposing the condition (19c) on the temperature T^(k) determined by the (38), we obtain the unknown temperature f_p (α₁, α₃) on the sleeve.

To obtain the oil velocity components (33), hydrodynamic pressure (36) and temperature (38) for Newtonian oil, and taking into account oil viscosity variations crosswise the gap height, we put k = 0 and hence: Si(0)=0,ΠNNi(0)=0 for i = 1.2.3. In this case we obtain oil velocity components vi(0), pressure p⁽⁰⁾, temperature T⁽⁰⁾. Assuming cylindrical coordinates: α₁ = ϕ, α₂ = r, α₃ = z, h₁ = R, h₂ = h₃ = 1, where R is the radius of the cylindrical journal, we put the temperature T⁽⁰⁾, determined from the reduced Eq. (38), into the reduced pressure Equation (36). Finally, the Eq. (36) determines the pressure p(ϕ, z), 0 ≤ ϕ ≤ 2π, –b≤ z ≤ b, and tends to the following form:

∂R∂ϕ∂p(0)R∂ϕKJ1+∂∂z∂p(0)∂zKJ1=ω∂KJ2∂ϕ,KJ1≡εT3ηoη1BE+δTεT3ηoη1BE{fc−q2κεT−15ωR2ηoη1BEκ−160ηoη1BEκεT4∂p(0)R∂ϕ2+∂p(0)∂z2},KJ2=6εT−6δTεT{16qκεT+124ωR2ηoη1BEκ+1144ηoη1BEκεT4∂p(0)R∂ϕ2+∂p(0)∂z2},(39)

0≤ ϕ ≤ 2π, –b ≤ z ≤ b, p = p(ϕ, z), η_1BE = exp(δ_BB + δ_EE).

For δ_T = 0, viscosity is constant in gap height direction and (39) tends to the well known classical form of Reynolds equation.

6 Conclusions and results

After semi-analytical solutions and initial numerical calculations it appears that the pressure obtained directly for the oil dynamic viscosity variations crosswise the film thickness caused by the temperature gradients in gap height directions are about 5 to 7 percent different in comparison with the pressure values calculated for constant oil viscosity across the film thickness i.e. in the case when the temperature is constant in gap height direction and varies only in bearing length and circumference direction.
The mutually direct interactions of influences between hydrodynamic pressure on temperature as well as temperature on pressure, are valid only if the temperature, hence oil dynamic viscosity, varies in the gap height direction.
It is worse to notice that constant temperature in the gap height direction does not correspond with 3D temperature field obtained directly from energy equation. The assumptions of constant temperature and viscosity crosswise the bearing gap are in contradiction with the contemporary achievements connected with new devices such as micro-bearing, nano-bearing, magnetic bearings, artificial joints in humanoid robots, micro-motors. Unfortunately, numerous authors in the field of hydrodynamic lubrication, avoid to assume the oil viscosity variations in the gap height direction.
The performed numerical calculations show that oil, presented in the calculation, increases the pressure and load capacity in the slide journal bearing for Newtonian and non-Newtonian oil, in the presence of magnetic induction field.

References

[1] Khonsari M.M., Wang S.H., On the Maximum Temperature in Double Layered Journal Bearing, Journal of Tribology, 1991, 113, 464-469.10.1115/1.2920647Suche in Google Scholar

[2] Pincus O., Wilcock D.J., Thermal Effects in Fluid Film Bearings, Proceeding of the 6th Leeds-Lyon Symposium on Tribology: Thermal effects on tribology, 1980, 3-23.Suche in Google Scholar

[3] Massoudi M., Local non-similarity solutions for the flow of a non-Newtnian fluid over wedge, International Journal of Non-Linear Mechanics, 2001, 36, 961-976.10.1016/S0020-7462(00)00061-5Suche in Google Scholar

[4] Makinde O.D., On thermal stability of a reactive third-grade fluid in a channel with convective coolingthe walls, Applied Mathematics and Computation, 2009, 213, 170-176.10.1016/j.amc.2009.03.003Suche in Google Scholar

[5] Schweitzer G., Traxler A., Bleuer H., Magnetlager, Springer Verlag, Berlin, Heidelberg, 199310.1007/978-3-662-08448-9Suche in Google Scholar

Received: 2017-11-02

Accepted: 2017-11-12

Published Online: 2018-05-30

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.

Artikel in diesem Heft

https://doi.org/10.1515/phys-2018-0040

Schlagwörter für diesen Artikel

3D-lubrication; viscosity variations crosswise film thickness; semi-analytical solutions; curvilinear coordinates

Creative Commons

BY-NC-ND 4.0