Prueba de la ley de Ampère a partir de la ley de Biot-Savart para distribuciones de corriente tridimensionales

Supongamos la validez de la ley de Biot-Savart para una distribución tridimensional de corriente:

B ( X ) = m 0 4 π V j ( y ) × X y X y 3 d 3 y

dónde j es la densidad de la corriente distribuida en la región V . ¿Cómo podríamos entonces probar la ley de Ampère, en la forma

γ B d X = m 0 I vinculado o × B = m 0 j ?
Preferiría una prueba que use solo cálculo de múltiples variables , sin hacer uso de Dirac. d , que es una herramienta matemática, bastante delicada de usar, que requeriría mucha explicación para que entienda las igualdades donde aparece (soy principiante, como puedes leer en mi perfil, pero solo quiero entender por qué las leyes físicas son derivada de una manera matemáticamente correcta y no me interesan los atajos a menos que conozca su demostración matemática) y cuya utilización a veces se encuentra en los libros de texto me causa enormes problemas cuando el texto no dice, por ejemplo, por qué conmuta signos integrales y de diferenciación, como y × . De todos modos, estaría muy agradecido a cualquiera que produzca alguna prueba, incluso usando el d , pero explicando claramente qué teoremas y resultados matemáticos justifican pasos delicados como aquellos donde los operadores integrales (definidos explícitamente como Riemann, Lebesgue, signos simbólicos para distribuciones u otros) y diferenciales conmutan. Agradezco de corazón a cualquier contestador. Sugerencias de pruebas con tales requisitos que se pueden encontrar en línea o en libros impresos (mejor si no se utiliza el d ) son muy bienvenidos! Agradezco de corazón a cualquier contestador.

PD: Una derivación de la ley de Ampère de la de Biot-Savart está aquí , pero abarca el caso de una distribución lineal de corriente , que es diferente , como pensé que era descaradamente obvio antes de ver votos cerrados "duplicados".


Ensayos : Aunque he estado buscando una prueba rigurosa y comprensible (para mí) durante más de un mes, no he podido encontrar ni producir una. La derivación más común de la ley de Ampère de Biot-Savart utiliza esencialmente los cálculos que se encuentran en el esquema de prueba de Wikipedia, que son incomprensibles para mí: conozco algunas propiedades de Dirac d como el hecho de que, si ( j 1 , j 2 , j 3 ) = j C ( R 3 ) es compatible de forma compacta, entonces

1 4 π R 3 yo 2 j i ( yo ) yo r d m yo = j i ( r ) =: d ( yo r ) j i ( yo ) d 3 yo
donde la integral de la izquierda es del tipo de Lebesgue, mientras que el signo integral de la derecha es solo una notación simbólica para un funcional lineal, pero tengo tremendos problemas para comprender los pasos que se encuentran en ese esquema de prueba, por las razones que explique aquí en detalle para los posibles respondedores que intentan modificarlo, agregando las explicaciones de los hechos matemáticos utilizados para justificar las conmutaciones entre y y del significado exacto de los signos integrales, para producir una prueba más detallada y completa. Solo para hacer un breve resumen de los problemas que encontré en el esquema de prueba de Wikipedia, como me lo solicitó un comentarista, no entiendo qué significan esos signos integrales (integrales de Lebesgue, signos simbólicos para funcionales o qué más), tampoco Entiendo cuáles son los componentes derivados de × B son: dado que los teoremas como el de Stokes generalmente se aplican al integrar × B , yo creo que son las derivadas ordinarias del cálculo multivariado elemental, pero entonces el d , que es una herramienta de la teoría de distribuciones, aparece en el bosquejo de la demostración , y en la teoría de distribuciones existen derivadas de distribuciones que son una cosa muy diferente, pero se toman, hasta donde yo sé, con respeto a las variables escritas como "variables de integración" en la notación integral de distribución, mientras que, aquí, el esquema de la prueba comienza con r × B con r , mientras que la integral aparece con d 3 yo ...

