{"id":"f1e9077b-cad7-4e35-90dd-84300bbce815","arxiv_id":"2411.18213","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Closed-form solutions for the radially symmetric extension of a circular cylinder are derived for the isotropic relaxed micromorphic model in polar coordinates, recovering classical elasticity in the limiting length-scale cases.","lead":"The paper derives the isotropic relaxed micromorphic model in polar coordinates and solves an axisymmetric extension problem in closed form using modified Bessel functions. It then shows that classical linear elasticity is recovered in the limits of vanishing and infinite characteristic length, and plots the influence of material parameters.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The plane-strain reduction to a 2D microdistortion is not justified: the full 3D relaxed micromorphic model has out-of-plane P components, so the claimed long-cylinder benchmark may solve a different model.","rationale":"The reader's conditional verdict is appropriate, and I do not recommend changing it. The reader identified the plane-strain assumption as a weak point, but the specific mechanism is stronger and more concrete than the reader stated: even if the axisymmetric ansatz and uθ=0 hold, the reduction from 3D to 2D is not a consequence of plane strain because the microdistortion P is kinematically independent of u. A nonzero, spatially constant Pzz is allowed, has zero CurlP energy, and is energetically favored through the trace couplings. This can change the in-plane displacement and microdistortion, so the closed-form benchmark may describe the 2D planar relaxed micromorphic model rather than the 3D long-cylinder problem advertised in §4.1. The concrete test above would settle this directly. A smaller internal issue also exists: Eq. (30) appears inconsistent with Eq. (13) under the axisymmetric reduction, but the out-of-plane reduction is the more fundamental concern and does not change the conditional verdict.","tokens_in":13208,"tokens_out":32095,"duration_ms":274631,"concrete_test":"Evaluate the stationarity condition for a constant out-of-plane component c=Pzz at a representative in-plane state, using the Table 1 parameter sets and Eq. (1). If the optimal c is nonzero, the Pzz=0 ansatz is inconsistent with the 3D energy. To settle the concern fully, run a 3D axisymmetric finite-element solve of Eqs. (3)-(4) on a long cylinder with the consistent-coupling lateral boundary conditions and compare ur(r), Prr(r), Pθθ(r) with Eqs. (68)-(70); if the full solution has Pzz≠0 or a different displacement profile, the paper's formulas are not the 3D plane-strain solution.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing gap is the passage from the 3D relaxed micromorphic model to the axisymmetric plane-strain problem in §4.1. The ansatz keeps only ur, Prr, Pθθ, Prθ, Pθr, and implicitly sets Pzz=Pzr=Prz=Pθz=Pzθ=0. The paper asserts that 'plane-strain conditions are warranted' for a long cylinder, but the relaxed micromorphic model has an independent microdistortion field: u_z=0 and ∂_z=0 do not force Pzz to vanish. In the 3D energy (1)-(4), Pzz enters through tr(Du−P) and trP in the λe and λmicro terms, and a constant Pzz produces no CurlP penalty. For fixed in-plane data with s=div u and p=Prr+Pθθ, the stationarity condition with respect to a constant c=Pzz is −λe(s−p−c)+2µmicro c+λmicro(p+c)=0, whose solution c=[λe(s−p)−λmicro p]/(λe+2µmicro+λmicro) is generically nonzero for the Table 1 parameters. Therefore the true 3D plane-strain minimizer has out-of-plane microdistortion, and the 2D formulas (68)-(70) are not the solution of the 3D long-cylinder problem unless a special parameter condition holds. The formulas may still be valid for the explicitly 2D model announced in §3, but the benchmark claim for a long circular cylinder is not supported without an additional reduction argument.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives the governing equations of the isotropic relaxed micromorphic model in two-dimensional polar coordinates and applies them to an axisymmetric extension problem for a long circular cylinder under uniform radial displacement. The authors assume plane strain, set uθ=0 and all fields independent of θ and z, and restrict the microdistortion to its in-plane components. Using the consistent coupling boundary conditions, they reduce the problem to linear ODEs and obtain closed-form expressions involving modified Bessel functions for the displacement and in-plane microdistortion components. They then study the limits Lc→0 and Lc→∞, which recover classical elasticity, and present numerical results for three parameter sets from the literature.","tokens_in":13511,"tokens_out":11403,"duration_ms":100787,"significance":"If the polar-coordinate equations and the closed-form solution are correct, the paper provides a useful analytical benchmark for validating finite element implementations of the planar relaxed micromorphic model. The derivation is self-contained, does not fit parameters to data, and the limit-case recovery of classical elasticity is a valuable internal check. However, the claimed application to a three-dimensional long cylinder is not justified by the presented reduction: the independent microdistortion field has out-of-plane components that are silently discarded. The significance of the paper as a 3D benchmark therefore depends on an additional reduction argument or on an explicit reframing as a planar 2D model.","major_comments":[{"comment":"The reduction from the three-dimensional model (1)-(4) to the plane-strain problem drops all out-of-plane microdistortion components without justification. For the axisymmetric ansatz u=(u_r(r),0,0) with ∂_z=0, the energy still depends on c=P_zz through the terms µe|sym(Du-P)|², λe tr²(Du-P), µm|sym P|² and λm tr²(P). The stationarity condition with respect to a constant mode c is λe(s-p-c) = [2(µe+µm)+λm]c + λm p, where s=u_r'+u_r/r and p=P_rr+P_θθ, so c = [λe(s-p)-λm p]/[λe+2(µe+µm)+λm]. For the parameters in Table 1 this is generically nonzero, and a constant P_zz produces no Curl P penalty. Hence the true 3D long-cylinder minimizer is not the solution of Eqs. (29)-(34) unless the constraint P_zz=0 is imposed and justified. The statement that plane-strain conditions are warranted concerns only the displacement field, not the independent microdistortion field, and the use of the 2D bulk moduli in Eq. (37) effectively adopts the planar model rather than proving a 3D reduction.","section":"§4.1, Eqs. (29)-(34) and Eq. (9)"},{"comment":"For a cylindrical lateral surface the tangent space is spanned by e_θ and e_z, but the consistent coupling boundary conditions are written only for τ=e_θ. Applying the general condition (7) with τ=e_z to the axisymmetric displacement gives P·e_z = Du·e_z = 0 at r=R, i.e., P_rz(R)=P_θz(R)=P_zz(R)=0. These additional conditions do not appear in the 2D formulation and cannot be recovered from the e_θ condition used to obtain Eq. (36). If the intended setting is the full 3D model, this is a missing boundary condition; if the intended setting is the planar 2D model, the paper should state this explicitly and should not present the solution as the long-cylinder limit of the 3D model.","section":"§3.2, Eqs. (23)-(28)"},{"comment":"The introduction claims that the general formulation of the relaxed micromorphic model in orthogonal curvilinear coordinates is absent from the literature, but Section 3 actually derives only the two-dimensional polar-coordinate formulation: Eqs. (8)-(9) restrict P to the r-θ plane and the Curl operator in Eq. (14) is the planar 2D curl. The title and abstract should be adjusted to state that the contribution is a planar polar-coordinate formulation, or the full 3D cylindrical-coordinate formulation should be derived. This is not merely a wording issue because the axisymmetric extension problem is framed for a 3D long cylinder rather than for a planar domain.","section":"§1 and §3"}],"minor_comments":[{"comment":"There is a typo: 'relaxed meromorphic model' should read 'relaxed micromorphic model'.","section":"§1, last paragraph"},{"comment":"The text says 'subtracting Eqs. (34) from (33)', but Eqs. (33) and (34) are the Pθr and Pθθ equations; the algebra appears to subtract different equations. Please check the equation numbering and clarify the manipulation.","section":"§4.1, Eq. (63)"},{"comment":"Several key algebraic reductions (e.g., from Eqs. (31)-(34) to Eqs. (40)-(47), and from Eqs. (54)-(60) to (61)-(62)) are described only as 'after simplification'. Given that the closed-form solution is the main result, a detailed derivation or a symbolic-verification appendix would substantially improve reproducibility.","section":"§4.1, Eqs. (29)-(68)"},{"comment":"The phrase 'Due to the length of the cylinder, plane-strain conditions are warranted' presents an assumption as a consequence; this should be flagged as an additional modeling assumption, especially because of the microdistortion field's independence.","section":"Figure 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The plane-strain reduction is the load-bearing gap. If the authors reframe the paper as an explicit planar 2D relaxed-micromorphic benchmark, the contribution is defensible and likely publishable after revision. I would not reject outright, but the current claim of solving the 3D long-cylinder problem needs either a rigorous reduction argument or a clearly narrowed claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper does two things: it derives the 2D relaxed micromorphic equations in polar coordinates and then solves an axisymmetric extension problem in closed form using modified Bessel functions. The derivation is long but plausible, and the final formulas are new as far as I can tell. The classical-elasticity limit cases come out right, which is a good internal check. If you work on the planar model, this gives you a concrete benchmark for finite element verification. That part is solid.\n\nThe main problem is the jump from the planar model to the long-cylinder problem. Section 3 is explicitly about the two-dimensional relaxed micromorphic model; P is a 2×2 tensor and u has only radial and angular components. Then Section 4 says that for a long cylinder, plane-strain conditions are warranted and drops all out-of-plane components of P, including Pzz. That does not follow. In the full 3D energy, Pzz enters through tr(Du−P) and trP, and a constant Pzz costs no CurlP energy. For fixed in-plane data, stationarity with respect to a constant Pzz gives a nonzero value for generic parameters, including the ones in Table 1. So the true 3D plane-strain minimizer has out-of-plane microdistortion, and Eqs. (68)–(70) are not the solution of the 3D long-cylinder problem. They may well be the solution of the explicitly 2D planar problem, but then the benchmark claim needs to be restated. This is a load-bearing gap, not a cosmetic one.\n\nThe novelty claim in Section 1 is also a bit stronger than the evidence. The paper says the general curvilinear-coordinate formulation is absent from the literature, but Refs. [30] and [31] already solve cylindrical torsion and bending problems in cylindrical coordinates. Those are special problems, not the full general formulation, so the gap is narrower than claimed. I would soften the wording.\n\nMinor issues: there is no independent verification of the closed form against finite elements or any other numerical method; the “numerical results” are just evaluations of the derived expressions. Some algebraic steps are skipped, especially in the reduction from (29)–(34) to the Bessel equation. That is not fatal for a paper like this, but a supplementary page of key manipulations would help referees trust the algebra.\n\nWho is this for? Researchers using the relaxed micromorphic model who want an analytical benchmark for planar polar problems. It deserves a serious referee, not a desk rejection. The derivation is substantial and the flaw is fixable: either add a rigorous reduction argument showing when Pzz vanishes, or explicitly frame the solution as a planar 2D benchmark. I would send it to review with the expectation of major revision.","headline":"A useful polar-coordinate formulation and closed-form solution for the planar relaxed micromorphic model, but the advertised 3D long-cylinder benchmark is not supported because the plane-strain reduction ignores out-of-plane microdistortion.","tokens_in":14020,"tokens_out":2943,"would_cite":false,"duration_ms":27314,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["74A35","74G05","33C10"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives the relaxed micromorphic model in polar coordinates and solves an axisymmetric extension problem exactly, obtaining closed-form displacement and microdistortion fields.","keywords":["relaxed micromorphic model","polar coordinates","axisymmetric extension","closed-form solution","modified Bessel functions","consistent coupling boundary condition","size-dependent elasticity","generalized continua"],"falsifier":"Measure the interior radial displacement of a long cylinder under uniform radial boundary displacement, with the model's parameters first calibrated independently from bending or wave-propagation tests; formula (68) predicts a specific profile controlled by the modified Bessel functions and by $R/L_c$, so a systematic mismatch across radii would falsify the model or the boundary condition. A cheaper check is to solve the full two-field variational problem numerically without the $u_\\theta=0$ ansatz and compare whether the computed minimizer indeed has vanishing angular displacement and the closed-form field.","tokens_in":13017,"feed_emoji":"🧮","tokens_out":6537,"duration_ms":55194,"temperature":0.7,"pith_summary":"The paper fills a gap: until now the relaxed micromorphic model had no general formulation in orthogonal curvilinear coordinates. The authors write the isotropic model's governing equations in polar coordinates and solve an elastostatic axisymmetric extension problem, a long circular cylinder under uniform radial boundary displacement. The solution is closed form, built from modified Bessel functions, and uses the consistent coupling boundary condition $P\\cdot\\tau = Du\\cdot\\tau$ on the boundary. They show that the classical linear elasticity solution reappears as the characteristic length $L_c$ tends to zero or infinity, and use the formulas to study how material parameters and $L_c$ shift the displacement profile. If correct, the result supplies an analytical benchmark for checking and calibrating finite element implementations of the relaxed micromorphic model.","feed_headline":"Polar-coordinate solution turns micromorphic model into a benchmark","feed_subtitle":"New polar-coordinate equations give explicit displacement and microdistortion fields, with classical elasticity as a limit.","key_machinery":"The machinery is the polar-coordinate representation of the relaxed micromorphic model's differential operators, together with the consistent coupling boundary condition. The paper writes the gradient, divergence, and curl of the relevant tensor fields in $(r,\\theta)$ coordinates, producing the governing equations (12)--(20). For the axisymmetric problem it introduces the combinations $X = du_r/dr - P_{rr}$, $Y = u_r/r - P_{\\theta\\theta}$, and $Z = P_{\\theta\\theta}+P_{rr}$, reducing the system to a single inhomogeneous Bessel-type equation $d^2Z/dr^2 + (1/r)\\,dZ/dr - aZ + b = 0$; the solution enters the displacement through $du_r/dr + u_r/r = C_1 A + D_1 B I_0(\\sqrt{a}r)$. The boundary condition $P\\cdot\\tau = Du\\cdot\\tau$ fixes the remaining constants and eliminates the off-diagonal microdistortion.","core_discovery":"The central claim is that the isotropic relaxed micromorphic model can be formulated explicitly in polar coordinates and that, for the axisymmetric extension of a long cylinder, its equilibrium equations reduce to a linear system of ODEs whose solution is closed form. Under the axisymmetric ansatz ($u_\\theta=0$, all fields depending only on $r$), the displacement $u_r(r)$ and the diagonal microdistortion components $P_{rr}$, $P_{\\theta\\theta}$ are expressed through $I_0$ and $I_1$ modified Bessel functions, while the off-diagonal components $P_{r\\theta}$ and $P_{\\theta r}$ vanish under the consistent coupling boundary condition. In the limits $L_c\\to 0$ and $L_c\\to\\infty$, the solution degenerates to the classical linear elastic field $u_r = U_0 r/R$, with explicit expressions for $P_{\\theta\\theta}$ and $P_{rr}$. The paper also shows that the Cosserat couple modulus $\\mu_c$ drops out of this axisymmetric problem, so the result isolates the size-dependent effects carried by the characteristic length and the micro-moduli.","pith_inferences":["A natural extension is to repeat the derivation in cylindrical or spherical coordinates, where the same operator identities are available, giving benchmarks for geometries closer to actual metamaterial specimens.","The solution implicitly predicts a size-dependent effective stiffness for the cylinder: for a fixed boundary displacement, the interior displacement profile changes with $R/L_c$, so measuring interior strains on cylinders of different radii could test the model without requiring full-field measurements.","Replacing the consistent coupling boundary condition by $P\\cdot\\tau = 0$ would produce a different solution, so the gap between the two closed forms is a quantitative measure of how much the boundary condition on $P$ matters.","The closed form assumes $u_\\theta=0$ and radial-only dependence; solving the full two-field variational problem numerically without that ansatz would certify whether the true minimizer is indeed axisymmetric."],"forward_implications":["If the solution is correct, it gives a benchmark for finite element codes: with the listed parameter sets, an implementation should reproduce the closed-form displacement and microdistortion profiles.","The limits $L_c\\to 0$ and $L_c\\to\\infty$ both recover classical linear elasticity, so code validation can target those endpoints as separate checks.","The Cosserat couple modulus does not affect the axisymmetric extension solution, which simplifies parameter identification for this geometry.","The explicit dependence on $\\lambda_m$ and $\\mu_m$ shows how micro-moduli bend the displacement profile relative to the classical field, with the model predicting smaller or larger displacements depending on the parameter ratio $\\beta_1$.","The closed-form expression makes the role of the characteristic length $L_c$ explicit, so it can be used to study size effects without solving the full boundary-value problem numerically."],"supporting_citations":[{"why":"defines the two-field energy functional whose Euler-Lagrange equations are solved.","marker":"[40]"},{"why":"provides the operator form of the balance equations (3) that the paper rewrites in polar coordinates.","marker":"[21]"},{"why":"introduces the consistent coupling boundary condition P·τ = Du·τ used to close the axisymmetric problem.","marker":"[7]"},{"why":"supplies the polar-coordinate expression for Curl P used in the moment tensor.","marker":"[29]"},{"why":"gives the polar-coordinate displacement gradient and the classical elasticity solution recovered as a limit.","marker":"[34]"},{"why":"gives the plane-strain relations between bulk micro-moduli and Lamé-type moduli used in the reduction.","marker":"[14]"},{"why":"gives the macro-micro modulus relations in plane strain used in the parameter definitions.","marker":"[23]"}],"fun_headline_variants":["Exact polar solution turns micromorphic model into a benchmark","Closed-form axisymmetric solution for relaxed micromorphic model","Polar coordinates yield exact micromorphic fields","Explicit polar solution provides micromorphic benchmark","Micromorphic axisymmetric problem solved exactly with Bessel functions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes the consistent coupling boundary condition $P\\cdot\\tau = Du\\cdot\\tau$ on the boundary and the axisymmetric plane-strain ansatz $u_\\theta=0$ with all fields depending only on $r$; if the physical boundary condition on $P$ differs, or the true minimizer develops angular or off-axis structure, the closed-form solution will not describe that setting.","fun_headline_variants_meta":{"raw":{"variants":["Exact polar solution turns micromorphic model into a benchmark","Closed-form axisymmetric solution for relaxed micromorphic model","Polar coordinates yield exact micromorphic fields","Explicit polar solution provides micromorphic benchmark","Micromorphic axisymmetric problem solved exactly with Bessel functions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000718,"raw_usage":{"total_tokens":3191,"prompt_tokens":876,"completion_tokens":2315,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":492,"completion_tokens_details":{"reasoning_tokens":2240}},"tokens_in":492,"tokens_out":2315,"duration_ms":15951,"temperature":1.0,"reasoning_tokens":2240,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:23:51.598773+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the interior radial displacement of a long cylinder under uniform radial boundary displacement, with the model's parameters first calibrated independently from bending or wave-propagation tests; formula (68) predicts a specific profile controlled by the modified Bessel functions and by $R/L_c$, so a systematic mismatch across radii would falsify the model or the boundary condition. A cheaper check is to solve the full two-field variational problem numerically without the $u_\\theta=0$ ansatz and compare whether the computed minimizer indeed has vanishing angular displacement and the closed-form field.","supporting_citations":[{"cited_title":"Primal and mixed finite element formulations for the relaxed micromorphic model","cited_arxiv_id":null,"evidence_quote":"defines the two-field energy functional whose Euler-Lagrange equations are solved."},{"cited_title":"Band gaps in the relaxed linear micromorphic continuum","cited_arxiv_id":null,"evidence_quote":"provides the operator form of the balance equations (3) that the paper rewrites in polar coordinates."},{"cited_title":"The consistent coupling boundary condition for the classical micromorphic model: existence, uniqueness and interpretation of parameters","cited_arxiv_id":null,"evidence_quote":"introduces the consistent coupling boundary condition P·τ = Du·τ used to close the axisymmetric problem."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the polar-coordinate expression for Curl P used in the moment tensor."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the polar-coordinate displacement gradient and the classical elasticity solution recovered as a limit."},{"cited_title":"Green’s functions for the isotropic planar relaxed micromorphic model—Concentrated force and concentrated couple","cited_arxiv_id":null,"evidence_quote":"gives the plane-strain relations between bulk micro-moduli and Lamé-type moduli used in the reduction."},{"cited_title":"Neff and S","cited_arxiv_id":null,"evidence_quote":"gives the macro-micro modulus relations in plane strain used in the parameter definitions."}],"review_version":1}