He movido todos los comentarios al chat para mantener la publicación limpia. Parecía que algunos de los comentarios habían quedado obsoletos (p. ej., los comentarios que ya se abordaron en las ediciones), pero no pude distinguir fácilmente cuáles estaban obsoletos, así que los moví a todos.
¿Hay algún motivo en particular por el que haya hecho esta pregunta aquí en lugar de en Math.SE? Siento que obtendrías una mejor respuesta allí.
@Javier Gracias por tu comentario y sugerencia!!! Bueno, debido a que esas leyes se estudian... en física, por lo tanto, estoy seguro de que estoy en el tema aquí... Ni siquiera estoy seguro de cuánto un estudiante o erudito de matemáticas puras estaría familiarizado con el tema. .
Si la percepción física fuera lo que quisieras, probablemente no estarías haciendo esta pregunta. Y la mayoría de los matemáticos probablemente hayan oído hablar de las ecuaciones de Maxwell, especialmente los que saben cómo responder a tu pregunta.
Solo como referencia futura, @Self-teachingworker, trate de evitar editar sus publicaciones demasiadas veces. En general, debe guardar sus cambios y, cuando desee editar la publicación, corrija todo lo que pueda encontrar para corregirlo todo a la vez. Idealmente, no debería tener que editar una sola publicación más de 3 o 4 veces.

Respuestas (2)

Creo que tengo una prueba para ti (aunque puede que no te resulte satisfactoria). En su mayor parte, estoy siguiendo la prueba en la página de wikipedia a la que se vincula. Sin embargo, evito la función delta de dirac. A partir de la ley de Biot-Savart:

B ( r ) = m 0 4 π V j ( yo ) × r yo r yo 3 d 3 yo

Una de sus críticas a la demostración fue que no se especifica el tipo de integral. Las funciones que necesitamos integrar son integrables de Riemann (y, por lo tanto, también integrables de Lebesgue), por lo que con el fin de elegir pruebas particulares, procederé según la definición de Riemann (aunque esta es una elección arbitraria de mi parte).

Un problema que ya encontramos es que la ley de Biot-Savart es una integral impropia. Resolvemos este problema diciendo que la integral es el valor del principio de Cauchy, que se define en términos de un límite (esto es relevante).

Hago la sustitución:

r yo r yo 3 = r ( 1 r yo )

B ( r ) = m 0 4 π V j ( yo ) × r ( 1 r yo ) d 3 yo

Usando una identidad de cálculo vectorial estándar:

B ( r ) = m 0 4 π V r × ( j ( yo ) r yo ) d 3 yo

A continuación, vamos a sacar el operador rotacional fuera de la integral. Dado que tanto el rotacional como la integral se definen en términos de límites, esto equivale a intercambiar el orden de los límites, lo que generalmente es aceptable dados ciertos criterios de convergencia (por ejemplo, convergencia uniforme . Hay algunos teoremas de convergencia diferentes que pueden ser apropiados).

B ( r ) = m 0 4 π r × V ( j ( yo ) r yo ) d 3 yo

Aplicando rotacional a ambos lados de la ecuación, y la identidad de cálculo vectorial para el rotacional de un rotacional:

r × B ( r ) = m 0 4 π r V r ( j ( yo ) r yo ) d 3 yo m 0 4 π V r 2 ( j ( yo ) r yo ) d 3 yo

Asumiré que está familiarizado con el argumento de por qué el primer término integral es 0. Esto deja:

r × B ( r ) = m 0 4 π V j ( yo ) r 2 ( 1 r yo ) d 3 yo

Al evaluar el operador de Laplace, estoy usando la expresión típica para coordenadas esféricas, pero desplazada por yo .

r 2 ( 1 r yo ) = 1 r yo 2 r ( r yo 2 r 1 r yo )

Notará que esto se evalúa como 0 en todas partes, excepto en r = yo , donde la expresión no está definida. Dado que la expresión dentro de la integral se evalúa como 0 en todas partes excepto en r = yo , podemos hacer esta sustitución en el término actual y moverlo fuera de la integral.

r × B ( r ) = m 0 4 π j ( r ) V r 2 ( 1 r yo ) d 3 yo

Para evaluar la integral final, aplicamos el teorema de la divergencia.

V r 2 ( 1 r yo ) d 3 yo = V norte ^ r ( 1 r yo ) d 2 yo = V norte ^ ( r yo r yo 3 ) d 2 yo

Debido a que el valor dentro de la integral de volumen es 0 en todas partes excepto en r = yo , podemos elegir que el volumen sea esféricamente simétrico con respecto a yo , y finalmente evaluar la integral.

V norte ^ ( r yo r yo 3 ) d 2 yo = 4 π

Esto nos deja con la ley de Ampere.

× B = m 0 j

Editar: como se señaló, si tuviéramos que usar nuestra definición de integral, la respuesta sería 0. ¿Recuerda que se permite intercambiar límites dados ciertos criterios de convergencia? Si las integrales o integrandos no convergen, entonces intercambiar los límites es un error basado en las formulaciones que he elegido. En este caso, el integrando no converge.

Sin embargo, hemos demostrado que excluyendo el punto divergente, la integral es 0. Podemos dividir nuestra integral original en dos dominios, V B ( r , ϵ ) , y B ( r , ϵ ) , que es una bola de radio ϵ centrado en r . Ya hemos demostrado que la integral sobre el primer dominio es 0. Con base en esto:

r × B ( r ) = m 0 4 π r × B ( r , ϵ ) r × ( j ( yo ) r yo ) d 3 yo

La densidad de corriente es constante en el límite de ϵ 0 , por lo que podemos sacarlo fuera de la integral. Además, aquí la integral/integrando converge, considerando el valor del principio de Cauchy, por lo que podemos mover todas las derivadas fuera de la integral. Los pasos son básicamente los mismos que antes, excepto que todos los operadores derivados están fuera de la integral. Esto conduce a la siguiente ecuación:

r × B ( r ) = m 0 4 π j ( r ) r 2 B ( r , ϵ ) ( 1 r yo ) d 3 yo

La única diferencia con respecto a nuestra prueba inicial es que el operador laplaciano está fuera de la integral. Aquí se puede evaluar la integral, seguida del operador laplaciano. Esto debería producir 4 π , terminando la prueba.

Una sutil dificultad para evaluar esta integral, aunque hemos escrito el dominio como B ( r , ϵ ) , esto es un abuso de notación. El dominio no es realmente una función de r , por lo que debe permanecer constante a efectos de evaluar r 2 . De lo contrario, obtienes el extraño resultado de que la integral no es una función de r , y toda la expresión se convierte en 0.

¡Te lo agradezco de corazón! Veo que mi pregunta era comprensible. Si no me equivoco, quieres decir B ( r ) = límite ϵ 0 V B ( r , ϵ ) r × ( j ( yo ) r yo ) d 3 yo : ¿tengo razón? El primer problema que encuentro es: ¿podrías enunciar el hecho matemático que nos permita ver que
límite ϵ 0 V B ( r , ϵ ) r × ( j ( yo ) r yo ) d 3 yo = r × ( límite ϵ 0 V B ( r , ϵ ) j ( yo ) r yo d 3 yo ) ?
veo eso
V B ( r , ϵ ) r × ( j ( yo ) r yo ) d 3 yo = r × V B ( r , ϵ ) j ( yo ) r yo d 3 yo porque si V R norte es compacto y F : V × [ a , b ] R , ( X , t ) F ( X , t ) es tal que F t C ( V × [ a , b ] ) , entonces
t [ a , b ] d d t V F ( X , t ) d norte X = V F ( X , t ) t d norte X
Entiende lo que quiero decir con valor principal. Me di cuenta de que cometí algunos errores descuidados. Dado que algunas de las cantidades con las que estoy tratando no convergen, es un error formal intercambiar ciertos límites, por lo que las integrales vuelven a estar mal definidas. Intentaré arreglar esto cuando tenga la oportunidad, y leeré sus comentarios con más detalle.
@Self-teachingDavide Acabo de terminar la prueba. No creo que haya cometido ningún error esta vez. Este problema es un ejemplo específico del uso de una función de Green para resolver una ecuación diferencial, por lo que podría interesarle leer sobre ellos. Este tipo de cosas también se conoce como una función de respuesta de impulso de un sistema lineal. Un resultado integral divergente similar en matemáticas complejas son las relaciones de Kramers-Kronig. Toda esta clase de problemas está estrechamente relacionada con las transformadas de Fourier, que también tienen interesantes problemas de convergencia. Suerte en tu objetivo de autodidacta en matemáticas y ciencias.
¡Gracias de nuevo! El primer paso donde tropiezo es que no sé qué resultado nos permite ver que
r × [ límite ϵ 0 m 0 4 π B ( r , ϵ ) r × ( j ( yo ) r yo ) d 3 yo ]
= límite ϵ 0 r × [ m 0 4 π B ( r , ϵ ) r × ( j ( yo ) r yo ) d 3 yo ]
¿Podría indicarlo, si es posible con un enlace a una prueba, o con una prueba?
Un teorema que nos permite cambiar esos límites es el teorema de convergencia uniforme. Wikipedia tiene una buena explicación en su página sobre la convergencia uniforme.
De hecho, puedes tomar el rotacional primero, luego dividir tu integral en sus dos dominios, en ese orden. Eso colocará naturalmente ese límite en el exterior, mostrando que los dos límites se conmutan.
Conozco el teorema de la convergencia dominada de Lebesgue, del cual se deriva que si F : X × [ C , d ] R , con X medible, y t [ C , d ] la función X F ( X , t ) es sumable (es decir F ( , t ) L 1 ( X ) ), y hay un barrio B ( t 0 , d ) de t 0 tal que, para casi todos X X y todo t B ( t 0 , d ) , | F ( X , t ) t | ϕ ( X ) dónde ϕ L 1 ( X ) , entonces
d d t X F ( X , t ) d m X = X F ( X , t ) t d m X .
En nuestro caso, por ejemplo con nuestro ilimitado r 2 [ j 1 r 2 yo 2 r F ] en lugar de F ( X , t ) t , que puede nuestro ϕ , tal que | F ( X , t ) t | ϕ ( yo ) L 1 ( V ) , ¿ser? ¡¡¡Gracias de nuevo!!!
Estoy publicando aquí lo que publiqué en el chat: me temo que las ecuaciones que escribí no son comprensibles sin LaTeX. No he entendido nada: ¿sabes qué resultado matemático permite conmutar y límite ϵ 0 V   B ( r , ϵ ) ¿O lo buscas como yo? En el primer caso, ¿qué es? De todos modos, una vez conmutados los signos nabla e integral, no entiendo por qué la densidad de corriente es constante en el límite de ϵ 0 , por lo que podemos tomarlo fuera de la integral : nunca me he encontrado con esa redacción, lo que supongo que significa que [...]
[...] hay un barrio de r tal que j es constante en ese vecindario, que no creo que esté en nuestras suposiciones (en ese caso sería una función constante en R 3 , sin desaparecer en el infinito). Entonces, ¿cómo se calcula el límite de r 2 V   B ( r , ϵ ) r yo 1 d 3 yo como 4 π ? [...]
[...] Aquí noto curiosamente que no somos libres de conmutar operadores de límite/integrales y diferenciales porque, como ϵ 0 , r 2 V   B ( r , ϵ ) r yo 1 d 3 yo 0 ( límite ϵ 0 0 = 0 )...
El argumento para mudarse j fuera de la integral es que el dominio de nuestra integral tiende a un solo punto. Aunque puede haber pequeñas variaciones en la corriente sobre este dominio, a medida que el dominio se hace más pequeño, j se aproxima a un solo valor (el valor en ese punto).
No necesita conmutar ningún límite en el paso final. ¿Intentaste evaluar la integral? ¿Qué resultado obtuviste para eso?
¿Qué resultado obtuviste para eso? Con la esperanza de que estoy entendiendo su notación (ha olvidado algunos límite ϵ 0 en algún lugar -que, de nuevo, no he entendido qué teorema nos permite conmutar con los operadores diferenciales en cualquier paso- y algunos " V " en el fondo de ?), veo que (ya que ϵ > 0 Creo que podemos aplicar la regla de Leibniz ) [...]
[..] r 2 V B ( r , ϵ ) 1 r yo d 3 yo = V B ( r , ϵ ) r 2 ( 1 r yo ) d 3 yo = 0 y por lo tanto
límite ϵ 0 r 2 V B ( r , ϵ ) 1 r yo d 3 yo = límite ϵ 0 V B ( r , ϵ ) r 2 ( 1 r yo ) d 3 yo = 0
y en este caso podemos escribir j ( r ) o cualquier otra cosa antes de la integral pero todavía es 0 (y su límite también).[..]
[...] si tuviste la intención límite ϵ 0 estar después r 2 , entonces V B ( r , ϵ ) 1 r yo d 3 yo (y su límite) dependería de la forma de V (y no estoy seguro de qué teorema nos permitiría escribir j ( r ) afuera). Perdóneme si le sugiero que revise su prueba e introduzca notaciones explícitas para los límites, pero creo que sería mucho más claro...
En cuanto a los muchos límite / conmutaciones que usa, si conoce una prueba de ellas (¿usted o está suponiendo que se cumplen y buscando una prueba como yo?), podría estar interesado en responder esto con algo como "si pretende que la integral sea la límite de una integral de Riemann límite ϵ 0 V B ( r , ϵ ) , entonces se prueba la legitimidad de la conmutación tal y tal". ¡Muchas gracias de nuevo!
Creo que has entendido muy mal la prueba, así que déjame repasar los pasos con palabras, para asegurarme de que estamos en la misma página. Primero, tomamos el rizo de la ley de Biot-Savart. Debido a problemas de convergencia, subdividimos el dominio de la integral resultante en V B ( r , ϵ ) y B ( r , ϵ ) , e intente evaluar cada integral por separado. La integral sobre el dominio V B ( r , ϵ ) evalúa a 0. Al evaluar la integral sobre el dominio B ( r , ϵ ) , podemos mover la corriente fuera de la integral. ¿Estos pasos tienen sentido para ti?
¡Gracias! ¿Qué tipo de integral se define en B ( r , ϵ ) ? No puede ser à la Riemann porque el integrando no está definido, ni acotado, si no puede ser Lebesgue porque, admitiendo que podemos conmutar todo a nuestro antojo, sería 0...
La integral sigue siendo Riemann. El integrando no está acotado, pero está definido, siempre que no tratemos de mover operadores diferenciales a la integral. Teniendo en cuenta la longitud de mi respuesta, y también este prolongado ir y venir en los comentarios, actualizaré mi respuesta para que sea menos ambigua/confusa, luego no responderé más preguntas en los comentarios. La filosofía de stackexchange es que las preguntas/comentarios sean ampliamente útiles, y estos comentarios no están logrando eso. Creo que mi respuesta es correcta y lo suficientemente clara como para que usted tenga la responsabilidad de intentar comprenderla.
Esa página aborda integrales múltiples de funciones acotadas específicamente. No dice que las integrales múltiples de funciones ilimitadas no estén definidas en el sentido de Riemann. Considere esta página sobre integrales singulares.
Una de las definiciones más utilizadas de integrales de Riemann es esta (por R 2 pero generalizable a R norte ) y, en ese caso, la integral no existe finita si F es ilimitado porque, en ese caso, un F ( X k , y k ) podría elegirse arbitrariamente grande. Las integrales de Riemann impropias se pueden definir como límites de tales integrales, que pensé que eran las integrales que usó al tomar los límites de V B ( r , ϵ ) como ϵ 0 [...]
[...] (con los derivados en algún lugar antes o después límite : No he entendido eso de tu respuesta). En cambio, ahora veo que calculas V = V B ( r , ϵ ) + B ( r , ϵ ) dónde B ( r , ϵ ) es una especie de integral de Riemann definida de una manera que no sé [...]
[...] (excluyo que tenga la intención B ( r , ϵ ) ser una integral de Riemann impropia límite d 0 B ( r , ϵ ) B ( r , d ) porque no le veo sentido a separarse B ( r , ϵ ) de V y calculando V B ( r , ϵ ) + límite d 0 ( B ( r , ϵ ) B ( r , d ) ) en lugar de límite d 0 V B ( r , d ) ya que son lo mismo)
tengo la intención de B ( r , ϵ ) ser una integral de Riemann impropia límite d 0 B ( r , ϵ ) B ( r , d ) . El propósito de esto es permitirnos mover la densidad de corriente fuera de la integral, mientras que la integral permanece bien definida en el sentido de Riemann. Tal vez puedas encontrar una forma más elegante de hacerlo sin el límite adicional.

Espero haber podido encontrar una prueba usando integrales de Lebesgue, donde las diferenciaciones bajo los signos integrales estén justificadas, si no me equivoco, por resultados matemáticos probados.

Dejar j : R 3 R 3 ser de clase C 2 ( R 3 ) y que su soporte esté contenido en el subconjunto compacto acotado V R 3 , m -medible según la medida habitual de Lebesgue m definido en R 3 .

Usaré el siguiente lema

Dejar φ : V R 3 R estar acotado y m y -medible, con m y como de costumbre 3 -medida de Lebesgue dimensional, donde V es acotado y medible (según la misma medida). Definamos, para todos X R 3 ,

Φ ( X ) := V φ ( y ) X y d m y
entonces Φ C 1 ( R 3 ) y para k = 1 , 2 , 3 ,
X R 3 Φ ( X ) X k = V X k [ φ ( y ) X y ] d m y = V φ ( y ) y k X k X y 3 d m y

lo cual se demuestra aquí .

Definamos

A ( X ) := m 0 4 π V j ( y ) X y d m y = m 0 4 π R 3 j ( y ) X y d m y
donde notamos que el punto y = X no impide, por cualquier X R 3 , la sumabilidad de y X y 1 en cualquier subconjunto medible acotado de R 3 , como podemos verificar usando coordenadas esféricas alrededor X y teniendo en cuenta las propiedades de la integral de Lebesgue y su relación con la integral de Riemann.

Por tanto, si definimos B ( X ) := m 0 4 π V j ( y ) × X y X y 3 d m y = m 0 4 π R 3 X × [ j ( y ) X y ] d m y , yo diría, aplicando el lema anterior con φ = m 0 4 π j i ( i = 1 , 2 , 3 ), eso

× A = B
y, mediante el uso de una identidad conocida para el rizo del rizo , que se puede aplicar desde este argumento , con F ( z ) = z 1 y gramo = m 0 4 π j i ( i = 1 , 2 , 3 ), demuestra que
A C 2 ( R 3 ) ,
Vemos, aplicando de nuevo este argumento , que
X × B ( X ) = m 0 4 π X [ X R 3 j ( y ) X y d m y ] m 0 4 π X 2 R 3 j ( y ) X y d m y
= m 0 4 π X R 3 y j ( y ) X y d m y m 0 4 π R 3 y 2 j ( y ) X y d m y
donde el primer apéndice es nulo si X R 3 X j ( X ) = 0 , que se cumple para una corriente estacionaria, y donde creo que se puede probar la segunda adición, análogamente a lo que se hace aquí para una C función de soporte compacto, para ser
m 0 4 π R 3 y 2 j ( y ) X y d m y = m 0 j ( X ) =: m 0 d ( y X ) j ( y ) d 3 y
que obviamente no es m 0 4 π R 3 X 2 [ X y 1 ] j ( y ) d m y 0 en general; Digo esto para aclarar la notación de un paso utilizado en el esquema de prueba de Wikipedia , que es una derivación muy común de la ley de Ampère de la ley de Biot-Savart, que también aparece en la Electrodinámica clásica de Jackson .

Notemos también, con referencia a ese esquema de prueba , que en el término X R 3 j ( y ) X y d m y podemos diferenciar, utilizando el lema citado anteriormente, con respecto a los componentes de X bajo el signo integral y por lo tanto

R 3 y j ( y ) X y d m y = X R 3 j ( y ) X y d m y
= R 3 j ( y ) X [ 1 X y ] d m y
así que supongo que una posible interpretación de las integrales de ese esquema de prueba es leerlas como integrales de Lebesgue excepto por la integral donde 2 ( 1 | r yo | ) aparece